# Exercício 8.2 require(geoR) sigma2<- 1 phi<-0.2 mi=0.5 n<-1000 X<-seq(n/2,(n^2-n/2),n) Y<-rep(0,length(X)) M<-cbind(X,Y) View(M) set.seed(171) dados<- grf(nrow(M),grid=M,cov.pars = c(sigma2,phi),nsim=1000, cov.model="exp", mean=mi) # Abaixo comandos para usar método da máxima verossimilhança ml1 <- likfit(dados, ini=c(sigma2, phi), cov.model="exp") ml1$beta # Comandos para visualizar a variação da média por histograma hist(apply(dados$data,2,mean))