[R-br] Problemas no ajuste de semivariogramas teóricos

[Éder Comunello]

Helder, bom dia!

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'.

Da forma como está estão ocorrendo muitas duplicatas e algumas com valores
muito distintos, o que impacta a estrutura de variação.

Sugiro verificar e substituir os dados se for o caso.


### <code r>
setwd("C:/LAB/RGIS/geostat")
require(geoR)

link <- "http://dl.dropbox.com/s/8blubgumcsss834/cobre.csv"
if (!file.exists(basename(link))) download.file(link, dest=basename(link),
mode='wb')

cobre <- read.table(file="cobre.csv", sep=",", header=T, dec=".")[,-1]
names(cobre) <- c('x','y','z'); head(cobre)
cobre <- as.geodata(cobre)
### ! as.geodata: 13 replicated data locations found.

attach(cobre)
cobre.dup <- dup.coords(cobre)
plot(coords); points(cobre.dup[,2:3], col=2, pch=20)

barplot(cobre.dup[,4], col=rep(2:7,each=2), pch=20)
### </code>

Éder Comunello <c <comunello.eder em gmail.com>omunello.eder em gmail.com>
Dourados, MS - [22 16.5'S, 54 49'W]


Em 1 de maio de 2014 23:37, Helder Gramacho <agrohelder em gmail.com> escreveu:

> 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,
>
> [image: 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 /
> * agrohelder em gmail.com <agrohelder em hotmail.com>*
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140502/751c15fb/attachment.html>