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] < [email protected]> 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-semivariogramas-teoricos-tp4662040p4662046.html > 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
_______________________________________________ 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.
