发帖
雷达卡
我写的这个程序不知道哪里出错了，大家帮忙看一下，多谢
invOC=function(alp,n,c,do,beta,q)
{
oc=function(d)
{
p=1-exp((d0/d)^beta*log(1-q))
pbinom(c,n,p,lower.tail=TRUE)-(1.0-alp)
}
uniroot(oc,interval=c(1E-50,1E50))$root
}
invOC(alp=0.05,n=106,c=0,do=0.5,beta=3,q=0.1)
错误于uniroot(oc, interval = c(1e-50, 1e+50)) : 'zeroin'函数值不对
此外: 警告信息：
1: In if (is.na(f.lower)) stop("f.lower = f(lower) is NA") :
条件的长度大于一，因此只能用其第一元素
2: In if (is.na(f.upper)) stop("f.upper = f(upper) is NA") :
条件的长度大于一，因此只能用其第一元素
3: In if (f.lower * f.upper > 0) stop("f() values at end points not of opposite sign") :
条件的长度大于一，因此只能用其第一元素
>
京公网安备 11010802022788号
