search for: powerb1

...ame error message. The bolded part is the part that doesn't work. Any help would be appreciated! reddmcmc=glsm.mcmc(geodata, coords=geodata$coords, data=geodata$data,model=list(family="poisson", cov.pars=c(50,8000),beta=c(-26,-.04,6.94),cov.model="spherical",nugget=5,trend=~powerb1+temp1),mcmc.input = mcmc.control(S.scale =.000001 ,thin=10)) prepareredd1=prepare.likfit.glsm(reddmcmc) lik.redd1.1 <- likfit.glsm(prepareredd1, ini.phi = 10,cov.model = "spherical",trend=~powerb1 + temp1) # predict at new locations trendpred=trend.spatial(~habitat2$powerb + habitat2$...