Julia通过RCall调用R的garchFit出现多变量输入需公式左侧报错

RCall传递Julia数组到R环境时，如果原Julia侧的y是二维数组（比如n×1的单列矩阵）而非一维向量，传到R后会被识别为具有行、列维度的矩阵，garchFit会将其判定为多变量输入，因此要求公式左侧指定响应变量，触发你遇到的报错。

你可以任选以下任意一种方案修正：

• 方案1：传递前先将Julia侧的y转换为一维向量，再执行@rput
  示例代码：

  using RCall
  # 若原y是矩阵，用vec打平为一维向量
  y = vec(y)
  @rput y
  R"""
  library("fGarch")
  library("rugarch")
  fit.arch <- garchFit(~garch(1,0),data=y,trace=F,include.mean=FALSE)
  """
  

• 方案2：在R侧显式将传入的y转换为向量/单列数据框再传入garchFit
  示例代码：

  using RCall
  @rput y
  R"""
  library("fGarch")
  library("rugarch")
  # 转为一维向量
  y_r <- as.vector(y)
  fit.arch <- garchFit(~garch(1,0),data=y_r,trace=F,include.mean=FALSE)
  # 也可以选择转为单列数据框
  # y_df <- data.frame(y = y)
  # fit.arch <- garchFit(~garch(1,0),data=y_df,trace=F,include.mean=FALSE)
  """
  

你可以在R代码块中加入str(y)打印传入数据的结构，确认维度属性：

R"""
str(y)
"""

如果输出显示为num [1:n, 1]这类二维矩阵结构，即可确认是维度问题导致的报错。

内容的提问来源于stack exchange，提问作者David Murphy