Abril 2012 - R-br - DINF UFPR

DUVIDA ANOVA
by Kalill V.Páscoa 04 Abr '12

04 Abr '12

Nos meus dados existem 13 tratamentos porém quando gero a ANOVA ela indica que existe apenas 1, onde estou errando? > dados$Trat [1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 [41] 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 5 5 5 5 5 5 5 5 [81] 5 5 5 5 5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 7 7 7 7 7 7 7 7 7 7 7 7 8 8 8 8 8 8 8 8 [121] 8 8 8 8 8 8 8 8 8 8 8 8 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 10 10 10 10 10 10 10 10 [161] 10 10 10 10 10 10 10 10 10 10 10 10 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 12 12 12 12 12 12 12 12 [201] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13> mDAP <- aov(DAP~Bloco+Trat, data=dados)> class(mDAP)[1] "aov" "lm" > anova(mDAP)Analysis of Variance Table Response: DAP Df Sum Sq Mean Sq F value Pr(>F) Bloco 1 0.0 0.022 0.0013 0.970882 Trat 1 114.0 113.992 6.7906 0.009754 ** Residuals 233 3911.3 16.787 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 -- *Kalill José Viana da Páscoa* *Engenheiro Florestal * *Mestrando em Engenharia Florestal Universidade Federal de Lavras* * * *Rua Geraldo dos Passos, 81 * *Água Limpa, **Lavras-MG* *CEP: 37200-000* *(35)3821-6833* *(35)9152-4828* *(35)9928-6574*

2 2

PROBLEMA NA ANOVA E NO SKOTT-KNOTT
by Kalill V.Páscoa 04 Abr '12

04 Abr '12

Boa noite, estou tentando fazer a Anova e o teste Sk porém no fim ele gera um erro, alguém faz ideia do que seja? #ANAVA# mDAP <- aov(DAP~Bloco+Trat, data=dados) class(mDAP) anova(mDAP) summary(mDAP) summary.lm(mDAP) #TESTE DE NORMALIDADE# shapiro.test(dados$DAP) shapiro.test(residuals(mDAP)) library(ScottKnott) #TESTE SCOTT-KNOTT# as.factor(dados$Bloco) as.factor(dados$Trat) sk.DAP=SK(x=dados, y=dados$DAP, model="dados$DAP~dados$Bloco+dados$Trat",which="dados$Trat",sig.level=0.05) summary(sk.DAP) Erro em SK.aov(x = av, which = which, id.trim = id.trim, sig.level = sig.level) : No factors in the fitted model!Além disso: Mensagens de aviso perdidas:1: In replications(paste("~", xx), data = mf) : non-factors ignored: dados$Bloco2: In replications(paste("~", xx), data = mf) : non-factors ignored: dados$Trat -- *Kalill José Viana da Páscoa* *Engenheiro Florestal * *Mestrando em Engenharia Florestal Universidade Federal de Lavras* * * *Rua Geraldo dos Passos, 81 * *Água Limpa, **Lavras-MG* *CEP: 37200-000* *(35)3821-6833* *(35)9152-4828* *(35)9928-6574*

1 0

(sem assunto)
by Estevão Prado 04 Abr '12

04 Abr '12

Olá pessoal, Preciso carregar o mapa de uma cidade no R, já baixei os arquivos (.dbf, .shp...), porém quando vou carregar o arquivo ele dá erro. require(maptools) getwd() setwd("/home/est/Área de Trabalho/maps") city <- readShapePoly("maps/4105805") # 4105805 é código da cidade e nome dos arquivos (.dbf, .shp...) Dei uma olhada no site: http://www.leg.ufpr.br/doku.php/pessoais:elias:cursos:seufop11, mas não consegui obter êxito. Se alguém puder me ajudar desde já agradeço.

5 6

Re: [R-br] abrir mapa shp
by Estevão Prado 04 Abr '12

04 Abr '12

Alexandro, Com o arquivo que você me enviou consegui visualizar. Obrigado. Date: Tue, 3 Apr 2012 17:33:49 -0700 From: vl.alexandro(a)yahoo.com.br Subject: Re: [R-br] abrir mapa shp To: estevao.prado(a)hotmail.com Segue anexo o shape file de spse vc abri-lo isto significa que o seu shape file está ... vamos testar esta hipotése!tudo bem! De: Estevão Prado <estevao.prado(a)hotmail.com> Para: Lista R-br <r-br(a)listas.c3sl.ufpr.br> Enviadas: Terça-feira, 3 de Abril de 2012 20:36 Assunto: Re: [R-br] abrir mapa shp Alexandro, Tentei o que você me sugeriu, mas o mesmo erro ainda esta ocorrendo. (Erro em getinfo.shape(filen) : Error opening SHP file). Mas obrigado mesmo assim. Date: Tue, 3 Apr 2012 15:23:57 -0700 From: vl.alexandro(a)yahoo.com.br To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] abrir mapa shp Tente o seguinte: library(maptools)cidade<- readShapePoly(file.choose()) #procure onde está o arquivo shp #espero ter ajudado!!!Alex De: Estevão Prado <estevao.prado(a)hotmail.com> Para: r-br(a)listas.c3sl.ufpr.br Enviadas: Terça-feira, 3 de Abril de 2012 18:49 Assunto: [R-br] (sem assunto) Olá pessoal, Preciso carregar o mapa de uma cidade no R, já baixei os arquivos (.dbf, .shp...), porém quando vou carregar o arquivo ele dá erro. require(maptools) getwd() setwd("/home/est/Área de Trabalho/maps") city <- readShapePoly("maps/4105805") # 4105805 é código da cidade e nome dos arquivos (.dbf, .shp...) Dei uma olhada no site: http://www.leg.ufpr.br/doku.php/pessoais:elias:cursos:seufop11, mas não consegui obter êxito. Se alguém puder me ajudar desde já agradeço. _______________________________________________ R-br mailing list R-br(a)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(a)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(a)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.

