- R-br - DINF UFPR

Módulo minguante em acf para série periódica
by Alexandre Castagna 14 Dez '11

14 Dez '11

Caros, Estudando séries temporais me deparei com a seguinte questão. Ao gerar um série periódica como por exemplo: a = rep(c(1,-1),50) esperava que a função de autocorrelação refletisse correlação 1 a cada dois lags temporais e -1 a cada lag temporal uma vez que a = {1, -1, 1, -1, 1, -1, 1, -1,...}, porém a função de autocorrelação do R me retorna correlações cujo módulo reduz com o aumento do lag: b = acf(a) onde b = {1.00, -0.99, 0.98, -0.97, 0.96, -0.95, 0.94, -0.93, 0.92, -0.91, ...} Alguém pode explicar porque? Obrigado, -- Alexandre Castagna Mourão e Lima Laboratório de Radioecologia e Mudanças Globais - UERJ Rua São Francisco Xavier, 524, Maracanã Pav. Haroldo Lisboa da Cunha, Sub-solo

1 0

MCMCglmm - Contraste
by Eder David Borges da Silva 14 Dez '11

14 Dez '11

Bom pessoal , Venho mais uma vez discutir uma duvida computacional e teorica. imaginem o seguinte modelo: require(MCMCglmm) dados <- expand.grid(LOCAL=factor(1:2), TRAT1=factor(1:2), TRAT2=LETTERS[1:4], DATA=factor(1:10), REP=1:3) head(dados) dados$RESP <- rpois(nrow(dados),5) model<-MCMCglmm(RESP~LOCAL/TRAT1/TRAT2, random=~DATA,family='poisson', data=dados, verbose=FALSE) summary(model) Dentro deste modelo meu interesse principal é o TRAT2 (4 níveis), pelo resultado mostrado a matriz de efeitos fixos X é montada contrastando os níveis B C e D com A, porem gostaria de outros constastes como por exemplo : ## Contraste (LOCAL1:TRAT11:TRAT2B) - (LOCAL1:TRAT11:TRAT2C) se tivesse a posteriori direta apenas de B e C bastaria fazer apenas uma menos a outra porem como tenho ela B-A e C-A, eu creio não poder fazer isso, ou estou equivocado? Uma das opções que tentei foi mudar o data.frame(dados) alterando a ordem de TRAT2 esperando que mudasse a configuração dos contraste mais isso não aconteceu, ou fiz errado, desta forma alguém tem uma dica de como mudar este matriz para que eu obtenhas os outros contrastes desejados? Ate mais Éder

2 2

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

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

Por que o array tr não recebe a matriz Exit?
by Felipe Buchbinder 08 Dez '11

08 Dez '11

1. Estou tendo um problema em um código que não sei por quê deveria dar problema… Essencialmente, o código é este: NSim<-5 #Number of simulations #Generates the trace array. 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) { #Aqui vem o cerne do código… #A cada iteração, o código faz vários cálculos #e armazena em uma matriz chamada Exit tr[,,count]<-Exit #A matriz i-ésima da iteração i=count é armazenada no array tr. #O problema é que isto não acontece. #Ou seja: quando eu digito Exit, a matriz aparece; #Mas quando digito tr, só recebo NA } Eu nao consigo entender porque a matriz Exit existe mas a matriz tr não a está recebendo. Alguém sabe o por quê? Desde já, obrigado pela ajuda! Felipe -

1 0

Mapas Brasil Município
by Vinicius Brito Rocha 08 Dez '11

08 Dez '11

Senhores, Estou precisando de um plot do mapa do brasil por municípios. A ideia é associar uma estatística e ver a distribuição de forma visual. Alguém pode me dar uma luz? 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.

3 2

Re: [R-br] R help
by Edmar Alino 07 Dez '11

07 Dez '11

