找回密码
 立即注册
查看: 7575|回复: 2

if语句中:需要TRUE/FALSE值的地方不可以用缺少值

[复制链接]
发表于 2016-3-5 21:28:11 | 显示全部楼层 |阅读模式
求大神指导一下,应该怎么修改代码呢(红色字体标明)?万分感激!!!
#vz
F2vcz<-function(AB,vcz,lch)
{alpha<-theta(AB,e3)
F2vczRD<-abs(1/2*r*d2*(vcz-lch)^2*((sin(alpha))^3*cn90+picf)*l0)
RL<-abs(1/2*r*d2*(vcz-lch)^2*(sin(alpha))^2*cos(alpha)*cn90*l0)
if (det(rbind(AB,e1,e2))==0)
  {if(vcz>=lch)
     {F2vczRD=F2vczRD
      F2vczRLx=0  
      F2vczRLy=0
      }
    else
     {F2vczRD=-F2vczRD
      F2vczRLx=0  
      F2vczRLy=0
      }
   }

else if (det(rbind(AB,e1,e3))==0)
{if (vcz>=lch)
   {if (alpha>=0&&alpha<=pi/2)
    {F2vczRD=F2vczRD
     F2vczRLx=RL
     F2vczRLy=0
     }
    else (alpha>pi/2&&alpha<=pi)
    {F2vczRD=F2vczRD
     F2vczRLx=-RL
     F2vczRLy=0
     }
    }
   else
   {if (alpha>=0&&alpha<=pi/2)
    {F2vczRD=-F2vczRD
     F2vczRLx=-RL
     F2vczRLy=0
     }
    else (alpha>pi/2&&alpha<=pi)
    {F2vczRD=-F2vczRD
     F2vczRLx=RL
     F2vczRLy=0
     }
    }
  }


Error in if (vcz>=lch) {: 需要TRUE/FALSE值的地方不可以用缺少值
QQ图片20160305212305.png
回复

使用道具 举报

发表于 2016-3-30 17:58:59 | 显示全部楼层
看看vcz,lch的哪个字段有缺失值,加个length(字段)是否等于0的if判断,希望能帮到你
回复

使用道具 举报

发表于 2016-3-31 11:25:02 | 显示全部楼层
na.rm=T行吗,我也是新手。。
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Archiver|手机版|小黑屋|R语言中文网

GMT+8, 2024-11-23 11:28 , Processed in 0.023329 second(s), 21 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表