Error in if (u < loga) { : 需要TRUEFALSE值的地方不可以用缺少值

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

各位大神，程序运行提示：Error in if (u < loga) { : 需要TRUE/FALSE值的地方不可以用缺少值，不明白，请指示
rt<-read.table("logistic.txt")#读取数据
rt
wais<-matrix(0,nrow=54,ncol=2)
wais[,1]=rt[,2]
wais[,2]=rt[,3] #数据读入wais中
y<-wais[,2]
x<-wais[,1]
m<-10000
beta0<-c(0,0)  #initial value
mu.beta<-c(0,0) #prior mean
s.beta<-c(0,0)  #prior sd
prop.s<-c(1.75,0.2)#sd of proposal normal
beta<-matrix(nrow=m,ncol=2)
    acc.prob<-c(0,0)
    current.beta<-beta0
    for (t in 1:m)
        {
        for(j in 1:2)
           {
            prop.beta<-current.beta
            prop.beta[j]<-rnorm(1,current.beta[j],prop.s[j])#先抽第一个分量
            cur.eta<-current.beta[1]+current.beta[2]*x
            prop.eta<-prop.beta[1]+prop.beta[2]*x           #做一次比较
            if(sum(prop.eta>700)>0){print(t);stop;}
            if (sum(cur.eta>700)>0){print(t);stop;}
            loga<-(sum(y*prop.eta-log(1+exp(prop.eta)))
                    -sum(y*cur.eta-log(1+exp(cur.eta)))
                    +sum(dnorm(prop.beta,mu.beta,s.beta,log=TRUE))
                     -sum(dnorm(current.beta,mu.beta,s.beta,log=TRUE)))
            u<-runif(1)
            u<-log(u)
            if(u<loga)
  {  
   current.beta<-prop.beta
   acc.prob<-acc.prob+1
   
  }
            
           }
            beta[t,]<-current.beta

            }

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：TRUEFALSE Error Log Rue err Error 程序

loga的值可能是NA  如下面显而易见的例子：


x<-NA
y<-3
z<-5
if (x<2)
{ x=y+z
}

运行结果就是：Error in if (x < 2) { : missing value where TRUE/FALSE needed

谢谢，我查一下