我试图在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函数,是否需要更新代码中的任何内容?