- R-br - DINF UFPR

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 03 Abr '12

03 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 02 Abr '12

02 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

[dúvida] colocar única legenda em 4 plots no mesmo frame
by Marcelo Claro de Souza 31 Mar '12

31 Mar '12

Plotei quatro figuras em um mesmo frame utilizando: layout(matrix(1:4,nrow=2,ncol=2)) plot(x1~y1,data=dados) plot(x2~y2,data=dados) plot(x3~y3,data=dados) plot(x4~y4,data=dados) Eu gostaria de colocar uma única legenda nos eixos x e y para representar os eixos das 4 figuras ao mesmo tempo, sem a necessidade de chamar uma legenda para cada figura evitando redundâncias. Alguém teria alguma sugestão? Muito obrigado. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

3 2

Somatório com matrizes
by Guilherme Biz 31 Mar '12

31 Mar '12

Bom dia Senhores, exite alguma função pronta para fazer somatório de matrizes? Estou tentando somar, dentro de um loop, muitas matrizes e não encontrei uma função para este somatório. Vou ter que fazer outro loop para estas somas? Ex Tenho as matrizes A1,......A100. e quero fazer A1+A2+....+A100. Att, Guilherme Biz

3 3

Re: [R-br] Tobit espacial
by Diego andre 31 Mar '12

31 Mar '12

Obrigado Leonard de Assis pela resposta, mas infelizmente não tem *tobit espacial* neste pacote... -- 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…>

2 1