- R-br - DINF UFPR

Dúvida acessar base de dados
by Marco Pollo 18 Abr '12

18 Abr '12

Caros, Algum direcionamento para eu conseguir acessar um arquivo com extensão .GDB (FireBird Database) direto de uma base de dados. Marco Pollo Almeida

2 1

Re: [R-br] Duvida String
by Sérgio Henrique almeida da silva ju 18 Abr '12

18 Abr '12

---------- Mensagem encaminhada ---------- De: Sérgio Henrique almeida da silva ju <sergio.edfisica(a)gmail.com> Data: 18 de abril de 2012 11:22 Assunto: Re: [R-br] Duvida String Para: r-br(a)listas.c3sl.ufpr.br Uma alternativa que encontrei foi f = function(v, dat, v1=sapply(gregexpr(dat,v),"[",i=1)){ n = ifelse(v1 > 1, nchar(v), NA) substr(v, v1,n) } dados1$causaa = f(dados1$linhaa, "A39") dados1$causab = f(dados1$linhab, "A39") dados1$causac = f(dados1$linhac, "A39") dados1$causad = f(dados1$linhad, "A39") dados1$causaii = f(dados1$linhaii, "A39") dados1$causa<- with(dados1, gsub('NA', '', paste(dados1$causaa,dados1$causab,dados1$causac,dados1$causad,dados1$causaii, sep=''))) Como faço para buscar mais de um código? Tem que fazer um monte de comando como esses? Existe uma forma mais simples de fazer isso? Segue o banco anexo Desculpe enviar o arquivo anexo, mas é uma tentativa para vocês me ajudarem, já que o dput não é uma boa saida! Abraços -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942 -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

2 3

Agregar colunas
by Sérgio Henrique almeida da silva ju 18 Abr '12

18 Abr '12

Olá Estou tentando agregar três colunas em uma, porém quando faço isso os NAs também se agregam, como faço para descosiderá-los Ex: # Dados tmp = data.frame(x=c('A1', NA, 'A2', NA, 'A3', NA, NA), y=c(NA, NA, NA, 'A4', NA, 'A5', NA), z=c(NA, 'A6', NA, NA, NA, NA, 'A7')) # Agregando com o PASTE tmp$W <-paste(tmp$x,tmp$y,tmp$z, sep = "") # Como os dados ficam dput(tmp) structure(list(x = structure(c(1L, NA, 2L, NA, 3L, NA, NA), .Label = c("A1", "A2", "A3"), class = "factor"), y = structure(c(NA, NA, NA, 1L, NA, 2L, NA), .Label = c("A4", "A5"), class = "factor"), z = structure(c(NA, 1L, NA, NA, NA, NA, 2L), .Label = c("A6", "A7"), class = "factor"), W = c("A1NANA", "NANAA6", "A2NANA", "NAA4NA", "A3NANA", "NAA5NA", "NANAA7")), .Names = c("x", "y", "z", "W"), row.names = c(NA, -7L), class = "data.frame") # Gostaria que minha coluna W ficasse assim e não como acima: W = c("A1", "A6", "A2", "A4", "A3", "A5", "A7") Como faço Abraços -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

2 3

Configuração de e-mails - Lista
by Danilo Cuzzuol Pedrini 17 Abr '12

17 Abr '12

Prezados, Gostaria de saber se é possível escolhermos de receber apenas um relatório diário sobre o que foi enviado à Lista? Pelo menos para mim, a caixa de e-mails fica cheia de assuntos os quais não posso ajudar e/ou tenho interesse. A coisa fica pior quando me ausento uns 4 dias do e-mail. Atenciosamente, Danilo Cuzzuol Pedrini www.danilocp.blogspot.com

2 1

Comparar duas regressões no mesmo plot
by Marcelo Claro de Souza 17 Abr '12

17 Abr '12

Plotei a seguinte regressão hipotética e gostaria de saber se existe alguma forma se verificar se as reg 1 e 2 são significativamente diferentes, talvez comparando ângulo de inclinação das retas? Tb gostaria de saber como faço p inserir na legenda o que significa cada reta, pois só consegui colocar o que significam os pontos. x1=c(1,2,3,4,5,6,7,8,9,10) x2=c(13,14,15,16,17,18,19) y1=c(30,40,50,60,70,80,90,100,110,120) y2=c(20,22,23,24,25,26,27) x=c(1,2,3,4,5,6,7,8,9,10,13,14,15,16,17,18,19) y=c(30,40,50,60,70,80,90,100,110,120,20,22,23,24,25,26,27) reg1=lm(x1~y1) reg2=lm(x2~y2) simbolos=c(rep(19,10),rep(2,7)) plot(x~y,pch=simbolos) abline(reg1,lty=2) abline(reg2,lty=1) legend("bottomright",c("D","W"),pch=c(19,2)) Muito obrigado pela ajuda. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

3 2

ajuda com sample
by DANIELA CUNHA DA SE 17 Abr '12

17 Abr '12