1 0

Re: [R-br] Tobit espacial
by Diego andre 04 Abr '12

04 Abr '12

Caro Leonard de Assis, nestes pacotes que você citou tem o tobit, mas não na sua forma espacial (SPATIAL TOBIT), que é o que eu estou procurando. Mais uma vez, obrigado pela disposição de ajudar. -- Diego André Doutorando em Economia - CAEN/UFC cel: 85-87576600 lattes: http://lattes.cnpq.br/6480130040427049<http://www.orkut.com.br/Interstitial?u=http://lattes.cnpq.br/64801300404270…>

3 2

Re: [R-br] Intervalo de Confiança da Predição
by Danilo Scorzoni Ré 03 Abr '12

03 Abr '12

Perfeito! Deu certo, mto obrigado Walmes! Danilo 2012/4/3 <r-br-request(a)listas.c3sl.ufpr.br> > Enviar submissões para a lista de discussão R-br para > r-br(a)listas.c3sl.ufpr.br > > Para se cadastrar ou descadastrar via WWW, visite o endereço > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > ou, via email, envie uma mensagem com a palavra 'help' no assunto ou > corpo da mensagem para > r-br-request(a)listas.c3sl.ufpr.br > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)listas.c3sl.ufpr.br > > Quando responder, por favor edite sua linha Assunto assim ela será > mais específica que "Re: Contents of R-br digest..." > > > Tópicos de Hoje: > > 1. Matriz de proximidade (Gilbert Queiroz) > 2. Re: Matriz de proximidade (Paulo J Ribeiro Jr) > 3. Matriz de proximidade (Gilbert Queiroz) > 4. Intervalo de Confiança da Predição (Danilo Scorzoni Ré) > 5. Re: Intervalo de Confiança da Predição (Walmes Zeviani) > 6. Re: Problemas com fixed breaks (Mauro Sznelwar) > 7. Re: Tobit espacial (Wagner bonat) > 8. Re: Matriz de proximidade (Wagner bonat) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Mon, 2 Apr 2012 12:07:08 -0700 (PDT) > From: Gilbert Queiroz <gilbert_queiroz(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Matriz de proximidade > Message-ID: > <1333393628.70157.YahooMailNeo(a)web160603.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Amigos, > É possível gerar a matriz de proximidade no R, sem ter que importar de um > outro software, como o GeoDA ou ArcGis? > Quais os comandos? > Abs. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120402/20f4ab01/atta… > > > > ------------------------------ > > Message: 2 > Date: Mon, 02 Apr 2012 16:26:47 -0300 > From: Paulo J Ribeiro Jr <paulojus(a)leg.ufpr.br> > To: r-br(a)listas.c3sl.ufpr.br, Gilbert Queiroz > <gilbert_queiroz(a)yahoo.com.br> > Subject: Re: [R-br] Matriz de proximidade > Message-ID: <1333394807.16481.2.camel@atikum> > Content-Type: text/plain; charset="UTF-8" > > de uma olhada no pacote spdep e em particular em > > poly2nb() > nb2listw() > > > > Em Seg, 2012-04-02 às 12:07 -0700, Gilbert Queiroz escreveu: > > Amigos, > > É possível gerar a matriz de proximidade no R, sem ter que importar de > > um outro software, como o GeoDA ou ArcGis? > > Quais os comandos? > > Abs. > > > > _______________________________________________ > > R-br mailing list > > R-br(a)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. > > -- > Paulo Justiniano Ribeiro Jr > LEG (Laboratorio de Estatistica e Geoinformacao) > Universidade Federal do Parana > Caixa Postal 19.081 > CEP 81.531-990 > Curitiba, PR - Brasil > Tel: (+55) 41 3361 3573 > VOIP: (+55) (41) (3361 3600) 1053 1066 > Fax: (+55) 41 3361 3141 > e-mail: paulojus AT ufpr br > http://www.leg.ufpr.br/~paulojus > > > > > ------------------------------ > > Message: 3 > Date: Mon, 2 Apr 2012 12:25:52 -0700 (PDT) > From: Gilbert Queiroz <gilbert_queiroz(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Matriz de proximidade > Message-ID: > <1333394752.69772.YahooMailNeo(a)web160604.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Amigos, > É possível gerar a matriz de proximidade no R, sem ter que importar de um > outro software, como o GeoDA ou ArcGis? > Quais os comandos? > Abs. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120402/b9fdfefa/atta… > > > > ------------------------------ > > Message: 4 > Date: Mon, 2 Apr 2012 17:37:35 -0300 > From: Danilo Scorzoni Ré <danilo.scorzoni(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Intervalo de Confiança da Predição > Message-ID: > <CA+bJVrd5TuHhwPPjEYWngLKG7b+VaxFj20MHge4MUrDky0js-w(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Prezados, > > Estou ajustando um conjunto de dados ao modelo não-linear de Gompertz (G ~ > a * exp(b * exp(c * Idade))), onde a, b e c são parâmetros. > Consegui realizar o ajuste pelo comando nls() sem problema algum. O que eu > gostaria de obter agora são os intervalos de confiança da predição. Sei que > nos modelos lineares, eu consigo obter através da função predict(), > passando a informação para o parâmetro interval="prediction", no entanto, > ao que me parece, essa opção ainda não está implementada para o comando > nls(), pois quando eu faço o comando > > ic = predict(model2, newdata=g, interval="prediction", level=0.95) > > Ele me retorna apenas a resposta média, e não os intervalos de confiança da > predição. Testei a mesma sequência de comandos com o lm() e funcionou > corretamente. > Alguém sabe se existe outra maneira de estimar esses intervalos? > > Segue o comando mínimo reproduzível: > > > -------------------------------------------------------------------------------------------------------------------------- > > # Código Mínimo Reproduzível > > x = seq(0,15,0.1) > erros = rnorm(151, 0, 3) > y = 22 * exp(-3 * exp(-0.3 * x)) + erros > plot(y ~ x) # Comportamento não-linear > > # Predição e Intervalos de Confiança no Modelo Linear > model.linear = lm(y ~ x + I(x^2)) > summary(model.linear) > xest = as.data.frame(seq(0,15,1)) > colnames(xest) = "x" > pred = as.data.frame(predict(model.linear, newdata=xest, > interval="prediction")) > lines(pred$fit ~ xest$x, lwd=2) > lines(pred$lwr ~ xest$x, lty=2, col=2) > lines(pred$upr ~ xest$x, lty=2, col=2) > > # Predição e Intervalos de Confiança no Modelo Não-Linear > model.nlinear = nls(y ~ a * exp(b * exp(c * x)), start=list(a=25, b=-4, > c=-0.3)) > summary(model1) > > plot(y ~ x) > xest = as.data.frame(seq(0,15,1)) > colnames(xest) = "x" > pred = predict(model.nlinear, newdata=xest, interval="prediction") > str(pred) > lines(pred ~ xest$x, lwd=2) > > -------------------------------------------------------------------------------------------------------------------------- > > Desde já agradeço o empenho! > > Abraços, > > -- > Danilo Scorzoni Ré > Engenheiro Florestal > Mestre em Ciência Florestal > (14) 8180-2494 > http://about.me/dscorzoni > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120402/98c19b32/atta… > > > > ------------------------------ > > Message: 5 > Date: Mon, 2 Apr 2012 19:07:26 -0300 > From: Walmes Zeviani <walmeszeviani(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Intervalo de Confiança da Predição > Message-ID: > <CAFU=EkZ+3nrxZjTuOmXLbkzsKuQJXxyzBXiR9Oc6GB6C=8r2Mw(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Veja se isso pode ser útil > > > http://ridiculas.wordpress.com/2011/05/19/bandas-de-confianca-para-modelo-d… > > À disposição. > Walmes. > > ========================================================================== > Walmes Marques Zeviani > LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) > Departamento de Estatística - Universidade Federal do Paraná > fone: (+55) 41 3361 3573 > VoIP: (3361 3600) 1053 1173 > e-mail: walmes(a)ufpr.br > twitter: @walmeszeviani > homepage: http://www.leg.ufpr.br/~walmes > linux user number: 531218 > ========================================================================== > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120402/20f5e4cd/atta… > > > > ------------------------------ > > Message: 6 > Date: Mon, 2 Apr 2012 22:08:08 -0300 > From: "Mauro Sznelwar" <sznelwar(a)uol.com.br> > To: <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Problemas com fixed breaks > Message-ID: <C90D1BEEB10647F8A610AB827EBF250D@ACER4520> > Content-Type: text/plain; charset="iso-8859-1" > > > Tentei rodar de novo, mas deu dois problemas: > > library(maptools) > Carregando pacotes exigidos: foreign > Carregando pacotes exigidos: lattice > Checking rgeos availability: FALSE > Note: when rgeos is not available, polygon geometry > computations in maptools depend on gpclib, > which has a restricted licence. It is disabled by default; > to enable gpclib, type gpclibPermit() > > e depois > > > mapa = readShapePoly("saudemapa.shp") > Erro em getinfo.shape(filen) : Error opening SHP file > > E ainda faltou disponibilizar o arquivo > logo_pb.png > > > Em 1 de abril de 2012 21:05, Mauro Sznelwar <sznelwar(a)uol.com.br> > escreveuo: > > Estava tentando rodar este código, e vi que falta este > arquivoc:/users/Julianna/Documents/SecSaude/saudemapa.shp > Daria para enviar?. > > Mensagem abaixo em nome da Julianna... (dados disponiveis via > datafilehost) > > > ==== > > > E quando fui conferir os breaks, os intervalos estavam certos,mas o > mapa resultante estava errado. > Porque quando eu não importo o csv e uso o dataframe resultante da > função, os breaks ficam ok(da maneira que eu especifico em FixedBreaks) e, > do contrario, não??? > > > segue o código: > > > ####################################################### > # CONSTRUINDO MAPAS > ####################################################### > ##Code (Comments are preceded by ##) > ## Load required packages > > > # Limpando memoria. > rm(list=ls(all=TRUE)) > > > library(maps) > library(sp) > library(maptools) > library(spdep) > library(classInt) > library(RColorBrewer) > library(shape) > library(SDMTools) > library(latticeExtra) > library(ReadImages) > library(png) > > > rosa_dos_ventos = function(loc=c(-34.5, -8), size = 0.15, bearing=0, > cols=1, cex=0.6,...) > { > cols <- rep(c("white","black"),8) > radii <- rep(size/c(1,4,2,4),4) > x <- radii[(0:15)+1]*cos((0:15)*pi/8+bearing)+loc[1] > y <- radii[(0:15)+1]*sin((0:15)*pi/8+bearing)+loc[2] > for (i in 1:15) > { > x1 <- c(x[i],x[i+1],loc[1]) > y1 <- c(y[i],y[i+1],loc[2]) > polygon(x1,y1,col=cols[i]) > } > polygon(c(x[16],x[1],loc[1]),c(y[16],y[1],loc[2]),col=cols[16]) > b <- c("O","S","L","N") > for (i in 0:3) > text((size+par("cxy")[1]-0.40)*cos(bearing+i*pi/2)+loc[1], > > (size+par("cxy")[2]-0.42)*sin(bearing+i*pi/2)+loc[2],b[i+1],cex=cex) > } > > > ## Lendo dados > ## http://www.datafilehost.com/download-e5cb7e86.html > dados = read.table("dadosJulianna.csv", sep = ",", header = TRUE) > > > # Lendo mapa. > mapa = > readShapePoly("c:/users/Julianna/Documents/SecSaude/saudemapa.shp") > proj4string(mapa) <- CRS("+init=epsg:4291") > scaleXY(mapa, 1836656) > > > > > # arquivo png > #img <- readPNG("c:/users/Julianna/Documents/SecSaude/logo_pb.png") # > seu arquivo png aqui > #set.seed(123) > #par(xpd = TRUE) > > > > > colors <- brewer.pal(7, "Reds") > brks <- classIntervals(dados$ano_2006, n = 7,style="fixed", > intervalClosure="right", > fixedBreaks=c(-2, -0.5, 0, 20, 40, 60, 80, > 10000), main="Fixed") > brks<- brks$brks > > > ID = > match(substr(mapa$CODIGO,1,6),dados$t_indicadores_municipio.f_municipio) > dados = dados[ID,] > > > png(filename="c:/mif2006.png", height=750, width=1500) > > > > > plot(mapa, col = colors[findInterval(dados$ano_2006, > brks,all.inside=TRUE)], axes = TRUE, cex = 1.5, mar=5) > ##, xlim=c(-39.7, -34.7), ylim=c(-8.5, -6) > > > ##invisible(text(getSpPPolygonsLabptSlots(mapa), > labels=as.character(mapa$NOMEMUN_1), cex=0.65)) > > > title(main = "PORCENTAGEM DE Ã BITOS INVESTIGADOS\ntipo: mulheres em > idade fÃ©rtil", > font.main = 2, cex.main = 1.8) > > > text(-35.7, -6.2, "ANO: 2006", cex = 1.8) > > > legend(-39.6,-7.6, c("Sem informaÃ§Ã£o", "0", "]0,20]", "]20,40]", > "]40,60]", "]60,80]","Acima de 80"), > fill = colors, bty="n", x.intersp = 1.1, y.intersp = 1.1, cex > = 1.6) > > > rosa_dos_ventos(loc=c(-34.5,-8),size=0.17,cex=1.20, bearing=0, cols=1) > > > SpatialPolygonsRescale(layout.scale.bar(), offset = c(-35,-8.30), > scale = 1, fill = c("transparent", "black"), > plot.grid= FALSE) > > > text(-34.5, -8.35, "110.3 KM", cex= 1.3) > > > #rasterImage(img, xleft = -39.5, ybottom = -6.25, xright = -38.9 , > ytop = -6) > > > dev.off() > > > > > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120402/e7f0095e/atta… > > > > ------------------------------ > > Message: 7 > Date: Tue, 3 Apr 2012 11:42:26 -0300 > From: Wagner bonat <wbonat(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Tobit espacial > Message-ID: > <CANt=4MietnxFQRVKX2FWTXnvNFVRt4NXrjtnBMLkPFMDFKb4=w(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Que tipo de dependência espacial vc quer modelar ? Onde encontrou este > modelo ? Pode não ter programado, mas com certeza pode ser programado. > > > > -- > Wagner Hugo Bonat > LEG - Laboratório de Estatística e Geoinformação > UFPR - Universidade Federal do Paraná > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120403/2b5e44e9/atta… > > > > ------------------------------ > > Message: 8 > Date: Tue, 3 Apr 2012 11:43:55 -0300 > From: Wagner bonat <wbonat(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br, Gilbert Queiroz > <gilbert_queiroz(a)yahoo.com.br> > Subject: Re: [R-br] Matriz de proximidade > Message-ID: > <CANt=4MhG4Gug6tLEJLX-7PKb+_VP0SsQsYDi4B6iP_SjPK0JiA(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Sim, o pacote spdep tem várias funções para montar matrizes de vizinhança > seguinte diveros critérios e de acordo com o tipo de dados que vc lê. > -- > Wagner Hugo Bonat > LEG - Laboratório de Estatística e Geoinformação > UFPR - Universidade Federal do Paraná > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120403/94bcdfd9/atta… > > > > ------------------------------ > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > Fim da Digest R-br, volume 14, assunto 3 > **************************************** > -- Danilo Scorzoni Ré Engenheiro Florestal Mestre em Ciência Florestal (14) 8180-2494 http://about.me/dscorzoni

1 0

Matriz de proximidade
by Gilbert Queiroz 03 Abr '12

03 Abr '12

Amigos, É possível gerar a matriz de proximidade no R, sem ter que importar de um outro software, como o GeoDA ou ArcGis? Quais os comandos? Abs.

2 1

Re: [R-br] Problemas com fixed breaks
by Julianna Trindade 03 Abr '12

03 Abr '12

Obs: as faixas de valores com "n" casas decimais foram colocadas em função do erro que mostrei no email anterior, no entanto, ao selecionar o dataframe brks, vcs verão que 2 valores das classes de intervalos vem com erro. 2012/3/30 Julianna Trindade <julianna(a)jubalitpb.com> > arquivos em anexo > > > > Em 29 de março de 2012 21:33, Benilton Carvalho < > beniltoncarvalho(a)gmail.com> escreveu: > > Julianna, >> >> vc precisa dar-nos algo que possamos reproduzir para que possamos ajudar >> (esse e' o motivo de vc nao ter obtido resposta ao enviar esta mesma msg >> dois dias atras). >> >> Assumindo que 'dados' nao tenha muito mais que 30-50 linhas, seria muito >> util que vc postasse o resultado de: >> >> dput(dados$ano_2006) >> >> e tbm nos dissesse qual eh o pacote que te da' a funcao classIntervals(), >> afinal nao e' algo do base. >> >> benilton >> >> 2012/3/29 Julianna Trindade <julianna(a)jubalitpb.com> >> >>> Boa Tarde, >>> >>> Gostaria de uma solução para o FixedBreaks obedecer os limites >>> superiores e inferiores das minhas classes de intervalos. No entanto, ele >>> deve associar algum método estatístico e burla minha configuração. >>> >>> Eis o código que eu utilizo: >>> >>> brks <- classIntervals(dados$ano_2006, n = 7,style="fixed", >>> intervalClosure="right", fixedBreaks=c(-2, -0.5, 0, 20, 40, 60, 80, 10000), >>> main="Fixed") >>> brks<- brks$brks >>> >>> Eis o breaks gerado pelo R: >>> >>> c(-11.1, -0.5, 16.665, 41.665, 58.335, 83.335, 110.1) >>> >>> Alguém pode me ajudar, por favor? >>> >>> >>> >>> >>> -- >>> Julianna Trindade >>> >>> "Exalta-te, Senhor, na Tua força! Nós cantaremos e louvaremos o Teu >>> poder." >>> >>> Salmos 21; 13 >>> >>> _______________________________________________ >>> R-br mailing list >>> R-br(a)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(a)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. >> > > > > -- > Julianna Trindade > > "Exalta-te, Senhor, na Tua força! Nós cantaremos e louvaremos o Teu poder." > > Salmos 21; 13 > -- Julianna Trindade "Exalta-te, Senhor, na Tua força! Nós cantaremos e louvaremos o Teu poder." Salmos 21; 13

4 10

Intervalo de Confiança da Predição
by Danilo Scorzoni Ré 02 Abr '12

02 Abr '12

Prezados, Estou ajustando um conjunto de dados ao modelo não-linear de Gompertz (G ~ a * exp(b * exp(c * Idade))), onde a, b e c são parâmetros. Consegui realizar o ajuste pelo comando nls() sem problema algum. O que eu gostaria de obter agora são os intervalos de confiança da predição. Sei que nos modelos lineares, eu consigo obter através da função predict(), passando a informação para o parâmetro interval="prediction", no entanto, ao que me parece, essa opção ainda não está implementada para o comando nls(), pois quando eu faço o comando ic = predict(model2, newdata=g, interval="prediction", level=0.95) Ele me retorna apenas a resposta média, e não os intervalos de confiança da predição. Testei a mesma sequência de comandos com o lm() e funcionou corretamente. Alguém sabe se existe outra maneira de estimar esses intervalos? Segue o comando mínimo reproduzível: -------------------------------------------------------------------------------------------------------------------------- # Código Mínimo Reproduzível x = seq(0,15,0.1) erros = rnorm(151, 0, 3) y = 22 * exp(-3 * exp(-0.3 * x)) + erros plot(y ~ x) # Comportamento não-linear # Predição e Intervalos de Confiança no Modelo Linear model.linear = lm(y ~ x + I(x^2)) summary(model.linear) xest = as.data.frame(seq(0,15,1)) colnames(xest) = "x" pred = as.data.frame(predict(model.linear, newdata=xest, interval="prediction")) lines(pred$fit ~ xest$x, lwd=2) lines(pred$lwr ~ xest$x, lty=2, col=2) lines(pred$upr ~ xest$x, lty=2, col=2) # Predição e Intervalos de Confiança no Modelo Não-Linear model.nlinear = nls(y ~ a * exp(b * exp(c * x)), start=list(a=25, b=-4, c=-0.3)) summary(model1) plot(y ~ x) xest = as.data.frame(seq(0,15,1)) colnames(xest) = "x" pred = predict(model.nlinear, newdata=xest, interval="prediction") str(pred) lines(pred ~ xest$x, lwd=2) -------------------------------------------------------------------------------------------------------------------------- Desde já agradeço o empenho! Abraços, -- Danilo Scorzoni Ré Engenheiro Florestal Mestre em Ciência Florestal (14) 8180-2494 http://about.me/dscorzoni

2 1

Matriz de proximidade
by Gilbert Queiroz 02 Abr '12

02 Abr '12

Amigos, É possível gerar a matriz de proximidade no R, sem ter que importar de um outro software, como o GeoDA ou ArcGis? Quais os comandos? Abs.

2 1