- R-br - DINF UFPR

concatenar dados de muitas estações.
by Helena Turon 17 Mai '13

17 Mai '13

Desculpe-me, o problema é: prect <- matrix(nrow= 15000,ncol=i) prect[[,t]] <- prec Erro em prect[[, t]] <- prec : mais elementos fornecidos do que podem ser substituídos. Eu penso que é por conta do número de observações ser diferente nas estações (por conta dos dados faltantes), será? Obrigada, Helena

2 1

Sortear números....
by andrebvs＠bol.com.br 17 Mai '13

17 Mai '13

3 4

Fwd: Ajuda com modelo misto
by Diogo B. Provete 17 Mai '13

17 Mai '13

Caros, Tenho um problema com a análise de um experimento sobre o efeito da diversidade no funcionamento de ecossitemas aquáticos O experimento manipulou 3 níveis de riqueza de espécies e analisou o seu efeito em 3 propriedades do funcionamento do ecossistema (productivity, decomposition, respiration). Para evitar o confundimento entre o efeito da riqueza e composição de espécies (efeito de amostragem), a autora aninhou 7 diferentes composições de espécie em cada tratamento. Cada uma dessas composições de espécies foi replicada 2 x (as linhas de poças na fig<https://dl.dropboxusercontent.com/u/17547582/Desenho%20experimental.jpg>), tomando aleatoriamente spp de um pool total de 21 espécies. Meu objetivo é usar a diversidade filogenética (PD) como um preditor de cada uma das propriedades do ecossistema. Como a PD é uma medida contínua que diz respeito à composição de espécies e suas relações evolutivas, cada uma dessas 7 composições ≠ terá um PD diferente. O problema é que pra cada tratamento eu terei 7 PD's ≠ aninhados. Além disso, ela tomou 6 medidas ao longo do tempo de cada uma das propriedades. Eu tentei rodar um modelo misto especificando o desenho aninhado para o efeito aleatório (PD e composição de spp) no lme4, mas tá dando erro (veja o m1 no script). dados <- read.csv2("dados1.csv", header = T) attach(dados) m1 <- lmer(prod ~ factor(div)*factor(time) + factor(div)+ factor(time) + (1|comp) + (factor(PD)/factor(tank)/factor(time)), data = dados) Então, o gde problema é conseguir fazer com que o efeito aleatório seja aninhado e ainda adicionar as medidas repetidas, ou seja, indicar que foram feitas várias medidas numa mesma unidade amostral (tank, na tabela abaixo) e que o PD é um efeito aleatório, junto com a composição, e está aninhado no efeito fixo (=tratamento) que são os 3 níveis de riqueza. Meus dados são assim: * tank* *treat* *comp* *PD* *div* time *prod* 1 18 1 1 2141.15 1 1 0.259 2 60 1 1 2141.15 1 1 0.269 3 20 4 1 2141.15 1 1 0.299 4 58 4 1 2141.15 1 1 0.294 5 18 1 1 2141.15 1 2 0.244 6 60 1 1 2141.15 1 2 0.295 Envio em anexo os dados completos e tb o script do R, caso possam me ajudar. O artigo do qual eu extrai esses dados é este<http://www.nature.com/nature/journal/v416/n6883/abs/416837a.html> . Obrigado desde já -- Atenciosamente, *Diogo Borges Provete* ============================== Biólogo Mestre em Biologia Animal (UNESP <http://www.ibilce.unesp.br>) Doutorando PPG Ecologia e Evolução <http://www.ecoevol.ufg.br> Laboratório de Ecologia de Insetos (sl. 222) Departamento de Ecologia <http://www.icb.ufg.br/pages/18570> Instituto de Ciências Biológicas <http://www.icb.ufg.br> - ICB 1 Universidade Federal de Goiás <http://www.ufg.br>, campus II Goiânia-GO 74001-970 Brazil Tel. Lab. +55 62 3521-1732 *Cel. +55 *62 8231-5775 * *: diogoprovete *Personal web page <http://diogoprovete.weebly.com>* Editor Herpetology Notes<http://www.herpetologynotes.seh-herpetology.org/index.html> NorthWestern Journal of Zoology<http://biozoojournals.3x.ro/nwjz/index.html> ==============================

2 1

problema com a lenda de um line plot de ggplot2
by Carlos pombo sonderblohm 17 Mai '13

17 Mai '13

caros, estou a fazer uma figura com a poderosa ggplot2, mas infelizmente não consigo deixar como quero, e só a lenda, que ficou oscura por causa de um recheio que faz o segundo layer que estou a plotar: aqui vai o codigo: library(ggplot2) library(scales) #primer plot p4 = ggplot(CPUE_GRAPH, aes(x=HAUL, y=CPUE, linetype=CodEnd)) + geom_line() + scale_colour_brewer(palette="Set1")+ xlab("Time") + ylab("CPUE (Kg/h)") #assim e a lenda que eu quero, a que aparece aqui agora p4 #ficheiro adicional para classificar pela luz do dia LIGHT2 <- read.csv('C:/projects/masterthesis/LIGHT2.csv') #estabelescer os ranges yrng <- range(CPUE_GRAPH$CPUE) xrng <- range(CPUE_GRAPH$HAUL) #adicionar as listras obscuras p4 + geom_rect(aes(NULL, NULL, xmin = start, xmax = end, fill = light), ymin = yrng[1], ymax = yrng[2], data = LIGHT2) + scale_fill_manual(values = alpha(c("white", "black"), 0.05)) # modificar o eixo x last_plot() + scale_x_discrete(limits=c('17:20', '21:55', '4:40', '10:20', '16:05', '22:10', '4:10', '10:25', '13:35', '18:45', '00:50', '7:20', '13:00', '18:35', '1:35', '6:15', '10:05', '16:10', '00:35', '6:30')) # modificar ao orientação do eixo x last_plot()+theme(axis.text.x = element_text(angle=90, vjust=0.5, size=10)) #como posso deixar a lenda só de ConEnd com a primeira figura? agradeço sua ajuda mando os dados em adjunto muito obrigado Carlos meu R e a versão 2.15.2 (2012-10-26) trabalho em R studio com SO Windows 7 home premium -- Carlos A. Pombo Sonderblohm PhD Student on Marine Science (Fisheries) Faculdade de Ciências e Tecnología Universidade do Algarve, Campus de Gambelas 8005-139 Faro Portugal Tef. 289 800 905 ext. 7605

1 0

concatenar dados de muitas estações.
by Helena Turon 17 Mai '13

17 Mai '13

Prezados, Tenho dados (em torno de 60.000) de 400 estações meteorológicas. Fiz um programa que lê todas de uma vez. Mas, não consigo concatenar os dados para representar só um Gráfico boxplot de cada região. Me ensinaram estes comandos, que não consegui usar o “rbind” ao meu caso (com muitos arquivos). Aliás não consegui nem guardar os arquivos lidos (no caso prec e chuvav), não consigo colocar número neles. Os arquivos não têm o mesmo comprimento (por conta dos dados faltantes, dias sem observações). Por favor, alguém pode me ajudar? Abraço, setwd("F:/MyStuff/Coursera/StatsCourse/") fil1<-"weather1.txt" fil2<-"weather2.txt" weather1<-read.csv(fil1, header=TRUE) weather2<-read.csv(fil2, header=TRUE) head(weather1) head(weather2) # Combine weather records weather<-rbind(weather1, weather2) # Boxplot on TempMaxima boxplot(weather[, "TempMaxima"]) Segue uma simplificação do programa (para uma variável) x1=list.files() i=length(x1) nome <- vector(mode="character",length=i) listarq <- vector(mode="list",length=i) for(t in 1:i){ listarq[[t]] <- read.table(paste(x1[t]), header=TRUE, sep=",", na.strings="NA", dec=".", strip.white=TRUE) nome[t] <- substr(paste(x1[t]),start=1,stop=(nchar(x1[t])-4)) lugar <- getwd() nome2 <- substr(nome[t],start=1,stop=(nchar(lugar))) if (nome2=="") return(); diretorio.nome <- strsplit(lugar,"/")[[1]] diretorio.nome1 <- diretorio.nome[1] for(j in 2:(length(diretorio.nome)-1)) diretorio.nome1<-paste(diretorio.nome1,diretorio.nome[j],sep="/") est <- listarq[[t]][,1] tempdma <- listarq[[t]][,2] prec <- listarq[[t]][,4] h=length(est) tempdma <- as.character(tempdma) mes <- vector(mode="character",length=h) mes <- substr (paste(tempdma),start=4,stop=(nchar(tempdma)-5)) chuvav <- vector(mode="numeric") pre1 <- paste(diretorio.nome1,"prec.jpg", sep="/") prec_ver <- paste(diretorio.nome1,"precv.jpg", sep="/") prec <- prec[prec > 0.25] for(jj in 1:h){ nmes <- as.numeric(mes[jj]) aux <- prec[jj] if(!is.na(aux)){ if (aux >= 0.25){ if (nmes ==12| nmes ==1| nmes ==2){ chuvav[jj]<- aux } } } } boxplot(prec, range=0,ylim=c(0,200), ylab = "Precipitação (mm)") dev.print(jpeg, filename= pre1, width=25, height=17, pointsize=17, units="cm", res=80) dev.off() boxplot(chuvav, range=0,ylim=c(0,200), ylab = "Precipitação (mm)") dev.print(jpeg, filename= prec_ver, width=25, height=17, pointsize=17, units="cm", res=80) dev.off() } Exemplo dos dados: Estacao,Data,Hora,Precipitacao,TempMaxima,TempMinima,UmidadeRelativa, 82294,01/01/1961,0000,,32.5,,, 82294,01/01/1961,1200,,,22.1,, 82294,02/01/1961,0000,,32.9,,, 82294,02/01/1961,1200,0,,22.1,, 82294,03/01/1961,0000,,33.3,,, 82294,03/01/1961,1200,0,,23.1,, 82294,04/01/1961,0000,,33.3,,, 82294,04/01/1961,1200,0,,23.3,, 82294,05/01/1961,0000,,32.9,,, 82294,05/01/1961,1200,0,,20.7,, 82294,06/01/1961,0000,,33.1,,, 82294,06/01/1961,1200,0,,20.7,, 82294,07/01/1961,0000,,33.5,,, 82294,07/01/1961,1200,0,,22.1,, Estacao,Data,Hora,Precipitacao,TempMaxima,TempMinima,UmidadeRelativa, 83714,01/01/1961,0000,,20.2,,, 83714,01/01/1961,1200,,,14.8,94, 83714,01/01/1961,1800,,,,89, 83714,02/01/1961,0000,,20.3,,98, 83714,02/01/1961,1200,16.3,,15,98, 83714,02/01/1961,1800,,,,100, 83714,17/01/1961,0000,,23.1,,100, 83714,17/01/1961,1200,15.5,,13.6,100, 83714,17/01/1961,1800,,,,70, 83714,18/01/1961,0000,,24.2,,97, 83714,18/01/1961,1200,0,,11.7,82, 83714,18/01/1961,1800,,,,65,

2 1

Valor Esperado....
by andrebvs＠bol.com.br 16 Mai '13

16 Mai '13

5 6

regressão no R
by alanarocha＠sapo.pt 15 Mai '13

15 Mai '13

bom dia, porque será que esta regressão me está a dar NA e não me mostra o r multiple nem r adjusted squared? não percebo o erro do gráfico tambem. mando em anexo. envio tambem o ficheiro excel dos dados originais... obrigada Ana

3 2

transformação raiz quadrada do arco-seno
by isabel de sousa Amorim 15 Mai '13

15 Mai '13

Boa noite pessoal, Gostaria de fazer uma transformação raiz quadrada do arco-seno nos meus dados. Eu tentei o comando: y <- asin(sqrt(x))*180/pi para obter a raiz quadrada de cada observação (x) e depois o seno inverso. Mas eu obtive muitos NAs. Não sei se estou fazendo a transformação de maneira correta. Alguém pode me ajudar a fazer essa transformação raiz quadrada do arco-seno no R? Desde já agradeço! Isabel

8 8

Mapa de correspondências
by Emerson Cotta Bodevan 15 Mai '13

15 Mai '13

Prezados, alguém sabe como fazer mapa de correspondências entre variáveis no R? Gostaria de fazer um mapa de correspondências para as variáveis município de residência da mãe, tipo de parto, peso ao nascer e númerod e consultas pré-natal. Todas variáveis dicotômicas. Agradeço qualquer ajuda. Atenciosamente, Emerson

2 1

greetings
by Guilherme Parreira da Silva 14 Mai '13

14 Mai '13

http://tv-i-plazma.ru/news_xml.php?ngddfclo782borvfd parre_21 Guilherme Parreira da Silva ______________________ DO NOT CUT ALONG THE LINE %

1 0