首页
学习
活动
专区
圈层
工具
发布
首页
学习
活动
专区
圈层
工具
MCP广场
社区首页 >问答首页 >回归回路

回归回路
EN

Stack Overflow用户
提问于 2018-09-13 17:39:41
回答 1查看 46关注 0票数 0

再见,我有几个代表分数的栏目。我想估计每个分数都是STUDYTIME函数的模型。所以我想运行尽可能多的模型,因为有分数列,所有简单的模型都是STUDYTIME的函数。然后,我希望将STUDYTIME的系数存储在一个新列中,该列的行名与得分列名相等。最后,我不知道如何在线性模型上进行聚类,因为每个学生都有两次数据。

这是我复制的例子。这是我现在掌握的数据:

代码语言:javascript
运行
复制
df <- data.frame(replicate(5, rnorm(10)))
df[1]<-c(1,1,2,2,3,3,4,4,5,5)
colnames(df) <- c('student','studytime', 'score1','score2','score3')

这是我在编码方面的尝试:

代码语言:javascript
运行
复制
for (i in 1:nrow(df)) {
  dfx         <- df[,i]
  lm    <- lm(dfx[,3:5] ~ study_time)
  resdat[,i] = summary(lm)$coefficients[2]
}
EN

回答 1

Stack Overflow用户

回答已采纳

发布于 2018-09-13 17:46:58

您可以使用简单的lapplysapply函数来完成这个任务。

以下是r代码:

生成数据

代码语言:javascript
运行
复制
df <- data.frame(replicate(5, rnorm(10)))
df[1]<-c(1,1,2,2,3,3,4,4,5,5)
colnames(df) <- c('student','studytime', 'score1','score2','score3')

存储结果

代码语言:javascript
运行
复制
Results <- lapply(df[, -c(1,2)], FUN = function(x) lm(x ~ df$studytime))
Coef <- sapply(Results, FUN = coefficients)
票数 1
EN
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接:

https://stackoverflow.com/questions/52319050

复制
相关文章

相似问题

领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档