読者です 読者をやめる 読者になる 読者になる

plot中に数式を表示させる

R

パラメータを変化させるとどうなるかをplotしつつ、数式も表示する時の叩き台。自分用のメモです。

n <- seq(1,30,length.out=1000)
a <- 2:7
plot(1,1,xlim=c(min(n),max(n)),ylim=c(0,1),type="n",xlab="n",ylab="y")
title(main = expression(paste(y == frac(1,n^plain(alpha)),"のplot")))
sapply(1:length(a),function(i){lines(n,sapply(n,function(n){1/n^(1-1/a[i])}),col=i)})
legend(15,0.8,
       parse(text=sapply(a,function(a){
         paste("paste(alpha,\"=\",",a,")",sep="")
       })),
       lwd=1,col=1:length(a))

f:id:syou6162:20161010163109p:plain
この辺を参考にした。Rで数式表示は使うようで使わないのでほとんど自分の中にノウハウがたまっていない。