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@gmail.com <agrohelder@hotmail.com>*
