Hi Ties.
pgls is in the package 'caper' not 'geiger'. You can fit the same models
using gls in the nlme package using correlation structures from ape.
It's possible that this implementation might better protect you against
regions of the likelihood space that are not defined. (Without a
Hello,
I have a question about an error message that occurs sporadically when I
run multiple PGLS models using the pgls function in geiger.
I am performing all possible models on 15 independent variables (~ 30, 000
models).
Whilst I am looping across the models I get the following error message:
