2015-05-04 83 views
0

在下面的代码,我遇到了一个错误:问题与逻辑变量和循环

运行时错误 错误而(x0tol){:缺失值,其中TRUE/FALSE需要 执行暂停

如果我注释掉循环,x0tol返回TRUE和类逻辑。我错过了什么?先谢谢你。

a = c(1,2,2,3,3,4,5,6,6,6,8,10) 
b = c(-890,1411,-1560,-2220,-2091,-2878,-3537,-3268,-3920,-4163,-5471,-5157) 
data = data.frame(a,b) 
tol = .0001 
m = dim(data)[1] 
x0 = 0 
x1 = 0 
x0temp = 1 
x1temp = 1 
alpha = 1 
x0tol = TRUE 
while(x0tol){ 
    x0sum = 0 
    x1sum = 0 
    for (i in 1:m){ 
     x0sum = x0 + x1*data[m,1] - data[m,2] 
     x1sum = (x0 + x1*data[m,1] - data[m,2])*data[m,1] 
    } 
    x0temp = x0 - alpha*1/m*x0sum 
    x1temp = x1 - alpha*1/m*x1sum 
    x0tol = (abs(x0 - x0temp) > tol) 
    x1tol = x1 - x1temp 
    x0 = x0temp 
    x1 = x1temp 
    class(x0tol) 
    x0tol 
} 
+1

'x0tol'在循环中变成NA(因为'x0'变成NaN且'abs(NaN)> 0'NA)。 '而(NA)'产生你得到的错误。你没有描述这个意图,所以你不清楚你期望的行为。 – MrFlick

+0

另外,x0会转到-Inf,这可能不是您想要的。 – goodtimeslim

+0

我试图为两个未知参数编写一个简单的线性回归模型。预期的结果是x0和x1基于提供的数据收敛到估计值。这应该小于一些容差值。 – derrick0285

回答

0

出于某种原因,假定在while循环的条件中存在NA值。只需更改语句并创建一个带有中断的If语句就可以让循环运行。干杯

while(TRUE){ 
    x0sum = 0 
    x1sum = 0 
    for (i in 1:m){ 
    x0sum = x0 + x1*data[m,1] - data[m,2] 
    x1sum = (x0 + x1*data[m,1] - data[m,2])*data[m,1] 
    } 
    x0temp = x0 - alpha*1/m*x0sum 
    x1temp = x1 - alpha*1/m*x1sum 
    if (abs(x0 - x0temp) > tol){break} 
    x1tol = x1 - x1temp 
    x0 = x0temp 
    x1 = x1temp 
    class(x0tol) 
    x0tol 
}