solution of boot.phylo()'s problem
0
0
Entering edit mode
Nora Muda ▴ 70
@nora-muda-2065
Last seen 10.2 years ago
Thank you to Martin and Earl for their response and answered. I got answered from Emmanuel Paradis. Below are the solutions: 'boot.phylo(phy,x,FUN,B=100,block=1)' x is a character matrix (not the distance matrix) that will be resampled for the bootstrap. My mistake is not to store the tree which is needed by boot.phylo; so tr <- nj(distallbacteria) Therefore; boot.phylo(tr,allbacteria,FUN=function(xx)nj(dist.dna(xx,pairwise.dele tion=TRUE) We also can refer to the book on using R for phylogenetic analyses, written by Emmanuel Paradis : http://www.springer.com/west/home/statistics/stats+life+sci?SGWID=4-10 134-22-156993712-0 --Nora. ______________________________________________________________________ ______________ TV dinner still cooling? Check out "Tonight's Picks" on Yahoo! TV.
• 921 views
ADD COMMENT
0
Entering edit mode

Dear Nora, I have a similar problem.

protein_seq <- readAAStringSet("./protein_seq.fasta")
msa <- msa(gene_seq, method = "Muscle")

msa_aln <- msaConvert(gene_seq, type="seqinr::alignment")

d <- dist.alignment(msa_aln, "identity")
tree <- nj(d)

bp <- boot.phylo(tree, msa, function(xx) 
  nj(dist.alignment(xx, "identity")))

I tried different inputs for the parameter x of boot.phylo (raw seq, msa object), but they did not work. Could you please give me an example of x's structure? Thank you very much for your help

ADD REPLY

Login before adding your answer.

Traffic: 356 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6