Vẽ trung bình nhiều biến trong chuỗi thời gian bằng ggplot

Question

Tôi có tệp chứa dữ liệu chuỗi thời gian cho nhiều biến từ a đến k.

Tôi muốn tạo biểu đồ vẽ mức trung bình của các biến a đến k theo thời gian và cao hơn và bên dưới dòng trung bình đó thêm một vùng được làm mịn biểu thị biến thể tối đa và tối thiểu trên mỗi ngày.

Vì vậy, điều gì đó giống như khoảng tin cậy nhưng trong một phiên bản được làm phẳng.

Dưới đây là số liệu: https://dl.dropbox.com/u/22681355/co.csv

và đây là đoạn code tôi có cho đến nay:

library(ggplot2) 
library(reshape2) 
meltdf <- melt(df,id="Year") 
ggplot(meltdf,aes(x=Year,y=value,colour=variable,group=variable)) + geom_line() 

Answer 1

này mô tả bootstrapped 95% khoảng tin cậy:

ggplot(meltdf,aes(x=Year,y=value,colour=variable,group=variable)) + 
    stat_summary(fun.data = "mean_cl_boot", geom = "smooth") 

Điều này mô tả t ông có nghĩa là tất cả các giá trị của tất cả các biến + -1SD:

ggplot(meltdf,aes(x=Year,y=value)) + 
    stat_summary(fun.data ="mean_sdl", mult=1, geom = "smooth") 

Bạn có thể muốn tính toán năm có nghĩa là trước khi tính toán các phương tiện và SD trên các biến, nhưng tôi lại đó cho bạn.

Tuy nhiên, tôi tin rằng khoảng tin cậy tăng cường sẽ hợp lý hơn vì phân phối rõ ràng không đối xứng. Nó cũng sẽ hẹp hơn. ;)

Và tất nhiên bạn có thể ghi nhật ký biến đổi giá trị của mình.