睡眠研究数据
研究人员在一项睡眠剥夺研究中招募了 18 名受试者。其观测到的 sleep_study 数据已加载到工作区中。数据包含每位 subject 在 day_0 的反应时以及经历 3 个睡眠剥夺之夜后的 day_3 反应时。
您将定义并探索 diff_3,即每位受试者反应时的"差值"。这需要使用 mutate() 和 summarize() 函数。例如,下面的代码会向 sleep_study 添加变量 day_0_s,表示以"秒"为单位的 day_0 反应时:
sleep_study <- sleep_study %>%
mutate(day_0_s = day_0 * 0.001)
随后,您可以使用 summarize() 汇总 day_0_s 的取值,这里以最小值和最大值为例:
sleep_study %>%
summarize(min(day_0_s), max(day_0_s))
本练习是课程的一部分
使用 RJAGS 的贝叶斯建模
练习说明
- 查看
sleep_study的前 6 行。 - 定义新的
sleep_study变量diff_3,其值为day_3减去day_0的反应时。 - 使用
ggplot()并添加geom_histogram()图层,绘制diff_3的直方图。 - 使用
summarize()汇总diff_3观测的均值和标准差。
交互式实操练习
通过完成这段示例代码来试试这个练习。
# Check out the first 6 rows of sleep_study
# Define diff_3
sleep_study <- sleep_study %>%
mutate(diff_3 = ___)
# Histogram of diff_3
ggplot(___, aes(x = ___)) +
geom_histogram(binwidth = 20, color = "white")
# Mean and standard deviation of diff_3
sleep_study %>%
summarize(___, ___)