Prezados A última coluna é resultado de um sistema tipo black box e nao tenho acesso aos modelos internos infelizmente. Ademais os modelos são não lineares e algumas vezes apresenta descontinuidades. Que tipo de inferência estatística posso almejar? Informações sobre a importância (peso) de cada variável, inter-relação entre as variáveis ou bloco de variáveis, sensibilidade, etc. Como o software R pode me auxiliar? Grato antecipadamente pela ajuda! Edmar Alino Em 07/12/2011 07:55, "FHRB Toledo" <fernandohtoledo(a)gmail.com>escreveu: Edmar, Se sua ultima coluna é a variável (resposta -- dependente) uma boa pedida é usar a função lm(), de ajuste de modelos lineares. Tente nos dizer mais sobre que tipo de relacionamento que as 20 primeiras colunas tem com a resposta, tenho certeza que você vai conseguir ajuda. abraço, FH 2011/12/7 Edmar Alino <edmaralino(a)gmail.com> > > > > Srs, > > > > Tenho uma matriz M(1000,21), onde a ultima coluna (21a) eh funcao das > demais colunas. E... > _______________________________________________ > 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. > _______________________________________________ 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 5

Ajuda com lapply
by Emmanuel Arnhold 07 Dez '11

07 Dez '11

Como usar a função lapply (ou então outra função se não for possível via lapply) de modo que a lista de parâmetros ‘a’ seja inclusa no procedimento abaixo? (~ai+bi*x????) a=c(1,2,3) b=c(0.2,1.6) x=c(1,2,3,4) # como incluir ai no modelo abaixo? exemplo1=lapply(b, function(bi) {bi*x}) exemplo1

2 1

GPU x CPU e aplicações reais
by Benilton Carvalho 07 Dez '11

07 Dez '11

http://wp.me/pbygO-gr

1 0

R help
by Edmar Alino 07 Dez '11

07 Dez '11

Srs, Tenho uma matriz M(1000,21), onde a ultima coluna (21a) eh funcao das demais colunas. Eh basicamente um sistema black box e preciso usar o R e suas ferramentas estatisticas para analisar a influencia de cada coluna (variavel), sua dominancia e o efeito na variavel de interesse (ou saida). Uma das variaveis de entrada eh o tempo. Como nao sou estatistico e ainda estou engatinhando no R, gostaria da vossa ajuda. Grato antecipadamente, Edmar

2 1

Criação de função fortran em pacote (base)
by Junior Beleti 07 Dez '11

07 Dez '11

Boa noite a todos.. Estou realizando estudos na função eigen.f, devido a seu elevado tempo de processamento em meu experimentos. Gostaria de saber como devo proceder para criar uma nova função fortran (por exemplo eigen2.f) e compilá-la juntamente no código R. Tentei criar tal função na pasta src/appl mas não obtive sucesso pois no momento da compilação gera um erro. Se alguem puder me auxiliar ficaria grato. Att, -- Carlos Roberto Beleti Junior Mestrado em Ciência da Computação Departamento de Informática Universidade Estadual de Maringá

3 2

Predição em regressão múltipla
by Edson Lira 06 Dez '11

06 Dez '11

Boa tarde a todos. Tenho um modelo de regressão múltipla com 4 variáveis preditoras, lp=1.489+1.41*10^(-2)elevacao+7.42*10^(-3)data-3.18*10^(-1)inudacao-4.41*10^(-5)esgoto Se eu quise fazer predição para elevação 9,10,14. como faço? []'s. Edson Lira Estatístico Manaus-Amazonas

3 2

Plotar uma lista de funções em um único gráfico
by Emmanuel Arnhold 06 Dez '11

06 Dez '11

Gostaria de algo semelhante ao operacionalizado pela função fabaixo. Mas gostaria de plotar uma lista de funções e não de diferentes y como é exemplificado, ou seja, as equações obtidas pelos diferentes y. x=c(1,2,3,4,5,6) y1=c(1,2,3,4,5,6) y2=c(10,11,12,13,14,15) y3=c(15,14,13,12,11,10) y4=c(6,5,4,3,2,1) f=function(x,...){x=x; y=data.frame(...); matplot(x,y, type="o")} f(x,y1,y2,y3,y4) Como fazer um processo semelhante com funções? Talvez partindo daqui? Usando lapply ou outro recurso?? a=list();a$a1=2;a$a2=3 x=c(1,2,3,4) l1=lapply(a, function(ai) {ai*x}) l1 a=c(1,2,3,4) x=c(1,2,3,4,5,6) l2=lapply(a, function(ai) {ai*x}) l2

5 7

regressão linear sem intercepto
by Daniele Portela de oliveira 06 Dez '11

06 Dez '11

