[统计软件] 求大神帮助，不知道为什么会有这个missing value where TRUE/FALSE needed [推广有奖]

0关注
0粉丝

高中生

95%

还不是VIP/贵宾

威望: 0 级
论坛币: 197 个
通用积分: 0
学术水平: 1 点
热心指数: 1 点
信用等级: 0 点
经验: 193 点
帖子: 9
精华: 0
在线时间: 51 小时
注册时间: 2018-6-9
最后登录: 2021-4-27

楼主

sjy3 发表于 2018-12-3 17:26:56 |AI写论文

是否 +2 论坛币

k人参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

j<-0;u0<-5;u10<-12;u11<-10;o0<-4;o10<-5;o11<-8;p<-0.1
repeat{
  j<-j+1
  q1<-0;q2<-0;q3<-0;q4<-0;q5<-0;q6<-0;q7<-0;q8<-0;q9<-0
  for(i in 1:m1){
f1<-function(a,b){ #qab
   v1<-(1/k)*(a*u10+b*u11+(k-a-b)*u0)#uab'
   c1<-(1/k^2)*(a*o10+b*o11+(k-a-b)*o0)#oab'
   s1<-dbinom(a,e[i],p)*dbinom(b,k-e[i],p)*(2*pi*c1)^(-1/2)*exp(-(d[i]-v1)^2/(2*c1))
   s2<-0
   for(a2 in 0:e[i]){
      for(b2 in 0:k-e[i]){
      v2<-(1/k)*(a2*u10+b2*u11+(k-a2-b2)*u0)#uab'
      c2<-(1/k^2)*(a2*o10+b2*o11+(k-a2-b2)*o0)#oab'
      s2<-s2+dbinom(a2,e[i],p)*dbinom(b2,k-e[i],p)*(2*pi*c2)^(-1/2)*exp(-(d[i]-v2)^2/(2*c2))
      }
   }
   return(s1/s2)
}
f2<-function(a,b){ #qab*(yi-uab')^2
   v1<-(1/k)*(a*u10+b*u11+(k-a-b)*u0)#uab'
   c1<-(1/k^2)*(a*o10+b*o11+(k-a-b)*o0)#oab'
   s3<-dbinom(a,e[i],p)*dbinom(b,k-e[i],p)*(2*pi*c1)^(-1/2)*exp(-(d[i]-v1)^2/(2*c1))*(d[i]-v1)^2
   s4<-0
   for(a2 in 0:e[i]){
      for(b2 in 0:k-e[i]){
      v2<-(1/k)*(a2*u10+b2*u11+(k-a2-b2)*u0)#uab'
      c2<-(1/k^2)*(a2*o10+b2*o11+(k-a2-b2)*o0)#oab'
      s4<-s4+dbinom(a2,e[i],p)*dbinom(b2,k-e[i],p)*(2*pi*c2)^(-1/2)*exp(-(d[i]-v2)^2/(2*c2))
      }
   }
   return(s3/s4)
}
q1<-q1+f1(0,0)*d[i]#q00*yi连加
q2<-q2+f1(0,0)#q00连加
q3<-q3+f1(0,1)*d[i]#q01*yi连加
q4<-q4+f1(0,1)#q01连加
q5<-q5+f1(1,0)*d[i]#q10*yi连加
q6<-q6+f1(1,0)#q10连加
q7<-q7+f2(0,0)#q00*(yi-u00')^2连加
q8<-q8+f2(0,1)#q01*(yi-u01')^2连加
q9<-q9+f2(1,0)#q10*(yi-u10')^2连加
  }
  oldu0<-u0;oldu11<-u11;oldu10<-u10;oldo0<-o0;oldo11<-o11;oldo10<-o10;oldp<-p#记录上一次的迭代值
  u0<-q1/q2
  u11<-k*q3/q4-(k-1)*u0
  u10<-k*q5/q6-(k-1)*u0
  o0<-k*q7/q2
  o11<-(k^2)*q8/q4-(k-1)*o0
  o10<-(k^2)*q9/q6-(k-1)*o0
  p<-1-(q2/m1)^(1/k)
  epsilo<-1
  if((u0-oldu0)^2<epsilo&
   (u11-oldu11)^2<epsilo&
   (u10-oldu10)^2<epsilo&
   (o0-oldo0)^2<epsilo&
   (o10-oldo10)^2<epsilo&
   (o11-oldo11)^2<epsilo&
   (p-oldp)^2<epsilo)break
}
Error in if ((u0 - oldu0)^2 < epsilo & (u11 - oldu11)^2 < epsilo & (u10 -  :
  missing value where TRUE/FALSE needed
> u0
[1] NaN
> u11
[1] NaN