Caros Hélder e Éder, Pelo visto a variável não tem dependencia espacial, é o que se percebe nos semirariogramas Hélio Em 3 de maio de 2014 01:18, Hélder [via R-br] < ml-node+s2285057n4662046h94@n4.nabble.com> escreveu:
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] <http://user/SendEmail.jtp?type=node&node=4662046&i=0> *
Em 2 de maio de 2014 12:22, Éder Comunello <[hidden email]<http://user/SendEmail.jtp?type=node&node=4662046&i=1>
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]<http://user/SendEmail.jtp?type=node&node=4662046&i=2>[hidden email] <http://user/SendEmail.jtp?type=node&node=4662046&i=3>> Dourados, MS - [22 16.5'S, 54 49'W]
Em 1 de maio de 2014 23:37, Helder Gramacho <[hidden email]<http://user/SendEmail.jtp?type=node&node=4662046&i=4>
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 / * [hidden email] <http://user/SendEmail.jtp?type=node&node=4662046&i=5>*
_______________________________________________ R-br mailing list [hidden email] <http://user/SendEmail.jtp?type=node&node=4662046&i=6> 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] <http://user/SendEmail.jtp?type=node&node=4662046&i=7> 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] <http://user/SendEmail.jtp?type=node&node=4662046&i=8> 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-semivariograma... To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=3357982&code=aGVsaW9nYWxsb3JvY2hhQGdtYWlsLmNvbXwzMzU3OTgyfC0xMzQ3NTkwMDY4> . NAML<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml>
-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho