- 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