multimode (version 1.4)

stamps: Stamps thickness

Description

This dataset, analysed in Izenman and Sommer (1988) and Ameijeiras--Alonso et al. (2018), consists of thickness measurements (in millimeters) of 485 unwatermarked used white wove stamps of the 1872 Hidalgo stamp issue of Mexico. All of them had an overprint with the year (1872 or either an 1873 or 1874) and some of them were watermarked (Papel Sellado or LA+-F).

Usage

data(stamps)
data(stamps1)
data(stamps2)
data(stampstable)

Arguments

Format

stamps includes the thickness (in millimeters) of the different stamps in the 1872 Hidalgo stamp issue. stampstable reproduces the original frequency table in Izenman and Sommer (1988), including the thickness, the overprinted years (1872 or 1873-1874) and the watermarks (Papel Sellado or LA+-F). stamps1 includes both the thickness and the overprinted year. stamps2 includes both the thickness and the watermark. Note that the stamp with a thickness equal to 0.118 and the watermark "Papel Sellado" was removed in stamps2 since it does not have a corresponding year, probably this watermark is associated with the observation of thickness 0.117 or 0.119.

Details

This is a classic example for determining the number of modes.

References

Jose Ameijeiras--Alonso, Rosa M. Crujeiras, Alberto Rodr<U+00ED>guez--Casal (2018). Mode testing, critical bandwidth and excess mass, to appear in Test.

Izenman, A. J., and Sommer, C. J. (1988). Philatelic mixtures and multimodal densities. Journal of the American Statistical association, 83, 941--953.

Examples

Run this code
# NOT RUN {
data("stamps")
# Histogram of Wilson (bin width 0.008)
seqx=seq(0.0585,0.1385,by=0.008)
hist(stamps,breaks=seqx)
# Histogram of Figure 1, Izenman and Sommer (bin width 0.002)
seqx=seq(0.0585,0.1385,by=0.002)
hist(stamps,breaks=seqx)

# Kernel density estimation of Izenman and Sommer (number of modes: 7)
locmodes(stamps,mod0=7,lowsup=0.04,uppsup=0.15,display=TRUE)
# Kernel density estimation of Ameijeiras-Alonso et al. (number of modes: 4)
locmodes(stamps,mod0=4,lowsup=0.04,uppsup=0.15,display=TRUE)

# }

Run the code above in your browser using DataCamp Workspace