- R-br - DINF UFPR

[Dúvida] Configuração de computador
by Pedro Rafael 04 Jan '12

04 Jan '12

Pessoal um processador intel Xeon E5620 é melhor que o core i7 para trabalhar com contas e grandes base de dados? Quero trabalhar com o R e algumas outras linguagens de programação neste computador. Mais detalhes segue abaixo: Processador: - Quantidade: 1 - Fabricante: Intel - Modelo: Xeon Nehalem E5520 - Freqüência de operação: 2,4GHz~2,66GHz - Cache L2: 12MB (6MB 6MB), compartilhado Memória: - Capacidade: 8GB (2 x 4GB) - Expanção Máxima Opcional: 36GB (9 x 4GB) - Voltagem: 1.5V - Arquitetura: DDR3 - Freqüência: PC3-1333Mhz - Número de pinos: 240 Disco Rígido (HD): - Capacidade: 1TB - Interface: Serial ATA 3.0Gbps - Velocidade: 7200 RPM Unidade Óptica: - Leitor/Gravador DVD-RW Placa Mãe: - Marca: Intel - Modelo: S5500HCV Chipset: - Intel® 5500 chipset with Intel® I/O Controller Hub ICH10R Soluções Termicas: - Intel® Thermal Solutions STS100 series -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

1 0

Adicionado índices do teste de médias no barplot
by ASANTOS 04 Jan '12

04 Jan '12

Boa tarde pessoal, Estou fazendo um barplot com a função barplot2 do pacote gplot, porém eu gostaria de colocar os índices do teste de média do Tukey sobre as barras do erro padrão, para isso fiz: ### require(gplot) # Criando o barplot múltiplo mp <- barplot2(TAB, beside = TRUE, axisnames = FALSE, plot.ci=T, ci.u=TAB+TAB2,ci.l=TAB-TAB2, ylab="Number of Bradysia difformis", legend.text=T) # Adicionando os meses mtext(1, at = colMeans(mp), text = rep(c("April", "May", "June", "July","August", "September","October"), 7),line = 1, cex = 0.75) Porém gostaria de saber se existe algum comando que consigo especificar um vetor contendo os índices e faze-los aparecer sobre as barras do erro padrão, ao invés de especificar índice por índice como text(). 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

5 5

Teste de Dunnett
by Adriano Carvalho Costa 04 Jan '12

04 Jan '12

Boa noite, Venho pedir ajuda de voces para perguntar se alguém conhece alguma rotina para realizar um teste de Dunnet. Pois estou com um experimento fatorial com um tratamento adicional e devo comparar o adicional versus os demais (um a um). Muito obrigado!

3 2

Calcular sobreposição de poligonos de um "SpatialPolygonsDataFrame" do pacote "sp"
by Augusto Ribas 03 Jan '12

03 Jan '12

Ola, primeiro feliz ano novo todos da lista. Eu gostaria de calcular a sobreposição entre poligonos, mas não consigo. Nas verdade eu tenho varios poligos gostaria de uma matriz de sobreposição, como a que demonstro ali no final do script. Nao achei como fazer isso ainda. Procurando por Overlap pelas listas de R não tenho tido muito sucesso então acho que o problema deva ser como estou procurando. Os poligonos são areas de vida de alguns bixinhos e a sobreposição da a informação se eles deixam outro individuos entrarem no territorio deles ou não, então é isso que estou tentando fazer e nao consigo. Segue um CMR pra vcs verem como são os dados e um exemplo de como gostaria do do resultado mais ou menos: #Pacote com dados de exemplo, ele necessita do pacote sp library(adehabitatHR) library(sp) #Dados de exemplo data(puechabonsp) rel<- puechabonsp$relocs #Estimando Area de Vida cp <- mcp(rel[,1]) #Olhando o os poligonos plot(cp) text(coordinates(cp)[,1],coordinates(cp)[,2],label=c("A","B","C","D")) #Como estao os dados class(cp) [1] "SpatialPolygonsDataFrame" attr(,"package") [1] "sp" #O output que eu gostaria # A B C D #A 100 0 0 0 #B 0 100 0 0 #C 0 0 100 X #D 0 0 X 100 -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

3 2

qcc por estrato
by Edson Lira 02 Jan '12

02 Jan '12

