Dezembro 2011 - R-br

Categorização de variáveis
by Edson Lira 14 Dez '11

14 Dez '11

Boa tarde, caros r-istas, estou com alguns probleminhas na categorização de variáveis. Tenho este banco. pad pas sexo pg 155 134 83 feminino 15 156 72 99 feminino 30 157 82 151 feminino 10 158 23 70 feminino 20 159 75 128 masculino 13 160 73 125 feminino 05 161 125 80 masculino 08 162 65 110 feminino 10 163 132 92 masculino 13 164 114 68 masculino 11 165 83 153 feminino 26 166 166 107 masculino 22 167 114 76 feminino 03 Preciso classificar da seguinte forma: se pad<85 e pas<130 é normal se 85<pad<89 e 130<pas<139 Normal limítrofe e assim por diante. a outra classificação é a seguinte: Para homens o percentual de gordura entre 15 e 20 é normal entre 20 e 25 levemente obeso entre 25 e 30 excessivamente obeso entre 10 e 14 gordura ideal 3 gordura essencial Para homens o percentual de gordura entre 25 e 30 é normal entre 30 e 35 levemente obeso entre 35 e 40 excessivamente obeso entre 14 e 18 gordura ideal 12 gordura essencial Como fazer estaS categorizações? []'s. Edson Lira Estatístico Manaus-Amazonas

5 7

Modificar valores de dataframe com base em outro dataframe
by João José de Simoni Gouveia 13 Dez '11

13 Dez '11

Prezados, Possuo dois dataframes (nomeados aqui de 1 e 2). Preciso modificar os valores do dataframe1 com base no dataframe2. Mais precisamente: Se os valores da coluna snp1 (dataframe1) forem iguais ao valor da linha snp1 (dataframe2), preciso converter estes valores para 0 (zero), caso contrário o valor deve ser convertido para 1 (um). Alguém tem alguma dica? Att JJ Abaixo uma ilustração do que eu tenho DataFrame 1 --> snp1 snp2 snp3 snp4 1 A A C T 2 A G T G 3 A G C G 4 A A C T 5 T G T T DataFrame 2 --> snp anc 1 snp1 A 2 snp2 G 3 snp3 C 4 snp4 T E o que eu preciso ter: Dataframe final --> snp1 snp2 snp3 snp4 1 0 1 0 0 2 0 0 1 1 3 0 0 0 1 4 0 1 0 0 5 1 0 1 0 -- ============================================ João José de Simoni Gouveia Médico Veterinário, Msc. Doutorando em Zootecnia - PDIZ (UFC/UFPB/UFRPE) Professor Assistente I - Colegiado Acadêmico de Zootecnia Universidade Federal do Vale do São Francisco (UNIVASF) ---------------------------------------------------------------------------------------- Endereço: UNIVASF (Campus Ciências Agrárias) Rodovia BR 407, Km 12, Lote 543 Projeto de Irrigação Senador Nilo Coelho, s/nº-C1 Petrolina, PE - CEP 56.300-000 Telefone: (87)3986-3800/3804-3801 Skype: joao_jose_de_simoni_gouveia E-mail alternativo: jjsgouveia(a)gmail.com ============================================

3 2

MQM
by FHRB Toledo 13 Dez '11

13 Dez '11

Senhores, Tenho uma dúvida que não necessáriamente se relaciona com o programa, mas que em muito tem a ver com teoria estatística e algebra de matrizes, espero que aqui ache pelo menos uma sugestão para a resposta ou uma recomendação de literatura. O problema é dado o sistema de equações na forma matricial: Xb = y, pelo MQM, dada a matriz X e o vetor y, obtenho o vetor de soluções por X^-1y, respeitando todas as pressuposições. Se caso conheço ao invés de X e y, conheço os vetores b e y, existe alguma forma de conhecer a matriz X, podendo apenas assegurar que ela é uma matriz quadrada simétrica e positiva definida? att, FH