Caros colegas gostaria de saber como faço para fazer uma analise de regressão linear no R sem considerar o intercepto. Abraços Daniele Portela de Oliveira Mestranda Zootecnia - UEM Cel: 44 9949 3525 67 9617 0965

2 2

Selecionar Contorno de região
by Luís Gustavo 06 Dez '11

06 Dez '11

Olá pessoal, Estou com o seguinte problema: Tenho os estados que compõem a região sudeste do Brasil, porém gostaria de selecionar apenas as coordenadas da região (o contorno da região sudeste). Ou seja criar um polígono da região sudeste. O objetivo é reduzir o tamanho do meu shape. Muito obrigado. -- Luís Gustavo Silva e Silva

3 2

Como categorizar variáveis
by Tiago Vieira 06 Dez '11

06 Dez '11

Boa noite, Minha dúvida é a seguinte: tenho um conjunto de variáveis dicretas e quero categorizalas. Como posso fazer isto no R de forma masi elegante? Já que, estou tendo que usar if's ou vários ifelse! Obrigado pela atenção! Taigo Vieira.

3 2

Ajuda
by Gilbert Queiroz 05 Dez '11

05 Dez '11

Pessoal, tenho a seguinte parte de um script que estou desenvolvendo: n<-c(87, 108, 80, 25) err.m<-matrix(NA, choose(nx,2),1) err.m nx for(i in 1:nx) { for(j in 2:nx) { if(i<j){ err<-sqrt((410.35/(n[i]+0.5))+(410.35/(n[j]+0.5))) print(err) } } } Informações úteis: nx = 4 (já foi definido em outra parte do script) err.m é uma matriz-coluna (6 x 1) O quero fazer é colocar os resultados de err dentro da matriz err.m Como faço isso??? O comando print me mostra que err está calculando o que quero. Mas me interessa colocar os resultados em err.m Abs.

5 6

Função com problema no Rprofile
by Emmanuel Arnhold 05 Dez '11

05 Dez '11

Criei uma função bem simples para simplificar a digitação da read.table("clipboard", h=T) que me é bastante útil. Mas a função não funciona via Rprofile. Alguém tem alguma sugestão? Testem a função do R e depois tentem inserir no Rprofile e usar novamente para ver o erro. A função é: clip=function(x=read.table("clipboard ", h=T)){x=read.table("clipboard", h=T); print(x)} clip()

2 1

Re: [R-br] teste Duncan
by Bruna Gabriela Wendpap 05 Dez '11

05 Dez '11

Há um pacote não oficial do R para análise de experimentos chamado ExpDes (Experimental Designs). A página abaixo contém link para download do pacote e tutorial em português. https://sites.google.com/site/ericbferreira/unifal/downloads-1 Att, Bruna Gabriela Wendpap > From: r-br-request(a)listas.c3sl.ufpr.br > Subject: Digest R-br, volume 10, assunto 4 > To: r-br(a)listas.c3sl.ufpr.br > Date: Sat, 3 Dec 2011 12:00:01 -0200 > > Enviar submissões para a lista de discussão R-br para > r-br(a)listas.c3sl.ufpr.br > > Para se cadastrar ou descadastrar via WWW, visite o endereço > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > ou, via email, envie uma mensagem com a palavra 'help' no assunto ou > corpo da mensagem para > r-br-request(a)listas.c3sl.ufpr.br > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)listas.c3sl.ufpr.br > > Quando responder, por favor edite sua linha Assunto assim ela será > mais específica que "Re: Contents of R-br digest..." > > > Tópicos de Hoje: > > 1. Re: teste Duncan (Andre Oliveira) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Fri, 2 Dec 2011 13:25:26 -0800 (PST) > From: Andre Oliveira <andreolsouza(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br>, > "crabak(a)acm.org" <crabak(a)acm.org> > Subject: Re: [R-br] teste Duncan > Message-ID: > <1322861126.39856.YahooMailNeo(a)web161805.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Uma alternativa é o pacote laercio. > require(laercio) > ?LDuncan > > > > André Oliveira Souza > > > ________________________________ > De: Cesar Rabak <cesar.rabak(a)gmail.com> > Para: r-br(a)listas.c3sl.ufpr.br > Enviadas: Sexta-feira, 2 de Dezembro de 2011 10:42 > Assunto: Re: [R-br] teste Duncan > > Em 2/12/2011 09:31, Luiz André Zardo escreveu: > > Prezados, > > Qual pacote posso utilizar para rodar o teste Duncan para comparação de > > medias no R?? > No prompt do seu terminal R digite o seguinte: > > > RSiteSearch("Duncan test") > > Na primeira página dos resultados já vi dois pacotes. . . > > > -- Cesar Rabak > GNU/Linux User 52247. > Get counted: http://counter.li.org/ > _______________________________________________ > 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. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111202/f2833cee/atta…> > > ------------------------------ > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > Fim da Digest R-br, volume 10, assunto 4 > ****************************************

