I have still been unable to enter site-level covariates that differ by year. I have tried using the code in a previous post, but replacing the 'rbind' function with 'cbind' and still have not been successful. Here is the code I used:
- Code: Select all
eps.cov<-data.frame(HabYear=cbind(habcov[,3],habcov[,4],habcov[,5],habcov[,6]))
mod<- occMod(model=list(psi~1, gamma~1,epsilon~HabYear, p~1),data= EAME.pao, type= "do.1",cov.list=list(eps.cov=eps.cov))
"Error in eval(predvars,data,env): object 'HabYear' not found
I have also tried creating the data.frame from the matrix as an alternative and have not had any luck.
Any advice would be appreciated.
Thank you,
Devin