Novembro 2011 - R-br

Gráficos no R ou Gnu Plot
by Pedro Rafael 29 Nov '11

29 Nov '11

Pessoal ultimamente venho precisando criar gráficos com bastantes "frescurinhas", com perdão da palavra. Eu pessoalmente não gosto de gráficos assim, prefiro os mais simples e informativos, contudo no lugar que eu trabalho o pessoal gostam dessas coisas porque precisa apresentar trabalhos para pessoas legas em estatística. Em geral o pessoal gosta dos gráficos simples do excel que eu me recuso a usar no Linux rodando via Wine. Qual seria uma boa biblioteca no R para construção de gráficos de linhas, setores em 3D? Estou meio frustado de ter que fazer essas coisas podendo me preocupar com coisas mais importantes mas tudo bem. rsrsrs Tava olhando também o GNU Plot. Ele é mais maleável que o R para construção de gráficos? Vale a pena estudar sua sintaxe ou melhor fazer tudo no R? -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

11 25

Dúvida
by Edson Lira 29 Nov '11

29 Nov '11

Caros R-istas, boa noite. Um amigo me pediu auxílio para fazer a seguinte operação: Seja um vetor: pegar o primeiro elemento e subtrair pelo segundo, depois pegar o primeiro e subtrair pelo terceiro, e assim sucessivamente, sendo que depois eu recomeço fazendo o segundo menos o terceiro, depois o segundo menos o quarto e assim até acabar todos os elementos do vetor. Como fazer isso??? Edson Lira Estatístico Manaus-Amazonas

5 5

Tratar código HTML no R
by Paulo Nogueira 29 Nov '11

29 Nov '11

Saudações a todos! Preciso aprender como tratar HTML no R para identificar alguns padrões na estrutura do código, extrair e guardar determinadas informações. Alguém costuma fazer algo do tipo e poderia informar os pacotes mais indicados para esse trabalho? Terei de acessar cerca de 300 páginas semanalmente, e extrais alguns parâmetros. O HTML abaixo é um trecho do código que terei de lidar, e preciso extrair as informações destacadas. Alguém pode me sugerir o início do caminho das pedras, como algum pacote e/ou funções específicas? <div id="individual-completo"> <div class="title-box"> <img src="http://address.jpg" gestor="foto" class="foto" /> <div> <div class="colright"> <span>Participações: 22</span> <a href="http://players" class="visite">Saiba Mais</a> </div> <h2>NomeElemento</h2> <p class="tflag"><img src="http://team" /> <a href="/refteam"> NomeGrupo</a></p> <p>NomeFunção</p> </div> </div> <div class="content-box"> <h2>Desempenho completo</h2> <table cellspacing="0" cellpadding="0"> <tbody> <tr> <th>Parâmetro1</th> <td><span gestor="param1">2</span></td> </tr> <tr class="divisor"> <td colspan="2"> </td> </tr> <tr> <th>Parâmetro2</th> <td><span gestor="param2">63</span></td> </tr> <tr class="odd"> <th>Parâmetro3</th> <td><span gestor="param3">87</span></td> </tr> <tr class="divisor"> <td colspan="2"> </td> </tr> <tr class="destaque"> </tbody> </table> </div> </div> Obrigado

3 3

Como ordenar com cruzamentos de variáveis.
by Edson Lira 28 Nov '11

28 Nov '11

Caros amigos, estou analisando um banco de dados de doações, usando o comando: rev(sort(table(tip_doa$BAIRRORESIDENCIA)))[1:20] Tenho os 20 bairros com maior número de doações de janeiro a outubro de 2011. Gostaria de cruzar o bairro com a aptidão na doação e ordenar pelos percentuais de aptidão. Como faço isto? table(tip_doa$BAIRRORESIDENCIA,tip_doa$APTO) Edson Lira Estatístico Manaus-Amazonas

2 1

correlação de Spearman e respectiva linha no grafico
by Pedro Emmanuel Alvarenga Americano do Brasil 28 Nov '11

28 Nov '11

