<div dir="ltr"><div>Boa noite pessoal,</div><div><br></div><div>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.</div>
<div>Será que se trata de algum problema com o código que estou executando, ou os dados é que não permitem mesmo o ajuste.</div><div>Reproduzi o código abaixo da figura e o arquivo está no Dropbox.</div><div><br></div><div>
Desde já agradeço qualquer ajuda,</div><div><br></div><div><img src="cid:ii_145bafc6d9121ac6" alt="Imagem inline 1" style="margin-right: 0px;" width="455" height="219"><br></div><div># Download de um arquivo do Dropbox<br>
</div><div><br></div><div>links <- c("<a href="https://www.dropbox.com/s/8blubgumcsss834/cobre.csv">https://www.dropbox.com/s/8blubgumcsss834/cobre.csv</a>")</div><div><br></div><div>tokens    <- gsub("^.*/s/","",dirname(links))</div>
<div>fileNames <- basename(links)</div><div>newLinks  <- file.path("<a href="http://dl.dropbox.com/s">http://dl.dropbox.com/s</a>", tokens, fileNames);</div><div>newLinks</div><div><br></div><div>for (a in newLinks) {</div>
<div>  tryCatch(download.file(a, dest=basename(a), mode='wb'),</div><div>           error=function(...) print("Falha no download!"))}</div><div><br></div><div>cobre<-read.table(file="cobre.csv", sep=",", header=T, dec=".")</div>
<div>attach(cobre)</div><div>cobre<-as.geodata(cbind(cobre$V1,cobre$V2,cobre$V3))</div><div><br></div><div># Semivariograma Experimental</div><div>par(mfrow=c(1,1),xpd=F)</div><div><a href="http://var.cu">var.cu</a><-variog(cobre,uvec=seq(0.1,3,l=10),pairs.min=30,estimador.type="classical", direction="omnidirectional",tolerance=pi/8)</div>
<div>plot(<a href="http://var.cu">var.cu</a>, xlab="distância", ylab="semivariância",main='Cobre',font.main = 3)</div><div><br></div><div>#--------------------------------------------------------------------------------------------------------</div>
<div>#      Ajustando Semivariograma Teórico pelo Método dos Mínimos Quadrados</div><div>#--------------------------------------------------------------------------------------------------------</div><div># Modelo exponencial</div>
<div><a href="http://exp.ols.cu">exp.ols.cu</a><-variofit(<a href="http://var.cu">var.cu</a>,ini=c(235, 0.41),weights="equal", cov.model="exp")</div><div><a href="http://exp.ols.cu">exp.ols.cu</a></div>
<div>plot(<a href="http://var.cu">var.cu</a>, xlab='Alcance', ylab='Semivariância', main='Semivariograma OLS - Exponencial')</div><div>lines(<a href="http://exp.ols.cu">exp.ols.cu</a>, col="blue")</div>
<div>summary(<a href="http://exp.ols.cu">exp.ols.cu</a>)</div><div><br></div><div>## Modelo esférico</div><div><a href="http://sph.ols.cu">sph.ols.cu</a><-variofit(<a href="http://var.cu">var.cu</a>,ini=c(235, 0.41),weights="equal", cov.model="sph")</div>
<div><a href="http://sph.ols.cu">sph.ols.cu</a></div><div>plot(<a href="http://var.cu">var.cu</a>, xlab='Alcance', ylab='Semivariância', main='Semivariograma OLS - Esférico')</div><div>lines(<a href="http://sph.ols.cu">sph.ols.cu</a>, col="blue")</div>
<div>summary(<a href="http://sph.ols.cu">sph.ols.cu</a>)</div><div><br></div><div>## Modelo gaussiano</div><div><a href="http://gaus.ols.cu">gaus.ols.cu</a><-variofit(<a href="http://var.cu">var.cu</a>,ini=c(235, 0.41),weights="equal", cov.model="gaus")</div>
<div><a href="http://gaus.ols.cu">gaus.ols.cu</a></div><div>plot(<a href="http://var.cu">var.cu</a>, xlab='Alcance', ylab='Semivariância', main='Semivariograma OLS - Gaussiano')</div><div>lines(<a href="http://gaus.ols.cu">gaus.ols.cu</a>, col="blue")</div>
<div>summary(<a href="http://gaus.ols.cu">gaus.ols.cu</a>)</div><div><div dir="ltr"><div style="text-align:center"><font size="1"><i><b>Hélder Gramacho </b></i></font></div><div style="text-align:center"><font size="1">Recife-PE / </font><i style="font-size:x-small;color:rgb(51,51,255)"><div style="display:inline!important">
<i><a href="mailto:agrohelder@hotmail.com" target="_blank">agrohelder@gmail.com</a></i></div></i></div><div style="text-align:center"><br></div></div></div>
</div>