3 2

Rotina order
by Edson Lira 13 Dez '11

13 Dez '11

Bom dia a todos! Com a rotina abaixo eu consigo ordenar a tabela. inap<-table(tip_doa$INAPTIDAO,tip_doa$mês) inap2 = inap[order(rowSums(inap), decreasing=TRUE),] ftable(addmargins(inap2, c(1,2), FUN=list(Total = sum))) Que me dá a saida abaixo(as primeiras linhas) jan fev mar abr mai jun jul ago set out nov Total ESTEVE EM ZONA DE MALARIA 465 443 387 362 343 347 343 385 550 397 409 4431 HEMATOCRITO IRREGULAR 299 387 478 473 576 393 266 307 362 381 333 4255 VARIAÇÃO DE PARCEIRO SEM PRESERVATIVO (MAIS DE 4 EM 12 MESES) 249 329 240 152 100 183 150 86 110 124 106 1829 VARIAÇÃO DE PARCEIROS SEM PRESERVATIVO (MAIS DE 1 EM 6 MESES) 190 206 147 83 73 104 87 58 103 56 102 1209 VARIAÇÃO DE PARCEIROS SEXUAIS SEM PRESERVATIVO 0 0 20 106 238 135 109 138 134 142 160 1182 GRIPE 119 112 125 98 105 112 75 96 90 84 102 1118 Gostaria de acrescentar o ano e ordenar. Quando acrescento o ano me dá o seguinte problema. > inap<-table(tip_doa$INAPTIDAO,tip_doa$ano,tip_doa$mês) > inap2 = inap[order(rowSums(inap),decreasing=TRUE),] Erro em inap[order(rowSums(inap), decreasing = TRUE), ] : número incorreto de dimensões > inap<-table(tip_doa$INAPTIDAO,tip_doa$ano,tip_doa$mês) > inap2 = inap[order(rowSums(inap),decreasing=TRUE),] Erro em inap[order(rowSums(inap), decreasing = TRUE), ] : número incorreto de dimensões > ftable(addmargins(inap2, c(1,2,3), + FUN=list(Total = sum))) Erro em newdimnames[[margin]] : índice fora de limites Qual o problema? Alguém pode me ajudar? []'s. Edson Lira Estatístico Manaus-Amazonas

2 1

arquivo xlsx
by Edson Lira 13 Dez '11

13 Dez '11

