22.16 迈克尔逊光速数据分析

1879 年迈克尔逊光速测量数据,记录了五次实验,每次试验测量 20 次光速,得到表格 22.2

reshape(
  data = morley, v.names = "Speed", idvar = "Expt",
  timevar = "Run", direction = "wide", sep = ""
) %>%
  knitr::kable(.,
    caption = "迈克尔逊光速数据",
    row.names = FALSE, col.names = gsub("(Speed)", "", names(.)),
    align = "c"
  )
表 22.2: 迈克尔逊光速数据
Expt 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
1 850 740 900 1070 930 850 950 980 980 880 1000 980 930 650 760 810 1000 1000 960 960
2 960 940 960 940 880 800 850 880 900 840 830 790 810 880 880 830 800 790 760 800
3 880 880 880 860 720 720 620 860 970 950 880 910 850 870 840 840 850 840 840 840
4 890 810 810 820 800 770 760 740 750 760 910 920 890 860 880 720 840 850 850 780
5 890 840 780 810 760 810 790 810 820 850 870 870 810 740 810 940 950 800 810 870

数据集 morley 中光速 Speed 已经编码过了,原始观测速度减去了 299000 (km/sec),为了展示方便

ggplot(data = morley, aes(x = Expt, y = Speed, group = Expt)) +
  geom_boxplot() +
  geom_jitter() +
  theme_minimal() +
  labs(x = "Expt", y = "Speed (km/sec)")
1879 年迈克尔逊光速实验数据

图 22.4: 1879 年迈克尔逊光速实验数据