Novembro 2012 - R-br

Ajuda com NA
by Sérgio Henrique almeida da silva ju 28 Nov '12

28 Nov '12

Olá Estou tentando excluir os NAs do meu banco, já fiz de tudo, inclusive mandinga, mas não consigo: structure(list(afia3 = c(0, 1, 0, 3, 0, 0, 0, 0, 0, 0, 0, 5, 0, 0, 0, 0, 2, 4, 0, 3, 0, 3, 0, 0, 0, 3, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 5, 0, 1, 0, 0, 0, 0, 0, 0, 3, 0, 0, 0, 0, 0, 0, 1, 2, 0, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 7, NA, 0, 0, 0, 0, 2, 2, 0, 3, 0, 0, 0, 0, 5, 0, 3, 1, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 4, 5, 0, 1, 2, 0, 0, 0, 2, 0, 0, 0, 0, 0, 4, 4, 0, 0, 0, 0, 0, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3, 0, 0, 3, 0, 3, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 2, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 3, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 3, 0, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 2, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3, 0, 0, 2, 0, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3, 0, 0, 3, 2, 5, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3, 5, 2, 1, 3, 3, 0, 0, 0, 0, 6, 0, 0, 0, 4, 0, 0, 0, 0, 0, 3, 4, 4, 3, 3, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 2, 1, 0, 0, 3, 2, 0, 0, 0, 0, 0, 0, 0, 7, 0, 0, 0, 0, 0, 0, 0, 0, 0, 5, 0, 0, 0, 0, 0, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3, 3, 0, 2, 0, 0, 3, 0, 0, 0, 0, 3, 0, 0, 0, 0, 3, 0, 0, 0, 0, 0, 0, 0, 5, 0, 4, 0, 0, 0, 2, 0, 0, 7, 0, 0, 0, 0, 0, 0, 5, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 3, NA, 0, 0, 0, 5, 0, 0, 4, 0, 0, 0, 0, 0, 0, 5, 6, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 6, 2, 3, 0, 3, 0, 3, 0, 0, 0, 0, 0, 0, 0, 3, 3, 3, 0, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 4, 6, 5, 3, 0, 0, NA, 0, 0, 5, 0, 0, 0, 4, 0, 0, 0, 6, 0, 0, 0, 0, 3, 0, 0, NA, 0, 0, 0, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3, 0, 0, 4, 0, 0, 0, NA, 0, 0, 0, 2, 0, 0, 0, 7, 0, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 0, 5, 3, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3, 0, 0, 0, 0, 0, 0, 0, 0, 1, 3, 0, 3, 5, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 1, 0, 0, 0, 0, 0, 2, 0, 0, 3, 0, 0, 0, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 5, 0, 0, 3, 0, 0, 4, 0, 0, 2, 3, 0, 0, 0, 0, 2, 0, 0, 0, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 3, 4, 0, 5, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3, 0, 2, 0, 0, 0, 0, 3, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 5, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 4, 0, 0) Como posso fazer isso? Abraços -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

6 11

ESPERANÇA DE QUADRADOS MÉDIOS
by tiago souza marçal 28 Nov '12

28 Nov '12

Boa noite prezados colegas, Vocês poderiam me indicar algum material onde eu tenha acesso a esperança de quadrados médios do delineamento em quadrado latino e arranjo de parcelas subsubdivididas. A baixo segue um exemplo do que estou querendo: Modelo: Yijk = m + Gi + Aj + GAij + Eijk FV GL E(QM) Ambiente a-1 σ2 + gσ2a Genótipo g-1 σ2 + arσ2a G x A (a-1)(g-1) σ2 + rσ2ga Resíduo ga(r-1) σ2 Em que: a - é o número de ambientes. g - é o número de genótipos. r - é o número de repetições Agradeço desde já pela ajuda. Att. Tiago.

3 4

Sobre o R Revolution - para esclarecimentos
by Diogo Ferrari 28 Nov '12

28 Nov '12

Caros, Eu tenho mencionado o Revolution, e vejo que algumas perguntas tem pipocado. Deixe eu tentar esclarecer algumas coisas que talvez seja util para as necessidades de alguns usuários. O Revolution faz tudo que o R faz. Ele é uma espécie de IDE, e a única diferença é que ele tem algumas funções a mais feitas em java pelos desenvolvedores. Essas funções a mais são para uso de Big-data e HPC, só isso. Tais funções fazem parte do pacote Revoscale (um ?Revoscale no R Revolution informa bastante), que já vem com o R Revolution e carregam automaticamente quando você abre o software. O Revolution é gratuito para uso acadêmico, e tem quase todas as funcionalidades da versão paga, exceto algumas coisas para HPC. Na prática, não faz falta. Basta pedir para baixar essa versão no site, fazer o mini cadastro que você recebe no e-mail o key para instalar e o link para baixar. A grande vantagem na minha opinião é o uso do big-data. Para isso, o Revolution tem um formato de arquivo que facilita o processamento. Assim como o SPSS tem o .sav, o Stata o .dta, o Revolution le todos, mas tem um arquivo próprio .xdf, que torna a análise mais eficiente (.xdf é um arquivo indexado, etc., etc.). Minha experiência: estou analisando 6 censos brazileiros (5x +/- 20.000.000 de casos em cada). Com um pouco de programação estatística + Revolution + .xdf, faço uma logit, de uma vez, nesses 5 censos (e mais 1 com +/-200.000 casos) em 10 segundos usando um note care2duo, 3GB e 160GB. Não encontrei desempenho assim em qualquer outro software. Para qualquer outra coisa que não envolva trabalhar com dados na hd ao invés da memória, ou que trabalhe com tamanho big-data, o R sem essas funcionalidades a mais funciona muito bem. Abs. -- Diogo Ferrari University of San Paulo (USP) Faculty of Philosophy, Languages, Literature and Human Sciences (FFLCH) Department of Political Science (DCP) San Paulo/SP - Brazil E-mail: diogo.ferrari(a)usp.br Open Source! Use R! Use Linux! "A vida é a arte de tirar conclusões suficientes de dados insuficientes"

4 5

Modelo aleatorio fixo ou variado
by Ze Henrique 28 Nov '12

28 Nov '12

Boa noite colegas. Gostaria que alguem me ajuda-se com algum material ou me explique quando utilizar modelo fixo e modelo aleatório. Obrigado desde já a todos.

2 2

Divulgação de R como software livre de estatística
by Richard Santos 28 Nov '12

28 Nov '12

Caros, Pelo que estou aprendendo do R, penso em alguns pontos importantes sobre como falar para as empresas e instituições do governo a respeito do software livre em Estatística, o único que existe que atende requisitos mínimos, o R. Peço que apontem quando discordar de algum ponto. Obrigado. - uma programação mais extensa do que os outros sotwares estatísticos viáveis (que são pagos, a maioria não tão caros - Stata: 3.290 reais, Minitab: 2990 reais, o próprio Revolution R: 2000 reais, SPSS: um pouco mais, SAS: mais caro); - demanda mais tempo de aprendizado do que os demais também, antes de começar a usar; - uma das soluções para grandes bases e com point in click do R é paga, então acaba não tendo a vantagem de ser livre; - sem ser pago, o R também tem soluções de grandes bases e de point in click, mas é preciso já conhecer R de intermediário para avançado para conseguir instalar e conseguir usas essas opções. A maioria das pessoas que falam que o R faz tudo e gratuito, defende mas sem na hora saber fazer, tendo apenas o nível básico do software. Ainda preciso ver alguém com razoável conhecimento do R usando grandes bases de dados, ou testar os pacotes sugeridos pela lista (vou fazer isso). Richard

3 2

código em R
by alanarocha＠sapo.pt 28 Nov '12

28 Nov '12

modC<-lm(pthi~dose+grupo) modC Call: lm(formula = pthi ~ dose + grupo) Coefficients: (Intercept) dose grupoH 752.8614 0.1791 -159.5864 summary(modC) Call: lm(formula = pthi ~ dose + grupo) Residuals: Min 1Q Median 3Q Max -862.97 -323.87 -60.89 222.13 1936.94 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 752.86139 84.56014 8.903 1.22e-15 *** dose 0.17913 0.08253 2.171 0.0315 * grupoH -159.58641 92.49382 -1.725 0.0864 . --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 540.3 on 158 degrees of freedom (1 observation deleted due to missingness) Multiple R-squared: 0.03447, Adjusted R-squared: 0.02225 F-statistic: 2.82 on 2 and 158 DF, p-value: 0.0626 eu fiz isto, percebo que o grupo C aparece por defeito do R. ams como esplico no relatorio qual o valor estimado para o grupo C e para o grupo H? E se no meu ficheiro de dados o H aparece primeiro porque é que o R considera o C primeiro? Ana

2 1

código em R
by alanarocha＠sapo.pt 28 Nov '12

28 Nov '12

modC<-lm(pthi~dose+grupo) modC Call: lm(formula = pthi ~ dose + grupo) Coefficients: (Intercept) dose grupoH 752.8614 0.1791 -159.5864 summary(modC) Call: lm(formula = pthi ~ dose + grupo) Residuals: Min 1Q Median 3Q Max -862.97 -323.87 -60.89 222.13 1936.94 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 752.86139 84.56014 8.903 1.22e-15 *** dose 0.17913 0.08253 2.171 0.0315 * grupoH -159.58641 92.49382 -1.725 0.0864 . --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 540.3 on 158 degrees of freedom (1 observation deleted due to missingness) Multiple R-squared: 0.03447, Adjusted R-squared: 0.02225 F-statistic: 2.82 on 2 and 158 DF, p-value: 0.0626 eu fiz isto, percebo que o grupo C aparece por defeito do R. ams como esplico no relatorio qual o valor estimado para o grupo C e para o grupo H? E se no meu ficheiro de dados o H aparece primeiro porque é que o R considera o C primeiro? Ana

1 0

Separando nomes com espaços em branco após o nome
by Fátima Lima Paula 28 Nov '12

28 Nov '12

Prezados, já aprendi a separar nomes, mas agora surgiu um problema: esse banco que estou usando, foi exportado do sihd2 e vem com os espaços do campo que é para ser preenchido. Assim, tenho nomes com esse aspecto: nome=c("Maria Tereza dos Santos ", "Jorge Roberto da Silva ", "José de Faria Almeida ") nome=factor(nome) idade=c(76,87,98) data=data.frame(nome, idade) Preciso que os nomes sejam divididos da seguinte forma: data$nome1=c("Maria", "Jorge", "José") data$nome2=c("Tereza dos", "Roberto da","de Faria") data$nome3=c("Santos", "Silva", "Almeida") Fui usar o script que guardei do Walmes, mas não deu certo. O problema está nos espaços que ficam depois do nome. Só para ilustrar, vejam o que retorna: separado=strsplit(as.character(nome),split=" ") > head(separado) [[1]] [1] "Maria" "Tereza" "dos" "Santos" "" "" "" "" [9] "" "" "" "" "" [[2]] [1] "Jorge" "Roberto" "da" "Silva" "" "" "" [8] "" "" "" "" "" "" [[3]] [1] "José" "de" "Faria" "Almeida" "" "" "" [8] "" "" "" "" "" "" "" [15] "" "" "" "" O pior é que não existe um número constante de espaços, nem de total de linhas. É exatamente assim como no exemplo. Alguém pode me ajudar, por favor?

2 2

dúvida em código
by alanarocha＠sapo.pt 28 Nov '12

28 Nov '12

sim eu já tinha removido essa linha e o problema é o mesmo.

1 0

dúvida em codigo
by alanarocha＠sapo.pt 28 Nov '12

28 Nov '12

Eu não estou a perceber porque é que o meu codigo le 3 níveis quando devia ser so 2 H e C. https://dl.dropbox.com/u/9337305/regressão.R Ana

2 1