Janeiro 2012 - R-br

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