<div dir="ltr">Olá Éder, boa noite!<div><br></div><div>Realmente haviam coordenadas duplicadas, corrigi aumentando o número de casas decimais do arquivo *.csv como você sugeriu.</div><div>Entretanto, o ajuste do modelo teórico continua apresentando o mesmo problema, segue o código dos dados já corrigidos:</div>
<div><br></div><div>### <code r><br></div><div><div><div>link <- "<a href="http://dl.dropbox.com/s/8blubgumcsss834/cobre.csv">http://dl.dropbox.com/s/8blubgumcsss834/cobre.csv</a>"</div><div>if (!file.exists(basename(link))) download.file(link, dest=basename(link), mode='wb')</div>
</div><div><br></div><div>cobre<-read.table(file="cobre.csv", sep=",", header=T, dec=".")[,-1]</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>### </code><br></div><div><br></div><div>obrigado,</div><div><br></div></div></div><div class="gmail_extra"><br clear="all"><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>
<br><br><div class="gmail_quote">Em 2 de maio de 2014 12:22, Éder Comunello <span dir="ltr"><<a href="mailto:comunello.eder@gmail.com" target="_blank">comunello.eder@gmail.com</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr"><div>Helder, bom dia!</div><div><br></div><div>Peguei pra dar uma olhada, mas houve um problema com dados. Talvez seja o caso de aumentar o número de casas decimais utilizadas pra gerar o '.csv'.</div>
<div><br></div><div>Da forma como está estão ocorrendo muitas duplicatas e algumas com valores muito distintos, o que impacta a estrutura de variação.</div><div><br></div><div>Sugiro verificar e substituir os dados se for o caso.</div>
<div><br></div><div><br></div><div><font face="courier new, monospace">### <code r></font></div><div><div><font face="courier new, monospace">setwd("C:/LAB/RGIS/geostat")</font></div><div><font face="courier new, monospace">require(geoR)</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">link <- "<a href="http://dl.dropbox.com/s/8blubgumcsss834/cobre.csv" target="_blank">http://dl.dropbox.com/s/8blubgumcsss834/cobre.csv</a>"</font></div>
<div><font face="courier new, monospace">if (!file.exists(basename(link))) download.file(link, dest=basename(link), mode='wb')</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">cobre <- read.table(file="cobre.csv", sep=",", header=T, dec=".")[,-1]</font></div>
<div><font face="courier new, monospace">names(cobre) <- c('x','y','z'); head(cobre)</font></div><div><font face="courier new, monospace">cobre <- as.geodata(cobre)</font></div><div><font face="courier new, monospace">### ! as.geodata: 13 replicated data locations found.</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">attach(cobre)</font></div><div><font face="courier new, monospace">cobre.dup <- dup.coords(cobre)</font></div><div><font face="courier new, monospace">plot(coords); points(cobre.dup[,2:3], col=2, pch=20)</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">barplot(cobre.dup[,4], col=rep(2:7,each=2), pch=20)</font></div></div><div><font face="courier new, monospace">### </code></font></div>
<div class="gmail_extra">
<br clear="all"><div><div dir="ltr">Éder Comunello <<a href="mailto:comunello.eder@gmail.com" target="_blank">c</a><a href="mailto:omunello.eder@gmail.com" target="_blank">omunello.eder@gmail.com</a>> <br>Dourados, MS - [22 16.5'S, 54 49'W]<br>
</div></div>
<br><br><div class="gmail_quote">Em 1 de maio de 2014 23:37, Helder Gramacho <span dir="ltr"><<a href="mailto:agrohelder@gmail.com" target="_blank">agrohelder@gmail.com</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<div><div class="h5">
<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 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" target="_blank">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" target="_blank">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" target="_blank">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" target="_blank">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" target="_blank">exp.ols.cu</a><-variofit(<a href="http://var.cu" target="_blank">var.cu</a>,ini=c(235, 0.41),weights="equal", cov.model="exp")</div><div><a href="http://exp.ols.cu" target="_blank">exp.ols.cu</a></div>
<div>plot(<a href="http://var.cu" target="_blank">var.cu</a>, xlab='Alcance', ylab='Semivariância', main='Semivariograma OLS - Exponencial')</div><div>lines(<a href="http://exp.ols.cu" target="_blank">exp.ols.cu</a>, col="blue")</div>
<div>summary(<a href="http://exp.ols.cu" target="_blank">exp.ols.cu</a>)</div><div><br></div><div>## Modelo esférico</div><div><a href="http://sph.ols.cu" target="_blank">sph.ols.cu</a><-variofit(<a href="http://var.cu" target="_blank">var.cu</a>,ini=c(235, 0.41),weights="equal", cov.model="sph")</div>
<div><a href="http://sph.ols.cu" target="_blank">sph.ols.cu</a></div><div>plot(<a href="http://var.cu" target="_blank">var.cu</a>, xlab='Alcance', ylab='Semivariância', main='Semivariograma OLS - Esférico')</div>
<div>lines(<a href="http://sph.ols.cu" target="_blank">sph.ols.cu</a>, col="blue")</div>
<div>summary(<a href="http://sph.ols.cu" target="_blank">sph.ols.cu</a>)</div><div><br></div><div>## Modelo gaussiano</div><div><a href="http://gaus.ols.cu" target="_blank">gaus.ols.cu</a><-variofit(<a href="http://var.cu" target="_blank">var.cu</a>,ini=c(235, 0.41),weights="equal", cov.model="gaus")</div>
<div><a href="http://gaus.ols.cu" target="_blank">gaus.ols.cu</a></div><div>plot(<a href="http://var.cu" target="_blank">var.cu</a>, xlab='Alcance', ylab='Semivariância', main='Semivariograma OLS - Gaussiano')</div>
<div>lines(<a href="http://gaus.ols.cu" target="_blank">gaus.ols.cu</a>, col="blue")</div>
<div>summary(<a href="http://gaus.ols.cu" target="_blank">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>
<br></div></div>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div>