Hatena::ブログ(Diary)

arupaka-_-arupakaの日記

2018-05-17

ggplot2で時系列的なデータを2つ並べた図を作る例。図を並べる

ggplot2で時系列的なデータを2つ並べた図を作る例。図を並べる

x1a<-cumsum(rnorm(1000,mean=0,sd=1))
x1b<-cumsum(rnorm(1000,mean=0,sd=1))
m1<-seq(0,2,length.out=length(x1a))
y1<-x1a*m1
y2<-x1b

z1<-x1a
z2<-x1b/m1
x1<-1:length(y1)
set1<-factor(c(rep("g(t)",length(x1)),rep("f(t)",length(x1))),levels=c("g(t)","f(t)"))
x0<-c(x1,x1)
y0<-c(y1,z1)
z0<-c(y2,z2)
d<-data.frame(x0,y0,z0,set1)

#d<-data.frame(y1,y2)


g3<- ggplot(data=d)
g3<-g3+theme_bw(base_size=32) 
g3<-g3+ geom_line(aes(x=x0,y=z0),colour=2,linetype=2)
g3<-g3+geom_line(aes(x=x0,y=y0))
g3<-g3+xlab("Date")+ylab("Counts")
#g3<-g3+geom_text(y=370,x=50,label="r: Random walk \n m: Scale",size=8,hjust = 0)

g3<-g3 + facet_grid(set1~.) 
print(g3)


2015-07-31

ggplot2で両対数グラフ その3. いくつかの系列をプロット

ggplot2で両対数グラフ その3. いくつかの系列をプロット

 d<-data.frame(rr_rd_mean_0[[1]],rr_rd_sd[[1]],rr_rd_mean_0[[42]],rr_rd_sd[[42]],rr_rd_mean_0[[79]],rr_rd_sd[[79]])
 colnames(d)<-c("x1","y1","x2","y2","x3","y3")
 g<-ggplot(data=d)+theme_bw(base_size=24)+scale_x_log10(labels=trans_format("log10",math_format(10^.x)))+scale_y_log10(labels=trans_format("log10",math_format(10^.x)))
 g<-g+geom_point(aes(x=x3,y=y3),size=I(1),colour=I(3),shape=I(3),alpha=0.5)
 g<-g+geom_point(aes(x=x2,y=y2),size=I(1),colour=I(2),shape=I(2),alpha=0.5)
 
 g<-g+geom_point(aes(x=x1,y=y1),size=I(1),colour=I(4),shape=I(1),alpha=0.5)
 
 g<-g+geom_line(aes(x=x1,y=sqrt(2*x1/len1[1]+0.05^2*x1^2)),size=1,linetype=2)
 g<-g+geom_line(aes(x=x1,y=sqrt(2*x1/len1[42]+0.05^2*x1^2)),colour=1,size=1,linetype=2)
 g<-g+geom_line(aes(x=x1,y=sqrt(2*x1/len1[79]+0.05^2*x1^2)),colour=1,size=1,linetype=2)
 g<-g+xlab("Mean")+ylab("S.D.")
 print(g)