Re: [R-sig-phylo] 3. partial correlation with gls residuals? (Tom Schoenemann)

2012-03-14 Thread Emmanuel Paradis
Hi Alejandro, I agree with almost all your points. A few additional comments below. Correlation gives the correlation among the estimated regression parameters: it is computed from the variance-covariance matrix extracted from the fitted object, model.fit$varBeta or vcov(model.fit) -- the

Re: [R-sig-phylo] 3. partial correlation with gls residuals? (Tom Schoenemann)

2012-03-13 Thread Alejandro Gonzalez
Hello Tom, I will answer some of your questions below, Emmanuel and others may be better placed to reply to some of your questions: On 12, Mar 2012, at 7:49 PM, Tom Schoenemann wrote: Thanks Rob and Alejandro, OK, I did as suggested and ran a PGLS with A ~ B + C. I was hoping for some

[R-sig-phylo] 3. partial correlation with gls residuals? (Tom Schoenemann)

2012-03-12 Thread Robert Barton
Dear Tom, There is no reason to assume that the residuals from your two PGLS analyses will be independent of phylogeny, so if you are going to do this you should correlate the residuals phylogenetically (i.e. run them through PGLS). General problems with using residuals as data have been