[R-sig-phylo] PIC or PGLS for genome-wide SNP screening

2012-05-23 Thread Mattia Prosperi
Dear all, I am working on a data set composed of bacterial genomic sequences (a few genes) associated to phenotypic values (in-vitro resistance to antibiotics, a numerical value discretised into a binary class). Of note, the bacterial isolates were sampled non-uniformly at different times and

Re: [R-sig-phylo] PIC or PGLS for genome-wide SNP screening

2012-05-23 Thread Theodore Garland Jr
Is this the way in which one decides for OLS vs. PGLS? If you have the same set of independent variables, then you just prefer the one (OLS or PGLS) with the higher likelihood. So far as I am told by Joe Felsenstein, you cannot do a ln maximum likelihood ratio test because the number of

[R-sig-phylo] General question about the impact of sampling bias in comparative and diversification analyses

2012-05-23 Thread Julien Lorion
Dear all, I have a 5-gene phylogenetic tree, dated, quite well resolved, including 50 species and covering 80~90% of the known diversity... I have also LASER diversification analyses and bayesian reconstruction of ancestral states for 5 characters , both taking into account the impact of

[R-sig-phylo] problem with fitcontinuos function

2012-05-23 Thread Agus Camacho
Dear all, fitcontinuous is giving me a warning generally related to non correspondence between tree tips and data. However name.check and the tree tips seems ok. Anybody has a hint on this? Thanks in advance. Agus name.check(tree, tp)[1] OK BMfit - fitContinuous(tree, tp$temperature,

Re: [R-sig-phylo] problem with fitcontinuos function

2012-05-23 Thread Renata Brandt
Hello Agus, I believe your problem will be solved if you set row names to the object tp$temperature. An example follows: names(temperature)-row.names(tp) Please let me know if it worked. Cheers, Renata -- Renata Brandt Departamento de Biologia - FFCLRP Universidade de São Paulo Ribeirão

Re: [R-sig-phylo] problem with fitcontinuos function

2012-05-23 Thread Agus Camacho
Thanks to all for your fast answer, this list is amazing! Graham's answer was perfectly satisfactory, once onnly that eliminated the warning. Thanks to all again. Agus 2012/5/23 Graham Slater gsla...@ucla.edu This doesn't indicate that they don't match - rather it's telling you that there