plot中に数式を表示させる

パラメータを変化させるとどうなるかを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で数式表示は使うようで使わないのでほとんど自分の中にノウハウがたまっていない。