|
# パッケージMCMCpack
library(MCMCpack)
# 逆Γ関数dinvgamma(n, alpha, beta)の作図
plot(dinvgamma(1:100,0.1,5),type="l",
xlab="x",ylab=expression(InvGamma(a,b)),
lwd=2,cex.axis=1.3,cex.lab=1.3)
# ベータ関数の作図
x <- seq(0,1,length=100)
plot(x,dbeta(x,0.5,0.5),type="l", xlab="x",ylab=expression(Beta(a,b)),
lwd=2,cex.axis=1.3,cex.lab=1.3)
|
