Question

我在运行R代码时遇到问题。我一直收到错误：while（Diff＆gt; 1e-06）{：缺少值需要TRUE / FALSE

时出错

以下是我正在使用的代码：

     Dose     with.tumor   total.animals
  1    0             3            50
  2  188             6            50
  3  375             6            50
  4  750            12            50

  X1=as.matrix(rep(1,4))
  X2a=as.matrix(c(0,188,375,750))
  X2=as.matrix(log(X2a))
  X=cbind(X1,X2)

a=1
b=1

B=matrix(c(a,b),nrow=2,ncol=1,byrow=F)

Diff=1
while(Diff>.000001){
P=(1-exp(-exp(a+b*X)))^-1
ni=c(rep(50,4))
V=diag(4)
for(i in 1:4)
{V[i,i]=diag(ni[i]%*%P[i]%*%(1-P[i]))}
R=c(3,6,6,12)
S=matrix(data=1,nrow=4,ncol=1)
for (i in 1:4)
{S[i]=R[i]-ni[i]*P[i]}

new.B=B+solve(t(X)%*%V%*%X)%*%(t(X)%*%S)
Diff=max(abs(new.B-B))
B=new.B}

任何见解都将受到赞赏。

Answer 1

计算NaN时您有new.B。通常，会弹出此错误，因为您要将NaN与数字进行比较。

您的t(X)有Inf作为其元素之一。这是因为您在log(0)然后X2a=as.matrix(c(0,188,375,750))

时执行X2=as.matrix(log(X2a))这一事实

具有Diff函数的R错误的Newton Raphson代码

1 个答案: