当前位置: 首页>>代码示例 >>用法及示例精选 >>正文


R morley 迈克尔逊光速数据


R语言 morley 位于 datasets 包(package)。

说明

迈克尔逊(但不是莫利的这个)1879 年测量光速的经典数据。该数据由五个实验组成,每个实验由 20 个连续的 ‘runs’ 组成。响应是光速测量值,经过适当编码(公里/秒,减去299000)。

用法

morley

格式

包含以下 3 个变量的 100 个观察值的 DataFrame 。

Expt

实验编号,从 1 到 5。

Run

每个实验中的运行编号。

Speed

Speed-of-light测量。

细节

这里的数据被视为以 ‘experiment’ 和 ‘run’ 作为因子的随机区组实验。 ‘run’ 也可以被视为定量变量,以解释单个实验过程中测量值的线性(或多项式)变化。

注意

这是与包 MASS 中的 michelson 相同的数据集。

例子

require(stats); require(graphics)
michelson <- transform(morley,
                       Expt = factor(Expt), Run = factor(Run))
xtabs(~ Expt + Run, data = michelson)  # 5 x 20 balanced (two-way)
plot(Speed ~ Expt, data = michelson,
     main = "Speed of Light Data", xlab = "Experiment No.")
fm <- aov(Speed ~ Run + Expt, data = michelson)
summary(fm)
fm0 <- update(fm, . ~ . - Run)
anova(fm0, fm)

来源

A. J. Weekes (1986) Genstat 入门。伦敦:爱德华·阿诺德。

S. M. Stigler (1977) 稳健的估计器是否适用于真实数据?统计年鉴 5,1055-1098。 (参见表 6。)

A. A. Michelson (1882) 在安纳波利斯美国海军学院进行的光速实验测定。天文论文 1 135-8。美国航海年鉴办公室。 (参见表 24。)

相关用法


注:本文由纯净天空筛选整理自R-devel大神的英文原创作品 Michelson Speed of Light Data。非经特殊声明,原始代码版权归原作者所有,本译文未经允许或授权,请勿转载或复制。