On this Picostat.com statistics page, you will find information about the morley data set which pertains to Michelson Speed of Light Data. The morley data set is found in the datasets R package. You can load the morley data set in R by issuing the following command at the console data("morley"). This will load the data into a variable called morley. If R says the morley data set is not found, you can try installing the package by issuing this command install.packages("datasets") and then attempt to reload the data. If you need to download R, you can go to the R project website. You can download a CSV (comma separated values) version of the morley R data set. The size of this file is about 880 bytes.
Michelson Speed of Light Data
A classical data of Michelson (but not this one with Morley) on
measurements done in 1879 on the speed of light. The data consists of
five experiments, each consisting of 20 consecutive ‘runs’.
The response is the speed of light measurement, suitably coded
A data frame with 100 observations on the following 3 variables.
The experiment number, from 1 to 5.
The run number within each experiment.
The data is here viewed as a randomized block experiment with
‘experiment’ and ‘run’ as the factors. ‘run’ may
also be considered a quantitative variate to account for linear (or
polynomial) changes in the measurement over the course of a single
This is the same dataset as
michelson in package
A. J. Weekes (1986)
A Genstat Primer.
London: Edward Arnold.
S. M. Stigler (1977)
Do robust estimators work with real data?
Annals of Statistics 5, 1055–1098. (See Table 6.)
A. A. Michelson (1882)
Experimental determination of the velocity of light made at the United
States Naval Academy, Annapolis.
Astronomic Papers 1 135–8.
U.S. Nautical Almanac Office. (See Table 24.)
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)
fm0 <- update(fm, . ~ . - Run)
Dataset imported from https://www.r-project.org.