R: update.formula 限制为 512 个字符

Question

我正在尝试将一个公式“粘贴”在一起以用于派对的暴民功能。这意味着公式的形式为 y ~ x1+ ...+xM | z1+..zN.

我正在对 y ~ x1+ ...+xM 进行一些初步拟合，然后想使用 update() 添加方程的条件部分。

测试代码如下:

var1 <- 1:78
x1 <- paste("x", var1, sep="")
f1 <- paste("f", var1[1:10], sep="")                                

# use first 77 variables
fmla <- as.formula( paste("y ~ ", paste(x1[1:77], collapse=" + ", sep=""), sep=""))
fmla2 <- update(fmla, paste(". ~ . | ", paste(f1, collapse= " + "), sep=""))

# CHANGE x to all 78 variables
fmla <- as.formula( paste("y ~ ", paste(x1, collapse=" + ", sep=""), sep=""))
fmla2 <- update(fmla, paste(". ~ . | ", paste(f1, collapse= " + "), sep=""))

我在 Windows 和 Linux(64 位)上运行过这个程序，但在使用所有 78 个术语(以及超过 78 个术语)时都失败了。错误消息包含 Error in parse(text = x) : :1:514: unexpected ')'。

更改 x 变量名称的长度会破坏变量较少的 update()，但总是会出现仅超过 512 个字符的错误。

有解决办法吗？

谢谢。

Answer 1

因为这看起来像是一个字符串长度问题，所以删除一些非必要的字符很简单，比如加号周围的空格，即 collapse= "+"。

这个有效:

fmla2 <- update(fmla, paste(". ~ . | ", paste(f1, collapse= "+"), sep=""))

fmla2
y ~ (x1 + x2 + x3 + x4 + x5 + x6 + x7 + x8 + x9 + x10 + x11 + 
    x12 + x13 + x14 + x15 + x16 + x17 + x18 + x19 + x20 + x21 + 
    x22 + x23 + x24 + x25 + x26 + x27 + x28 + x29 + x30 + x31 + 
    x32 + x33 + x34 + x35 + x36 + x37 + x38 + x39 + x40 + x41 + 
    x42 + x43 + x44 + x45 + x46 + x47 + x48 + x49 + x50 + x51 + 
    x52 + x53 + x54 + x55 + x56 + x57 + x58 + x59 + x60 + x61 + 
    x62 + x63 + x64 + x65 + x66 + x67 + x68 + x69 + x70 + x71 + 
    x72 + x73 + x74 + x75 + x76 + x77 | f1 + f2 + f3 + f4 + f5 + 
    f6 + f7 + f8 + f9 + f10)