循环通过Dataframe for中的列以构建重复测量线性模型lme in R

install.packages("nlme")

library(nlme)

# 假设我们有一个名为df的Dataframe，包含了多个重复测量的变量
# 变量之间的重复测量关系由一个名为Subject的变量指定
# 我们想要循环通过每一列构建重复测量线性模型

# 首先，创建一个空的list，用于存储模型结果
model_results <- list()

# 循环遍历df的每一列
for (col in colnames(df)) {
  if (col != "Subject") {
    # 提取当前列和Subject列的数据
    data <- cbind(df[[col]], df$Subject)
    
    # 构建重复测量线性模型
    model <- lme(fixed = col ~ Time, random = ~1 | Subject, data = data)
    
    # 将模型结果添加到list中
    model_results[[col]] <- model
  }
}

# 输出每个列的模型结果
for (col in colnames(df)) {
  if (col != "Subject") {
    print(paste("模型结果 -", col))
    print(model_results[[col]])
  }
}