powered by
rvbeta
rvnbeta(n, a, b) ("neutral" Beta distribution) is equivalent to rvbeta(n, 1/3+a, 1/3+b).
rvnbeta(n, a, b)
rvbeta(n, 1/3+a, 1/3+b)
rvbeta(n=1, shape1, shape2) rvnbeta(n=1, shape1, shape2)
See also vignette("rv").
vignette("rv")
n <- 12 # sample size y <- (0:(n-1)) # observations a <- b <- 1/3 # the neutral beta prior rvbeta(1, shape1=a+y, shape2=b+n-y) rvnbeta(1, shape1=y, shape2=n-y)
Run the code above in your browser using DataLab