Hello, I ran a PGLS analysis like this: model = gls(y ~ x, correlation = corBrownian(phy = tree, form=~species), data = data_df, method = "ML")
I was able to extract the coefficients (intercept + b1 in this case), as well as the relevant p-values. What I couldn't figure out is how to get a correlation coefficient between y and x, something like R-squared or pearson's coefficient. Can you please help with that? Thanks! [[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/