1 0

teste Duncan
by Luiz André Zardo 04 Dez '11

04 Dez '11

Prezados, Qual pacote posso utilizar para rodar o teste Duncan para comparação de medias no R?? Obrigado.

4 3

Distribuição de Weibull
by Rodrigo Muniz 04 Dez '11

04 Dez '11

Boa noite pessoal. Estou com o seguinte script: require(MASS) getwd() setwd("C:/Users/Rodrigo/Documents/Pos_graduacao/DoutoradoUSP/Geoestatistica") Vento <- read.table("Weibull_full.txt", head=T) WVento50m <- fitdistr(Vento[[1]], densfun="weibull") WVento30m <- fitdistr(Vento[[2]], densfun="weibull") WVento15m <- fitdistr(Vento[[3]], densfun="weibull") WVento10m <- fitdistr(Vento[[4]], densfun="weibull") WVento3m <- fitdistr(Vento[[5]], densfun="weibull") O arquivo Weibull_full é um data frame com 5 coluna sendo a variável velocidade do vento em cinco diferentes níveis de altura. Estou tentanto calcular os parâmetros de forma e escala da distribuição de Weibull para cada nível de altura. No entanto, só consigo calcular para o nível de altura de 50 m (WVento50m), quando mando rodar os outros níveis dá o seguinte erro: Error in fitdistr(Vento[[2]], densfun = "weibull") : Weibull values must be > 0 Já conferi o arquivo dos dado para ver se há algo errado com os dados, mais não há. Tentei fazer a análise separada mais deu o mesmo erro. Um observação que eu faço é que embora eu não consiga rodar a função "fitdistr", eu consigo plotar os gráficos de histograma normalmente para todos os níveis de altura. Se alguém puder me ajudar eu agradeço. Abraço a todos. -- *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)

3 4

Algoritmo para simular VA GAMA
by Victor Eduardo 02 Dez '11

02 Dez '11

Gente, onde posso encontrar um algoritmo para simular VA Gama?? Eu uso o livro Simulation do Sheldon Ross, e ao que parece, lá não tem esse algoritmo. Posso simular uma gama quando o alfa for um numero inteiro por soma de exponenciais, e quando o mesmo não for, ou seja, e quando, por exemplo, o alfa for < 1, como poderia fazer? OBS: As restrições já foram adicionadas na minha função, ou seja, o alpha tem que ser positivo e o beta também. Abraços!

2 1

Estimar vários modelos não lineares de uma só vez
by Emmanuel Arnhold 02 Dez '11

02 Dez '11

Gostaria de saber se é possível fazer (e como fazer se for possível) via nls algo semelhante ao que segue abaixo via lm, ou seja, estimar vários modelos de uma só vez. No caso o modelo tem a mesma variável explicativa e os mesmos starts, mudam as variáveis respostas. x=c(0:5);y1=c(50:55);y2=c(500:505) modelo=lm(cbind(y1,y2)~x);modelo

3 3

Re: [R-br] Versão R-2.14.0 no Ubunto - RESOLVIDO
by marlene marchena 01 Dez '11

01 Dez '11

