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 / 

[hidden email]

Em 2 de maio de 2014 12:22, Éder Comunello <[hidden email]> 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 <[hidden email][hidden email]>
Dourados, MS - [22 16.5'S, 54 49'W]

Em 1 de maio de 2014 23:37, Helder Gramacho <[hidden email]> 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,

# 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 / 

[hidden email]

_______________________________________________
R-br mailing list
[hidden email]
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
[hidden email]
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
[hidden email]
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.

---

If you reply to this email, your message will be added to the discussion below:

http://r-br.2285057.n4.nabble.com/R-br-Problemas-no-ajuste-de-semivariogramas-teoricos-tp4662040p4662046.html

To unsubscribe from R-br, click here.
NAML