Amigos de R, Andei procurando RSiteSearch funçoes sobre a correlação de spearman. Eu achei diversas funçoes que fazem o rho de spearman e algumas extensões e correções, mas não consegui achar nenhuma que colocasse a linha da correlação em um grafico, a semelhança da reta na correlação de pearson. Alguem conhece alguma função ou tem algum script que faz isso? v1 v2 1 26 0.000 2 73 0.000 3 76 143.072 4 72 172.587 5 59 123.325 6 78 142.367 7 73 190.804 8 77 14.497 9 63 120.499 10 72 165.570 11 73 152.220 12 75 170.482 13 73 137.436 14 71 169.079 15 72 170.482 16 65 160.654 17 76 0.000 18 28 0.000 19 61 107.763 20 66 8.643 21 68 133.911 22 65 113.427 23 27 0.000 24 71 104.929 25 73 116.257 26 68 0.000 27 75 142.367 28 63 0.000 29 77 7.168 30 66 127.561 31 66 0.000 32 60 144.480 33 74 136.731 34 65 32.527 35 74 145.184 36 65 112.012 37 67 145.184 38 81 103.511 39 82 3.448 40 72 130.384 41 66 244.553 42 73 0.000 43 65 0.000 44 57 0.000 45 63 105.637 46 60 98.546 47 76 137.436 48 68 1.558 49 26 0.000 50 72 0.000 51 64 30.378 52 76 0.000 53 25 3.072 54 58 0.000 55 28 0.000 56 69 72.592 57 26 0.000 58 27 2.318 59 26 0.790 60 63 0.000 61 71 121.912 62 60 41.808 63 33 0.000 64 31 0.000 65 71 0.000 66 24 0.000 67 69 0.000 68 69 22.468 69 35 0.000 70 65 33.243 71 78 7.906 72 31 0.000 73 74 0.000 74 71 4.943 75 65 38.244 76 34 3.448 77 19 11.944 78 35 19.577 79 26 0.000 80 35 197.415 81 13 0.000 82 35 0.000 83 35 0.000 84 28 0.000 85 17 30.378 86 37 0.000 library('Hmisc') spearman2(v1 ~ v2, data=z, p=2) Spearman rho^2 Response variable:v1 rho2 F df1 df2 P Adjusted rho2 n v2 0.181 9.16 2 83 0.0003 0.161 86 Abraço forte e que a força esteja com você, Dr. Pedro Emmanuel A. A. do Brasil Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365, CEP 21040-360, Tel 55 21 3865-9648 email: pedro.brasil(a)ipec.fiocruz.br email: emmanuel.brasil(a)gmail.com ---Apoio aos softwares livres www.zotero.org - gerenciamento de referências bibliográficas. www.broffice.org ou www.libreoffice.org - textos, planilhas ou apresentações. www.epidata.dk - entrada de dados. www.r-project.org - análise de dados. www.ubuntu.com - sistema operacional

4 4

Inversão de matriz diagonal gigante
by Marcelo Cardoso mello 28 Nov '11

28 Nov '11

Olá pessoal. Gostaria de inverter uma matriz diagonal com algum processo que use menos memória que o do exemplo abaixo. (tenho uma matriz enorme e não esta invertendo via solve()). d<-c(1,1,0.5,0.75,0.5,0.5) d<-diag(d) i<-solve(d) i diag(i) # resultado que eu quero, mas de uma forma a gastar menos memória para diagonais gigantes!!!!

3 2

Versão R-2.14.0 no Ubunto
by Marcelo Cardoso mello 27 Nov '11

27 Nov '11

Pessoal, executei no terminal: $ sudo apt-get install r-base $ sudo apt-get install r-recommendad $ sudo apt-get etc etc que me foi indicado por aqui Tudo ok. Mas o R instalado é da versão 2.10 Tudo bem até aí. O problema é que não instala pacotes recentes. usei: instal.packages("vegan") tudo ok ... instalou mas não instala o lme4 por exemplo. O problema deve ser da versão do R que tenho? Como faço para instalar a versão R-2.14.0 no Ubunto? De forma rápida e simples. Obrigado por qualquer dica.

