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

[Helder Gramacho]

Olá Éder, boa noite!

Realmente haviam coordenadas duplicadas, corrigi aumentando o número de
casas decimais do arquivo *.csv como você sugeriu.
Entretanto, o ajuste do modelo teórico continua apresentando o mesmo
problema, segue o código dos dados já corrigidos:

### <code r>
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]
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)
### </code>

obrigado,


*Hélder Gramacho *
Recife-PE /
*agrohelder em gmail.com <agrohelder em hotmail.com>*



Em 2 de maio de 2014 12:22, Éder Comunello <comunello.eder em gmail.com>escreveu:

> 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.
>>
>
>
> _______________________________________________
> 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/20140503/a9b8527a/attachment.html>