# \dontshow{
set.seed(25)
data(dataca20)
d1=dataca20[1:158,]
######covariance matern: kappa=0.3 prior:reference
res=tsroba(calcont~altitude+area, kappa = 0.3, data=d1,
ini.pars=c(10,3,10),iter=2,thin=1,burn=0)
datapred=dataca20[159:178,]
formula=calcont~altitude+area
xpred=model.matrix(formula,data=datapred)
tsrobapred(res,xpred=xpred,coordspred=dataca20[159:178,1:2])
# }
# \donttest{
set.seed(25)
data(dataca20)
d1=dataca20[1:158,]
######covariance matern: kappa=0.3 prior:reference
res=tsroba(calcont~altitude+area, kappa = 0.3, data=d1,
ini.pars=c(10,3,10),iter=50,thin=1,burn=5)
datapred=dataca20[159:178,]
formula=calcont~altitude+area
xpred=model.matrix(formula,data=datapred)
tsrobapred(res,xpred=xpred,coordspred=dataca20[159:178,1:2])
# }
Run the code above in your browser using DataLab