<div dir="ltr">Prof. Paulo Justiniano, boa tarde!<div><br></div><div>Os dados são do dataset JURA mesmo, compreendi as suas explicações.</div><div>Com relação ao ajuste por verossimilhança tenho duas dúvidas:</div><div>Não foi especificado o argumento cov.model neste caso, em qual modelo ele ajusta?</div>
<div>O valor do lambda deveria ser lam=-0.2, com o valor positivo não estava conseguindo plotar o resultado.</div><div><br></div><div>Obrigado pela ajuda, esse grupo tem sido muito importante para minha evolução com o R.</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 5 de maio de 2014 10:08, Paulo Justiniano <span dir="ltr"><<a href="mailto:paulojus@leg.ufpr.br" target="_blank">paulojus@leg.ufpr.br</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
os dados (dataset JURA?) são EXTREMAMENTE assimétricos e o variograma usual nao vai ter bom comportamento e não é surpresa que<br>
ajuste o modelo com efeito pepita puro<br>
<br>
seguem algums comandos com algumas ideias explorando uma possivel trensformação de variáveis<br>
<br>
2 comentarios adicionais:<br>
1. se o variograma é omnidirecional (default) o argumento tolerance nao é usada<br>
<br>
2. no eixo dis do veriograma (ao final) o xlab deve ser distancia (e nao alcance)<br>
<br>
<br>
<br>
<br>
##<br>
plot(cobre)<br>
summary(cobre)<br>
##<br>
require(MASS)<br>
boxcox(cobre, lam=seq(-1, 1, len=50))<br>
<br>
##<br>
var.cu0 <- variog(cobre,uvec=seq(0.1,3,l=<u></u>10), estimador.type="classical",<br>
lam=-0.2, max.dist=2)<br>
plot(var.cu0, xlab="distância", ylab="semivariância",main='<u></u>Cobre',font.main= 3)<br>
<br>
## ajuste de variogramas.... talvez manos adequado....<br>
exp.ols.cu0<-variofit(var.cu0,<u></u>ini=c(0.2, 0.5),weights="equal",<br>
cov.model="exp")<br>
exp.ols.cu0<br>
plot(var.cu0, xlab='Distância', ylab='Semivariância', main='Semivariograma OLS<br>
- Exponencial')<br>
lines(exp.ols.cu0, col="blue")<br>
summary(exp.ols.cu0)<br>
<br>
<br>
## ajuste por verossimilhança (eu considero preferivel)<br>
<a href="http://cu0.ml" target="_blank">cu0.ml</a> <- likfit(cobre, ini=c(0.5, 0.1), lam=0.2)<br>
<a href="http://cu0.ml" target="_blank">cu0.ml</a><span class="HOEnZb"><font color="#888888"><br>
<br>
PJ</font></span><div class=""><br>
<br>
On Fri, 2 May 2014, Helder Gramacho wrote:<br>
<br>
</div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="">
Boa noite pessoal,<br>
<br>
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<br>
possível realizar o ajuste, ficando como na figura abaixo.<br>
Será que se trata de algum problema com o código que estou executando, ou os dados é que não permitem mesmo o ajuste.<br>
Reproduzi o código abaixo da figura e o arquivo está no Dropbox.<br>
<br>
Desde já agradeço qualquer ajuda,<br>
<br></div><div><div class="h5">
Imagem inline 1<br>
# Download de um arquivo do Dropbox<br>
<br>
links <- c("<a href="https://www.dropbox.com/s/8blubgumcsss834/cobre.csv" target="_blank">https://www.dropbox.com/s/<u></u>8blubgumcsss834/cobre.csv</a>")<br>
<br>
tokens <- gsub("^.*/s/","",dirname(<u></u>links))<br>
fileNames <- basename(links)<br>
newLinks <- file.path("<a href="http://dl.dropbox.com/s" target="_blank">http://dl.dropbox.<u></u>com/s</a>", tokens, fileNames);<br>
newLinks<br>
<br>
for (a in newLinks) {<br>
tryCatch(download.file(a, dest=basename(a), mode='wb'),<br>
error=function(...) print("Falha no download!"))}<br>
<br>
cobre<-read.table(file="cobre.<u></u>csv", sep=",", header=T, dec=".")<br>
attach(cobre)<br>
cobre<-as.geodata(cbind(cobre$<u></u>V1,cobre$V2,cobre$V3))<br>
<br>
# Semivariograma Experimental<br>
par(mfrow=c(1,1),xpd=F)<br>
<a href="http://var.cu" target="_blank">var.cu</a><-variog(cobre,uvec=seq(<u></u>0.1,3,l=10),pairs.min=30,<u></u>estimador.type="classical", direction="omnidirectional",<u></u>tolerance=pi/8)<br>
plot(<a href="http://var.cu" target="_blank">var.cu</a>, xlab="distância", ylab="semivariância",main='<u></u>Cobre',font.main = 3)<br>
<br>
#-----------------------------<u></u>------------------------------<u></u>------------------------------<u></u>---------------<br>
# Ajustando Semivariograma Teórico pelo Método dos Mínimos Quadrados<br>
#-----------------------------<u></u>------------------------------<u></u>------------------------------<u></u>---------------<br>
# Modelo exponencial<br>
<a href="http://exp.ols.cu" target="_blank">exp.ols.cu</a><-variofit(<a href="http://var.cu" target="_blank">var.cu</a>,<u></u>ini=c(235, 0.41),weights="equal", cov.model="exp")<br>
<a href="http://exp.ols.cu" target="_blank">exp.ols.cu</a><br>
plot(<a href="http://var.cu" target="_blank">var.cu</a>, xlab='Alcance', ylab='Semivariância', main='Semivariograma OLS - Exponencial')<br>
lines(<a href="http://exp.ols.cu" target="_blank">exp.ols.cu</a>, col="blue")<br>
summary(<a href="http://exp.ols.cu" target="_blank">exp.ols.cu</a>)<br>
<br>
## Modelo esférico<br>
<a href="http://sph.ols.cu" target="_blank">sph.ols.cu</a><-variofit(<a href="http://var.cu" target="_blank">var.cu</a>,<u></u>ini=c(235, 0.41),weights="equal", cov.model="sph")<br>
<a href="http://sph.ols.cu" target="_blank">sph.ols.cu</a><br>
plot(<a href="http://var.cu" target="_blank">var.cu</a>, xlab='Alcance', ylab='Semivariância', main='Semivariograma OLS - Esférico')<br>
lines(<a href="http://sph.ols.cu" target="_blank">sph.ols.cu</a>, col="blue")<br>
summary(<a href="http://sph.ols.cu" target="_blank">sph.ols.cu</a>)<br>
<br>
## Modelo gaussiano<br>
<a href="http://gaus.ols.cu" target="_blank">gaus.ols.cu</a><-variofit(<a href="http://var.cu" target="_blank">var.cu</a>,<u></u>ini=c(235, 0.41),weights="equal", cov.model="gaus")<br>
<a href="http://gaus.ols.cu" target="_blank">gaus.ols.cu</a><br>
plot(<a href="http://var.cu" target="_blank">var.cu</a>, xlab='Alcance', ylab='Semivariância', main='Semivariograma OLS - Gaussiano')<br>
lines(<a href="http://gaus.ols.cu" target="_blank">gaus.ols.cu</a>, col="blue")<br>
summary(<a href="http://gaus.ols.cu" target="_blank">gaus.ols.cu</a>)<br>
Hélder Gramacho <br>
Recife-PE / <a href="mailto:agrohelder@gmail.com" target="_blank">agrohelder@gmail.com</a><br>
<br>
<br>
</div></div></blockquote>
<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></div>