Prof. Paulo Justiniano, boa tarde! Os dados são do dataset JURA mesmo, compreendi as suas explicações. Com relação ao ajuste por verossimilhança tenho duas dúvidas: Não foi especificado o argumento cov.model neste caso, em qual modelo ele ajusta? O valor do lambda deveria ser lam=-0.2, com o valor positivo não estava conseguindo plotar o resultado.
Obrigado pela ajuda, esse grupo tem sido muito importante para minha evolução com o R. *Hélder Gramacho * Recife-PE / *[email protected] <[email protected]>* Em 5 de maio de 2014 10:08, Paulo Justiniano <[email protected]>escreveu: > os dados (dataset JURA?) são EXTREMAMENTE assimétricos e o variograma > usual nao vai ter bom comportamento e não é surpresa que > ajuste o modelo com efeito pepita puro > > seguem algums comandos com algumas ideias explorando uma possivel > trensformação de variáveis > > 2 comentarios adicionais: > 1. se o variograma é omnidirecional (default) o argumento tolerance nao é > usada > > 2. no eixo dis do veriograma (ao final) o xlab deve ser distancia (e nao > alcance) > > > > > ## > plot(cobre) > summary(cobre) > ## > require(MASS) > boxcox(cobre, lam=seq(-1, 1, len=50)) > > ## > var.cu0 <- variog(cobre,uvec=seq(0.1,3,l=10), estimador.type="classical", > lam=-0.2, max.dist=2) > plot(var.cu0, xlab="distância", ylab="semivariância",main='Cobre',font.main= > 3) > > ## ajuste de variogramas.... talvez manos adequado.... > exp.ols.cu0<-variofit(var.cu0,ini=c(0.2, 0.5),weights="equal", > cov.model="exp") > exp.ols.cu0 > plot(var.cu0, xlab='Distância', ylab='Semivariância', main='Semivariograma > OLS > - Exponencial') > lines(exp.ols.cu0, col="blue") > summary(exp.ols.cu0) > > > ## ajuste por verossimilhança (eu considero preferivel) > cu0.ml <- likfit(cobre, ini=c(0.5, 0.1), lam=0.2) > cu0.ml > > PJ > > > On Fri, 2 May 2014, Helder Gramacho wrote: > > 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, >> >> 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 / [email protected] >> >> >> > _______________________________________________ > R-br mailing list > [email protected] > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. >
_______________________________________________ R-br mailing list [email protected] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
