1

我正在尝试使用线性混合模型编写一个闪亮的应用程序。在应用程序中，我想根据用户的输入计算 MM 的残差以用于以后的绘图。像这样的东西：

datos0=reactive({
    trait=input$trait
    a=droplevels(subset(a.0, (YEAR %in% input$year)&(SEAS %in% input$season)&(LOC %in% input$location)))
    # fit a mixed linear model with COL & ROW as fixed
    # GENO as random
    MModel=as.formula(paste(trait, "~ROW+COL"))
    modelo=lme(MModel, random=~1|GENO,data=a, na.action="na.omit")
    a$res[which(!is.na(a[,trait]))]=residuals(modelo)
    a[order(a$COL, a$ROW),]
     })

从 Shiny 外部运行它时，我没有问题。

我可以提供完整的代码，但一切正常，除非我添加混合模型的拟合（library(nlme)正在调用global.r）

有什么问题？

4

0 回答 0