facet_grid を使用して 8 年間 (2004 ～ 2011 年) のシリーズの月間降水量データをプロットして、各年の月間値を 1 つのプロットにプロットしようとしています。しかし、多くの場合、プロット時に問題が発生します (おそらく私のせいです)。各プロットの 2004 年から 2011 年までの x 軸で取得したプロットを見ることができますが、年/ファセットごとに 1 月から 12 月までプロットを実行する必要があります。どうすればそのエラーを取り除くことができますか?
どんな助けでも大歓迎です。
ここでは、dput(pcp.mensual) http://ubuntuone.com/0EZvmwXGnSkqNVwXZOnx8Eの出力を取得できます。以下は、私が使用するコードです (より長いスクリプトの一部)。
library(ggplot2)
library(plyr)
library(lubridate)
library(scales)
#---------------------------------------------------------------------
system("awk '{ print $1 \";\" $2 \";\" $3 \";\" $5 \";\" $10 \";\" $11 \";\" $13 }' visbel.cor > kk.dat",intern=T)
# Lectura de datos
datos=read.csv("kk.dat",sep=";",header=T,na.strings="-99.900")
# Inclusión de una columna con formato temporal
# (dato original en formato %y/%m/%d %H:%M:%S)
datos=within(datos, datetime <- as.POSIXct(paste(FECHA,H_SOLAR),format = "%y/%m/%d %H:%M:%S"))
# Creación del objeto datetime
#datetime <- as.POSIXct(paste(datos$FECHA,datos$H_SOLAR),format = "%y/%m/%d %H:%M:%S")
# Eliminación de valores NA en Precipitación para barplot
datos$Precipitacion[is.na(datos$Precipitacion)]=0
#---------------------------------------------------------------------
pcp.diaria=aggregate(cbind(Precip.diaria=Precipitacion) ~ FECHA, datos, sum)
pcp.diaria=within(pcp.diaria, date <- as.POSIXct(paste(pcp.diaria$FECHA),format = "%y/%m/%d"))
# Agregación mensual
pcp.diaria$mes <- floor_date(pcp.diaria$date, "month")
pcp.mensual=ddply(pcp.diaria, "mes", summarise, x = sum(Precip.diaria,na.rm=TRUE))
pcp.mensual$ano=year(pcp.mensual$mes)
colnames(pcp.mensual) <- c("mes", "Precip.mensual","ano")
#---------------------------------------------------------------------
ggplot(data=pcp.mensual,aes(x=mes, y=Precip.mensual)) +
facet_grid(ano ~. ) +
geom_bar(colour="blue",stat="identity",fill="blue") +