如何在GLS for Shiny应用程序中传递反应模型声明

Question

我试图在GLS中传递反应性模型语句，以便用户可以使用自己选择的独立变量和因变量。

covs <- c("TRT01P*AVISIT")
    model <- reactive({
              gls(as.formula(paste(input$sub1,"~",input$sub2,"+",paste(covs,collapse = "+"))),
                  correlation = corSymm(form = ~ 1 | USUBJID),
                  weights = varIdent(form = ~ 1 | AVISIT),
                  data = data3(),
                 method = "REML")
        })

但是我收到错误消息“找不到输入”，尽管这在LM功能中正常工作

model <- reactive({
      lm(as.formula(paste(input$sub1,"~",input$sub2,"+",paste(covs,collapse = "+"))),
          data = data3())
    })

但是对于我来说，我必须使用GLS函数，是否需要更新代码中的任何内容？