为什么没有人回答啊？
经过几天的研究，用尝试法做了个计算。缺点是当不良率很小的时候，计算量大。

ssize <- function (p, conf.level = 0.95)
{
p.u=1
n=0
while (p.u-p > 0)
{n=n+1
z <- qnorm((1 + conf.level)/2)
z22n <- z^2/(2 * n)
p.u=(0 + z22n + z * sqrt(z22n/(2 *
            n)))/(1 + 2 * z22n)
}

list(n)
}

我收藏了，努力学习啊！！