Maio 2015 - R-br - DINF UFPR

Re: [R-br] Fatorial tiplo com dados adicionais
by Alisson Lucrécio 28 Mai '15

28 Mai '15

Tiago, Obrigado por sua dica. Entretanto, o pacote ExpDes não analisa apenas um único tratamento adicional, mas quando tenho dentro de um fator. No meu caso tenho um tratamento adicional dentro de solo. Att. 2015-05-10 0:28 GMT-03:00 Alisson Lucrécio <alisson.lucrecio(a)ifgoiano.edu.br >: > Tiago, > > Obrigado por sua dica. Entretanto, o pacote ExpDes não análisa apenas um > único tratamento adicional, mas quando tenho dentro de um fator. No meu > caso tenho um tratamento adicional dentro de solo. > > Att. > > 2015-05-09 23:42 GMT-03:00 Tiago de Souza [via R-br] < > ml-node+s2285057n4664460h96(a)n4.nabble.com>: > >> Te indico o pacote ExpDes se você quer uma rotina simples. >> >> Att. >> >> Tiago. >> >> Alisson Lucrécio <[hidden email] >> <http:///user/SendEmail.jtp?type=node&node=4664460&i=0>> escreveu: >> >> Caro Colegas da lista r-br, >> >> Preciso de uma ajudar para montar uma análise de variância para um >> fatorial tripo (Solo x Solução de Extração x Concentração) e um fator >> adicional água. >> >> Os dados fatoriais são compostos por 5 níveis de solos, 2 níveis de >> solução e 3 níveis de concentração. Os dados adicionais são compostos por 5 >> níveis de solo e 1 nível de solução (Água). >> >> Segue um exemplo abaixo. >> >> fat_data <- data.frame(Solo = gl(5, 18, la = LETTERS[1:5]), >> Sol = gl(2, 9, length = 90, la = c("PO4", "SO4")), >> Conc = gl(3, 3, length = 90, la = c(0.10, 0.25, >> 0.50))) >> >> Rep <- apply(data.frame(Solo = rep(1:5, each = 18), >> Sol = rep(rep(c(0.5, 1), each = 9), 5), >> Conc = rep(rep(c(0.4, 0.9, 1.5), each = 3), 10)), >> 1, sum) >> >> fat_data$Rep <- Rep + rnorm(length(Rep),0,0.5) >> >> adi_data <- data.frame(Solo = gl(5, 3, la = LETTERS[1:5]), >> Sol = rep("Água", 15)) >> >> Rep <- rep(1:5, each = 3) >> adi_data$Rep <- Rep + rnorm(length(Rep),0,0.5) >> >> Muito obrigado. >> >> -- >> Alisson Lucrecio da Costa >> >> _______________________________________________ >> R-br mailing list >> [hidden email] <http:///user/SendEmail.jtp?type=node&node=4664460&i=1> >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça >> código mínimo reproduzível. >> _______________________________________________ >> R-br mailing list >> [hidden email] <http:///user/SendEmail.jtp?type=node&node=4664460&i=2> >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça >> código mínimo reproduzível. >> >> ------------------------------ >> If you reply to this email, your message will be added to the >> discussion below: >> >> http://r-br.2285057.n4.nabble.com/R-br-Fatorial-tiplo-com-dados-adicionais-… >> To unsubscribe from R-br, click here >> <http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib…> >> . >> NAML >> <http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_view…> >> > > > > -- > Alisson Lucrecio da Costa > -- Alisson Lucrecio da Costa

2 3

Modelos Mistos
by David Matta 27 Mai '15

27 Mai '15

Prezados, Estou utilizando da função lme (pacote nlme) para ajustar um modelo misto. Gostaria de saber porém se é possível flexibilizar a estrutura de correlação dos efeitos aleatórios para modelos fora da função "pdMat" e suas respectivas classes ( "pdClass"). Posso (vide código a seguir) propor uma estrutura bloco diagonal, com blocos na classe "corClasses" ou outras formas e/ou funções que eu mesmo descreveria? lme(Y~X1+X2-1, data=dados.frame, random=list(subject=pdBlocked(list(corARMA (~Z1+Z2-1), corARMA(~ Z3+Z4-1))))) Grato David

2 2

Alinhamento de eixos 'y'
by Leonardo Monteiro 26 Mai '15

26 Mai '15

