Abril 2011 - R-br - DINF UFPR

Interação significativa em regressão logística
by Marcelo Cardoso mello 18 Abr '11

18 Abr '11

Prezados, tenho um modelo logístico do tipo: Prenhez (+/_) -> (Método de inseminação (0,1); Época do ano (0,1); interação) A interação foi significativa. Então, avaliei os métodos em cada época, ou seja, fiz dois modelos mais simples, um para cada época: Prenhez (+/_) -> (Método de inseminação (0,1)) (considerando a época 1 e teve efeito do método) Prenhez (+/_) -> (Método de inseminação (0,1)) (considerando a época 0 e nãp teve efeito do método) O procedimento é esse mesmo? Fiz via glm. Tem uma maneira diferente de fazer essa avaliação? Tenho mesmo que partir para modelos separados dentro de cada época? Agradeço a opinião dos colegas. Marcelo

3 7

Arvore Geradora Minima
by Emerson Cotta Bodevan 18 Abr '11

18 Abr '11

Prezados, existe algum pacote no R para construção de árvores geradoras mínimas? Abs, Emerson.

2 2

Fwd: [Dúvida] Limpando memória RAM
by Pedro Rafael 18 Abr '11

18 Abr '11

Realmente Jakson apesar de tá usando o WIndows que não tem um gerenciamento de memória que preste eu fiz o teste e não vi grandes melhoras. Consegui otimizar reescrevendo o programa de uma forma mais eficiente. Pedro ---------- Mensagem encaminhada ---------- De: Pedro Rafael <pedro.rafael.marinho(a)gmail.com> Data: 18 de abril de 2011 11:01 Assunto: Re: [R-br] [Dúvida] Limpando memória RAM Para: "Jakson A. Aquino [via R-br]" < ml-node+3454967-637477867-223914(a)n4.nabble.com> Realmente Jakson apesar de tá usando o WIndows que não tem um gerenciamento de memória que preste eu fiz o teste e não vi grandes melhoras. Consegui otimizar reescrevendo o programa de uma forma mais eficiente. Pedro Em 17 de abril de 2011 01:13, Jakson A. Aquino [via R-br] < ml-node+3454967-637477867-223914(a)n4.nabble.com> escreveu: On Sat, Apr 16, 2011 at 03:03:49PM -0300, Henrique Dallazuanna wrote: > > remova o objeto com o rm e após execute gc(reset = TRUE) > > Eu não conhecia a opção reset da função gc(), mas num teste que fiz > aqui ela não parece ter feito muita diferença. Carreguei o banco de > dados do World Values Survey com read.dta (o arquivo está disponível > em três formatos, sendo um deles o do STATA). O comando top do Linux > mostrou que o R estava ocupando em torno de 2,7 GB da RAM. > > Depois dos comandos > > rm(list=ls()) > gc(reset = TRUE) > > a memória ocupada continuou um pouco acima de 1 GB. > > Existe alguma forma mais eficiente do que gc() de mandar o R > reorganizar (desfragmentar?) a memória? > > Atualmente, tenho feito o seguinte em situações em que preciso dessa > reorganização da memória: > > 1) Salvo a área de trabalho: save.image() > 2) Fecho o R: quit(save = no) > 3) Abro o R novamente e carrego a área de trabalho salva. > > [É claro que poderia ter usado quit(save = "yes"), mas prefiro > escolher um nome para o RData.] > > -- > Jakson > _______________________________________________ > R-br mailing list > [hidden email]<http://user/SendEmail.jtp?type=node&node=3454967&i=0&by-user=t> > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > ------------------------------ > If you reply to this email, your message will be added to the discussion > below: > > http://r-br.2285057.n4.nabble.com/R-br-Duvida-Limpando-memoria-RAM-tp345434… > To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib…>. > > -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6 -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

1 0

Movimentar título nos eixos x y de um gráfico
by Marcelo Claro de Souza 18 Abr '11

18 Abr '11

Além de utilizar xlab="nome da variável" e ylab="nome da variável" para determinar o título dos eixos, como faço para movimentar os títulos para mais perto ou mais longe da figura? Muito obrigado.

3 3

Re: [R-br] Movimentar título nos eixos x y de um gráfico
by Marcelo Claro de Souza 18 Abr '11

18 Abr '11

José, estou utilizando o script abaixo: cor=as.data.frame(cor) simbolos=c(19,19,19,19,19,19,19,2,2,2,2,2,2,2,3,3,3,3,3,3,3) N.P.cor=lm(N~P,data=cor) require(graphics) X11() jpeg('NPcor.jpeg') plot(N~P,data=cor,pch=simbolos,ylab=expression(N~(g~Kg^{-1})),xlab=expression(P~(g~Kg^{-1}))xlim=c(0,1.6),ylim=c(0,35)) abline(N.P.cor) legend("topleft",c("Ep","Es","Et"),pch=c(19,2,3)) dev.off() Eu gostaria de mudar o título do eixo y p mais próximo da figura (p o lado direito) pq os valores sobrescritos estão ficando muito colados na borda da figura enão estão aparecendo direito. Segue uma figura em anexo marcando o problema. Muito obrigado. --- Em seg, 18/4/11, José Cláudio Faria <joseclaudio.faria(a)gmail.com> escreveu: De: José Cláudio Faria <joseclaudio.faria(a)gmail.com> Assunto: Re: [R-br] Movimentar título nos eixos x y de um gráfico Para: r-br(a)listas.c3sl.ufpr.br Data: Segunda-feira, 18 de Abril de 2011, 8:44 Marcelo, Isso irá depender muito do ambiete gráfico que vc está usando (graphics, grig, rgl) e do pacote/função (graphics, lattice, ggplot2). Se você mandar um script de exemplo a ajuda poderá ser facilitada e mais certeira. Abs, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - prof. Titular UESC/DCET/Brasil joseclaudio.faria(a)gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

