Tem este arquivo 'dado' para rodar? Ele não está mais disponível no dropbox!
Bom dia,
Eu estou com um problema para suavizar meu gráfico. Eu tenho dados horários e desejo plotar individualmente cada data (no período das 06h as 18h). Entretanto o gráfico tem ficado com muito ruído devido as observações serem horárias e existirem falhas (NA) nas observações. E por isso eu gostaria de suavizar o gráfico com as médias horárias, porém não estou conseguindo implementar isso.
Espero que alguém possa me ajudar a resolver.
dado utilizado:
https://www.dropbox.com/s/b80zzzeaesnh8qp/dado.txt?dl=0
Eu estou com um problema para suavizar meu gráfico. Eu tenho dados horários e desejo plotar individualmente cada data (no período das 06h as 18h). Entretanto o gráfico tem ficado com muito ruído devido as observações serem horárias e existirem falhas (NA) nas observações. E por isso eu gostaria de suavizar o gráfico com as médias horárias, porém não estou conseguindo implementar isso.
Espero que alguém possa me ajudar a resolver.
dado utilizado:
https://www.dropbox.com/s/b80zzzeaesnh8qp/dado.txt?dl=0
library(openair)
mydata<-read.table("dado.txt", h = T, sep = ",")
mydata$direta <- as.numeric(as.character(mydata$direta))
mydata$date <- as.POSIXct(strptime(mydata$date, format = "%d/%m/%Y %H:%M:%S",tz = "Etc/GMT-1"))
means <- aggregate(mydata["direta"], format(mydata["date"],"%d-%m-%Y %H:%M:%S"), mean, na.rm = TRUE)
means[ ,"date"] <- seq(min(mydata[ , "date"]), max(mydata[ , "date"]), length = nrow(means))
Figura1<-plot(mydata$date[2:721], mydata$direta[2:721], type = "l", col="red",xlab = "minutes",main = "01 de Janeiro de 2007", ylab = "Direct Solar Radiation (W/m²)")
Figura1
dev.copy(device=jpeg,filename="Brasilia_01dejaneiro2007.jpg", width = 800, height = 500)
dev.off ()
desde já, agradeço.
desde já, agradeço.
_______________________________________________________________________________________
Postgraduate Program in Meteorology - PGMET - CPTEC / INPE
Center for Weather Forecasting and Climate Studies / National Institute for Space Research
Cachoeira Paulista - SP - Brasil