[R-sig-phylo] p-value of PGLS models in caper

Hi all, I'm having a hard time figuring out how to extract the p-value from a PGLS model fitted in caper. When I run summary(mod) and print it, I get, among other things: >summ=summary(mod) >summ Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept)0.0179553

Re: [R-sig-phylo] p-value of PGLS models in caper

It looks like the p-value is calculated as follows: pf(summ$fstatistic[1], summ$fstatistic[2], summ$fstatistic[3], lower.tail = FALSE) You can format it like it is in the summary by wrapping that in format.pval(). Best, Will William Gearty williamgearty.com On Tue, Dec 10, 2019 at 1:33 PM