2 1

[Dúvida] Limpando memória RAM
by Pedro Rafael 18 Abr '11

18 Abr '11

Pessoal eu tenho duas base de dados de 600 MB totalizando 1200 MB e carro ela no R mas não necessito da base toda carregada na memória RAM. As variáveis que me interessam eu jogo em um vetor e com o comando rm() eu removo os dois data.frames. Esse comando retira os objetos completamente da memória RAM ou deixa lixo? Existe uma forma mais eficiente? -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

7 10

Link para grupo no Site do R
by Augusto Ribas 18 Abr '11

18 Abr '11

Ola pessoal. Estava navegando pelo site do R (http://www.r-project.org/) e no menu lateral tem um link uma documantação com os user groups do mundo. Ai ja tinha o Brasil la mas cadastrada a lista R-Stat. Eu tomei a liberdade de atualizar para a lista r-br e coloquei o Leg como responsavel la. Espero que esteja tudo bem, se alguem quiser dar uma conferida la: http://rwiki.sciviews.org/doku.php?id=rugs:r_user_groups Assim, quem vier do site oficial do R ja vai cair nessa lista e nao no R-stat. Julguei nao necessario deixar o link para a lista R-stat la ja que ela vai ser desativa. Bem é isso Abraços -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org/ http://dhost.info/augusto/ Lattes: http://lattes.cnpq.br/7355685961127056

1 0

Re: [R-br] Superficie de resposta publicavel no R?!
by José Cláudio Faria 18 Abr '11

18 Abr '11

Pessoal, Parece que o anexo passou nos servidores! Na parte da mensagem aqui para mim aparece assim: -------------- Próxima Parte ---------- Um anexo não-texto foi limpo... Nome: srp.r Tipo: application/x-wine-extension-r Tamanho: 2317 bytes Descrição: não disponível URL:<http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110417/cfe9351f/atta…> Basta fazer o download do link URL acima e mudar o nome de attachement.bin para srp.r. Abs, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - prof. Titular UESC/DCET/Brasil joseclaudio.faria(a)gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

1 0

Superficie de resposta publicavel no R?!
by José Cláudio Faria 18 Abr '11

18 Abr '11

Pessoal, Sempre que penso em fazer uma superfície de resposta publicável no R me dá um frio na barriga... A alguns anos fiz (com a ajuda do J.Fox e D.Murdoch) uma função (usando o rgl) muito flexível (disponível em: http://rwiki.sciviews.org/doku.php?id=graph_gallery:new-graphics&s=plotlm3d <http://rwiki.sciviews.org/doku.php?id=graph_gallery:new-graphics&s=plotlm3d> O problema dela é que labels do eixos não admitem expressões matemáticas, geralmente necessárias para expressar a escala da variável (mm^2, kg h^-1, etc). A limitação não é da função, mas do rgl. Dá para contornar editando o figura gerada ou inserindo bitmaps usando os recursos do rgl, mas ambas alternativas trabalhooosas. Pois bem, estou diante da necessidade de refazer várias superfícies de resposta, feitas no passado (9 anos) com o Statiscs for Windows (nunca mais usei, mas ainda sinto saudades de alguns recursos...). Com um frio na barriga e ralando no fim de semana, decidi tentar simplificar esta tarefa para o presente e futuro. Com a ajuda de alguma coisa que achei nas listas de discussão, a juda do Walmes e do Henrique (vcs não descansam no domingo né? ;-)) rascunhei umas mal traçadas linhas no R, que embora não perfeita, está resolvendo o básico. Decidi socializar, pois é assunto recorrente nas listas de discussão, e não fiz só. Possivelmente seja útil para outros usuários, melhor ainda se receber alguém fizer melhoramentos. No momento estou sem tempo, mas se tiver alguém interessado, podemos pensar em juntar a "plolm3d" e a atual "sr" em um pacote e publicar na CRAN. Penso que seria de muita utilidade! A função "sr", assim como um pequeno script demo estão em anexo (vou testar para ver se a lista recebe, se não for possível envio outra mensagem com o script). Abs, (obrigado novamente so Walmes e Henrique pela disposição de ajudar) -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - prof. Titular UESC/DCET/Brasil joseclaudio.faria(a)gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

1 0

Re: [R-br] Ajuda: duvida passando uma string como modelo para uma função
by Jose Claudio Faria 17 Abr '11

17 Abr '11

Walmes, Embora ache que essa função ainda possa ser melhorada, sua ajuda foi valiosa! Sua sugestão pode ser mais simples, e próximo do que realmente estou precisando, veja: mF <- function(x=seq(1, 5, by=1), y=seq(0, 10, by=5), model){ z <- outer(x, y, function(x, y){ eval(model) }) z } mF(model=quote(2 - 0.4*y - 0.5*x + 0.8*sqrt(x))) > mF(model=quote(2 - 0.4*y - 0.5*x + 0.8*sqrt(x))) [,1] [,2] [,3] [1,] 2.300000 0.3000000 -1.700000 [2,] 2.131371 0.1313708 -1.868629 [3,] 1.885641 -0.1143594 -2.114359 [4,] 1.600000 -0.4000000 -2.400000 [5,] 1.288854 -0.7111456 -2.711146 Talvez na segunda-feira outras pessoas possam dar alguma contribuição... Abs e muito obrigado! -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - prof. Titular UESC/DCET/Brasil joseclaudio.faria(a)gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

3 3