Oi! Somente para fechar este post. Depois de pedir ajuda em outros forums descubrí que tinha a versão R (2.13.01) instalada manualmente (sem usar a infraestructura apt-get) no meu diretório HOME. Por isso não conseguia desinstalá-lo. Bom, a única coisa que fiz foi deletar a pasta R 2.13.01 do diretório HOME e reinstalar a versão 2.14, isso resolviu tudo. Obrigada pela ajuda! Marlene Em 27 de novembro de 2011 20:53, Benilton Carvalho < beniltoncarvalho(a)gmail.com> escreveu: > R --vanilla > > No terminal > On Nov 27, 2011 4:45 PM, "marlene marchena" <marchenamarlene(a)gmail.com> > wrote: > >> Gente estou precisando de ajuda! >> >> Eu instalei o R 2.14 seguindo o link acima e os passos indicados pelo >> Pedro no ubuntu 11.10. Eu rodei o R no terminal e apareceu q a versão era >> a 2.14. Porém agora pouco eu percebi que após reiniciado o computador meu R >> volta para a versão que timha antes a 2.13.1 >> >> Não entendo o que fiz de errado? >> >> Alguém pode me ajudar a entender o que acontece. >> >> Grata, >> >> Marlene. >> >> >> Em 25 de novembro de 2011 13:57, Daniel C Bezerra < >> danielcbezerra(a)gmail.com> escreveu: >> >>> Olá Marlene (e quem estiver iniciando Ubuntu), >>> >>> Já que vc está usando Ubuntu 10.04 seguem mais duas sugestões que podem >>> te ajudar: >>> >>> Agora que vc instalou o repositório do R, vale à pena vc entrar em >>> Sistema>Administração>Gerenciador de Pacotes Synaptic e procurar por todo >>> os programas que tenham "r-cran" (sem aspas). Vc vai ver diversos pacotes >>> do R disponíveis. >>> >>> Vale à pena instalar todos os pacotes que vc usa por ali (pelo menos da >>> primeira vez) pois isto aumentará a probabilidade de que todas as >>> dependências que os pacotes precisam para funcionar também sejam >>> instaladas. >>> >>> Depois é só rodar no R: >>> update.packages(ask=F,dep=T,checkBuilt=T) >>> para que tudo seja atualizado. >>> >>> Por fim, mesmo assim pode acontecer que ao tentar instalar um pacote vc >>> não consiga por falta de dependências. Por exemplo imagine que vc não >>> esteja conseguindo instalar o pacote rgl. Para solucionar esta questão >>> basta rodar no terminal: >>> >>> sudo apt-get build-dep r-cran-rgl >>> >>> Assim todas as dependências serão instaladas. >>> >>> Abs, >>> >>> D >>> >>> 2011/11/25 Pedro Rafael <pedro.rafael.marinho(a)gmail.com> >>> >>>> Ok, que bom que conseguiu Marlene. Qualquer coisa, estamos as ordens. >>>> >>>> Em 24 de novembro de 2011 21:39, marlene marchena [via R-br] < >>>> ml-node+s2285057n4106046h3(a)n4.nabble.com> escreveu: >>>> >>>>> Oi Pedro, >>>>> >>>>> Obrigada pelo roteiro. Graças a teu post finalmente conseguí instalar >>>>> o R 2.14. Espero que o Marcelo tenha conseguido tb. >>>>> >>>>> Um abraço, >>>>> >>>>> Marlene. >>>>> >>>>> >>>>> >>>>> Em 22 de novembro de 2011 23:50, Pedro Rafael <[hidden email]<http://user/SendEmail.jtp?type=node&node=4106046&i=0> >>>>> > escreveu: >>>>> >>>>>> Mas o R está meio antigo, já estamos na versão 2.14. Qual a versão do >>>>>> UbuntU que você está utilizando? Lembrando que é Ubuntu e não Ubunto. Faz >>>>>> assim: >>>>>> >>>>>> sudo su >>>>>> cd /etc/apt/ >>>>>> gedit sources.list >>>>>> >>>>>> O último comando abrirá um arquivo de texto. No final do arquivo >>>>>> coloque: deb http://cran-r.c3sl.ufpr.br/bin/linux/ubuntu oneiric/ >>>>>> Após copiar salve o arquivo e feche. Posteriormente ainda no terminal >>>>>> rode: apt-get update para reconhecer o novo respositorio. >>>>>> >>>>>> Após isso rode no terminal os dois comandos abaixo: >>>>>> >>>>>> gpg --keyserver keyserver.ubuntu.com --recv-key E084DAB9 >>>>>> gpg -a --export E084DAB9 | sudo apt-key add - >>>>>> >>>>>> Novamente rode: apt-get update. Posteriormente remova a versão antiga >>>>>> do R com: apt-get remove r-base. Após a remoção do R instale novamente >>>>>> usando apt-get install r-base e você terá a versão 2.14 do R. Depois >>>>>> instale o pacote do R para ver se dar certo. >>>>>> >>>>>> >>>>>> >>>>>> Em 22 de novembro de 2011 19:40, FHRB Toledo [via R-br] <[hidden >>>>>> email] <http://user/SendEmail.jtp?type=node&node=4106046&i=1>>escreveu: >>>>>> >>>>>> Você deve adicionar o repositório no seu arquivo >>>>>>> /et/apt/sources.list ... veja a explicação em: >>>>>>> http://cran.r-project.org/bin/linux/ubuntu/README >>>>>>> >>>>>>> att, >>>>>>> FH >>>>>>> >>>>>>> 2011/11/22 Marcelo Cardoso mello <[hidden email]<http://user/SendEmail.jtp?type=node&node=4097703&i=0> >>>>>>> > >>>>>>> >>>>>>> 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. >>>>>>>> >>>>>>>> >>>>>>>> _______________________________________________ >>>>>>>> R-br mailing list >>>>>>>> [hidden email]<http://user/SendEmail.jtp?type=node&node=4097703&i=1> >>>>>>>> >>>>>>>> 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. >>>>>>>> >>>>>>> >>>>>>> >>>>>>> _______________________________________________ >>>>>>> R-br mailing list >>>>>>> [hidden email]<http://user/SendEmail.jtp?type=node&node=4097703&i=2> >>>>>>> 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. >>>>>>> >>>>>>> ------------------------------ >>>>>>> If you reply to this email, your message will be added to the >>>>>>> discussion below: >>>>>>> >>>>>>> http://r-br.2285057.n4.nabble.com/R-br-Versao-R-2-14-0-no-Ubunto-tp4097586p… >>>>>>> To unsubscribe from R-br, click here. >>>>>>> NAML<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_view…> >>>>>>> >>>>>> >>>>>> >>>>>> >>>>>> -- >>>>>> Saudações, >>>>>> Pedro Rafael Diniz Marinho. >>>>>> Estatístico - Secretaria de Estado da Saúde - PB. >>>>>> >>>>>> >>>>>> >>>>>> _______________________________________________ >>>>>> R-br mailing list >>>>>> [hidden email] <http://user/SendEmail.jtp?type=node&node=4106046&i=2> >>>>>> 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. >>>>>> >>>>> >>>>> >>>>> _______________________________________________ >>>>> R-br mailing list >>>>> [hidden email] <http://user/SendEmail.jtp?type=node&node=4106046&i=3> >>>>> 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. >>>>> >>>>> ------------------------------ >>>>> If you reply to this email, your message will be added to the >>>>> discussion below: >>>>> >>>>> http://r-br.2285057.n4.nabble.com/R-br-Versao-R-2-14-0-no-Ubunto-tp4097586p… >>>>> To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib…> >>>>> . >>>>> NAML<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_view…> >>>>> >>>> >>>> >>>> >>>> -- >>>> Saudações, >>>> Pedro Rafael Diniz Marinho. >>>> Estatístico - Secretaria de Estado da Saúde - PB. >>>> >>>> >>>> >>>> _______________________________________________ >>>> 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. >>>> >>> >>> >>> _______________________________________________ >>> 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. >>> >> >> >> _______________________________________________ >> 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. >> > > _______________________________________________ > 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. >

