Não foi possível enviar o arquivo. Será algum problema com as permissões?
Essa é uma revisão anterior do documento!
Simulation study
Test beta bias
gs <- expand.grid((0:10)/10, (0:10)/10)
niter <- 100
arr <- array(NA, dim=c(niter, 2, 2), dimnames=list(niter=1:niter, stat=c("beta","beta.var"), dep=c(1,0)))
# gau
set.seed(111)
gd <- grf(grid=gs, cov.pars=c(1,0.2), mean=2, nsim=niter)
for(i in 1:niter){
cat(i)
ggdd <- gd
ggdd$data <- gd$data[,i]
lf <- likfit(ggdd, ini.cov.pars=c(1,0.2), messages=FALSE)
arr[i,,1] <- c(lf$beta, lf$beta.var)
}
# gau independente
set.seed(111)
gd <- grf(grid=gs, cov.pars=c(0,0), mean=2, nsim=niter, nugget=1)
for(i in 1:niter){
cat(i)
ggdd <- gd
ggdd$data <- gd$data[,i]
lf <- likfit(ggdd, ini.cov.pars=c(0.5,0.5), messages=F)
arr[i,,2] <- c(lf$beta, lf$beta.var)
}
pdf("betabias.pdf")
par(mfrow=c(2,2))
hist(arr[,1,1], main="beta for mu=2 s2=1 phi=0.2 t2=0")
hist(arr[,1,2], main="beta for mu=2 s2=0 phi=0 t2=1")
hist(arr[,2,1], main="beta.var")
hist(arr[,2,2], main="beta.var")
dev.off()
O problema está na variância do beta que não é um estimador da variância do processo. Logo quando fazemos a backtrans há uma subestimação da média do processo.
