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/