2 1

lquadmath
by FHRB Toledo 30 Nov '11

30 Nov '11

Senhores, Uso Debian com R 2.14, desde o upgrade instalar (e/ou atualizar) pacotes com implementações acusa o seguinte "erro /usr/bin/ld: cannot find -lquadmath", vi alguma coisa na internet desse problema com R/Debian e as bibliotecas do gcc! Referência: http://groups.google.com/group/rattle-users/browse_thread/thread/f4d004a094… . Como proceder, pois seguindo o recomendado pela referência nada acontece! att, FH

2 4

Perguntas dentro de function
by Guilherme Biz 30 Nov '11

30 Nov '11

Bom dia, Gostaria de saber como fazer uma pergunta dentro de uma function e assim continuar esta mesma função. Por exemplo: Dentro da função é perguntada o peso da pessoa e com a resposta dada continuar a função. Note que a continuação desta função vai depender da resposta. Att, Guilherme Biz

4 4

setwidth: configuração automática de options("width") em terminal
by Jakson Alves de Aquino 30 Nov '11

30 Nov '11

Olá! Problema: Quando o R está sendo usado em emulador de terminal, a opção "width" não é atualizada automaticamente quando o terminal é redimensionado. Solução: Capturar o sinal SIGWINCH (que avisa ao aplicativo que o terminal foi redimensionado) e modificar o valor de options("width"). Estou lançando a versão para teste do pacote setwidth com a implementação da solução mencionada: http://www.lepem.ufc.br/jaa/setwidth_0.9-0.tar.gz Para testar, basta instalar o pacote R CMD INSTALL setwidth_0.9-0.tar.gz e, no R, digitar library(setwidth) Ao ser carregado o pacote, fica habilitada a captura do sinal SIGWINCH. Não há nenhuma função em R visível para o usuário. Agradeço relatos de bug e comentários. -- Jakson Alves de Aquino Universidade Federal do Ceará Departamento de Ciências Sociais www.lepem.ufc.br/aquino.php

