【发布时间】:2016-11-08 21:58:09
【问题描述】:
我对许多科目都有一组观察结果,我想为每个科目拟合一个模型。
我正在使用包data.table 和fitdistrplus,但也可以尝试使用dlpyr。
假设我的数据是这种形式:
#subject_id #observation
1 35
1 38
2 44
2 49
这是我迄今为止尝试过的:
subject_models <- dt[,fitdist(observation, "norm", method = "mme"), by=subject_id]
我认为这会导致错误,因为对fitdist 的调用返回了一个无法存储在数据表/数据帧中的fitdist 对象。
有没有使用data.table 或dplyr 的直观方法来做到这一点?
编辑:提供了一个 dplyr 答案,但我也希望有一个 data.table 答案,我会尝试针对这两者运行一些基准测试。
【问题讨论】:
-
将模型粘贴到列表列中,直到可以简化为止。 Here's a video about doing it in the tidyverse.
标签: r data.table dplyr model-fitting