Hello!

I'm running single-species, single-season models using RPresence and would like to incorporate quadratic terms into my models. I've attempted to run models using the following structure for a quadratic term:

mod1<-occMod(model = list(psi~distcov+(distcov^2), p~1), data=data, type="so")

However, the resulting model estimates seem to discard the quadratic term I've tried to force and only report the linear term. Is there another way I can code for a quadratic term or do I have to manually generate the quadratic term and add it to my dataframe of covariates? Thank you!