Displaying 1 result from an estimated 1 matches for "prepareredd1".
2011 Oct 28
0
problem with glsm.krige: trendd and trend l must have similar specifications error
...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$avgtemp)
test2 <- glsm.krige(reddmcmc,coords,trend.l=trendpred) # tren...