<div dir="ltr">Caros Hélder e Éder,<div><br></div><div>Pelo visto a variável não tem dependencia espacial, é o que se percebe nos semirariogramas</div><div><br></div><div>Hélio</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">
Em 3 de maio de 2014 01:18, Hélder [via R-br] <span dir="ltr"><<a href="mailto:ml-node+s2285057n4662046h94@n4.nabble.com" target="_blank">ml-node+s2285057n4662046h94@n4.nabble.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">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 class=""><div><div>link <- "<a href="http://dl.dropbox.com/s/8blubgumcsss834/cobre.csv" rel="nofollow" link="external" target="_blank">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><div>cobre<-read.table(file="cobre.csv", sep=",", header=T, dec=".")[,-1]</div><div class=""><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" rel="nofollow" link="external" 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" rel="nofollow" link="external" target="_blank">var.cu</a>, xlab="distância", ylab="semivariância",main='Cobre',font.main = 3)</div><div><br></div></div><div>#---------------------------------------------------------------------------------------------------------</div>
<div class="">
<div># Ajustando Semivariograma Teórico pelo Método dos Mínimos Quadrados</div></div><div>#----------------------------------------------------------------------------------------------------------</div><div class="">
<div># Modelo exponencial</div>
<div><a href="http://exp.ols.cu" rel="nofollow" link="external" target="_blank">exp.ols.cu</a><-variofit(<a href="http://var.cu" rel="nofollow" link="external" target="_blank">var.cu</a>,ini=c(235, 0.41),weights="equal", cov.model="exp")</div>
<div><a href="http://exp.ols.cu" rel="nofollow" link="external" target="_blank">exp.ols.cu</a></div>
<div>plot(<a href="http://var.cu" rel="nofollow" link="external" target="_blank">var.cu</a>, xlab='Alcance', ylab='Semivariância', main='Semivariograma OLS - Exponencial')</div><div>lines(<a href="http://exp.ols.cu" rel="nofollow" link="external" target="_blank">exp.ols.cu</a>, col="blue")</div>
<div>summary(<a href="http://exp.ols.cu" rel="nofollow" link="external" target="_blank">exp.ols.cu</a>)</div><div><br></div><div>## Modelo esférico</div><div><a href="http://sph.ols.cu" rel="nofollow" link="external" target="_blank">sph.ols.cu</a><-variofit(<a href="http://var.cu" rel="nofollow" link="external" target="_blank">var.cu</a>,ini=c(235, 0.41),weights="equal", cov.model="sph")</div>
<div><a href="http://sph.ols.cu" rel="nofollow" link="external" target="_blank">sph.ols.cu</a></div><div>plot(<a href="http://var.cu" rel="nofollow" link="external" target="_blank">var.cu</a>, xlab='Alcance', ylab='Semivariância', main='Semivariograma OLS - Esférico')</div>
<div>lines(<a href="http://sph.ols.cu" rel="nofollow" link="external" target="_blank">sph.ols.cu</a>, col="blue")</div>
<div>summary(<a href="http://sph.ols.cu" rel="nofollow" link="external" target="_blank">sph.ols.cu</a>)</div><div><br></div><div>## Modelo gaussiano</div><div><a href="http://gaus.ols.cu" rel="nofollow" link="external" target="_blank">gaus.ols.cu</a><-variofit(<a href="http://var.cu" rel="nofollow" link="external" target="_blank">var.cu</a>,ini=c(235, 0.41),weights="equal", cov.model="gaus")</div>
<div><a href="http://gaus.ols.cu" rel="nofollow" link="external" target="_blank">gaus.ols.cu</a></div><div>plot(<a href="http://var.cu" rel="nofollow" link="external" target="_blank">var.cu</a>, xlab='Alcance', ylab='Semivariância', main='Semivariograma OLS - Gaussiano')</div>
<div>lines(<a href="http://gaus.ols.cu" rel="nofollow" link="external" target="_blank">gaus.ols.cu</a>, col="blue")</div>
<div>summary(<a href="http://gaus.ols.cu" rel="nofollow" link="external" target="_blank">gaus.ols.cu</a>)</div></div><div>### </code><br></div><div><br></div><div>obrigado,</div><div><br></div></div></div><div class="gmail_extra">
<div class=""><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="http://user/SendEmail.jtp?type=node&node=4662046&i=0" rel="nofollow" link="external" target="_blank">[hidden email]</a></i></div>
</i></div><div style="text-align:center"><br></div></div></div>
<br><br></div><div class="gmail_quote">Em 2 de maio de 2014 12:22, Éder Comunello <span dir="ltr"><<a href="http://user/SendEmail.jtp?type=node&node=4662046&i=1" rel="nofollow" link="external" target="_blank">[hidden email]</a>></span> escreveu:<div>
<div class="h5"><br><blockquote style="border-left:2px solid #cccccc;padding:0 1em" class="gmail_quote">
<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" rel="nofollow" link="external" 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="http://user/SendEmail.jtp?type=node&node=4662046&i=2" rel="nofollow" link="external" target="_blank">[hidden email]</a><a href="http://user/SendEmail.jtp?type=node&node=4662046&i=3" rel="nofollow" link="external" target="_blank">[hidden email]</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="http://user/SendEmail.jtp?type=node&node=4662046&i=4" rel="nofollow" link="external" target="_blank">[hidden email]</a>></span> escreveu:<br>
<blockquote style="border-left:2px solid #cccccc;padding:0 1em" class="gmail_quote">
<div><div>
<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" rel="nofollow" link="external" 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" rel="nofollow" link="external" 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" rel="nofollow" link="external" 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" rel="nofollow" link="external" 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" rel="nofollow" link="external" target="_blank">exp.ols.cu</a><-variofit(<a href="http://var.cu" rel="nofollow" link="external" target="_blank">var.cu</a>,ini=c(235, 0.41),weights="equal", cov.model="exp")</div>
<div><a href="http://exp.ols.cu" rel="nofollow" link="external" target="_blank">exp.ols.cu</a></div>
<div>plot(<a href="http://var.cu" rel="nofollow" link="external" target="_blank">var.cu</a>, xlab='Alcance', ylab='Semivariância', main='Semivariograma OLS - Exponencial')</div><div>lines(<a href="http://exp.ols.cu" rel="nofollow" link="external" target="_blank">exp.ols.cu</a>, col="blue")</div>
<div>summary(<a href="http://exp.ols.cu" rel="nofollow" link="external" target="_blank">exp.ols.cu</a>)</div><div><br></div><div>## Modelo esférico</div><div><a href="http://sph.ols.cu" rel="nofollow" link="external" target="_blank">sph.ols.cu</a><-variofit(<a href="http://var.cu" rel="nofollow" link="external" target="_blank">var.cu</a>,ini=c(235, 0.41),weights="equal", cov.model="sph")</div>
<div><a href="http://sph.ols.cu" rel="nofollow" link="external" target="_blank">sph.ols.cu</a></div><div>plot(<a href="http://var.cu" rel="nofollow" link="external" target="_blank">var.cu</a>, xlab='Alcance', ylab='Semivariância', main='Semivariograma OLS - Esférico')</div>
<div>lines(<a href="http://sph.ols.cu" rel="nofollow" link="external" target="_blank">sph.ols.cu</a>, col="blue")</div>
<div>summary(<a href="http://sph.ols.cu" rel="nofollow" link="external" target="_blank">sph.ols.cu</a>)</div><div><br></div><div>## Modelo gaussiano</div><div><a href="http://gaus.ols.cu" rel="nofollow" link="external" target="_blank">gaus.ols.cu</a><-variofit(<a href="http://var.cu" rel="nofollow" link="external" target="_blank">var.cu</a>,ini=c(235, 0.41),weights="equal", cov.model="gaus")</div>
<div><a href="http://gaus.ols.cu" rel="nofollow" link="external" target="_blank">gaus.ols.cu</a></div><div>plot(<a href="http://var.cu" rel="nofollow" link="external" target="_blank">var.cu</a>, xlab='Alcance', ylab='Semivariância', main='Semivariograma OLS - Gaussiano')</div>
<div>lines(<a href="http://gaus.ols.cu" rel="nofollow" link="external" target="_blank">gaus.ols.cu</a>, col="blue")</div>
<div>summary(<a href="http://gaus.ols.cu" rel="nofollow" link="external" 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="http://user/SendEmail.jtp?type=node&node=4662046&i=5" rel="nofollow" link="external" target="_blank">[hidden email]</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="http://user/SendEmail.jtp?type=node&node=4662046&i=6" rel="nofollow" link="external" target="_blank">[hidden email]</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" link="external" 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" rel="nofollow" link="external" 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="http://user/SendEmail.jtp?type=node&node=4662046&i=7" rel="nofollow" link="external" target="_blank">[hidden email]</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" link="external" 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" rel="nofollow" link="external" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div></div></div>
<br></div><div><div class="h5">
<br>_______________________________________________
<br>R-br mailing list
<br><a href="http://user/SendEmail.jtp?type=node&node=4662046&i=8" rel="nofollow" link="external" target="_blank">[hidden email]</a>
<br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" link="external" 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" rel="nofollow" link="external" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.
<br>
<br>
<hr noshade size="1" color="#cccccc">
</div></div><div style="color:#444;font:12px tahoma,geneva,helvetica,arial,sans-serif"><div><div class="h5">
<div style="font-weight:bold">If you reply to this email, your message will be added to the discussion below:</div>
</div></div><a href="http://r-br.2285057.n4.nabble.com/R-br-Problemas-no-ajuste-de-semivariogramas-teoricos-tp4662040p4662046.html" target="_blank">http://r-br.2285057.n4.nabble.com/R-br-Problemas-no-ajuste-de-semivariogramas-teoricos-tp4662040p4662046.html</a>
</div><div class="HOEnZb"><div class="h5">
<div style="color:#666;font:11px tahoma,geneva,helvetica,arial,sans-serif;margin-top:.4em;line-height:1.5em">
To unsubscribe from R-br, <a href="http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=3357982&code=aGVsaW9nYWxsb3JvY2hhQGdtYWlsLmNvbXwzMzU3OTgyfC0xMzQ3NTkwMDY4" target="_blank">click here</a>.<br>
<a href="http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml" rel="nofollow" style="font:9px serif" target="_blank">NAML</a>
</div></div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>Hélio Gallo Rocha<br>IFSULDEMINAS - Câmpus Muzambinho<br>
</div>