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)