Olá Pessoal, Gostaria de saber se alguém pode me ajudar: Quero sortear uma amostra dentro de um data.frame com varias colunas.....porém quero resgatar da amostra selecionada as informações das outras colunas. Exemplo a1<-sample(dados$parc,10) [1] 21 51 59 3 30 65 10 2 13 47 exemplo a primeira unidade selecionada é a parcela 21 (da coluna dados$parc do data.frame), como relacionar a essas parcelas os valores das outras colunas? Obrigada pela atenção.

2 1

Ajuda com mapas
by Renato Braga 17 Abr '12

17 Abr '12

Prezados, Sou iniciante com mapas: vocês podem me encaminhar/indicar material de estudo? Como plotar as bordas de um mapa, definir limites (municípios, por exemplo)? Obrigado,

2 1

modelos mistos - funçao para nested ancova
by Catarina Jakovac 17 Abr '12

17 Abr '12

Olá, Alguém sabe qual a funçao para rodar uma ancova hierárquica? O modelo que quero testar é o crescimento de plantulas (RGR) em funçao do tipo de capoeira (Standtype), espécies de plantulas (spp), tratamento de abertura do dossel (treatment) e tendo como covariáveis a idade da capoeira (age) e os niveis de luz que atingem a plantula (light). Como o tratamento de abertura do dossel é pareado dentro de cada capoeira, preciso incluí-lo com um random effect. Desejo testar também todas as interaçoes possíveis. Pelo que eu estudei até agora, usando o pacote lme4, o modelo ficaria algo parecido com: ancovaRGR<-lme(RGR~Standtype + treat + spp + age + light + Standtype* treat *spp *age + Standtype* treat *spp *light, data = seedlings, random = ~ 1|num/treat/seedlingID) Mas esse modelo não está rodando corretamente. Alguém sabe o que está errado, e se essa é mesmo a função que eu devo usar? Obrigada, Catarina End -- MSc. Ana Catarina C. Jakovac PhD. Candidate Forest Ecology and Forest Management Group Centre for Ecosystem Studies Wageningen University and Research Centre P.O. Box 47 - Lumen Building room 1.212 6700 AA Wageningen, The Netherlands -- MSc. Ana Catarina C. Jakovac PhD. Candidate Forest Ecology and Forest Management Group Centre for Ecosystem Studies Wageningen University and Research Centre P.O. Box 47 - Lumen Building room 1.212 6700 AA Wageningen, The Netherlands -- MSc. Ana Catarina C. Jakovac PhD. Candidate Forest Ecology and Forest Management Group Centre for Ecosystem Studies Wageningen University and Research Centre P.O. Box 47 - Lumen Building room 1.212 6700 AA Wageningen, The Netherlands

4 9

(sem assunto)
by sznelwar 16 Abr '12

16 Abr '12

3 2

Problema com barplot múltiplo
by ASANTOS 15 Abr '12

15 Abr '12

Boa noite pessoal, Tenho 3 tratamentos hipotéticos avaliados em 3 tempos diferentes, porém no tempo dois não foi avaliado o tratamento 2. Para fazer a representação deles em um barplot múltiplo, vai ficar faltando uma barra, portanto vai ficar um espaço em branco entre o tratamento 1 e o tratamento 3 e gostaria de remover esse espaço (por motivo de uma força maior chamada reviewer do periódico) e deixar lado a lado as barras dos tratamentos 1 e 3 no tempo 2. Meu CMR seria: require(gplots) #tratmentos no tempo1 trat1<-rep(1:3,5) resp1<-rpois(15,25) med1<-tapply(resp1, trat1,mean)##Media dos tratamentos ep1<-tapply(resp1,trat1,sd)/sqrt(tapply(resp1,trat1,length))###Erros padrao da media #tratmentos no tempo2 trat2<-c(1,1,1,1,1,2,2,2,2,2,3,3,3,3,3) resp2<-c(7,7,11,12,8,NA,NA,NA,NA,NA,8,8,10,9,6)##Não foi avaliado o tratamento 2 med2<-tapply(resp2, trat2,mean, na.rm = TRUE)##Media dos tratamentos ep2<-tapply(resp2,trat2,sd, na.rm = TRUE)/sqrt(tapply(resp2,trat2,length))###Erros padrao da media #tratmentos no tempo3 trat3<-rep(1:3,5) resp3<-rpois(15,45) med3<-tapply(resp3, trat3,mean)##Media dos tratamentos ep3<-tapply(resp3,trat3,sd)/sqrt(tapply(resp3,trat3,length))###Erros padrao da media ### # Montando as tabelas para cada mêss: a<- med1 b<- med2 cc<-med3 a2<-ep1 b2<-ep2 c2<-ep3 # Unindo as tabelas TAB <- cbind(a,b,cc) TAB2<-cbind(a2,b2,c2) # 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 insects", legend.text=T, ylim=c(0, 60),col =c("grey25","grey50","white")) # Adicionando os meses mtext(1, at = colMeans(mp), text = rep(c("April", "May", "June"), 7),line = 1, cex = 0.8) # Então eu pergunto se é possível fazer isso, pois a função barplot2() não me permite alterar isso uma vez que tenho no tempo 2 o tratamento 2 representado por NA, por outro lado se retiro o tratamento 2 na hora de unir as tabelas não posso fazer um cbind(), pois tenho vetores de média e erro padrão com comprimentos diferentes. Se algum puder me dar uma luz agradeço, -- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304

4 9