Michelle, bom dia!Encaminho uma ideia de resolução com base no que pude entender e aproveitando a sugestão do colega Rubem.### <code r>setwd("c:\\lab\\temp")AA <- read.table("AXP_NYS_USD_1min_EST5EDT_200407s12e.dat", header=FALSE,col.names=c("data","hora","preco","volume"), as.is=T)str(AA)AA$data <- as.Date(AA$data, "%d.%m.%Y")AA$hora <- format(strptime(AA$hora, "%H:%M:%S"), "%H:%M")AA$horaid <- as.numeric(as.factor(AA$hora)) ### valores utilizados no seu código originalhead(AA, 3)# data hora preco volume horaid# 1 01.07.2004 09:31 51.37 27900 1# 2 01.07.2004 09:32 51.32 100 2# 3 01.07.2004 09:33 51.45 700 3# write.table(AA,"testeteste.txt")### seleção de uma data simplesAA[AA$hora=="09:31",]# AA[AA$horaid==1,] ### resultado equivalente### reestruturando data.framerequire(reshape2)AA2 <- dcast(AA, data~hora, value.var='volume') ### Utilizando texto das horasAA3 <- dcast(AA, data~horaid, value.var='volume') ### Utilizando índice das horashead(AA2[,1:8], 3)# data 09:31 09:32 09:33 09:34 09:35 09:36 09:37# 1 2004-07-01 27900 100 700 1000 800 3000 1600# 2 2004-07-02 145400 500 NA 500 6900 NA 4400# 3 2004-07-06 3200 2300 3800 1500 NA 1100 3700head(AA3[,1:8], 3)# data 1 2 3 4 5 6 7# 1 2004-07-01 27900 100 700 1000 800 3000 1600# 2 2004-07-02 145400 500 NA 500 6900 NA 4400# 3 2004-07-06 3200 2300 3800 1500 NA 1100 3700### extradObs <- unique(AA2$data); length(dObs) ### dias observadosdEsp <- seq(min(dObs), max(dObs), by="day"); length(dEsp) ### dias esperadosdEsp[!dEsp %in% dObs] ### dias faltantes### </code>
_______________________________________________
R-br mailing list
R-br@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.