powered by
b <- sw.beta(S, t, p, is.theta=FALSE);
# 0.72088e-3 (from p1964 of McDougall 1987) b <- sw.beta(40, 10, 4000, is.theta=TRUE)
Run the code above in your browser using DataLab