Dear All,
Does anyone know how to combine MCMC results from BiSSE analyses (used 100
trees to address phylogenetic uncertainty)?
testsamples <- lapply(lik.lambda, mcmc, c(0.0380925074, 0.0274189791,
0.0128575318, 0.0009811349), nsteps=100, prior=prior.exp, w=0.1,
print.every=10)
> testsamplesSub *(post-burnin)*
[[1]]
i lambda0 lambda1 mu0 q10 p
21 21 0.03470906 0.02580640 0.011623783 0.0007190859 -3638.357
22 22 0.03645545 0.02369784 0.011018956 0.0011557427 -3638.152
23 23 0.03604715 0.02193379 0.009545690 0.0007547081 -3638.385
...
[[2]]
i lambda0 lambda1 mu0 q10 p
21 21 0.03470906 0.02580640 0.011623783 0.0007190859 -3638.357
22 22 0.03645545 0.02369784 0.011018956 0.0011557427 -3638.152
23 23 0.03604715 0.02193379 0.009545690 0.0007547081 -3638.385
...
...
[[100]]
i lambda0 lambda1 mu0 q10 p
21 21 0.03470906 0.02580640 0.011623783 0.0007190859 -3638.357
22 22 0.03645545 0.02369784 0.011018956 0.0011557427 -3638.152
23 23 0.03604715 0.02193379 0.009545690 0.0007547081 -3638.385
...
Sincerely, LY
[[alternative HTML version deleted]]
_______________________________________________
R-sig-phylo mailing list - R-sig-phylo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
Searchable archive at http://www.mail-archive.com/r-sig-phylo@r-project.org/
