- R-br - DINF UFPR

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

Plot - Séries Temporais
by Juliana Freitas de Mello e Silva 10 Dez '11

10 Dez '11

Bom dia! Eu estou tentando plotar uma série temporal, de 10 anos, através do comando plot(z) onde z é a série. E dá o seguinte erro: Erro em plotts(x = x, y = y, plot.type = plot.type, xy.labels = xy.labels, : cannot plot more than 10 series as "multiple" Qual seria o problema? Obrigada Abaixo segue o início da série> z Jan Fev Mar Abr Mai Jun Jul Ago Set Out Nov Dez1951 14.3 14.8 32.0 306.2 92.7 202.3 11.5 0.0 3.8 16.2 68.9 15.01952 83.1 90.8 566.0 390.5 134.1 27.1 23.4 9.6 7.5 6.0 0.9 7.21953 47.3 73.2 201.8 404.3 92.5 80.2 2.0 0.0 16.4 5.6 0.7 6.8

3 3

Re: [R-br] análise de sarampo - como proceder com o teste de hipóteses?
by Pedro Emmanuel Alvarenga Americano do Brasil 09 Dez '11

09 Dez '11

Uma serie temporal para diferentes faixas etarias? St nao e a minha praia apesar ja ter estudado ha muito tempo, talvez arima ou sarima voce consegue fazer uma funcao com uma banda de confianca para diferetes faixas etarias. Series temorais sao bastante intuitivas para os amigos que trabalham com vigilancia de doencas transmissiveis. Pedro Brasil via android (:)= Em 08/12/2011 11:37, "Walmes Zeviani" <walmeszeviani(a)gmail.com>escreveu: Daniela, A princípio, tratando-se do número de eventos por ano (que caracteriza uma Poisson sob certas suposições), faria análise com um modelo glm(..., family=poisson), incluindo o efeito do ano como preditora. Pode-se entrar com a faixa de idade também, e até mesmo, a interação entre elas, dado que você mecionou que ao longo do tempo diminui as crianças e aumentam os adultos com a doença (interação presente). Ainda é possível usar uma versão mais simples, fazendo uma tabela de contingência faixa_de_idade x ano e usar o teste chi-quadrado. Embora eu (e não só eu) goste de dar palpites sobre análises, recomendo fortemente que consulte um estatístico local. Falta de alguma informação importante e/ou minha falta de experiência, pode me fazer recomendar o errado ou menos correto para o caso. Mensagens desse caráter devem ser enviadas com prefixo "OFF-topic". À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

4 4

Por que o array tr não recebe a matriz Exit?
by Ivan Bezerra Allaman 09 Dez '11

09 Dez '11

Felipe, Desde quando "Exit" é uma matriz? Segundo seu código, "Exit" não é nada. Talvez sua intenção fosse: Exit <- matrix(ncol=X,nrow=Y) tr[,,count]<- Exit ????? (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 5

Categoria de referência (glm)
by Tiago Vieira 09 Dez '11

09 Dez '11

Boa tarde a todos! Hoje minha dúvida é: Como posso fazer para alterar a categoria de referênicias dos fatores em um modelo linear generalizado? Meu modelo é log-linear, e uso a função glm(resposta ~ factor1*factor2, poisson(link = "log"), dados). Desde já agradeço! Tiago Vieira.

2 2

análise de sarampo - como proceder com o teste de hipóteses?
by Daniela Recchia 09 Dez '11

09 Dez '11

Olá pessoal, a pergunta é mais teórica. Estou analisando o comportamento do sarampo ao longo dos anos. Aqui, em 1999 foi criada a lei que diz que as criancas devem ser vacinadas contra sarampo, cachumba e rubéola em 2 doses de vacinas, para a garantia de aproximadamente 99% de protecao contra as doencas. Mas a princípio, parece que ao decorrer do tempo (analiso de 2001 a 2011) o numero de criancas com sarampo realmente diminui, após essa campanha. Existem sim ainda algumas, devido a varios fatores (pais nao querem vacinar os filhos, as criancas só tomam 1 vacina...) mas o numero de adultos que apresentam sarampo tem aumentado, se comparado ao inicio da campanha. Por exemplo, se olhar os dados de doentes no peridio de 2003-2005 e 2009-2011 os adultos infectados nesses ultimos anos sao bem mais do que nos primeiros. (Na verdade nas idades de 15-39 anos) Entao a pincipio tenho dados para os doentes ao longo dos anos e por idade, algo como: > dados2 Altersgruppe 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 1 00-00 107 156 39 14 38 156 24 28 50 29 66 2 01-01 524 405 81 22 57 157 35 50 63 41 96 3 02-02 420 313 42 9 30 68 20 31 15 25 49 4 03-03 359 320 71 10 40 73 25 42 20 16 47 5 04-04 397 375 56 5 42 82 21 46 17 19 42 6 05-09 1594 1527262 17 250 446 123 210 83 127 293 7 10-14 1034 803 112 8 136 515 116 217 72 179 340 8 15-19 786 356 48 11 56 375 68 139 74 133 218 9 20-24 376 147 26 2 37 135 39 44 55 64 126 10 25-29 175 85 12 8 31 110 32 41 49 45 99 11 30-39 168 124 14 6 45 122 49 46 51 67 124 12 40-49 49 29 11 8 13 39 11 13 17 28 70 13 50-59 17 7 3 2 3 19 3 6 3 4 8 14 60-69 9 3 0 1 3 8 0 1 0 3 3 15 70 1 3 0 0 0 3 0 0 2 0 3 A minha idéia entao é rodar umt este estatistico, para a hipótese nula H0: a proporcao(?) ou o numero de adultos doentes no periodo 2003-2005 nao se diferencia significativamente dos doentes no periodo de 2009-2011. Tudo bem, meu próximo passo entao foi agrupar os dados, de talforma que tenho agora grupo de pessoas de 0 a 4 anos, 5 a 14 anos, 15 a 39 anos e acima de 40 anos. Para analisar por pediordos, simplesmente somei os casos para os anos de 2003-2005 e 2009-2011. Ok, agora tenho um numero por periodo e por grupo. Como compara-los efetivamente? Pois do jeito que retrabalhei os dados, nao consigo aplicar o teste T para 2 populacoes, nao tenho como medir normalidade e homocedasticidade. Na verdade, nao sei se o jeito que formatei os dados para a análise está correta para analisar a hipótese desejada. Como vcs fariam? Abracos, -- Daniela Rodrigues Recchia Master Student of Statistics - Technische Universität Dortmund. "It is better to have an inexact answer to the right question than a precise one to the wrong question!" John Tuckey

5 7

Re: [R-br] Plot - Séries Temporais
by Henrique Ewbank de Miranda Vieira 09 Dez '11

09 Dez '11

z é um data.frame com mais de 10 colunas? experimente usar o comando abaixo, onde colunadaSérie é o número da coluna do data.frame z que você quer plotar. ts.plot(z[,colunadaSérie]) Abs

1 0

Re: [R-br] OFF-TOPIC] análise de sarampo - como proceder com o teste de hipóteses?
by Fernando Colugnati 09 Dez '11

09 Dez '11

Legal Leonard, mas as coisas no Brasil estão indo na boa direção, e essas análises cruas não são prática comum somente aqui. Há 2 semanas atrás aconteceu o 3o ESAMP (escola de amostragem e metodologia de pesquisa), em Juiz de Fora, e estas questões foram muito discutidas, e saí de lá bem impressionado e motivado. Acho que em breve passaremos a ter boas análises feitas por bons bons analistas. Abs Em 9 de dezembro de 2011 11:31, Leonard de Assis <assis.leonard(a)gmail.com>escreveu: > Fernando, > > foi exatamente isso que quis dizer quando sugeri processos estocásticos > (Modelos de fluxo migratório, estudantes na escola, dentre outros são bem > representados por meio de processos estocásticos) . Um estudo demográfico > caberia bem no problema, principalmente porque neste problema, craramente > vemos que 10 casos em um ano e 10 no ano seguinte não são grandezas > comparáveis, a primeira providência a ser tomada é torná-la comparável, na > mesma ordem de grandeza. > > Quanto aos modelos lineares generalizados, creio ser possível utilizar, > pois existem modelos de efeito misto (Mixed effects models) que lidam com > estruturas de erro correlacionadas de maneira temporal. Nunca utilizei > estes modelos neste contexto, pode ser possível, teria que pensar mais em > cima deste problema. > > Este "erro" de se comparar grandezas incomparáveis vem de muito tempo no > Brasil, basta uma busca simples no Google ou em nossas mentes mesmo para > acharmos divulgação de "Número de acidentes nas rodovias" divulgado por > estado ou rodovia. Ora vejamos: Além de cada estado possuir malhas viárias > de tamanhos diferentes e com circulação de veículos diferente, temos > diversos outros fatores "esquecidos", como condições da rodovia, tipo de > relevo, etc. > > > []s > Leonard de Assis > assis <dot> leonard <at> gmail <dot> com > > > Em 09/12/2011 10:23, Fernando Colugnati escreveu: > > Pessoal, este é um problema clássico de epidemiologia e saúde pública, e a > comparação deve ser feita estimando-se as incidênicas (casos novos/ > ano*10.000habitantes, por exemplo, as unidades de tempo e população podem > mudar), que podem ser via modelos de poissson, Binom.-negativo ou algum > processo de contagem estocástico, como já sugerido. Mas não devem ser > esquecidas correções para taxas de crescimento populacional, etc...enfim, > como disse, é bem clássico. Este aumento pode ser um efeito de coorte > (age-cohort effect), por exemplo. > > Lembre-se de onde vêm estes dados, são nacionais, regionais? Vieram de > alguma base do MS (datasus)? É apenas a cobertura do SUS? A modelagem > estatística só pode ser bem definida se vc souber o que estes dados > representam, quais correções utilizar, etc...senão estará estimado apenas > viéses, que até é importante, mas pra outros objetivos. > > Abs > > Em 8 de dezembro de 2011 23:20, Leonard de Assis <assis.leonard(a)gmail.com>escreveu: > >> Dani, >> >> Existem várias formas de se analisar estes dados. Eu particularmente >> tentaria unir estas informações fazendo uma estratégia de processos >> estocásticos (Acho teorico demais, mas é mais desafiante). O desafio >> encontra-se na forte correlação temporal envolvida em tudo. >> >> Se partir para modelos lineares, creio que Poisson não seja uma boa >> idéia, pois ficva meio "fora de escala'. O ideal seria smnalisar a >> proporç~~ao que estes casos representam na população em estudo. >> >> Mas tudo isso é só palpite, cada um tem seu jeito de pensar e abordar >> determinados problemas. >> >> >> []s >> Leonard de Assis >> assis <dot> leonard <at> gmail <dot> com >> >> >> Em 08/12/2011 11:24, Daniela Recchia escreveu: >> >> >> Olá pessoal, >> >> a pergunta é mais teórica. Estou analisando o comportamento do sarampo >> ao longo dos anos. Aqui, em 1999 foi criada a lei que diz que as criancas >> devem ser vacinadas contra sarampo, cachumba e rubéola em 2 doses de >> vacinas, para a garantia de aproximadamente 99% de protecao contra as >> doencas. >> >> Mas a princípio, parece que ao decorrer do tempo (analiso de 2001 a >> 2011) o numero de criancas com sarampo realmente diminui, após essa >> campanha. Existem sim ainda algumas, devido a varios fatores (pais nao >> querem vacinar os filhos, as criancas só tomam 1 vacina...) mas o numero de >> adultos que apresentam sarampo tem aumentado, se comparado ao inicio da >> campanha. Por exemplo, se olhar os dados de doentes no peridio de 2003-2005 >> e 2009-2011 os adultos infectados nesses ultimos anos sao bem mais do que >> nos primeiros. (Na verdade nas idades de 15-39 anos) >> >> Entao a pincipio tenho dados para os doentes ao longo dos anos e por >> idade, algo como: >> >> > dados2 >> Altersgruppe 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 >> 1 00-00 107 156 39 14 38 156 24 28 50 29 66 >> 2 01-01 524 405 81 22 57 157 35 50 63 41 96 >> 3 02-02 420 313 42 9 30 68 20 31 15 25 >> 49 >> 4 03-03 359 320 71 10 40 73 25 42 20 16 47 >> 5 04-04 397 375 56 5 42 82 21 46 17 19 42 >> 6 05-09 1594 1527 262 17 250 446 123 210 83 127 >> 293 >> 7 10-14 1034 803 112 8 136 515 116 217 72 179 >> 340 >> 8 15-19 786 356 48 11 56 375 68 139 74 133 >> 218 >> 9 20-24 376 147 26 2 37 135 39 44 55 64 126 >> 10 25-29 175 85 12 8 31 110 32 41 49 45 99 >> 11 30-39 168 124 14 6 45 122 49 46 51 67 124 >> 12 40-49 49 29 11 8 13 39 11 13 17 28 70 >> 13 50-59 17 7 3 2 3 19 3 6 3 4 >> 8 >> 14 60-69 9 3 0 1 3 8 0 1 0 3 >> 3 >> 15 70 1 3 0 0 0 3 0 0 2 0 >> 3 >> >> A minha idéia entao é rodar umt este estatistico, para a hipótese nula >> H0: a proporcao(?) ou o numero de adultos doentes no periodo 2003-2005 nao >> se diferencia significativamente dos doentes no periodo de 2009-2011. >> >> Tudo bem, meu próximo passo entao foi agrupar os dados, de talforma que >> tenho agora grupo de pessoas de 0 a 4 anos, 5 a 14 anos, 15 a 39 anos e >> acima de 40 anos. Para analisar por pediordos, simplesmente somei os casos >> para os anos de 2003-2005 e 2009-2011. Ok, agora tenho um numero por >> periodo e por grupo. Como compara-los efetivamente? Pois do jeito que >> retrabalhei os dados, nao consigo aplicar o teste T para 2 populacoes, nao >> tenho como medir normalidade e homocedasticidade. Na verdade, nao sei se o >> jeito que formatei os dados para a análise está correta para analisar a >> hipótese desejada. >> >> Como vcs fariam? >> >> Abracos, >> >> -- >> Daniela Rodrigues Recchia >> >> Master Student of Statistics - Technische Universität Dortmund. >> >> "It is better to have an inexact answer to the right question than a >> precise one to the wrong question!" >> John Tuckey >> >> __._,_.___ >> >> <daniela_recchia@yahoo.com.br?subject=Res%3A%20an%E1lise%20de%20sarampo%20-%20como%20proceder%20com%20o%20teste%20de%20hip%F3teses%3F>| através >> de email<STAT-MATH@yahoogrupos.com.br?subject=Res%3A%20an%E1lise%20de%20sarampo%20-%20como%20proceder%20com%20o%20teste%20de%20hip%F3teses%3F>| Responder >> através da web<http://br.groups.yahoo.com/group/STAT-MATH/post;_ylc=X3oDMTJxOXVnYzhyBF9TAz…>| Adicionar >> um novo tópico<http://br.groups.yahoo.com/group/STAT-MATH/post;_ylc=X3oDMTJlZnJnMnAyBF9TAz…> >> Mensagens neste tópico<http://br.groups.yahoo.com/group/STAT-MATH/message/23373;_ylc=X3oDMTM2OWhkN…>( >> 1) >> Atividade nos últimos dias: >> >> - Novos usuários<http://br.groups.yahoo.com/group/STAT-MATH/members;_ylc=X3oDMTJmczk1dDY1BF9…> >> 5 >> >> Visite seu Grupo<http://br.groups.yahoo.com/group/STAT-MATH;_ylc=X3oDMTJldm9ycDRlBF9TAzk3NDk…> >> O que é melhor do que pagar pouco por algo que custa muito? Acesse Já.<http://global.ard.yahoo.com/SIG=15m93faov/M=804485.14750442.14520128.129601…> >> ------------------------------ >> Y! Encontros. É hora de dar uma chance a quem quer te conhecer!<http://global.ard.yahoo.com/SIG=15mnse4m1/M=758712.14532720.14632738.129601…> >> [image: Yahoo! Grupos]<http://br.groups.yahoo.com/;_ylc=X3oDMTJkZXY0MmYzBF9TAzk3NDkwNDM1BGdycElkAz…> >> Trocar para: SÃ³ Texto<STAT-MATH-traditional@yahoogrupos.com.br?subject=Mudar%20Formato%20de%20Envio:%20Tradicional>, >> Resenha DiÃ¡ria<STAT-MATH-digest@yahoogrupos.com.br?subject=Envio%20de%20email:%20Resenha>• Sair >> do grupo<STAT-MATH-unsubscribe(a)yahoogrupos.com.br?subject=Sair%20do%20grupo>• Termos >> de uso <http://br.yahoo.com/info/utos.html> >> . >> >> __,_._,___ >> >> >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça >> código mínimo reproduzível. >> > > > > -- > Fernando A.B. Colugnati > > > -- Fernando A.B. Colugnati

1 0

Statistics with R
by ASANTOS 09 Dez '11

09 Dez '11

Bom dia Pessoal, Só repassando uma novidade que recebi do R-Bloggers, um bom e completo site sobre estatística com R, em Statistics with R <http://zoonek2.free.fr/UNIX/48_R/all.html> em http://zoonek2.free.fr/UNIX/48_R/all.html Abraço, -- 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

1 0