Displaying 1 result from an estimated 1 matches for "trendpr".
Did you mean:
trendp
2011 Oct 28
0
problem with glsm.krige: trendd and trend l must have similar specifications error
....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) # trend.d and trend.l must have similar specifications ?????
test2 <- glsm.krige(reddmcmc,coords,trend.l=habitat2$powerb + habitat2$avgtemp) # same error message...
[[alternative HTM...