ベータ分布を円周に描く

# (0,pi/2)の角座標 theta <- seq(from=0,to=pi/2,length=100) theta <- theta[-c(1,length(theta))] # それに対応する y1=cos(theta), p1 = y1^2なる「成功確率変数」 p <- cos(theta)^2 # 円周上での一様分布 x.theta <- rep(1,length(theta)) # それに対…