Boa noite a todos, estou tentando ler um arquivo em excel (xlsx) com a rotina abaixo. require(xlsx) ag<-read.xlsx("/Users/edsonlira/IBGE/Censo_2010/Agregado/Agre2010_sinopse_AM.xlsx",1,encoding="macroman",row.names=1) Que está mostrando o log abaixo. Erro em .jcall("RJavaTools", "Ljava/lang/Object;", "invokeMethod", cl, : java.lang.OutOfMemoryError: Java heap space Será que é memória ou coisa parecida? Minha máquina é um MacBook Pro com 4 Gb de RAM (DDR3). Processador I7, 2,7 Ghz. O que pode estar errado? []'s Edson Lira Estatístico Manaus-Amazonas

4 3

Notação
by Edson Lira 12 Dez '11

12 Dez '11

Bom dia a todos, estou com um probleminha de notação. Tenho uma saída do R com a notação abaixo: 1.59578362e+00 1.52990264e+00 1.31029939e+00 1.17121733e+00 1.33225972e+00 1.31761950e+00 1.52185052e+01 Como mudar para a forma padrão 1.5957 1.529 ... 15.2 []'S. Edson Lira Estatístico Manaus-Amazonas

3 3

Extrair p valor de uma matriz de correlação
by ASANTOS 12 Dez '11

12 Dez '11

**Bom dia pessoal, Gostaria de extrair o p valor de uma matriz de correlação de modo similar a: x<-rnorm(100) y<-rnorm(100) cor.test(x, y, method='pearson')$p.value mais se faço, x<-rnorm(100) y<-rnorm(100) z<-rnorm(100) dados<-cbind(x,y,z) cor(dados, method='pearson')$p.value o cor.test() me permite obter esse valor, mas o cor(), não faz o mesmo, Alguém teria alguma sugestão, Obrigado -- Alexandre DOS SANTOS Engenheiro Florestal, Msc. Laboratório de Entomologia Florestal Departamento de Entomologia Universidade Federal de Lavras Caixa Postal 3037 37200-000 - Lavras - Minas Gerais - Brasil Tel: +55 35 92230304

3 6

Exclusão da lista
by David orais 12 Dez '11

12 Dez '11

Srs, Não estou conseguindo fazer a exclusão do meu email desta lista. Seria possível ao coordenador da confer~encia agilizar essa socilitação? Desde já agradeço. Atenciosamente, David Morais Universidade Candido mendes - RJ - RJ

1 0

Dúvida com arrays em for
by Felipe Buchbinder 12 Dez '11

12 Dez '11

Estou fazendo um programa de simulação cujo objetivo é comparar Teoria de Resposta ao Item com outros modelos de análise de variáveis latentes, como análise fatorial confirmatória com correlações de Pearson ou policórica. Em cada simulação, um array Exit é gerada e atribuído a um array tr (traço) que guarda todos os arrays "Exit" para posterior cálculo de média e desvio-padrão dos resultados (que serão, por sua vez, colocados no array Final). O problema é que, embora o array "Exit" seja gerado, os arrays tr e Final não o são. Não consigo entender o porquê....O código segue abaixo (retirei toda a parte de análise por TRI e CFA porque se não ficaria absurdametne grande) > > NSim<-5 #Number of simulations > > #Generates the trace matrix. trace[,,i] contains the i-th simulation of > the Exit matrix > tr<-array(rep(NA,7*14*NSim),dim=c(7,14,NSim)) > > > for(count in 1:NSim) > { > #%%%%%%%%%%%%%%%%%%%%%%%%%% > #THE ENTIRE CODE COMES HERE > #%%%%%%%%%%%%%%%%%%%%%%%%%% > #O código faz uma simulação de TRI. > #O código gera questões de um teste > #e respostas dadas por vários indivíduos com diferentes habilidades > #a essas questões. Essas respostas serão analisadas por TRI e análise > fatorial confirmatória. > #Os resultados destas análises são colocados nas várias células do array > Exit. > #Este código é rodado NSim vezes e os vários arrays Exit são guardados no > array tr > > > Exit<-array(rep(NA,7*14),dim=c(7,14)) > #O array Exit recebe os resultados de cálculos anteriores. > #Por exemplo, as linhas abaixo colocam na célula 1,1 do array o alpha de > Cronbach de um subteste chamado E1. > > Exit[1,1]<-alphaE1 > > #O array Exit existe. Aí, eu tenho passá-lo para o array tr. > > tr[,,count]<-Exit > } > Desde já, obrigado, Felipe > >

1 0

RStudio Server
by Pedro Rafael 11 Dez '11

11 Dez '11

Pessoal como configurar o RStudio Server? Estou querendo acessar um servidor da amazon. Segui as instruções no site do programa e o instalei no ubuntu com os comandos abaixo: $ sudo apt-get install libssl0.9.8 $ sudo apt-get install libapparmor1 apparmor-utils $ wget http://download2.rstudio.org/rstudio-server-0.94.110-i386.deb $ sudo dpkg -i rstudio-server-0.94.110-i386.deb Pelo que entendi posso agora na minha máquina rodar os scripts e ações na máquina da amazon onde está meu banco de dados. Pensei que apos instalar o RStudio Server seria instalado um novo programa RStudio pois na minha máquina já tinha o RStudio Desktop. Eu configuro o acesso no servidor via terminal ou na própria interface gráfica do RStudio Desktop? Obrigado à todos desde já. -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

3 3