6 12

Re: [R-br] valor maximo e minimo
by Ivan Bezerra Allaman 27 Nov '11

27 Nov '11

Gilson, Busque por modelos que podem lhe dar respostas mais interessantes. Veja, A.Det=c(-0.595750,0.050250,0.270000,0.711500,1.608750,3.353500,4.937652,4.573750,4.955000) luz=c(0,25,50,75,100,500,1000,1500,2000) ajuste <- nls((A.Det~(a + 0.595750) * (1 - exp (-b * luz/(a-0.595750)))+0.595750), start=list(a=500, b=5)) plot(luz,A.Det) curve(3.9669 + 0.595750 * (1 - exp (-0.0061 * x/(3.9669-0.595750)))+0.595750,add=TRUE,col='blue') curve(SSasympOff(x,4.8292,-5.7353,23.452),add=TRUE,col='red') legend('bottomright', xjust=1,bty="n", c('Alternativo','Gilson'), lwd=2,col=c('red','blue')) Sobre este modelo 'alternativo', consulte "Mixed Effects Models in S and S-Plus, Pinheiro & Bates (2000)" OBS: Estes modelos não tem ponto de inflexão, pois quanto x tende ao infinito, y tende a zero, veja: D(expression(theta1*(1-exp(-exp(theta2)*(x-theta3)))),'x') curve(4.8292 * (exp(-exp(-5.7353) * (x - 23.452)) * exp(-5.7353)),to=0,from=2000,ylab='y') (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

1 0

valor maximo e minimo
by Gilson Sanchez 26 Nov '11

26 Nov '11

Pessoal mais uma vez solicitando ajuda de vc's O seguinte estou fazendo ajuste de curva com uma exponencial, então quero saber se da pra saber o valor do ponto maximo e minimo dos valores ajustados tanto no eixo X quanto no eixo Y? espe1 A.Det=c(-0.595750,0.050250,0.270000,0.711500,1.608750,3.353500,4.937652,4.573750,4.955000) luz=c(0,25,50,75,100,500,1000,1500,2000) ajuste = nls((A.Det~(a + 0.595750) * (1 - exp (-b * luz/(a-0.595750)))+ 0.595750), start=list(a=500, b=5)) summary(ajuste) obrigado! -- Bach. Agr.Gilson Sánchez Chia, MSc. Laboratório de Fisiologia Vegetal Embrapa Amazônia Ocidental Fone: (92) 3303-7841 http://lattes.cnpq.br/3684553005529634 --- Antes de imprimir, pense na sua responsabilidade com o Meio Ambiente. Please consider the environment before printing this email. Renew, Reduce and Recycle. The Planet Knows.*** Salva un Árbol.* * *** * PIENSA**** EN EL MUNDO PARA TUS HIJOS** *

2 3

Distribuição densidade de probabilidade de Weibull
by Rodrigo Muniz 25 Nov '11

25 Nov '11

Olá Pessoal, Sou usuário recente do R e como tal ainda tenho muitas dúvidas básicas na utilização do programa. Tenho uma série de dados de velocidade de vento na qual quero ajustar os parâmetros c e K da distribuição de Wuibull, no entanto quando uso as funções abaixo, da um erro que x deve ser um vetor numérico, no entanto se eu colocar um valor numérico ele gera uma série de números pseudoaleatórios, O que eu quero abrir os meus dados dentro da função. Se alguém puder me ajudar eu agradeço. Vento50m <- dweibull(Vento, shape=2, scale = 1) fitdistr(Vento50m, densfun="weibull") ATT. -- *Rodrigo A. Muniz* Eng. Agrônomo. Ms Produção Vegetal (UENF) Doutorando em Engenharia de Sistemas Agrícolas (ESALQ/USP) Cel (19) 8300-4333 (Pessoal) Tel (19) 3375-1196 (Residêncial) Tel (19) 3447-8537 (Trabalho)

4 6