Caros amigos, tenho essa estrutura de dados abaixo, aonde os exames são HIV1 HIV2. Exame Total Reativo 141 HIV 1 4270 14 142 HIV 2 4270 18 143 HIV 1 3558 7 144 HIV 2 3558 9 145 HIV 1 3530 4 146 HIV 2 3530 10 147 HIV 1 3241 6 148 HIV 2 3241 9 149 HIV 1 3388 2 150 HIV 2 3388 13 151 HIV 1 4249 1 152 HIV 2 4249 13 153 HIV 1 3763 11 154 HIV 2 3763 30 155 HIV 1 3821 6 156 HIV 2 3821 13 157 HIV 1 3972 0 158 HIV 2 3972 18 159 HIV 1 3227 3 160 HIV 2 3227 16 161 HIV 1 3576 4 162 HIV 2 3576 10 163 HIV 1 2941 4 164 HIV 2 2941 13 Com a rotina abaixo, require(qcc) qcc(h1_05$Reativo,type="np",nsigmas=2,sizes=h1_05$Total) Me dá o gráfico para os dois exames conjuntamente. Como fazer separando HIV1 e HIV2? Edson Lira Estatístico Manaus-Amazonas

2 2

Re: [R-br] qcc por estrato
by Ivan Bezerra Allaman 02 Jan '12

02 Jan '12

?subset (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}

1 0

Dúvida: Erro em shapiro test
by Flávia Machado 31 Dez '11

31 Dez '11

Olá! Estou tendo problemas com o teste shapiro wilk. Inseri o arquivo txt (com pontos ao invés de vírgulas) dos meus dados com o comando "dados<-read.table(file.choose(), header=T). Os dados foram carregados e ao digitar dados todos os objetos apareceram normalmente, mas quando fui fazer o teste de shapiro apareceu " Erro em sort(x[complete.cases(x)]) : objeto 'FRUTO' não encontrado". Para o teste de shapiro, digitei: "shapiro.test(FRUTO)" Estou achando bem estranho porque a pouco tempo fiz o mesmo teste, da mesma forma, com outra série de dados e deu certo. Hoje usei estes mesmos dados apra testar e deu o mesmo erro. Vocês saberiam me dizer como solucionar esse problema? Obs.: sou bem iniciante no R, acho que deu para perceber. Atenciosamente, -- Flávia Machado

2 3

modelo multinomial logistico
by Luiz André Zardo 30 Dez '11

30 Dez '11

Prezados, estou usando a função multinom para ajustar um modelo logistico. Porém não sei como calcular os valores da estatistica da razão de verossimilhança, do pseudo-R2 e da estatistica de Wald. Alguém saberia me indicar como conseguir esses calculos??? ou indicar algum outro pacote que possa disponibilisar essas estatisticas?? Obrigado Luiz Zardo

2 1

aggregate
by Edson Lira 30 Dez '11

30 Dez '11

Caros amigos do R, na rotina abaixo gostaria de calcular além da média, o desvio padrão e o CV=(sd/mean). Como fazer estes cálculos? ## example with character variables and NAs testDF <- data.frame(v1 = c(1,3,5,7,8,3,5,NA,4,5,7,9), v2 = c(11,33,55,77,88,33,55,NA,44,55,77,99) ) by1 <- c("red","blue",1,2,NA,"big",1,2,"red",1,NA,12) by2 <- c("wet","dry",99,95,NA,"damp",95,99,"red",99,NA,NA) aggregate(x = testDF, by = list(by1, by2), FUN = "mean") # and if you want to treat NAs as a group fby1 <- factor(by1, exclude = "") fby2 <- factor(by2, exclude = "") aggregate(x = testDF, by = list(fby1, fby2), FUN = "mean") FELIZ 2012 a todos os amigos e familiares da lista. Edson Lira Estatístico Manaus-Amazonas

2 1

Retirar elemento de uma lista
by Rodrigo Coster 29 Dez '11

29 Dez '11

Povo, como faço para retirar apenas um elemento de uma lista? No caso, eu tenho essa lista: list(fruta='melancia',x=matrix(rpois(100,5),10),y=matrix(rpois(200,7),10)) Eu quero uma lista só com x e y (no caso a minha lista tem um numero variado de 'termos', não apenas 3). É o equivalente a índice negativo em vetor, mas não funciona em list :( []'s

3 3