2015-11-04 83 views
0
我有麻烦

开始使用分段包河[R分段包“可变长度不同”

运行简单的例子当下面我得到的错误:

Error in model.frame.default(formula = y ~ x + U1.x + psi1.x, data = mfExt, : variable lengths differ (found for 'x')

我期待分段以返回具有2段的分段线性模型。我明显在我的电话中犯了一个错误,但是我无法从错误信息和文档中找出我的错误所在。帮助将不胜感激。

require(segmented) 

test.df = data.frame(x = c(1:100), 
        y = c(c(1:50),seq(from = 52, by = 2, length = 50))) 

test.mod = lm(y ~ x, 
       test.df) 

segmented(test.mod, 
      seg.Z = ~ x, 
      psi = list(x = 40)) 
+0

我无法重现您的示例中的错误。只能得到一个“完全适合的警告:总结可能不可靠” – 2015-11-04 04:19:01

回答

0

事实证明,我在我的工作空间中有一个名为'x'的对象。删除这个对象后,分段调用给出了预期的结果。

无论对象是列表还是矢量,我都可以在任何时候复制错误,我有一个名为x的长度为1的对象。

如果对象的长度大于1,则错误消失并且分段行为与预期相同。

奇怪。感谢@Pascal的输入。