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