Dezembro 2011 - R-br

Ajuda para plotar vários gráficos em uma única janela
by Maria Papa 27 Dez '11

27 Dez '11

Olá pessoal, boa tarde, Eu estou tentanto, de uma forma muito deselegante, plotar vários gráficos na mesma janela, mas somente o último gráfico é impresso. Por favor, alguém poderia me dizer o que está errado? par(mfrow = c(2,5)) for(i in length(temp)) { p = temp[,i] ht = hist(p, freq = F, col = "gray86", main = "", xlab = "Temperaturas", ylab = "Densidade") dn = density(p) lines(dn, col = 4, lwd = 2) } temp é uma matriz (9X900) com valores de temperatura. desde já agradeço a todos. Célia

2 1

Comparações Múltiplas em modelo linear com co-variável.
by Danilo Scorzoni Ré 26 Dez '11

26 Dez '11

Prezados, Estou realizando uma análise estatística através de um modelo de análise de covariância: y ~ sítio + bloco + tratamento + densidade + sitio:bloco + sitio:tratamento onde: sitio, bloco e tratametno são qualitativos densidade é quantitativo, foi inserido no modelo para reduzir a variação dos dados Os resultados da análise apontam efeitos significativos de tratamento e densidade. O desafio, agora, é fazer comparações múltiplas entre os tratamentos, no entanto, fixando-se densidades médias que variam entre os tratamentos, ou seja, descontando os efeitos de sítios e blocos que não são significativos. Tenho densidades médias para os quatro tratamentos: (400, 1200, 1150 e 900 para os tratamentos A, B, C e D), desta forma, utilizando os parâmetros estimados para cada tratamento + efeitos da densidade, gostaria de saber como eu poderia comparar as estimativas e verificar se os tratamentos diferem entre si. Eu poderia usar o tukey? Lembrando que a variável resposta é contínua... Caso não tenha sido claro, trago alguns resultados para exemplificar o problema. Desde já agradeço o empenho de todos e Feliz Natal! Att. -- Danilo Scorzoni Ré Engenheiro Florestal Mestrando em Ciência Florestal FCA / UNESP - Botucatu (14) 8180-2494

2 1

Transformação de variáveis.
by Juliana Freitas de Mello e Silva 25 Dez '11

25 Dez '11

Eu achei que esse tópico já existia, porém não o encontrei. Quero codificar observações de uma variável no meu banco de dados.

6 13

Função que faz o código voltar até determinada linha
by Felipe Buchbinder 23 Dez '11

23 Dez '11

Prezados r-istas, Existe em R alguma função que, caso uma determinada condição seja satisfeita, faça o código continuar a correr a partir de uma linha anterior? Por exemplo, algo do tipo: if(........) go to line 126 Sei que na programação da HP isso existe. Queria saber se existe também em R. Obrigado, Felipe Buchbinder

4 3

[Dúvida] Automatizar um modelo de previsão.
by Pedro Rafael 23 Dez '11

23 Dez '11

Pessoal a função auto.arima do pacote forecast "presta" realmente para ajustar modelos AR, MA, ARIMA, SARIMA? Existe algum pacote que faz de forma automática e eficiente do ponto de vista estatístico as análises de independências e de auto-correlações? Fico com o pé atrás de usar pacotes de prometem ajustar modelos auto-regressivos de forma automática. No meu caso tenho vários dados de mortalidade por doenças de aparelho circulatório, neoplasias, doenças do aparelho respiratório, etc. Todos os meses aqui na secretaria de saúde tenho mais informações e minha série vai aumentando. No meu computador tenho vários dados de milhares de variáveis distribuídos mensalmente e gostaria de criar algum programa para ajustar essas séries de forma automática. Não necessito do melhor modelo possível, apenas queria um bom modelo. É possível automatizar de tal forma a "dispensar" uma análise cuidadosa de um estatístico? O único problema é tempo para se ajustar série por série, ficar fazendo gráficos de auto correlações, verificando independências etc. Hoje possuo informações mensais desde 1999 a novembro de 2011, 143 pontos para cada uma das variáveis que preciso estimar. -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

4 4

Transformação Box-Cox
by Alexandre Santos 23 Dez '11

23 Dez '11

Boa tarde pessoal, Tenho uma variável resposta constituída de medidas de peso, porém alguns valores são representados por zeros, estava tentando fazer uma transformação de box-cox para tentar normalizar os dados, porém quando fiz: boxcox(peso ~ trat, data = dados, plotit = T) Erro em boxcox.default(peso ~ trat, data = dados, plotit = T) : response variable must be positive Apareceu o seguinte erro, então gostaria de saber se são os zeros que estão causando o erro? Obrigado, Alexandre dos Santos Engenheiro Florestal, MSc. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 3829-5122

4 6

duvida cat ou paste
by Vinicius Brito Rocha 22 Dez '11

22 Dez '11

Senhores, estou tentando transformar essas linhas em um único objeto para realizar uma contagem, ou seja, preciso que as 6 colunas tornem-se apenas 1. tentei o cat e o paste combinados com o apply, mas não tive sucesso dado original: [,1] [,2] [,3] [,4] [,5] [,6] [1,] 4 5 30 33 41 52 [2,] 9 37 39 41 43 49 [3,] 10 11 29 30 36 47 [4,] 1 5 6 27 42 59 [5,] 1 2 6 16 19 46 objetivo: [1] "4 5 30 33 41 52" [2] "9 37 39 41 43 49" ... Abs -- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

4 4

Boas festas de fim de ano
by Rodrigo Muniz 22 Dez '11

22 Dez '11

Desejo a todos os companheiros de R, um feliz natal, feliz ano novo. Gostaria de agradeçer a todos pela ajuda dada, começei a alguns meses a ser usuário de R e pretendo continuar usando, mais não teria persistido se não fosse a ajuda de vcs. Mais uma vez agradeço a todos. 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)

1 0

Re: [R-br] Ajuda para combinar Levels
by Ivan Bezerra Allaman 22 Dez '11

22 Dez '11

Usando a imaginação! rep <- combn(3,2) results <- paste(rep,'vs',combn(3,2),sep='')[-c(1,3,5)] (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}

3 3

Pacote Experimental Designs
by Eric Ferreira 21 Dez '11

21 Dez '11

Prezados colegas, Com satisfação, informo que o pacote ExpDes (Experimental Designs), dedicado a análise de variância de experimentos de modelo fixo balanceados, após 2 anos de desenvolvimento e divulgação, se encontra no CRAN. Sugestões de melhorias e indicação de "bugs" são muito bem vindas! Atenciosamente, -- Prof Dr Eric B Ferreira Exact Sciences Department Federal University of Alfenas Brazil

3 2