5 11

Uso do ifelse para variáveis qualitativas
by Edson Lira 30 Nov '11

30 Nov '11

Como faço para usar o ifelse em variáveis qualitativas. Na rotina abaixo estou tendo um problema. tip_doa$bairro_resid<-ifelse(tip_doa$BAIRRORESIDENCIA=="NA", tip_doa$BAIRRO_NAO_TABELADO,tip_doa$BAIRRORESIDENCIA) Não indica erro, mas vejam table. > table(tip_doa$bairro_resid) 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 372 557 179 1956 531 108 537 252 60 417 57 861 2 328 512 3870 173 3 123 395 332 2046 995 141 240 580 306 1285 41 161 242 1072 261 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 883 581 254 365 114 54 117 245 256 154 561 479 146 3 2 173 681 92 246 828 140 403 1166 592 109 328 183 12 271 1032 117 490 34 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 193 354 616 364 177 713 443 1204 366 1 309 35 513 216 117 20 1 292 765 > Se eu usar NA ao invés de "NA" me dá a saída: > table(tip_doa$bairro_resid) character(0) [].s Edson Lira Estatístico Manaus-Amazonas

4 3

excluir variável específica de uma tabela de dados
by Rodrigo Plei 30 Nov '11

30 Nov '11

Prezados, Eu tenho uma tabela com 32 variáveis dispostas em colunas, todas com o mesmo número de linhas. Gostaria de saber como posso fazer uma segunda tabela com essa tabela original, mas excluindo somente a primeira variável. Agradeço desde logo pela ajuda. Rodrigo -- =8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8--> Rodrigo Silvestre Martins, PhD Bolsista FAPESP Pós-Doutorado Laboratório de Ecossistemas Pesqueiros (LabPesq) Universidade de São Paulo, Instituto Oceanográfico Praça do Oceanográfico, 191. Cidade Universitaria (sala 107-A/B) Butantã - São Paulo/SP, Brasil 05508-900 Tel: +55 11 3091 6549 Email: rodrigo.plei(a)gmail.com ; ocersm(a)lycos.com; rsmartins(a)usp.br CV Lattes: http://lattes.cnpq.br/5350064124902777

2 2

Dúvida em simulação de processos não homogêneo
by Victor Eduardo 29 Nov '11

29 Nov '11

Pessoal, estou com dúvida no script abaixo, pois as vezes fica me dando como resposta NULL e não estou conseguindo identificar o erro. O meu objetivo é Implementar uma funçao que, para um dado T , retorna os instantes das ocorrências de um processo de Poisson nao-homogeneo com funcao de intensidade lambda(t) = 3/(t + 1) lambda_t<-function(t){ if(t<0){ print("tempo negativo") break() } lambdat<-3/(t+1) return(lambdat) } pois_naohomog<-function(){ lambda<-3 i<-1 T<-0 saida<-NULL T<-T+exponencial(lambda) u<-runif(1) while(u<lambda_t(T)/lambda){ saida[i]<-T i<-i+1 T<-T+exponencial(lambda) } return(saida) }

2 1

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