-
Notifications
You must be signed in to change notification settings - Fork 3
Open
Description
I'm SO excited to use this package and thank you so much for making it available.
In trying to run:
JRfilt = genefilter_sample(JRphylo, filterfun_sample(function(x) x > 5), A = 0.1 * nsamples(JRphylo))
JR1 = prune_taxa(JRfilt, JRphylo)
ravel.mcmc = DirFactor( t(otu_table(JR1)), hyper, step = 50, thinning = 10 )
ravel.all.res = lapply( paste( ravel.mcmc$save.path, seq( 20, 50, 10 ), sep = "_" ), readRDS )
ConvDiagnosis(ravel.all.res)
I get the error
[1] "yes"
Error: $ operator is invalid for atomic vectors
I've tried to understand the function using the argument
ConvDiagnosis(ravel.all.res[[1]][1])
But this yields the same error. I see the package uses CODA which requires mcmc.list and the first argument for ConvDiagnosis is a list - are these the same? Where is it located?
Please advise how to properly use this function!
Reactions are currently unavailable
Metadata
Metadata
Assignees
Labels
No labels