[R-br] Problemas no ajuste de semivariogramas teóricos

[Helder Gramacho]

Boa noite pessoal,

Gostaria de uma ajuda, estou tentando fazer o ajuste de semivariogramas
teóricos pelo método dos mínimos quadrados, entretanto, aparentemente não
está sendo possível realizar o ajuste, ficando como na figura abaixo.
Será que se trata de algum problema com o código que estou executando, ou
os dados é que não permitem mesmo o ajuste.
Reproduzi o código abaixo da figura e o arquivo está no Dropbox.

Desde já agradeço qualquer ajuda,

[image: Imagem inline 1]
# Download de um arquivo do Dropbox

links <- c("https://www.dropbox.com/s/8blubgumcsss834/cobre.csv")

tokens    <- gsub("^.*/s/","",dirname(links))
fileNames <- basename(links)
newLinks  <- file.path("http://dl.dropbox.com/s", tokens, fileNames);
newLinks

for (a in newLinks) {
  tryCatch(download.file(a, dest=basename(a), mode='wb'),
           error=function(...) print("Falha no download!"))}

cobre<-read.table(file="cobre.csv", sep=",", header=T, dec=".")
attach(cobre)
cobre<-as.geodata(cbind(cobre$V1,cobre$V2,cobre$V3))

# Semivariograma Experimental
par(mfrow=c(1,1),xpd=F)
var.cu<-variog(cobre,uvec=seq(0.1,3,l=10),pairs.min=30,estimador.type="classical",
direction="omnidirectional",tolerance=pi/8)
plot(var.cu, xlab="distância", ylab="semivariância",main='Cobre',font.main
= 3)

#--------------------------------------------------------------------------------------------------------
#      Ajustando Semivariograma Teórico pelo Método dos Mínimos Quadrados
#--------------------------------------------------------------------------------------------------------
# Modelo exponencial
exp.ols.cu<-variofit(var.cu,ini=c(235, 0.41),weights="equal",
cov.model="exp")
exp.ols.cu
plot(var.cu, xlab='Alcance', ylab='Semivariância', main='Semivariograma OLS
- Exponencial')
lines(exp.ols.cu, col="blue")
summary(exp.ols.cu)

## Modelo esférico
sph.ols.cu<-variofit(var.cu,ini=c(235, 0.41),weights="equal",
cov.model="sph")
sph.ols.cu
plot(var.cu, xlab='Alcance', ylab='Semivariância', main='Semivariograma OLS
- Esférico')
lines(sph.ols.cu, col="blue")
summary(sph.ols.cu)

## Modelo gaussiano
gaus.ols.cu<-variofit(var.cu,ini=c(235, 0.41),weights="equal",
cov.model="gaus")
gaus.ols.cu
plot(var.cu, xlab='Alcance', ylab='Semivariância', main='Semivariograma OLS
- Gaussiano')
lines(gaus.ols.cu, col="blue")
summary(gaus.ols.cu)
*Hélder Gramacho *
Recife-PE /
*agrohelder em gmail.com <agrohelder em hotmail.com>*
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140502/5a92a9c3/attachment.html>
-------------- Próxima Parte ----------
Um anexo não-texto foi limpo...
Nome: image.png
Tipo: image/png
Tamanho: 19445 bytes
Descrição: não disponível
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140502/5a92a9c3/attachment.png>