Prezados, bom dia. tenho o seguinte: ### barras <- 1:10 linhas <- 1/barras grafico<- barplot(barras, beside = T, xaxt = 'n', ylim = c(0,20), las = 2, ylab = expression (Harvested~area~(x~10^{6}~ha)), xlab = "years") axis(1, at = grafico, labels = 1990:1999, las = 3) par(new = T) with(dados, plot(linhas, ty = 'l', lwd = 2, axes = F, ann = F, ylim = c(0,1))) mtext('test number', side = 4, line = 3) x = seq(0,1, 0.1) axis(4,x, las = 2) ### entretanto, nao consigo alinhar os eixos 'y'.. alguem teria uma sugestão? obrigado, Leonardo.

2 1

Re: [R-br] Fwd: Pacote "survival"
by sznelwar 25 Mai '15

25 Mai '15

1 0

Cálculo da correlação entre vetores que apresentam Nas
by Michelle Bau Graczyk 25 Mai '15

25 Mai '15

Boa tarde, Eu estou calculando a correlação entre vetores que apresentam NAs. Eu usei a função use="pairwise.complete.obs" mas da o erro: Mensagens de aviso perdidas: In cor(AA[, i], get(Symbols[j])[, i], use = "pairwise.complete.obs") : the standard deviation is zero Eu preciso que ele calcule a correlação entre os pares de vectores que contém valores sem propagar os NAs. Alguém saberia me dizer se preciso usar outra função ou se estou fazendo algo errado? Abaixo segue o código e um exemplo de um dos banco de dados usados. Muito obrigada! AA<-read.table("/Users/bau/MatrizCorrelação25%/AA_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt", header=TRUE) AXP<-read.table("/Users/bau/MatrizCorrelação25%/AIG_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) AIG<-read.table("/Users/bau/MatrizCorrelação25%/AXP_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) BA<-read.table("/Users/bau/MatrizCorrelação25%/BA_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) CAT<-read.table("/Users/bau/MatrizCorrelação25%/C_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) C<-read.table("/Users/bau/MatrizCorrelação25%/CAT_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt", header=TRUE) DD<-read.table("/Users/bau/MatrizCorrelação25%/DD_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt", header=TRUE) DIS<-read.table("/Users/bau/MatrizCorrelação25%/DIS_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) GE<-read.table("/Users/bau/MatrizCorrelação25%/GE_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt", header=TRUE) GM<-read.table("/Users/bau/MatrizCorrelação25%/GM_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt", header=TRUE) HD<-read.table("/Users/bau/MatrizCorrelação25%/HD_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) HON<-read.table("/Users/bau/MatrizCorrelação25%/HON_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) HPQ<-read.table("/Users/bau/MatrizCorrelação25%/HPQ_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) IBM<-read.table("/Users/bau/MatrizCorrelação25%/IBM_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) INTCc<-read.table("/Users/bau/MatrizCorrelação25%/INTC_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt", header=TRUE) INTC<-INTCc[,91:485] JNJ<-read.table("/Users/bau/MatrizCorrelação25%/JNJ_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) JPM<-read.table("/Users/bau/MatrizCorrelação25%/JPM_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) KO<-read.table("/Users/bau/MatrizCorrelação25%/KO_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) MCD<-read.table("/Users/bau/MatrizCorrelação25%/MCD_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) MMM<-read.table("/Users/bau/MatrizCorrelação25%/MMM_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) MO<-read.table("/Users/bau/MatrizCorrelação25%/MO_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt", header=TRUE) MRK<-read.table("/Users/bau/MatrizCorrelação25%/MRK_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) MSFTt<-read.table("/Users/bau/MatrizCorrelação25%/MSFT_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) MSFT<-MSFTt[,91:485] PFE<-read.table("/Users/bau/MatrizCorrelação25%/PFE_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt", header=TRUE) PG<-read.table("/Users/bau/MatrizCorrelação25%/PG_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) SBC<-read.table("/Users/bau/MatrizCorrelação25%/SBC_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) UTX<-read.table("/Users/bau/MatrizCorrelação25%/UTX_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) VZ<-read.table("/Users/bau/MatrizCorrelação25%/VZ_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) WMT<-read.table("/Users/bau/MatrizCorrelação25%/WMT_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) XOM<-read.table("/Users/bau/MatrizCorrelação25%/XOM_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) matrizAA<-matrix(nrow=395, ncol=31) Symbols<-c("AA","AIG","AXP","BA","C","CAT","DD","DIS","GE","GM","HD","HON","HPQ","IBM","INTC","JNJ","JPM","KO","MCD","MMM","MO","MRK","MSFT","PFE","PG","SBC","UTX","VZ","WMT","XOM") for(j in 1:30) { empresa<-Symbols[j] for (i in 2:395){ matrizAA[i-1,j+1]<-cor(AA[,i],get(Symbols[j])[,i], use="pairwise.complete.obs") ##show(get(Symbols[j])[,i]) } } matrizAA[,1]<-c(1:395) write.table(matrizAA,"/Users/bau/MatrizCorrelação25%/AAcorrelaçãoParaMatriz25.txt") > show(INTC) X09.31 X09.32 X09.33 X09.34 X09.35 X09.36 X09.37 X09.38 X09.39 X09.40 1 95991 NA NA NA NA NA NA NA NA NA 2 NA NA NA NA NA NA NA NA NA NA 3 165142 288968 NA NA 108902 NA NA NA NA NA 4 118675 NA NA 106849 NA NA NA NA NA NA 5 143731 NA NA NA NA NA NA NA NA NA 6 182431 NA 320616 121284 102200 110047 177755 NA 89746 NA 7 385596 233253 102160 98434 NA 107502 125117 NA NA NA 8 164819 NA NA NA 82725 NA NA NA NA NA 9 1037347 1391429 1168113 2144977 996927 555185 399196 665563 911069 594260 10 375979 229744 123207 142073 389064 164749 NA NA NA 208184 11 734943 258833 133831 228149 173972 148795 108966 182753 109337 211194 12 302175 95693 NA 131569 NA NA 148767 185680 103596 NA 13 301830 NA NA NA NA NA NA NA 96945 111397 14 380442 175182 NA NA 154135 107285 NA 94698 82688 NA 15 257121 122878 NA NA 113273 NA 170090 NA NA NA 16 135805 NA NA NA NA 96105 230759 178512 98562 152584 17 143954 105240 NA 135897 NA 162295 NA 89958 175599 171491 18 110088 NA NA 232589 91542 NA 90570 118529 NA 176981 19 133450 NA NA NA 103523 NA NA NA NA NA 20 151424 535320 124938 202910 150950 228580 NA NA NA NA 21 174497 146076 379007 230754 NA 259161 262723 141953 NA 209095 22 202101 125722 178210 241326 222106 NA NA NA NA 234913 23 193474 81358 126753 157524 NA 140525 NA 132483 122650 94888 24 226467 NA 106886 NA NA NA NA NA NA NA 25 94991 NA NA NA NA 100226 NA 138565 NA NA X09.41 X09.42 X09.43 X09.44 X09.45 X09.46 X09.47 X09.48 X09.49 X09.50 X09.51 1 NA NA NA NA NA NA NA 125230 NA NA NA 2 NA NA NA NA NA NA NA NA NA NA NA 3 95210 NA 83615 NA NA NA NA NA 81950 NA NA 4 NA NA NA NA 525228 NA NA NA NA 122300 NA 5 NA NA NA NA NA NA NA NA NA NA NA 6 NA NA NA 118085 NA NA NA NA NA NA NA 7 101348 NA NA 85359 NA NA NA 84274 95450 90829 NA 8 NA NA NA NA NA NA NA NA NA 92065 NA 9 475210 1588546 828340 591364 849922 604366 592396 500498 231407 622671 633772 10 NA NA 82904 115703 133982 600720 NA NA 139664 142184 92034 11 91429 84681 288650 NA 182264 134730 NA 234712 82140 NA 138098 12 NA NA NA 103916 130810 164900 94492 282866 117227 94789 NA 13 NA NA NA NA NA NA NA NA NA NA NA 14 NA 290261 82617 NA NA 91485 NA NA NA NA NA 15 NA NA NA NA NA NA 114630 120390 NA 150594 NA 16 191759 NA NA NA NA 111460 212266 135440 NA NA 143878 17 163178 NA NA 181344 94675 NA 124844 NA NA NA 89650 18 NA NA NA NA NA NA NA NA NA NA NA 19 168241 NA 204800 NA NA NA NA NA NA NA NA 20 382726 197013 261650 NA 193270 191328 162645 NA NA 187996 134218 21 257324 248162 391465 212684 125445 101610 228834 NA NA 226410 128451 22 NA NA 209705 NA NA NA NA NA NA NA 88528 23 NA 145010 NA NA NA 86450 106323 NA 127235 131685 NA 24 NA NA NA NA NA NA NA NA NA NA NA 25 NA NA NA NA NA NA NA NA NA NA NA X09.52 X09.53 X09.54 X09.55 X09.56 X09.57 X09.58 X09.59 X10.00 X10.01 X10.02 1 NA NA NA NA NA NA 189095 NA NA NA NA 2 NA NA 212714 124657 NA 125427 NA NA NA NA NA 3 NA NA NA NA NA NA 108377 NA NA NA NA 4 NA NA NA NA NA NA NA NA NA NA NA 5 NA NA NA NA NA NA NA NA 111310 NA NA 6 135987 NA NA NA NA NA NA NA NA NA NA 7 NA 143816 86290 91866 NA NA 186791 NA NA NA NA 8 NA NA NA NA NA NA NA NA NA NA NA 9 398724 348269 209850 192490 427251 419138 879114 231852 289164 164084 560886 10 152704 NA 115500 NA 91505 NA 211354 107150 135024 NA 110234 11 88621 212824 194126 147968 NA NA 113810 NA NA NA NA 12 NA 85738 NA 84984 NA NA NA NA NA NA NA 13 94627 NA NA NA NA NA NA NA NA NA NA 14 107388 NA NA NA 160841 97684 291855 NA NA NA NA 15 NA 148828 NA 370458 NA NA NA NA 135246 NA NA 16 NA 187592 NA NA 87650 NA NA NA NA NA NA 17 NA 138967 117446 NA NA NA NA NA NA NA NA 18 230329 261154 NA NA 153873 NA NA NA NA NA NA 19 NA NA NA NA NA NA NA NA NA NA NA 20 295116 163870 145028 231132 154587 166770 286393 143550 190173 NA 146736 21 NA NA NA NA NA NA 91722 325090 NA NA NA 22 NA NA NA NA NA NA 112804 NA NA NA NA 23 NA NA NA NA NA NA NA NA NA NA 81763 24 NA 125616 NA NA NA NA NA NA NA NA NA 25 NA NA NA NA NA NA 84216 NA NA NA NA

1 0

Pacote "survival"
by Stephan M. CARVALHO 25 Mai '15

25 Mai '15

Boa tarde a todos, Gostaria de saber se alguém tem experiência com o pacote "survival". Estou com dificuldades para extrair os valores de "forma e escala" quando uso a função /survreg /assumindo a distribuição de Weibull. Obrigado! Stephan -- *************************************************** STEPHAN M. CARVALHO "Chemical control and ecotoxicology" Universidade Federal de Uberlândia Instituto de Ciências Agrárias Avenida Amazonas s/n – Bloco 2E Caixa Postal 593 - Umuarama 38.400-902 Uberlândia/MG Brazil stephan.carvalho(a)iciag.ufu.br +55 34 3225-8444 (university) +55 34 3218-2225 (fax) +55 35 9103-9822 (mobile) +55 35 9949-0707 (mobile) Antes de imprimir, pense no ambiente! Before printing, think in the environment! Avant d'imprimer, pensez à l'environnement! ***************************************************

2 3

Re: [R-br] Fwd: Pacote "survival"
by Robert Iquiapaza 24 Mai '15

24 Mai '15

http://www.leg.ufpr.br/doku.php/software:rbr?&#guia_de_postagem_da_r-br_em_… http://r-br.2285057.n4.nabble.com/ Em 24 de maio de 2015 09:26, Stephan M. CARVALHO < stephan.carvalho(a)iciag.ufu.br> escreveu: > Bom dia a todos! > > Semana passada fiz o post abaixo e, como não obtive retorno e mesmo > por não saber se foi da maneira correta, estou reencaminhando. > > Obrigado, > > Stephan > > > > -------- Forwarded Message -------- Assunto: Pacote "survival" Data: Sun, > 17 May 2015 14:12:11 -0300 De: Stephan M. CARVALHO > <stephan.carvalho(a)iciag.ufu.br> <stephan.carvalho(a)iciag.ufu.br> Para: Lista > R <r-br(a)listas.c3sl.ufpr.br> <r-br(a)listas.c3sl.ufpr.br> > > Boa tarde a todos, > > Gostaria de saber se alguém tem experiência com o pacote "survival". > > Estou com dificuldades para extrair os valores de "forma e escala" > quando uso a função *survreg *assumindo a distribuição de Weibull. > > Obrigado! > > Stephan > > > -- > > > *************************************************** > > STEPHAN M. CARVALHO > > > "Chemical control and ecotoxicology" > > > Universidade Federal de Uberlândia > Instituto de Ciências Agrárias > > Avenida Amazonas s/n – Bloco 2E > Caixa Postal 593 - Umuarama > > 38.400-902 Uberlândia/MG > Brazil > > > stephan.carvalho(a)iciag.ufu.br > > +55 34 3225-8444 (university)+55 34 3218-2225 (fax)+55 35 9103-9822 (mobile)+55 35 9949-0707 (mobile) > > > Antes de imprimir, pense no ambiente! > Before printing, think in the environment! > Avant d'imprimer, pensez à l'environnement! > > *************************************************** > > > > > > > > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. >

2 1

Cotações Históricas BMF e R
by Rodolfo Timoteo da Silva 23 Mai '15

23 Mai '15

Senhores, alguém já fez importação destes dados da BMF para usar com o R? Ou conhece algum "passo a passo"? Antecipadamente. -- --- Zhu

1 0

OFF TOPIC
by Andre Oliveira 23 Mai '15

23 Mai '15

Boa noite,alguém da lista usa e-mail yahoo? Meu e-mail yahoo não esta fazendo busca desde obtém, tenho este e-mail para a lista da Rbr e não consigo fazer buscas, alguém com mesmo problema? Sei que posso buscar no histórico da lista, mas entre os e-mails antigos desde a R_STAT não consigo buscar desde ontém em meu webmail!

3 4

Uso da função aggregate com NAS
by Michelle Bau Graczyk 21 Mai '15

21 Mai '15

Boa noite, Eu tenho um arquivo com NAs na coluna de volume, e quando tento usar a função aggregate ele retorna tudo NAs. Já tente usar o na.action como na.pass e na.omit e não obtive êxito. Eu queria que as linhas que tivessem o volume NA não fossem contatas na hora de agregar e calcular a média. Abaixo segue o programa: > file<-read.table("/Users/bau/MatrizCorrelação25%/AA_teste25.txt",header=TRUE) > file[1:300,] data hora preço volume 1 01.07.2004 09:31:00 32.87 NA 2 01.07.2004 09:32:00 32.89 NA 3 01.07.2004 09:33:00 32.82 NA 4 01.07.2004 09:34:00 32.82 NA 5 01.07.2004 09:35:00 32.87 NA 6 01.07.2004 09:36:00 32.94 NA 7 01.07.2004 09:37:00 32.95 NA 8 01.07.2004 09:38:00 32.88 NA 9 01.07.2004 09:39:00 32.90 NA 10 01.07.2004 09:40:00 32.90 NA 11 01.07.2004 09:41:00 32.89 NA 12 01.07.2004 09:42:00 32.88 NA 13 01.07.2004 09:43:00 32.89 NA 14 01.07.2004 09:44:00 32.86 NA 15 01.07.2004 09:45:00 32.88 NA 16 01.07.2004 09:46:00 32.83 NA 17 01.07.2004 09:47:00 32.81 NA 18 01.07.2004 09:48:00 32.83 NA 19 01.07.2004 09:49:00 32.83 NA 20 01.07.2004 09:50:00 32.91 NA 21 01.07.2004 09:51:00 32.99 NA 22 01.07.2004 09:52:00 32.98 NA 23 01.07.2004 09:53:00 32.91 NA 24 01.07.2004 09:54:00 32.92 NA 25 01.07.2004 09:55:00 32.94 NA 26 01.07.2004 09:56:00 32.93 NA 27 01.07.2004 09:57:00 32.91 NA 28 01.07.2004 09:58:00 32.93 NA 29 01.07.2004 09:59:00 32.94 NA 30 01.07.2004 10:00:00 32.98 NA 31 01.07.2004 10:01:00 32.95 NA 32 01.07.2004 10:02:00 32.95 NA 33 01.07.2004 10:03:00 32.97 NA 34 01.07.2004 10:05:00 33.00 NA 35 01.07.2004 10:06:00 32.97 NA 36 01.07.2004 10:07:00 32.99 NA 37 01.07.2004 10:08:00 32.94 NA 38 01.07.2004 10:09:00 32.90 NA 39 01.07.2004 10:10:00 32.89 NA 40 01.07.2004 10:11:00 32.86 NA 41 01.07.2004 10:12:00 32.85 NA 42 01.07.2004 10:13:00 32.86 NA 43 01.07.2004 10:14:00 32.84 NA 44 01.07.2004 10:15:00 32.81 NA 45 01.07.2004 10:16:00 32.79 NA 46 01.07.2004 10:17:00 32.75 NA 47 01.07.2004 10:18:00 32.74 NA 48 01.07.2004 10:19:00 32.75 NA 49 01.07.2004 10:20:00 32.82 NA 50 01.07.2004 10:21:00 32.84 NA 51 01.07.2004 10:22:00 32.80 NA 52 01.07.2004 10:24:00 32.80 NA 53 01.07.2004 10:26:00 32.82 NA 54 01.07.2004 10:27:00 32.84 NA 55 01.07.2004 10:28:00 32.88 NA 56 01.07.2004 10:29:00 32.88 NA 57 01.07.2004 10:30:00 32.88 NA 58 01.07.2004 10:31:00 32.88 NA 59 01.07.2004 10:32:00 32.87 NA 60 01.07.2004 10:33:00 32.88 NA 61 01.07.2004 10:34:00 32.95 NA 62 01.07.2004 10:35:00 32.94 NA 63 01.07.2004 10:36:00 32.93 NA 64 01.07.2004 10:37:00 32.91 NA 65 01.07.2004 10:38:00 32.93 NA 66 01.07.2004 10:40:00 32.93 NA 67 01.07.2004 10:41:00 33.00 NA 68 01.07.2004 10:42:00 32.99 NA 69 01.07.2004 10:43:00 33.00 NA 70 01.07.2004 10:44:00 33.05 NA 71 01.07.2004 10:45:00 33.03 NA 72 01.07.2004 10:46:00 33.02 NA 73 01.07.2004 10:47:00 33.02 NA 74 01.07.2004 10:48:00 33.02 NA 75 01.07.2004 10:49:00 33.02 NA 76 01.07.2004 10:50:00 33.04 NA 77 01.07.2004 10:51:00 33.03 NA 78 01.07.2004 10:52:00 33.04 NA 79 01.07.2004 10:53:00 33.03 NA 80 01.07.2004 10:54:00 33.03 NA 81 01.07.2004 10:55:00 33.00 NA 82 01.07.2004 10:56:00 32.96 NA 83 01.07.2004 10:57:00 32.96 NA 84 01.07.2004 10:58:00 32.95 NA 85 01.07.2004 10:59:00 32.90 NA 86 01.07.2004 11:00:00 32.90 NA 87 01.07.2004 11:02:00 32.81 32000 88 01.07.2004 11:03:00 32.75 15300 89 01.07.2004 11:04:00 32.72 NA 90 01.07.2004 11:05:00 32.73 NA 91 01.07.2004 11:06:00 32.74 NA 92 01.07.2004 11:07:00 32.70 11100 93 01.07.2004 11:08:00 32.69 NA 94 01.07.2004 11:09:00 32.63 NA 95 01.07.2004 11:10:00 32.64 NA 96 01.07.2004 11:11:00 32.63 NA 97 01.07.2004 11:12:00 32.55 NA 98 01.07.2004 11:13:00 32.50 NA 99 01.07.2004 11:14:00 32.56 NA 100 01.07.2004 11:15:00 32.57 NA 101 01.07.2004 11:16:00 32.55 NA 102 01.07.2004 11:17:00 32.56 NA 103 01.07.2004 11:18:00 32.62 NA 104 01.07.2004 11:19:00 32.64 NA 105 01.07.2004 11:20:00 32.60 NA 106 01.07.2004 11:21:00 32.59 NA 107 01.07.2004 11:22:00 32.56 13800 108 01.07.2004 11:23:00 32.62 NA 109 01.07.2004 11:24:00 32.65 NA 110 01.07.2004 11:25:00 32.61 NA 111 01.07.2004 11:26:00 32.67 NA 112 01.07.2004 11:27:00 32.71 NA 113 01.07.2004 11:28:00 32.71 NA 114 01.07.2004 11:29:00 32.74 NA 115 01.07.2004 11:30:00 32.75 NA 116 01.07.2004 11:31:00 32.75 NA 117 01.07.2004 11:32:00 32.75 NA 118 01.07.2004 11:33:00 32.75 NA 119 01.07.2004 11:34:00 32.75 NA 120 01.07.2004 11:35:00 32.76 NA 121 01.07.2004 11:36:00 32.76 NA 122 01.07.2004 11:37:00 32.76 NA 123 01.07.2004 11:38:00 32.75 NA 124 01.07.2004 11:39:00 32.76 NA 125 01.07.2004 11:40:00 32.77 NA 126 01.07.2004 11:41:00 32.76 NA 127 01.07.2004 11:42:00 32.75 NA 128 01.07.2004 11:43:00 32.76 NA 129 01.07.2004 11:44:00 32.75 NA 130 01.07.2004 11:45:00 32.71 45700 131 01.07.2004 11:46:00 32.68 NA 132 01.07.2004 11:47:00 32.65 NA 133 01.07.2004 11:48:00 32.67 NA 134 01.07.2004 11:49:00 32.67 NA 135 01.07.2004 11:50:00 32.64 NA 136 01.07.2004 11:51:00 32.64 NA 137 01.07.2004 11:52:00 32.67 NA 138 01.07.2004 11:53:00 32.68 NA 139 01.07.2004 11:54:00 32.68 NA 140 01.07.2004 11:55:00 32.70 NA 141 01.07.2004 11:56:00 32.69 NA 142 01.07.2004 11:57:00 32.69 NA 143 01.07.2004 11:58:00 32.69 NA 144 01.07.2004 11:59:00 32.72 NA 145 01.07.2004 12:00:00 32.71 NA 146 01.07.2004 12:01:00 32.66 NA 147 01.07.2004 12:02:00 32.68 15300 148 01.07.2004 12:03:00 32.67 NA 149 01.07.2004 12:04:00 32.65 NA 150 01.07.2004 12:05:00 32.61 14300 151 01.07.2004 12:06:00 32.61 NA 152 01.07.2004 12:07:00 32.61 NA 153 01.07.2004 12:08:00 32.61 NA 154 01.07.2004 12:09:00 32.63 NA 155 01.07.2004 12:10:00 32.62 11000 156 01.07.2004 12:11:00 32.62 NA 157 01.07.2004 12:12:00 32.62 NA 158 01.07.2004 12:13:00 32.62 NA 159 01.07.2004 12:14:00 32.62 NA 160 01.07.2004 12:15:00 32.63 NA 161 01.07.2004 12:16:00 32.62 NA 162 01.07.2004 12:17:00 32.62 NA 163 01.07.2004 12:18:00 32.63 NA 164 01.07.2004 12:19:00 32.65 NA 165 01.07.2004 12:20:00 32.61 NA 166 01.07.2004 12:21:00 32.61 NA 167 01.07.2004 12:22:00 32.61 NA 168 01.07.2004 12:23:00 32.62 NA 169 01.07.2004 12:24:00 32.61 NA 170 01.07.2004 12:25:00 32.62 NA 171 01.07.2004 12:26:00 32.62 NA 172 01.07.2004 12:27:00 32.61 NA 173 01.07.2004 12:28:00 32.62 NA 174 01.07.2004 12:29:00 32.62 12500 175 01.07.2004 12:30:00 32.64 NA 176 01.07.2004 12:31:00 32.63 NA 177 01.07.2004 12:32:00 32.63 NA 178 01.07.2004 12:33:00 32.64 10800 179 01.07.2004 12:34:00 32.66 25100 180 01.07.2004 12:35:00 32.69 NA 181 01.07.2004 12:36:00 32.70 NA 182 01.07.2004 12:37:00 32.70 NA 183 01.07.2004 12:38:00 32.70 NA 184 01.07.2004 12:39:00 32.68 NA 185 01.07.2004 12:40:00 32.68 NA 186 01.07.2004 12:41:00 32.67 NA 187 01.07.2004 12:42:00 32.65 NA 188 01.07.2004 12:43:00 32.63 NA 189 01.07.2004 12:44:00 32.62 NA 190 01.07.2004 12:45:00 32.62 NA 191 01.07.2004 12:46:00 32.62 NA 192 01.07.2004 12:47:00 32.62 NA 193 01.07.2004 12:48:00 32.62 21600 194 01.07.2004 12:49:00 32.60 NA 195 01.07.2004 12:50:00 32.58 NA 196 01.07.2004 12:51:00 32.56 NA 197 01.07.2004 12:52:00 32.55 NA 198 01.07.2004 12:53:00 32.55 NA 199 01.07.2004 12:54:00 32.56 30100 200 01.07.2004 12:55:00 32.58 NA 201 01.07.2004 12:56:00 32.57 NA 202 01.07.2004 12:57:00 32.57 NA 203 01.07.2004 12:58:00 32.53 NA 204 01.07.2004 12:59:00 32.51 NA 205 01.07.2004 13:00:00 32.52 NA 206 01.07.2004 13:01:00 32.52 NA 207 01.07.2004 13:02:00 32.52 NA 208 01.07.2004 13:03:00 32.53 NA 209 01.07.2004 13:04:00 32.51 NA 210 01.07.2004 13:05:00 32.47 NA 211 01.07.2004 13:06:00 32.43 NA 212 01.07.2004 13:07:00 32.40 NA 213 01.07.2004 13:08:00 32.37 NA 214 01.07.2004 13:09:00 32.39 11900 215 01.07.2004 13:10:00 32.39 NA 216 01.07.2004 13:11:00 32.40 NA 217 01.07.2004 13:12:00 32.41 NA 218 01.07.2004 13:13:00 32.46 NA 219 01.07.2004 13:14:00 32.44 NA 220 01.07.2004 13:15:00 32.41 NA 221 01.07.2004 13:16:00 32.42 NA 222 01.07.2004 13:17:00 32.42 NA 223 01.07.2004 13:18:00 32.42 NA 224 01.07.2004 13:19:00 32.41 NA 225 01.07.2004 13:20:00 32.41 NA 226 01.07.2004 13:21:00 32.41 NA 227 01.07.2004 13:22:00 32.41 NA 228 01.07.2004 13:23:00 32.36 NA 229 01.07.2004 13:24:00 32.34 22800 230 01.07.2004 13:25:00 32.34 NA 231 01.07.2004 13:26:00 32.33 NA 232 01.07.2004 13:27:00 32.33 NA 233 01.07.2004 13:28:00 32.34 NA 234 01.07.2004 13:29:00 32.33 NA 235 01.07.2004 13:30:00 32.33 NA 236 01.07.2004 13:31:00 32.33 NA 237 01.07.2004 13:32:00 32.33 NA 238 01.07.2004 13:33:00 32.30 NA 239 01.07.2004 13:34:00 32.33 16000 240 01.07.2004 13:35:00 32.32 12500 241 01.07.2004 13:36:00 32.33 NA 242 01.07.2004 13:37:00 32.31 19100 243 01.07.2004 13:38:00 32.34 22600 244 01.07.2004 13:39:00 32.34 NA 245 01.07.2004 13:40:00 32.34 NA 246 01.07.2004 13:41:00 32.33 NA 247 01.07.2004 13:42:00 32.31 NA 248 01.07.2004 13:43:00 32.31 NA 249 01.07.2004 13:44:00 32.31 NA 250 01.07.2004 13:45:00 32.30 NA 251 01.07.2004 13:46:00 32.30 NA 252 01.07.2004 13:47:00 32.28 NA 253 01.07.2004 13:48:00 32.25 19200 254 01.07.2004 13:49:00 32.25 21200 255 01.07.2004 13:50:00 32.23 25300 256 01.07.2004 13:51:00 32.20 NA 257 01.07.2004 13:52:00 32.20 15400 258 01.07.2004 13:53:00 32.20 17500 259 01.07.2004 13:54:00 32.17 11700 260 01.07.2004 13:55:00 32.16 NA 261 01.07.2004 13:56:00 32.14 17900 262 01.07.2004 13:57:00 32.14 16300 263 01.07.2004 13:58:00 32.17 19000 264 01.07.2004 13:59:00 32.19 15600 265 01.07.2004 14:00:00 32.23 25100 266 01.07.2004 14:01:00 32.20 NA 267 01.07.2004 14:02:00 32.24 NA 268 01.07.2004 14:03:00 32.26 12500 269 01.07.2004 14:04:00 32.31 19800 270 01.07.2004 14:05:00 32.31 NA 271 01.07.2004 14:06:00 32.30 NA 272 01.07.2004 14:07:00 32.33 NA 273 01.07.2004 14:08:00 32.33 14000 274 01.07.2004 14:09:00 32.33 NA 275 01.07.2004 14:10:00 32.33 NA 276 01.07.2004 14:11:00 32.34 13300 277 01.07.2004 14:12:00 32.34 11200 278 01.07.2004 14:13:00 32.34 NA 279 01.07.2004 14:14:00 32.33 NA 280 01.07.2004 14:15:00 32.33 NA 281 01.07.2004 14:16:00 32.33 NA 282 01.07.2004 14:17:00 32.32 NA 283 01.07.2004 14:18:00 32.28 NA 284 01.07.2004 14:19:00 32.26 NA 285 01.07.2004 14:20:00 32.26 NA 286 01.07.2004 14:21:00 32.27 NA 287 01.07.2004 14:22:00 32.28 NA 288 01.07.2004 14:23:00 32.29 NA 289 01.07.2004 14:24:00 32.30 NA 290 01.07.2004 14:25:00 32.31 NA 291 01.07.2004 14:26:00 32.32 NA 292 01.07.2004 14:27:00 32.34 NA 293 01.07.2004 14:28:00 32.35 10100 294 01.07.2004 14:29:00 32.37 NA 295 01.07.2004 14:30:00 32.36 NA 296 01.07.2004 14:31:00 32.37 NA 297 01.07.2004 14:32:00 32.35 NA 298 01.07.2004 14:33:00 32.37 14100 299 01.07.2004 14:34:00 32.31 NA 300 01.07.2004 14:35:00 32.30 NA > agg <- aggregate.data.frame(x=file, by=list(file$hora), FUN=mean) Muito obrigada! Michelle

2 2