| Version: | 1.0.9 | 
| Date: | 2025-02-07 | 
| Title: | Datasets and Functions for Books by Julian Faraway | 
| Description: | Books are "Linear Models with R" published 1st Ed. August 2004, 2nd Ed. July 2014, 3rd Ed. February 2025 by CRC press, ISBN 9781439887332, and "Extending the Linear Model with R" published by CRC press in 1st Ed. December 2005 and 2nd Ed. March 2016, ISBN 9781584884248 and "Practical Regression and ANOVA in R" contributed documentation on CRAN (now very dated). | 
| Depends: | R (≥ 3.5.0) | 
| License: | GPL-2 | GPL-3 [expanded from: GPL] | 
| URL: | https://github.com/julianfaraway/faraway | 
| LazyData: | yes | 
| Imports: | methods, lme4, nlme | 
| Suggests: | leaps | 
| Encoding: | UTF-8 | 
| RoxygenNote: | 7.3.2 | 
| NeedsCompilation: | no | 
| Packaged: | 2025-02-07 11:30:48 UTC; jjf23 | 
| Author: | Julian Faraway [aut, cre] | 
| Maintainer: | Julian Faraway <jjf23@bath.ac.uk> | 
| Repository: | CRAN | 
| Date/Publication: | 2025-02-07 13:00:02 UTC | 
Cp plot
Description
Makes a Cp plot
Usage
Cpplot(cp)
Arguments
| cp | A leaps object returned from leaps() | 
Details
Requires leaps package
Value
none
Author(s)
Julian Faraway
See Also
leaps()
Helper function to convert VarCorr objects
Description
Used internally when showing summary for lmer objects
Usage
VarCorrmat(varc, useScale, digits)
Arguments
| varc | a numeric vector | 
| useScale | a numeric vector | 
| digits | a numeric vector | 
Value
a matrix
Author(s)
copied from the arm package
Annual mean temperatures in Ann Arbor, Michigan
Description
The data comes from the U.S. Historical Climatology Network.
Format
A data frame with 115 observations on the following 2 variables.
- year
- year from 1854 to 2000 
- temp
- annual mean temperatures in degrees F in Ann Arbor 
Source
United States Historical Climatology Network: https://www.ncei.noaa.gov/products/land-based-station/us-historical-climatology-network
Wear on materials according to type, run and position
Description
The abrasion data frame has 16 rows and 4 columns.  Four materials
were fed into a wear testing machine and the amount of wear recorded. Four
samples could be processed at the same time and the position of these
samples may be important.  A Latin square design was used.
Format
This data frame contains the following columns:
- run
- The run number 1-4 
- position
- The position number 1-4 
- material
- The material A-D 
- wear
- The wear measured loss of weight in 0.1mm over testing period 
Source
The Design and Analysis of Industrial Experiments by O. Davies, 1954, published by Wiley
aflatoxin dosage and liver cancer in lab animals
Description
Aflatoxin B1 was fed to lab animals at vary doses and the number responding with liver cancer recorded.
Format
A data frame with 6 observations on the following 3 variables.
- dose
- dose in ppb 
- total
- number of test animals 
- tumor
- number with liver cancer 
Source
Gaylor DW (1987) "Linear nonparametric upper limits for low dose extrapolation" ASA Proceedings of the Biopharmaceutical Section.
Examples
data(aflatoxin)
miltary coups and politics in sub-Saharan Africa
Description
Data is a subset of a larger study on factors affecting regime stability in Sub-Saharan Africa
Format
A data frame with 47 observations on the following 9 variables.
- miltcoup
- number of successful military coups from independence to 1989 
- oligarchy
- number years country ruled by military oligarchy from independence to 1989 
- pollib
- Political liberalization - 0 = no civil rights for political expression, 1 = limited civil rights for expression but right to form political parties, 2 = full civil rights 
- parties
- Number of legal political parties in 1993 
- pctvote
- Percent voting in last election 
- popn
- Population in millions in 1989 
- size
- Area in 1000 square km 
- numelec
- Total number of legislative and presidential elections 
- numregim
- Number of regime types 
Source
Bratton, Michael, and Nicholas Van De Walle. 1997. “Political Regimes and Regime Transitions in Africa, 1910-1994.” Study Number I06996. Ann Arbor: Inter-University Consortium for Political and Social Research.
References
"Bayesian Methods: A Social and Behavioral Sciences Approach" by Jeff Gill 2002.
Airline passengers
Description
Monthly totals of airline passengers from 1949 to 1951
Format
A data frame with 144 observations on the following 2 variables.
- pass
- number of passengers in thousands 
- year
- the date as a decimal 
Details
Well known time series example dataset
Source
Brown, R.G.(1962) Smoothing, Forecasting and Prediction of Discrete Time Series. Englewood Cliffs, N.J.: Prentice-Hall.
References
Box, G.E.P., Jenkins, G.M. and Reinsel, G.C. (1994) Time Series Analysis, Forecasting and Control, 3rd edn. Englewood Cliffs, N.J.: Prentice-Hall.
Examples
data(airpass)
## maybe str(airpass) ; plot(airpass) ...
Effects of seed inoculum, irrigation and shade on alfalfa yield
Description
The alfalfa data frame has 25 rows and 4 columns. Data comes from an
experiment to test the effects of seed inoculum, irrigation and shade on
alfalfa yield. A latin square design has been used.
Format
This data frame contains the following columns:
- shade
- Distance of location from tree line divided into 5 shade areas 
- irrigation
- Irrigation effect divided into 5 levels 
- inoculum
- Four types of seed incolum, A-D with E as control. 
- yield
- Dry matter yield of alfalfa 
Source
Petersen, R.G. 1994. Agricultural Field Experiments, Design and Analysis. Marcel Dekker, Inc., New York. Pages 70-74. 1994
Match pair study for AML and Xray link
Description
A matched case control study carried out to investigate the connection between X-ray usage and acute myeloid leukemia in childhood. The pairs are matched by age, race and county of residence.
Format
A data frame with 238 observations on the following 11 variables.
- ID
- a factor denoting the matched pairs 
- disease
- 0=control, 1=case 
- Sex
-  ForM
- downs
- Presence of Downs syndrome: - noor- yes
- age
- Age in years 
- Mray
- Did the mother ever have an Xray: - noor- yes
- MupRay
- Did the mother have an Xray of the upper body during pregnancy: - noor- yes
- MlowRay
- Did the mother have an Xray of the lower body during pregnancy: - noor- yes
- Fray
- Did the father ever have an Xray: - noor- yes
- Cray
- Did the child ever have an Xray: - noor- yes
- CnRay
- Total number of Xrays of the child - 1=none <- 2=1 or 2 <- 3=3 or 4 <- 4= 5 or more
Source
Chap T. Le (1998) "Applied Categorical Data Analysis" Wiley.
Time in minutes to eye opening after reversal of anaesthetic.
Description
A doctor at major London hospital compared the effects of 4 anaesthetics used in major operations. 80 patients were divided into groups of 20.
Format
A data frame with 80 observations on the following 2 variables.
- breath
- time in minutes to start breathing unassisted 
- tgrp
- Four treatment groups - A- B- C- D
Source
Chatfield C. (1995) Problem Solving: A Statistician's Guide, 2ed Chapman Hall.
Examples
data(anaesthetic)
## maybe str(anaesthetic) ; plot(anaesthetic) ...
Respiratory disease rates of babies fed in different ways
Description
Study on infant respiratory disease, namely the proportions of children developing bronchitis or pneumonia in their first year of life by type of feeding and sex.
Format
A data frame with 6 observations on the following 4 variables.
- disease
- number with disease 
- nondisease
- number without disease 
- sex
- a factor with levels - Boy- Girl
- food
- a factor with levels - Bottle- Breast- Suppl
Source
Payne, C. (1987). The GLIM System Release 3.77 Manual (2 ed.). Oxford: Nu- merical Algorithms Group.
Examples
data(babyfood)
## maybe str(babyfood) ; plot(babyfood) ...
Beetles exposed to fumigant
Description
Grain beetles were exposed to ethylene oxide
Format
A data frame with 10 observations on the following 3 variables.
- conc
- concentration of ethylene oxide in mg/l 
- affected
- number affected 
- exposed
- number exposed 
Source
Busvine (1938)
References
Collet D. "Modelling Binary Data"
Examples
data(beetle)
## maybe str(beetle) ; plot(beetle) ...
Insect mortality due to insecticide
Description
An experiment measuring death rates for insects, with 30 insects at each of five treatment levels.
Format
A data frame with 5 observations on the following 3 variables.
- dead
- number dead 
- alive
- number alive 
- conc
- concentration of insecticide 
Source
Bliss (1935). The calculation of the dosage-mortality curve. Annals of Applied Biology 22, 134-167.
Examples
data(bliss)
## maybe str(bliss) ; plot(bliss) ...
Breaking strength of materials
Description
An experiment was conducted to select the supplier of raw materials for production of a component. The breaking strength of the component was the objective of interest. Four suppliers were considered. The four operators can only produce one component each per day. A Latin square design was used.
Format
A data frame with 16 observations on the following 4 variables.
- y
- The breaking strength of the component 
- operator
- the operator - a factor with levels - op1- op2- op3- op4
- day
- the day of production - a factor with levels - day1- day2- day3- day4
- supplier
- the supplier of the raw material - a factor with levels - A- B- C- D
Source
Lentner M. and Bishop T. (1986) Experimental Design and Analysis, Valley Book Company
Broccoli weight variation
Description
A number of growers supply broccoli to a food processing plant. The plant instructs the growers to pack the broccoli into standard size boxes. There should be 18 clusters of broccoli per box and each cluster should weigh between 1.33 and 1.5 pounds. Because the growers use different varieties, methods of cultivation etc, there is some variation in the cluster weights. The plant manager selected 3 growers at random and then 4 boxes at random supplied by these growers. 3 clusters were selected from each box.
Format
A data frame with 36 observations on the following 4 variables.
- wt
- weight of broccoli 
- grower
- the grower - a factor with levels - 1- 2- 3
- box
- the box - a factor with levels - 1- 2- 3- 4
- cluster
- the cluster - a factor with levels - 1- 2- 3
Source
Lentner M. and Bishop T. (1986) Experimental Design and Analysis, Valley Book Company
Butterfat content of milk by breed
Description
Average butterfat content (percentages) of milk for random samples of twenty cows (ten two-year old and ten mature (greater than four years old)) from each of five breeds. The data are from Canadian records of pure-bred dairy cattle.
Format
A data frame with 100 observations on the following 3 variables.
- Butterfat
- butter fat content by percentage 
- Breed
- a factor with levels - Ayrshire- Canadian- Guernsey- Holstein-Fresian- Jersey
- Age
- a factor with levels - 2year- Mature
Source
Sokal, R. R. and Rohlf, F. J. (1994) Biometry. W. H. Freeman, New York, third edition.
Examples
data(butterfat)
## maybe str(butterfat) ; plot(butterfat) ...
Cathedral nave heights and lengths in England
Description
Example Dataset from "Practical Regression and Anova"
Format
A dataset with 25 cases
- style
- of the cathedral - romanesque or gothic 
- height
- in feet 
- width
- in feet 
Source
Weisberg, S. (2005). Applied Linear Regression, 3rd edition. New York: Wiley
References
Reference details may be found in "Practical Regression and Anova" by Julian Faraway
Taste of Cheddar cheese
Description
In a study of cheddar cheese from the LaTrobe Valley of Victoria, Australia, samples of cheese were analyzed for their chemical composition and were subjected to taste tests. Overall taste scores were obtained by combining the scores from several tasters.
Format
A data frame with 30 observations on the following 4 variables.
- taste
- a subjective taste score 
- Acetic
- concentration of acetic acid (log scale) 
- H2S
- concentration of hydrogen sulfide (log scale) 
- Lactic
- concentration of lactic acid 
Source
David S. Moore and George P. McCabe (1993) Introduction to the Practice of Statistics, W. H. Freeman and company, second edition.
Examples
data(cheddar)
## maybe str(cheddar) ; plot(cheddar) ...
Chicago insurance redlining
Description
Data from a 1970's study on the relationship between insurance redlining in Chicago and racial composition, fire and theft rates, age of housing and income in 47 zip codes.
Format
This dataframe contains the following columns
- race
- racial composition in percent minority 
- fire
- fires per 100 housing units 
- theft
- 
theft per 1000 population 
- age
- percent of housing units built before 1939 
- involact
- new FAIR plan policies and renewals per 100 housing units 
- income
- median family income in thousands of dollars 
- side
- North or South side of Chicago 
Source
Adapted from "Data : A Collection of Problems from Many Fields for the Student and Research Worker" by D. Andrews and A. Herzberg published by Springer-Verlag, in 1985
Chicago zip codes north-south
Description
Complements the chicago and chmiss datasets by dividing the zip codes into north and south
Format
- chiczip
- takes the values "n" (north) and "s" south 
References
Reference details may be found in "Practical Regression and Anova" by Julian Faraway
See Also
chicago
Chicago insurance redlining
Description
Data from a 1970's study on the relationship between insurance redlining in Chicago and racial composition, fire and theft rates, age of housing and income in 47 zip codes. Missing values have been randomly added.
Format
This dataframe contains the following columns
- race
- racial composition in percent minority 
- fire
- fires per 100 housing units 
- theft
- 
theft per 1000 population 
- age
- percent of housing units built before 1939 
- involact
- new FAIR plan policies and renewals per 100 housing units 
- income
- median family income in thousands of dollars 
- side
- North or South side of Chicago 
Source
Adapted from "Data : A Collection of Problems from Many Fields for the Student and Research Worker" by D. Andrews and A. Herzberg published by Springer-Verlag, in 1985
Chocolate cake experiment with split plot design
Description
An experiment was conducted to determine the effect of recipe and baking temperature on chocolate cake quality. 15 batches of cake mix for each recipe were prepared. Each batch was sufficient for six cakes. Each of the six cakes was baked at a different temperature which was randomly assigned. Several measures of cake quality were recorded of which breaking angle was just one.
Format
A data frame with 270 observations on the following 4 variables.
- recipe
- Chocolate for recipe 1 was added at 40C, Chocolate for recipe 2 was added at 60C and recipe 3 had extra sugar 
- batch
- batch number from 1 to 15 
- temp
- temperature at which cake was baked: - 175C- 185C- 195C- 205C- 215C- 225C
- breakang
- the breaking angle of the cake 
Source
Cochran W. and Cox G. (1992) Experimental Designs, 2nd Edition Wiley
Chicago insurance redlining
Description
Data from a 1970's study on the relationship between insurance redlining in Chicago and racial composition, fire and theft rates, age of housing and income in 47 zip codes
Format
This dataframe contains the following columns
- race
- racial composition in percent minority 
- fire
- fires per 100 housing units 
- theft
- 
theft per 1000 population 
- age
- percent of housing units built before 1939 
- involact
- new FAIR plan policies and renewals per 100 housing units 
- income
- median family income in thousands of dollars 
- side
- North or South side of Chicago 
Source
Adapted from "Data : A Collection of Problems from Many Fields for the Student and Research Worker" by D. Andrews and A. Herzberg published by Springer-Verlag, in 1985
Blood clotting times
Description
The clotting times of blood for plasma diluted with nine different percentage concentrations with prothrombin-free plasma
Format
This data frame contains the following columns:
- time
- time in seconds to clot 
- conc
- concentration in percent 
- lot
- lot number - either one or two 
Source
Hurn et al (1945)
References
Nelder & McCullagh (1989) Generalized Linear Models (2ed)
Social class mobility from 1971 to 1981 in the UK
Description
Social class mobility from 1971 to 1981 for 42425 men from the United Kingdom census. Subjects were aged 45-64.
Format
A data frame with 36 observations on the following 3 variables.
- y
- Frequency of observation 
- class71
- social class in 1971 - a factor with levels - I, professionals,- IIsemi-professionals,- IIINskilled non-manual,- IIIMskilled manual,- IVsemi-skilled,- Vunskilled
- class81
- social class in 1971 - a factor with levels - I- II- IIIN- IIIM- IV- Vwith same classification
Source
D. Blane and S. Harding and M. Rosato (1999) "Does social mobility affect the size of the socioeconomic mortality differential?: Evidence from the Office for National Statistics Longitudinal Study" JRSS-A, 162 59-70.
Malformations of the central nervous system
Description
Frequencies of various malformations of the central nervous system recorded on live births in South Wales, UK. Study was designed to determine the effect of water hardness on the incidence of such malformations.
Format
A data frame with 16 observations on the following 7 variables.
- Area
- a factor with levels - Cardiff- GlamorganC- GlamorganE- GlamorganW- MonmouthOther- MonmouthV- Newport- Swanseabeing areas of South Wales
- NoCNS
- count of births with no CNS problem 
- An
- count of Anencephalus births 
- Sp
- count of Spina Bifida births 
- Other
- count of other CNS births 
- Water
- water hardeness 
- Work
- a factor with levels - Manual- NonManualbeing the type of work done by the parents
Source
C. Lowe and C. Roberts and S. Lloyd, (1971) Malformations of the central nervous system and softness of local water supplies, British Medical Journal, 15,357-361.
References
P. McCullagh and J. Nelder (1989), Generalized Linear Models, Chapman and Hall, 2nd Ed.
Blood coagulation times by diet
Description
Dataset comes from a study of blood coagulation times. 24 animals were randomly assigned to four different diets and the samples were taken in a random order.
Format
This dataframe contains the following columns
- coag
- coagulation time in seconds 
- diet
- diet type - A,B,C or D 
Source
"Statistics for Experimenters" by G. P. Box, W. G. Hunter and J. S. Hunter, Wiley, 1978
Strength of a thermoplastic composite depending on two factors
Description
The composite data frame has 9 rows and 3 columns. Data comes from an
experiment to test the strength of a thermoplastic composite depending on
the power of a laser and speed of a tape.
Format
This data frame contains the following columns:
- strength
- interply bond strength of the composite 
- laser
- laser power at 40, 50 or 60W 
- tape
- 
tape speed, slow=6.42 m/s, medium=13m/s and fast=27m/s 
Source
Mazumdar, S and Hoa S (1995) "Application of a Taguchi Method for Process enhancement of an online consolidation technique" Composites 26, 669-673
Corn yields from nitrogen application
Description
The relationship between corn yield (bushels per acre) and nitrogen (pounds per acre) fertilizer application were studied in Wisconsin.
Format
A data frame with 44 observations on the following 2 variables.
- yield
- corn yield in bushels per acre 
- nitrogen
- pounds per acre 
Source
Unknown
Corrosion loss in Cu-Ni alloys
Description
Data consist of thirteen specimens of 90/10 Cu-Ni alloys with varying iron content in percent. The specimens were submerged in sea water for 60 days and the weight loss due to corrosion was recorded in units of milligrams per square decimeter per day.
Format
This dataframe contains the following columns
- Fe
- Iron content in percent 
- loss
- Weight loss in mg per square decimeter per day 
Source
"Applied Regression Analysis" by N. Draper and H. Smith, Wiley, 1998
Projected and actual sales of 20 consumer products
Description
Projected and actual sales of 20 consumer products. Data have been disguised from original form.
Format
A data frame with 20 observations on the following 2 variables.
- projected
- projected sales in dollars 
- actual
- actual sales in dollars 
Source
G. Whitmore (1986) "Inverse Gaussian Ratio Estimation" Applied Statistics, 35, 8-15.
Crawling babies by month
Description
A study investigated whether babies take longer to learn to crawl in cold months when they are often bundled in clothes that restrict their movement, than in warmer months. The study sought an association between babies' first crawling age and the average temperature during the month they first try to crawl (about 6 months after birth). Parents brought their babies into the University of Denver Infant Study Center between 1988-1991 for the study. The parents reported the birth month and age at which their child was first able to creep or crawl a distance of four feet in one minute. Data were collected on 208 boys and 206 girls (40 pairs of which were twins)
Format
A data frame with 12 observations on the following 4 variables.
- crawling
- average crawling age in weeks 
- SD
- standard deviation of crawling age 
- n
- sample size 
- temperature
- average temperature(F) six months after birth 
Source
Benson, Janette. (1993). Infant Behavior and Development
Examples
data(crawl)
## maybe str(crawl) ; plot(crawl) ...
Effects of surface and vision on balance.
Description
An experiment was conducted to study the effects of surface and vision on balance. The balance of subjects were observed for two different surfaces and for restricted and unrestricted vision. Balance was assessed qualitatively on an ordinal four-point scale based on observation by the experimenter. Forty subjects were studied, twenty males and twenty females ranging in age from 18 to 38, with heights given in cm and weights in kg. The subjects were tested while standing on foam or a normal surface and with their eyes closed or open or with a dome placed over their head. Each subject was tested twice in each of the surface and eye combinations for a total of 12 measures per subject.
Format
A data frame with 480 observations on the following 8 variables.
- Subject
- an indicator 
- Sex
- a factor with levels - female- male
- Age
- in years 
- Height
- in cm 
- Weight
- in kg 
- Surface
- a factor with levels - foam- norm
- Vision
- a factor with levels - closed- dome- open
- CTSIB
- a four point scale measuring balance 
Source
Steele, R. (1998). Effect of surface and vision on balance. Ph. D. thesis, Depart- ment of Physiotherapy, University of Queensland.
References
OzDasl
Examples
data(ctsib)
## maybe str(ctsib) ; plot(ctsib) ...
Death penalty in Florida 1977
Description
Data on 326 defendents in homicide indictments in 20 Florida counties during 1976-77.
Format
A data frame with 8 observations on the following 4 variables.
- y
- a numeric vector 
- penalty
- Did the subject recieve the death penalty? - noor- yes
- victim
- Was the victim - black or- white?
- defend
- Was the defendent - black or- white?
Source
Radelet M. (1981) Racial characteristics and the imposition of the death penalty. Amer. Sociol. Rev. 46 918-927.
References
Agresti A. (1990) Categorical Data Analysis, Wiley.
psychology of debt
Description
The data arise from a large postal survey on the psychology of debt.
Format
A data frame with 464 observations on the following 13 variables.
- incomegp
- income group (1=lowest, 5=highest) 
- house
- security of housing tenure (1=rent, 2=mortgage, 3=owned outright) 
- children
- number of children in household 
- singpar
- is the respondent a single parent? 
- agegp
- age group (1=youngest) 
- bankacc
- does the respondent have a bank account? 
- bsocacc
- does the respondent have a building society account? 
- manage
- self-rating of money management skill (high values=high skill) 
- ccarduse
- how often did s/he use credit cards (1=never... 3=regularly) 
- cigbuy
- does s/he buy cigarettes? 
- xmasbuy
- does s/he buy Christmas presents for children? 
- locintrn
- score on a locus of control scale (high values=internal) 
- prodebt
- score on a scale of attitudes to debt (high values=favourable to debt 
Details
All yes/no questions are coded 0=no, 1=yes. Locus of control is a personality measure introduced by Rotter, which claims to differentiate people according to how much they feel things that happen to them are as a result of processes within themselves (internal locus of control) or outside events (external locus of control).
Source
Lea, Webley & Walker, 1995, Journal of Economic Psychology, 16, 181-201 Data obtained from http://au.exeter.ac.uk/SEGLea/.
Denim wastage by supplier
Description
Five suppliers cut denim material for a jeans manufacturer. An algorithm is used to estimate how much material will be wasted given the dimensions of the material supplied. Typically, a supplier wastes more material than the target based on the algorithm although occasionally they waste less. The percentage of waste relative to target was collected weekly for the 5 suppliers. In all, 95 observations were recorded.
Format
A data frame with 95 observations on the following 2 variables.
- waste
- percentage wastage 
- supplier
- a factor with levels - 1- 2- 3- 4- 5
Source
Unknown
Examples
data(denim)
## maybe str(denim) ; plot(denim) ...
Diabetes and obesity, cardiovascular risk factors
Description
403 African Americans were interviewed in a study to understand the prevalence of obesity, diabetes, and other cardiovascular risk factors in central Virginia.
Format
A data frame with 403 observations on the following 19 variables.
- id
- Subject ID 
- chol
- Total Cholesterol 
- stab.glu
- Stabilized Glucose 
- hdl
- High Density Lipoprotein 
- ratio
- Cholesterol/HDL Ratio 
- glyhb
- Glycosolated Hemoglobin 
- location
- County - a factor with levels - Buckingham- Louisa
- age
- age in years 
- gender
- a factor with levels - male- female
- height
- height in inches 
- weight
- weight in pounds 
- frame
- a factor with levels - small- medium- large
- bp.1s
- First Systolic Blood Pressure 
- bp.1d
- First Diastolic Blood Pressure 
- bp.2s
- Second Systolic Blood Pressure 
- bp.2d
- Second Diastolic Blood Pressure 
- waist
- waist in inches 
- hip
- hip in inches 
- time.ppn
- Postprandial Time (in minutes) when Labs were Drawn 
Details
Glycosolated hemoglobin greater than 7.0 is usually taken as a positive diagnosis of diabetes
Source
Willems JP, Saunders JT, DE Hunt, JB Schorling: Prevalence of coronary heart disease risk factors among rural blacks: A community-based study. Southern Medical Journal 90:814-820; 1997
References
Schorling JB, Roach J, Siegel M, Baturka N, Hunt DE, Guterbock TM, Stewart HL: A trial of church-based smoking cessation interventions for rural African Americans. Preventive Medicine 26:92-101; 1997
Radiation dose effects on chromosomal abnormality
Description
An experiment was conducted to determine the effect of gamma radiation on the numbers of chromosomal abnormalities observed
Format
A data frame with 27 observations on the following 4 variables.
- cells
- Number of cells in hundreds 
- ca
- Number of chromosomal abnormalities 
- doseamt
- amount of dose in Grays 
- doserate
- rate of dose in Grays/hour 
Source
Purott R. and Reeder E. (1976) The effect of changes in dose rate on the yield of chromosome aberrations in human lymphocytes exposed to gamma radiation. Mutation Research. 35, 437-444.
References
Frome E. and DuFrain R. (1986) Maximum Likelihood Estimation for Cytogenic Dose-Response Curves. Biometrics. 42, 73-84.
Divorce in the USA 1920-1996
Description
Divorce rates in the USA from 1920-1996
Format
A data frame with 77 observations on the following 7 variables.
- year
- the year from 1920-1996 
- divorce
- divorce per 1000 women aged 15 or more 
- unemployed
- unemployment rate 
- femlab
- percent female participation in labor force aged 16+ 
- marriage
- marriages per 1000 unmarried women aged 16+ 
- birth
- births per 1000 women aged 15-44 
- military
- military personnel per 1000 population 
Source
Unknown
Choice of drug treatment for psychiatry patients
Description
A sample of psychiatry patients were cross-classified by their diagnosis and whether a drug treatment was prescribed.
Format
A data frame with 10 observations on the following 3 variables.
- y
- the number of patients 
- diagnosis
- a factor with levels - Affective.Disorder- Neurosis- Personality.Disorder- Schizophrenia- Special.Symptoms
- drug
- a factor with levels - no- yes
Source
Helmes E. and Fekken G. (1986) Effects of psychotropic drugs and psychiatric illness on vocational aptitude and interest assessment. J. Clin. Psychol. 42 569-576
References
Agresti A. (1990) "Categorical Data Analysis" Wiley
Doctor visits in Australia
Description
The data come from the Australian Health Survey of 1977-78 and consist of 5190 single adults where young and old have been oversampled.
Format
A data frame with 5190 observations on the following 19 variables.
- sex
- 1 if female, 0 if male 
- age
- Age in years divided by 100 (measured as mid-point of 10 age groups from 15-19 years to 65-69 with 70 or more coded treated as 72) 
- agesq
- age squared 
- income
- Annual income in Australian dollars divided by 1000 (measured as mid-point of coded ranges Nil, less than 200, 200-1000, 1001-, 2001-, 3001-, 4001-, 5001-, 6001-, 7001-, 8001-10000, 10001-12000, 12001-14000, with 14001- treated as 15000 
- levyplus
- 1 if covered by private health insurance fund for private patient in public hospital (with doctor of choice), 0 otherwise 
- freepoor
- 1 if covered by government because low income, recent immigrant, unemployed, 0 otherwise 
- freerepa
- 1 if covered free by government because of old-age or disability pension, or because invalid veteran or family of deceased veteran, 0 otherwise 
- illness
- Number of illnesses in past 2 weeks with 5 or more coded as 5 
- actdays
- Number of days of reduced activity in past two weeks due to illness or injury 
- hscore
- General health questionnaire score using Goldberg's method. High score indicates bad health 
- chcond1
- 1 if chronic condition(s) but not limited in activity, 0 otherwise 
- chcond2
- 1 if chronic condition(s) and limited in activity, 0 otherwise 
- doctorco
- Number of consultations with a doctor or specian the past 2 weeks 
- nondocco
- Number of consultations with non-doctor health professionals (chemist, optician, physiotherapist, social worker, district community nurse, chiropodist or chiropractor) in the past 2 weeks 
- hospadmi
- Number of admissions to a hospital, psychiatric hospital, nursing or convalescent home in the past 12 months (up to 5 or more admissions which is coded as 5) 
- hospdays
- Number of nights in a hospital, etc. during most recent admission: taken, where appropriate, as the mid-point of the intervals 1, 2, 3, 4, 5, 6, 7, 8-14, 15-30, 31-60, 61-79 with 80 or more admissions coded as 80. If no admission in past 12 months then equals zero 
- medicine
- Total number of prescribed and nonprescribed medications used in past 2 days 
- prescrib
- Total number of prescribed medications used in past 2 days 
- nonpresc
- Total number of nonprescribed medications used in past 2 days 
Source
Cameron A, Trivedi P, Milne F and Piggot J (1988) A Microeconometric model of the demand for health care and health insurance in Australia, Review of Economic Studies 55, 85-106
Ecological regression example
Description
Relationship between 1998 per capita income dollars from all sources and the proportion of legal state residents born in the United States in 1990 for each of the 50 states plus the District of Columbia
Format
This dataframe contains the following columns
- usborn
- Percentage of population born in the United States 
- income
- Per capita annual income in dollars 
- home
- Percentage born in state 
- pop
- 
Population of state 
Source
US Bureau of the Census
Treatment and block effects on egg production
Description
The eggprod data frame has 12 rows and 3 columns.  Six pullets were
placed into each of 12 pens. Four blocks were formed from groups of 3 pens
based on location. Three treatments were applied. The number of eggs
produced was recorded
Format
This data frame contains the following columns:
- treat
- Three treatments: O, E or F 
- block
- 
Four blocks labeled 1-4 
- eggs
- Number of eggs produced 
Source
Mead, R., R.N. Curnow, and A.M. Hasted. 1993. Statistical Methods in Agriculture and Experimental Biology. Chapman and Hall, London, p. 64. 1993
Nested data on lab testing of eggs
Description
Consistency between laboratory tests is important and yet the results may depend on who did the test and where the test was performed. In an experiment to test levels of consistency, a large jar of dried egg powder was divided up into a number of samples. Because the powder was homogenized, the fat content of the samples is the same, but this fact is withheld from the laboratories. Four samples were sent to each of six laboratories. Two of the samples were labeled as G and two as H, although in fact they were identical. The laboratories were instructed to give two samples to two different technicians. The technicians were then instructed to divide their samples into two parts and measure the fat content of each. So each laboratory reported eight measures, each technician four measures, that is, two replicated measures on each of two samples.
Format
A data frame with 48 observations on the following 4 variables.
- Fat
- a numeric vector 
- Lab
- a factor with levels - I- II- III- IV- V- VI
- Technician
- a factor with levels - one- two
- Sample
- a factor with levels - G- H
Source
Bliss, C. I. (1967). Statistics in Biology. New York: McGraw Hill.
Examples
data(eggs)
## maybe str(eggs) ; plot(eggs) ...
Epileptic seizures in clinical trial of drug
Description
Data from a clinical trial of 59 epileptics. For a baseline, patients were observed for 8 weeks and the number of seizures recorded. The patients were then randomized to treatment by the drug Progabide (31 patients) or to the placebo group (28 patients). They were observed for four 2-week periods and the number of seizures recorded.
Format
A data frame with 295 observations on the following 6 variables.
- seizures
- number of seizures 
- id
- identifying number 
- treat
- 1=treated, 0=not 
- expind
- 0=baseline period, 1=treatment period 
- timeadj
- weeks of period 
- age
- in years 
Source
Thall, P. F. and S. C. Vail (1990). Some covariance models for longitudinal count data with overdispersion. Biometrics 46, 657-671.
References
Breslow, N. E. and D. G. Clayton (1993). Approximate inference in generalized linear mixed models. Journal of the American Statistical Association 88, 9-25. Diggle, P. J., P. Heagerty, K. Y. Liang, and S. L. Zeger (2002). Analysis of Longitudinal Data (2 ed.). Oxford: Oxford University Press.
Examples
data(epilepsy)
## maybe str(epilepsy) ; plot(epilepsy) ...
Complaints about emergency room doctors
Description
Data was recorded on 44 doctors working in an emergency service at a hospital to study the factors affecting the number of complaints received.
Format
A data frame with 44 observations on the following 6 variables.
- visits
- the number of patient visits 
- complaints
- the number of complaints 
- residency
- is the doctor in residency training - Nor- Y
- gender
- gender of doctor - For- M
- revenue
- dollars per hour earned by the doctor 
- hours
- total number of hours worked 
Source
Chap T. Le (1998) "Applied Categorical Data Analysis" Wiley
Simulated non-parametric regression data
Description
True function is f(x)=sin^3(2pi x^3).
Format
A data frame with 256 observations on the following 3 variables.
- x
- input 
- y
- response 
- m
- true value 
Source
Haerdle, W. (1991). Smoothing Techniques with Implementation in S. New York:Springer.
Examples
data(exa)
## maybe str(exa) ; plot(exa) ...
Simulated non-parametric regression data
Description
True function is f(x)=0
Format
A data frame with 256 observations on the following 3 variables.
- x
- input 
- y
- response 
- m
- true value 
Source
Haerdle, W. (1991). Smoothing Techniques with Implementation in S. New York:Springer.
Examples
data(exa)
## maybe str(exa) ; plot(exa) ...
grading of eye pairs for distance vision
Description
A sample of women are rated for the performance of distance vision in each eye.
Format
A data frame with 16 observations on the following 3 variables.
- y
- the observed count 
- right
- rated vision in the right eye - a factor with levels - best- second- third- worst
- left
- rated vision in the left eye - a factor with levels - best- second- third- worst
Source
A. Stuart (1955) A test for homogeneity of the marginal distributions in a two-way classification, Biometrika, 42, 412-416.
Percentage of Body Fat and Body Measurements in Men
Description
Age, weight, height, and 10 body circumference measurements are recorded for 252 men. Each man's percentage of body fat was accurately estimated by an underwater weighing technique.
Format
A data frame with 252 observations on the following 18 variables.
- brozek
- Percent body fat using Brozek's equation, 457/Density - 414.2 
- siri
- Percent body fat using Siri's equation, 495/Density - 450 
- density
- Density (gm/$cm^3$) 
- age
- Age (yrs) 
- weight
- Weight (lbs) 
- height
- Height (inches) 
- adipos
- Adiposity index = Weight/Height$^2$ (kg/$m^2$) 
- free
- Fat Free Weight = (1 - fraction of body fat) * Weight, using Brozek's formula (lbs) 
- neck
- Neck circumference (cm) 
- chest
- Chest circumference (cm) 
- abdom
- Abdomen circumference (cm) at the umbilicus and level with the iliac crest 
- hip
- Hip circumference (cm) 
- thigh
- Thigh circumference (cm) 
- knee
- Knee circumference (cm) 
- ankle
- Ankle circumference (cm) 
- biceps
- Extended biceps circumference (cm) 
- forearm
- Forearm circumference (cm) 
- wrist
- Wrist circumference (cm) distal to the styloid processes 
Source
Johnson R. Journal of Statistics Education v.4, n.1 (1996)
Mortality due to smoking according age group in women
Description
In 1972-74, a survey of one in six residents of Whickham, near Newcastle, England was made. Twenty years later, this data recorded in a follow-up study. Only women who are current smokers or who have never smoked are included.
Format
A data frame with 28 observations on the following 4 variables.
- y
- observed count for given combination 
- smoker
- a factor with levels - yes- no
- dead
- a factor with levels - yes- no
- age
- a factor with agegroup levels - 18-24- 25-34- 35-44- 45-54- 55-64- 65-74- 75+
Source
D. Appleton, J. French, M. Vanderpump (1996) "Ignoring a Covariate: An Example of Simpson's Paradox" American Statistician, 50, 340-341
Billionaires' wealth and age
Description
Fortune magazine publishes a f the world's billionaires each year. The 1992 list includes 233 individuals. Their wealth, age, and geographic location (Asia, Europe, Middle East, United States, and Other) are reported.
Format
A data frame with 232 observations on the following 3 variables.
- wealth
- Billions of dollars 
- age
- age in years 
- region
- a factor with levels - A, Asia,- E, Europe,- M, Middle East,- OOther,- UUSA
Source
Fortune magazine
Examples
data(fortune)
## maybe str(fortune) ; plot(fortune) ...
1981 French Presidential Election
Description
Elections for the French presidency proceed in two rounds. In 1981, there were 10 candidates in the first round. The top two candidates then went on to the second round, which was won by Francois Mitterand over Valery Giscard-d'Estaing. The losers in the first round can gain political favors by urging their supporters to vote for one of the two fina Since voting is private, we cannot know how these votes were transferred, we might hope to infer from the published vote totals how this might have happened. Data is given for vote totals in every fourth department of France:
Format
This dataframe contains the following columns (vote totals are in thousands)
- list("EI
- Electeur Inscrits (registered voters) 
- A
- Voters for Mitterand in the first round 
- B
- 
Voters for Giscard in the first round 
- C
- Voters for Chirac in the first round 
- D
- Voters for Communists in the first round 
- E
- Voters for Ecology party in the first round 
- F
- Voters for party F in the first round 
- G
- 
Voters for party G in the first round 
- H
- Voters for party H in the first round 
- I
- Voters for party I in the first round 
- J
- Voters for party J in the first round 
- K
- 
Voters for party K in the first round 
- A2
- Voters for Mitterand in the second round 
- B2
- Voters for party Giscard in the second round 
- N
- Difference between the number of voters in the second round and in the first round 
Source
"The Teaching of Practical Statistics" by C.W. Anderson and R.M. Loynes, Wiley,1987
Formating the Rounding of Numbers
Description
fround rounds the values in its first argument to the specified
number of decimal places with surrounding quotes.
Usage
fround(x, digits)
Arguments
| x | a numeric vector. | 
| digits | integer indicating the precision to be used. | 
Details
pfround rounds the values in its first argument to the specified
number of decimal places without surrounding quotes.
Author(s)
Andrew Gelman; Yu-Sung Su
References
Copied from the arm package
See Also
Examples
    x <- 3.1415926
    fround(x, digits=2)
    pfround(x, digits=2)
Longevity of fruitflies depending on sexual activity and thorax length
Description
The fruitfly data frame has 9 rows and 3 columns.  125 fruitflies
were divided randomly into 5 groups of 25 each. The response was the
longevity of the fruitfly in days. One group was kept solitary, while
another was kept individually with a virgin female each day. Another group
was given 8 virgin females per day. As an additional control the fourth and
fifth groups were kept with one or eight pregnant females per day.  Pregnant
fruitflies will not mate. The thorax length of each male was measured as
this was known to affect longevity. One observation in the many group has
been lost.
Format
This data frame contains the following columns:
- thorax
- Thorax length 
- longevity
- Lifetime in days 
- activity
- The group: isolated = fly kept solitary, one = fly kept with one pregnant fruitfly, many = fly kept with eight pregnant fruitflies, low= fly kept with one virgin fruitfly, high = fly kept with eight virgin fruitflies. 
Source
"Sexual Activity and the Lifespan of Male Fruitflies" by L. Partridge and M. Farquhar, Nature, 1981, 580-581
Species diversity on the Galapagos Islands
Description
There are 30 Galapagos islands and 7 variables in the dataset. The
relationship between the number of plant species and several geographic
variables is of interest. The original dataset contained several missing
values which have been filled for convenience. See the galamiss
dataset for the original version.
Format
The dataset contains the following variables
- Species
- the number of plant species found on the island 
- Endemics
- the number of endemic species 
- Area
- 
the area of the island (km$^2$) 
- Elevation
- the highest elevation of the island (m) 
- Nearest
- the distance from the nearest island (km) 
- Scruz
- the distance from Santa Cruz island (km) 
- Adjacent
- the area of the adjacent island (square km) 
Source
M. P. Johnson and P. H. Raven (1973) "Species number and endemism: The Galapagos Archipelago revisited" Science, 179, 893-895
Species diversity on the Galapagos Islands
Description
There are 30 Galapagos islands and 7 variables in the dataset. The relationship between the number of plant species and several geographic variables is of interest. This is the original version of the dataset containing missing values.
Format
The dataset contains the following variables
- Species
- the number of plant species found on the island 
- Endemics
- the number of endemic species 
- Area
- 
the area of the island (km$^2$) 
- Elevation
- the highest elevation of the island (m) 
- Nearest
- the distance from the nearest island (km) 
- Scruz
- the distance from Santa Cruz island (km) 
- Adjacent
- the area of the adjacent island (square km) 
Source
M. P. Johnson and P. H. Raven (1973) "Species number and endemism: The Galapagos Archipelago revisited" Science, 179, 893-895
Xray decay from a gamma ray burst
Description
The X-ray decay light curve of Gamma ray burst 050525a obtained with the X-Ray Telescope (XRT) on board the Swift satellite. The dataset has 63 brightness measurements in the 0.4-4.5 keV spectral band at times ranging from 2 minutes to 5 days after the burst.
Format
A data frame with 63 observations on the following 3 variables.
- time
- in seconds since burst 
- flux
- X-ray flux in units of 10^-11 erg/cm2/s, 2-10 keV 
- error
- measurement error of the flux based on detector signal-to-noise values 
Source
A. J. Blustin and 64 coauthors, Astrophys. J. 637, 901-913 2006. Available at http://arxiv.org/abs/astro-ph/0507515.
Examples
data(gammaray)
## maybe str(gammaray) ; plot(gammaray) ...
Undercounted votes in Georgia in 2000 presidential election
Description
The data comes from the US presidential election in the state of Georgia. The undercount is the difference between the number of ballots cast and votes recorded. Voters may have chosen not to vote for president, voted for more than one candidate (disqualified) or the equipment may have failed to register their choice.
Format
A data frame with 159 observations on the following 10 variables. Each case represents a county in Georgia.
- equip
- The voting equipment used: - LEVER,- OS-CC(optical, central count),- OS-PC(optical, precinct count)- PAPER,- PUNCH
- econ
- economic status of county: - middle- poor- rich
- perAA
- percent of African Americans in county 
- rural
- indicator of whether county is - ruralor- urban
- atlanta
- indicator of whether county is in - Atlantaor not:- notAtlanta
- gore
- number of votes for Gore 
- bush
- number of votes for Bush 
- other
- number of votes for other candidates 
- votes
- number of votes 
- ballots
- number of ballots 
Source
Meyer M. (2002) Uncounted Votes: Does Voting Equipment Matter? Chance, 15(4), 33-38
Northern Hemisphere temperatures and climate proxies in the last millenia
Description
Average Northen Hemisphere Temperature from 1856-2000 and eight climate proxies from 1000-2000AD. Data can be used to predict temperatures prior to 1856.
Format
A data frame with 1001 observations on the following 10 variables.
- nhtemp
- Northern hemisphere average temperature (C) provided by the UK Met Office (known as HadCRUT2) 
- wusa
- Tree ring proxy information from the Western USA. 
- jasper
- Tree ring proxy information from Canada. 
- westgreen
- Ice core proxy information from west Greenland 
- chesapeake
- Sea shell proxy information from Chesapeake Bay 
- tornetrask
- Tree ring proxy information from Sweden 
- urals
- Tree ring proxy information from the Urals 
- mongolia
- Tree ring proxy information from Mongolia 
- tasman
- Tree ring proxy information from Tasmania 
- year
- Year 1000-2000AD 
Details
See the source and references below for the original data. Only some proxies have been included here. Some missing values have been imputed. The proxy data have been smoothed. This version of the data is intended only for demonstration purposes. If you are specifically interested in the subject matter, use the original data.
Source
P.D. Jones and M.E. Mann (2004) "Climate Over Past Millennia" Reviews of Geophysics, Vol. 42, No. 2, RG2002, doi:10.1029/2003RG000143
References
www.ncdc.noaa.gov/paleo/pubs/jones2004/jones2004.html
Examples
data(globwarm)
## maybe str(globwarm) ; plot(globwarm) ...
Hair and eye color
Description
Data collected from 592 students in an introductory statistics class
Format
A data frame with 16 observations on the following 3 variables.
- y
- count of the number of student with given hair/eye combination 
- eye
- a factor with levels - green- hazel- blue- brown
- hair
- a factor with levels - BLACK- BROWN- RED- BLOND
Source
Snee R. (1974) Graphical display of two-way contingency tables. American Statistician, 28, 9-12
Half Normal Plot
Description
Makes a half-normal plot
Usage
halfnorm(
  x,
  nlab = 2,
  labs = as.character(1:length(x)),
  ylab = "Sorted Data",
  ...
)
Arguments
| x | a numeric vector | 
| nlab | number of points to label | 
| labs | labels for points | 
| ylab | label for Y-axis | 
| ... | arguments passed to plot() | 
Value
none
Author(s)
Julian Faraway
See Also
qqnorm
Examples
halfnorm(runif(10))
MBA students experience love, sex, work and happiness
Description
Data were collected from 39 students in a University of Chicago MBA class
Usage
mba
Format
A data frame with 39 observations on the following 5 variables.
- happy
- Happiness on a 10 point scale where 10 is most happy 
- money
- family income in thousands of dollars 
- sex
- 1 = satisfactory sexual activity, 0 = not 
- love
- 1 = lonely, 2 = secure relationships, 3 = deep feeling of belonging and caring 
- work
- 5 point scale where 1 = no job, 3 = OK job, 5 = great job 
An object of class data.frame with 39 rows and 5 columns.
Source
George and McCulloch (1993) "Variable Selection via Gibbs Sampling" JASA, 88, 881-889
Treatment of insulin dependent diabetic children
Description
16 insulin-dependent diabetic children were enrolled in a study involving a new treatment. 8 children received the new treatment(N) while the other 8 received the standard treatment(S). The age and sex of the child was recorded along with the measured value of gycosolated hemoglobin both before and after treatment.
Format
A data frame with 16 observations on the following 5 variables.
- age
- age in years 
- sex
- a factor with levels - F- M
- treatment
- a factor with levels - N- S
- pre
- measured value of hemoglobin before treatment 
- post
- measured value of hemoglobin after treatment 
Source
Unknown
Examples
data(hemoglobin)
## maybe str(hemoglobin) ; plot(hemoglobin) ...
Ankylosing Spondylitis
Description
Data from Royal Mineral Hospital in Bath. AS is a chronic form of arthritis. A study conducted to determine whether daily stretching of the hip tissues would improve mobility. 39 “typical” AS patients were randomly allocated to control (standard treatment) group or the treatment group in a 1:2 ratio. Responses were flexion and rotation angles at the hip measured in degrees. Larger numbers indicate more flexibility.
Format
A data frame with 78 observations on the following 7 variables.
- fbef
- flexion angle before 
- faft
- flexion angle after 
- rbef
- rotation angle before 
- raft
- rotation angle after 
- grp
- treatment group - a factor with levels - control- treat
- side
- side of the body - a factor with levels - right- left
- person
- id for the individual 
Source
Chatfield C. (1995) Problem Solving: A Statistician's Guide, 2ed Chapman Hall.
Examples
data(hips)
## maybe str(hips) ; plot(hips) ...
Hormone concentrations in gay and straight men
Description
Urinary androsterone (androgen) and etiocholanolone (estrogen) values were recorded from 26 healthy males.
Format
A data frame with 26 observations on the following 3 variables.
- androgen
- concentration 
- estrogen
- concentration 
- orientation
- sexual orientation with levels - g- s
Source
Margolese, M. (1970). Homosexuality: A new endocrine correlate. Hormones and Behavior 1, 151-155.
References
Hand, D. (1981). Discrimination and Classification. Chichester, UK: Wiley.
Examples
data(hormone)
## maybe str(hormone) ; plot(hormone) ...
Housing prices in US cities 86-94
Description
Data on housing prices in 36 US metropolitan statistical areas (MSAs) over 9 years from 1986-1994 were collected.
Format
A data frame with 324 observations on the following 8 variables.
- narsp
- natural log average sale price in thousands of dollars 
- ypc
- average per capita income 
- perypc
- percentage growth in per capita income 
- regtest
- Regulatory environment index (high values = more regulations) 
- rcdum
- Rent control - a factor with levels - 0=no- 1=yes
- ajwtr
- Adjacent to a coastline - a factor with levels - 0=no- 1=yes
- msa
- indicator for the MSA 
- time
- Year 1=1986 to 9=1994 
Source
Longitudinal and Panel Data: Analysis and Applications in the Social Sciences, by Edward W. Frees, Cambridge University Press, August 2004.
Career choice of high school students
Description
Data was collected as a subset of the "High School and Beyond" study conducted by the National Education Longitudinal Studies (NELS) program of the National Center for Education Statistics (NCES).
Format
A data frame with 200 observations on the following 11 variables.
- id
- ID of student 
- gender
- a factor with levels - female- male
- race
- a factor with levels - african-amer- asian- hispanic- white
- ses
- socioeconomic class - a factor with levels - high- low- middle
- schtyp
- school type - a factor with levels - private- public
- prog
- choice of high school program - a factor with levels - academic- general- vocation
- read
- reading score 
- write
- writing score 
- math
- math score 
- science
- science score 
- socst
- social science score 
Details
One purpose of the study was to determine which factors are related to the choice of the type of program, academic, vocational or general, that the students pursue in high school.
Source
National Education Longitudinal Studies (NELS) program of the National Center for Education Statistics (NCES).
Inverse Logit Transformation
Description
Computes the inverse logit transformation
Usage
ilogit(x)
Arguments
| x | a numeric vector | 
Value
exp(x)/(1+exp(x))
Author(s)
Julian Faraway
See Also
logit
Examples
ilogit(1:3)
#[1] 0.7310586 0.8807971 0.9525741
Infant mortality according to income and region
Description
The infmort data frame has 105 rows and 4 columns.  The infant
mortality in regions of the world may be related to per capita income and
whether oil is exported. The dataset is not recent.
Format
This data frame contains the following columns:
- region
- Region of the world, Africa, Europe, Asia or the Americas 
- income
- Per capita annual income in dollars 
- mortality
- Infant mortality in deaths per 1000 births 
- oil
- Does the country export oil or not? 
Source
Unknown
Effects of insulation on gas consumption
Description
Data on natural gas usage in a house. The weekly gas consumption (in 1000 cubic feet) and the average outside temperature (in degrees Celsius) was recorded for 26 weeks before and 30 weeks after cavity-wall insulation had been installed. The house thermostat was set at 20C throughout.
Format
A data frame with 44 observations on the following 3 variables.
- Insulate
- a factor with levels - After- Before
- Temp
- Outside temperature 
- Gas
- Weekly consumption in 1000 cubic feet 
Source
MASS package as whiteside
Examples
data(insulgas)
## maybe str(insulgas) ; plot(insulgas) ...
Irrigation methods in an agricultural field trial
Description
In an agricultural field trial, the objective was to determine the effects of two crop varieties and four different irrigation methods. Eight fields were available, but only one type of irrigation may be applied to each field. The fields may be divided into two parts with a different variety planted in each half. The whole plot factor is the method of irrigation, which should be randomly assigned to the fields. Within each field, the variety is randomly assigned.
Format
A data frame with 16 observations on the following 4 variables.
- field
- a factor with levels - f1- f2- f3- f4- f5- f6- f7- f8
- irrigation
- a factor with levels - i1- i2- i3- i4
- variety
- a factor with levels - v1- v2
- yield
- a numeric vector 
Source
Found online but source not recorded.
Examples
data(irrigation)
## maybe str(irrigation) ; plot(irrigation) ...
Junior School Project
Description
Junior School Project collected from primary (U.S. term is elementary) schools in inner London.
Format
A data frame with 3236 observations on the following 9 variables.
- school
- 50 schools code 1-50 
- class
- a factor with levels - 1- 2- 3- 4
- gender
- a factor with levels - boy- girl
- social
- class of the father I=1; II=2; III nonmanual=3; III manual=4; IV=5; V=6; Long-term unemployed=7; Not currently employed=8; Father absent=9 
- raven
- test score 
- id
- student id coded 1-1402 
- english
- score on English 
- math
- score on Maths 
- year
- year of school 
Source
Mortimore, P., P. Sammons, L. Stoll, D. Lewis, and R. Ecob (1988). School Matters. Wells, UK: Open Books.
References
Goldstein, H. (1995). Multilevel Statistical Models (2 ed.). London: Arnold.
Examples
data(jsp)
## maybe str(jsp) ; plot(jsp) ...
Kangaroo skull measurements
Description
Sex and species of an specimens of kangaroo.
Format
A data frame with 148 observations on the following 20 variables.
- species
- a factor with levels - fuliginosus- giganteus- melanops
- sex
- a factor with levels - Female- Male
- basilar.length
- a numeric vector 
- occipitonasal.length
- a numeric vector 
- palate.length
- a numeric vector 
- palate.width
- a numeric vector 
- nasal.length
- a numeric vector 
- nasal.width
- a numeric vector 
- squamosal.depth
- a numeric vector 
- lacrymal.width
- a numeric vector 
- zygomatic.width
- a numeric vector 
- orbital.width
- a numeric vector 
- rostral.width
- a numeric vector 
- occipital.depth
- a numeric vector 
- crest.width
- a numeric vector 
- foramina.length
- a numeric vector 
- mandible.length
- a numeric vector 
- mandible.width
- a numeric vector 
- mandible.depth
- a numeric vector 
- ramus.height
- a numeric vector 
Source
Andrews and Herzberg (1985) Chapter 53.
References
Andrews, D. F. and Herzberg, A. M. (1985). Data. Springer-Verlag, New York.
Examples
data(kanga)
## maybe str(kanga) ; plot(kanga) ...
Cut-off times of lawnmowers
Description
Data on the cut-off times of lawnmowers was collected. 3 machines were randomly selected from those produced by manufacturers A and B. Each machine was tested twice at low speed and high speed.
Format
A data frame with 24 observations on the following 4 variables.
- manufact
- Manufacturer - a factor with levels - A- B
- machine
- Lawn mower - a factor with levels - m1- m2- m3- m4- m5- m6
- speed
- Speed of testing - a factor with levels - H- L
- time
- cut-off time 
Source
Unknown.
Leaf blotch on barley
Description
The data gives the proportion of leaf area affected by leaf blotch on 10 varieties of barley at 9 different sites.
Format
A data frame with 90 observations on the following 3 variables.
- blotch
- proportion of the barley leaf affected by blotch 
- site
- the physical location - a factor with levels - 1- 2- 3- 4- 5- 6- 7- 8- 9
- variety
- variety of barley - a factor with levels - 1- 2- 3- 4- 5- 6- 7- 8- 9- 10
Source
R. W. M. Wedderburn (1974) "Quasilikelihood functions, generalized linear models and the Gauss-Newton method" Biometrika, 61, 439-447.
References
P. McCullagh and J. Nelder (1989) "Generalized Linear Models" Chapman and Hall, 2nd ed.
Data on the burning time of samples of tobacco leaves
Description
Data on the burning time of samples of tobacco leaves
Format
A data frame with 30 observations on the following 4 variables.
- nitrogen
- nitrogen content by percentage weight 
- chlorine
- chlorine content by percentage weight 
- potassium
- potassium content by percentage weight 
- burntime
- burn time in seconds 
Source
Steel, R. G. D. and Torrie, J. H. (1980), Principles and Procedures of Statistics, Second Edition, New York: McGraw-Hill
Logit transformation
Description
Computes the logit transformation
Usage
logit(x)
Arguments
| x | a numeric vector | 
Details
x <=0 or >=1 will return NA
Value
log(x/(1-x))
Author(s)
Julian Faraway
See Also
ilogit
Examples
logit(c(0.1,0.5,1.0,1.1))
#[1] -2.197225  0.000000        NA        NA
Sleep in Mammals: Ecological and Constitutional Correlates
Description
The mammalsleep data frame has 62 rows and 10 columns. Sleep in
Mammals: Ecological and Constitutional Correlates
Format
This data frame contains the following columns:
- body
- body weight in kg 
- brain
- brain weight in g 
- nondream
- slow wave ("nondreaming") sleep (hrs/day) 
- dream
- paradoxical ("dreaming") sleep (hrs/day) 
- sleep
- total sleep (hrs/day) (sum of slow wave and paradoxical sleep) 
- lifespan
- maximum life span (years) 
- gestation
- gestation time (days) 
- predation
- 
predation index (1-5) 1 = minimum (least likely to be preyed upon) to 5 = maximum (most likely to be preyed upon) 
- exposure
- sleep exposure index (1-5) 1 = least exposed (e.g. animal sleeps in a well-protected den) 5 = most exposed 
- danger
- overall danger index (1-5) (based on the above two indices and other information) 1 = least danger (from other animals) 5 = most danger (from other animals) 
Source
"Sleep in Mammals: Ecological and Constitutional Correlates" by Allison, T. and Cicchetti, D. (1976), Science, November 12, vol. 194, pp. 732-734.
Mayer's 1750 data on the Manilius crater on the moon
Description
In 1750, Tobias Mayer collected data on various landmarks on the moon in order to determine its orbit. The data involving the position of the Manilius crater resulted in a least squares like problem. The example is discussed in Steven Stigler's History of Statistics.
Format
A data frame with 27 observations on the following 4 variables.
- arc
- an angle known as h in Stigler's notation 
- sinang
- the sin(g-k) where g and k are two angles in Stigler 
- cosang
- the cos(g-k) where g and k are two angles in Stigler 
- group
- one of three groups determined by Mayer 
Details
See Stigler for a detailed description.
Source
Stigler, S. (1986) History of Statistics. Belknap Press, Harvard.
References
Mayer, T. (1750) Abhandlung uber die Umwaltzung des Monds um seine Axe und die scheinbare Bewegung der Mondsflecken published in the Kosmographische Nachrichten und Sammlungen auf das Jahr 1748. 52-183
Examples
data(manilius)
Maximum Adjusted R-squared
Description
Displays the best models from a leaps object
Usage
maxadjr(l, best = 3)
Arguments
| l | A leaps object returned from leaps() | 
| best | An optional argument specify the number of models to be returned taking the default value of 3 | 
Details
Requires leaps package
Value
A list of the best models
Author(s)
Julian Faraway
See Also
leaps()
Meat spectrometry to determine fat content
Description
A Tecator Infratec Food and Feed Analyzer working in the wavelength range 850 - 1050 nm by the Near Infrared Transmission (NIT) principle was used to collect data on samples of finely chopped pure meat. 215 samples were measured. For each sample, the fat content was measured along with a 100 channel spectrum of absorbances. Since determining the fat content via analytical chemistry is time consuming we would like to build a model to predict the fat content of new samples using the 100 absorbances which can be measured more easily.
Format
Dataset contains the following variables
- V1-V100
- absorbances across a range of 100 wavelengths 
- fat
- fat content 
Source
H. H. Thodberg (1993) "Ace of Bayes: Application of Neural Networks With Pruning", report no. 1132E, Maglegaardvej 2, DK-4000 Roskilde, Danmark
Melanoma by type and location
Description
Data comes from a study of Malignant Melanoma involving 400 subjects.
Format
A data frame with 12 observations on the following 3 variables.
- count
- number of cases 
- tumor
- type of tumor - a factor with levels - freckle- indeterminate- nodular- superficial
- site
- location of tumor on the body - a factor with levels - extremity- head- trunk
Source
Dobson A. (2002) An introduction to generalized linear models, Chapman Hall.
Third party motor insurance claims in Sweden in 1977
Description
In Sweden all motor insurance companies apply identical risk arguments to classify customers, and thus their portfolios and their claims statistics can be combined. The data were compiled by a Swedish Committee on the Analysis of Risk Premium in Motor Insurance. The Committee was asked to look into the problem of analyzing the real influence on claims of the risk arguments and to compare this structure with the actual tariff.
Format
A data frame with 1797 observations on the following 8 variables.
- Kilometres
- an ordered factor representing kilomoters per year with levels 1: < 1000, 2: 1000-15000, 3: 15000-20000, 4: 20000-25000, 5: > 25000 
- Zone
- a factor representing geographical area with levels 1: Stockholm, Goteborg, Malmo with surroundings 2: Other large cities with surroundings 3: Smaller cities with surroundings in southern Sweden 4: Rural areas in southern Sweden 5: Smaller cities with surroundings in northern Sweden 6: Rural areas in northern Sweden 7: Gotland 
- Bonus
- No claims bonus. Equal to the number of years, plus one, since last claim 
- Make
- A factor representing eight different common car models. All other models are combined in class 9 
- Insured
- Number of insured in policy-years 
- Claims
- Number of claims 
- Payment
- Total value of payments in Skr 
- perd
- payment per claim 
Source
http://www.statsci.org/data/general/motorins.html
References
Hallin, M., and Ingenbleek, J.-F. (1983). The Swedish automobile portfolio in 1977. A statistical study. Scandinavian Actuarial Journal, 49-64.
Questionnaire study of neighborly help
Description
Subjects were asked questions in a study of neighborly help. Questions below are a subset of the full study.
Format
A data frame with 181 observations on the following 8 variables.
- longlive
- About how long have you lived where you do now? Ans is a factor with levels - <6mos- 6-12mos- 1-3yrs- 3-10yrs- 10yrs
- wherebfr
- Where were you living before you moved to your present house? Ans is a factor with levels - same- Exeter- Devon- Britain- Abroad
- hownbly
- How neighborly do you think the area where you now live is? Ans is a factor with levels - Vunfriendly- NVfriendly- Average- FFriendly- VFriendly
- knowname
- Roughly how many people in your street, or in the streets just near you, do you know the names of? Ans is a factor with levels - none- 1-5- 6-20- 20+
- callname
- How many of those people (not counting children) would you call by their first names? Ans is a factor with levels - none- 1-5- 6-20- 20+
- age
- a factor with levels - -18- 18-30- 31-50- 51-65- 65+
- district
- a factor with levels - 1- 2- 3- 4
- sex
- a factor with levels - female- male
Details
Exeter is a city in the county of Devon which is in Britain. The four districts can be briefly described as follows. District 1 was a long-established residential area near the city centre, with housing dating from the late nineteenth century. Originally working class, it now has a considerable middle class population with some student and other temporary accommodation. District 2 was a working-class housing estate dating from the 1930s, with mainly rented accommodation but some owner occupation. District 3 was the oldest part of a more recently developed, mainly middle-class, almost exclusively owner-occupied estate, dating from the 1960s. District 4 was the most recently developed part of a more sought-after middle-class residential area, with smaller but almost entirely owner-occupied properties dating from the 1970s and 1980s.
Source
P. Webley & S. Lea 1993, Human Relations 46, 65-76.
National Education Longitudinal Study of 1988
Description
A subset of the National Education Longitudinal Study of 1988
Format
A data frame with 260 observations on the following 5 variables.
- sex
- a factor with levels - Female- Male
- race
- a factor with levels - White- Asian- Black- Hispanic
- ses
- a numeric vector 
- paredu
- a factor with levels - ba- college- hs- lesshs- ma- phd
- math
- a numeric vector 
Source
http://www.icpsr.umich.edu/icpsrweb/ICPSR/series/107
Examples
data(nels88)
## maybe str(nels88) ; plot(nels88) ...
Nepali child heath study
Description
The data are a subset from public health study on Nepalese children.
Format
A data frame with 1000 observations on the following 9 variables.
- id
- There is a six digit code for the child's ID: 2 digits for the panchayat number; 2 digits for the ward within panchayat; 1 digits for the household; 1 digit for child within household. 
- sex
- 1 = male; 2 = female 
- wt
- Child's weight measured in kilograms 
- ht
- Child's height measured in centimeters 
- mage
- Mother's age in years 
- lit
- Indicator of mother's literacy: 0 = no; 1 = yes 
- died
- The number of children the mother has had that died. 
- alive
- The number of children the mother has ever had born alive 
- age
- age of child 
Source
West KP, Jr., LeClerq SC, Shrestha SR, Wu LS, Pradhan EK, Khatry SK, Katz J, Adhikari R, Sommer A. Effects of vitamin A on growth of vitamin A deficient children: field studies in Nepal. J Nutr 1997;10:1957-1965.
US 1996 national election study
Description
10 variable subset of the 1996 American National Election Study. Missing values and "don't know" responses have been se deleted. Respondents expressing a voting preference other than Clinton or Dole have been removed.
Format
A data frame with 944 observations on the following 10 variables.
- popul
- population of respondent's location in 1000s of people 
- TVnews
- days in the past week spent watching news on TV 
- selfLR
- Left-Right self-placement of respondent: an ordered factor with levels extremely liberal, - extLib< liberal,- Lib< slightly liberal,- sliLib< moderate,- Mod< slightly conservative,- sliCon< conservative,- Con< extremely conservative,- extCon
- ClinLR
- Left-Right placement of Bill Clinton (same scale as selfLR): an ordered factor with levels - extLib<- Lib<- sliLib<- Mod<- sliCon<- Con<- extCon
- DoleLR
- Left-Right placement of Bob Dole (same scale as selfLR): an ordered factor with levels - extLib<- Lib<- sliLib<- Mod<- sliCon<- Con<- extCon
- PID
- Party identification: an ordered factor with levels strong Democrat, - strDem< weak Democrat,- weakDem< independent Democrat,- indDem< independent independent- indind< indepedent Republican,- indRep< waek Republican,- weakRep< strong Republican,- strRep
- age
- Respondent's age in years 
- educ
- Respondent's education: an ordered factor with levels 8 years or less, - MS< high school dropout,- HSdrop< high school diploma or GED,- HS< some College,- Coll< Community or junior College degree,- CCdeg< BA degree,- BAdeg< postgraduate degree,- MAdeg
- income
- Respondent's family income: an ordered factor with levels - $3Kminus<- $3K-$5K<- $5K-$7K<- $7K-$9K<- $9K-$10K<- $10K-$11K<- $11K-$12K<- $12K-$13K<- $13K-$14K<- $14K-$15K<- $15K-$17K<- $17K-$20K<- $20K-$22K<- $22K-$25K<- $25K-$30K<- $30K-$35K<- $35K-$40K<- $40K-$45K<- $45K-$50K<- $50K-$60K<- $60K-$75K<- $75K-$90K<- $90K-$105K<- $105Kplus
- vote
- Expected vote in 1996 presidential election: a factor with levels - Clintonand- Dole
Source
Sapiro, Virginia, Steven J. Rosenstone, Donald R. Kinder, Warren E. Miller, and the National Election Studies. AMERICAN NATIONAL ELECTION STUDIES, 1992-1997: COMBINED FILE [Computer file]. 2nd ICPSR version. Ann Arbor, MI: University of Michigan, Center for Political Studies [producer], 1999. Ann Arbor, MI: Inter-university Consortium for Political and Social Research [distributor], 1999.
References
Found at http://www.stat.washington.edu/
New Hampshire Democratic Party Primary 2008
Description
Votes and other demographic information from 276 wards in the 2008 Democratic Party presidential primary.
Format
A data frame with 276 observations on the following 12 variables.
- votesys
- The voting system used where H is counted by hand and D is counted by machine. 
- Obama
- The number of votes for Barack Obama. 
- Clinton
- The number of votes for Hillary Clinton. 
- dem
- The total number of votes cast in the Democratic primary (there were other candidates besides Clinton and Obama). 
- povrate
- The poverty rate as a proportion as determined by the 2000 census. 
- pci
- Per capita annual income in USD in 1999. 
- Dean
- The proportion of voters for Howard Dean in the 2004 Democratic primary. 
- Kerry
- The proportion of voters for John Kerry in the 2004 Democratic primary. 
- white
- The proportion of non-Hispanic whites according to the 2000 census. 
- absentee
- The proportion voting by absentee ballot. 
- population
- An estimate of the population from 2002. 
- pObama
- Proportion voting for Obama 
Details
On the 8th January 2008, primaries to select US presidential candidates were held in New Hampshire. In the Democratic party primary, Hillary Clinton defeated Barack Obama contrary to the expectations pre-election opinion polls. Essentially two different voting technologies were used in New Hampshire. Some wards used paper ballots, counted by hand while others used optically scanned ballots, counted by machine. Among the paper ballots, Obama had more votes than Clinton while Clinton defeated Obama on just the machine counted ballots. Since the method of voting should make no causal difference to the outcome, suspicions have been raised regarding the integrity of the election.
Source
Herron, M., W. M. Jr, and J. Wand (2008). Voting Technology and the 2008 New Hampshire Primary. Wm. & Mary Bill Rts. J. 17, 351-374.
Yields of oat varieties planted in blocks
Description
Data from an experiment to compare 8 varieties of oats. The growing area was heterogeneous and so was grouped into 5 blocks. Each variety was sown once within each block and the yield in grams per 16ft row was recorded.
Format
The dataset contains the following variables
- yield
- Yield in grams per 16ft row 
- block
- 
Blocks I to V 
- variety
- Variety 1 to 8 
Source
"Statistical Theory in Research" by R. Anderson and T. Bancroft, McGraw Hill,1952
Odor of chemical by production settings
Description
Data from an experiment to determine the effects of column temperature, gas/liquid ratio and packing height in reducing unpleasant odor of chemical product that was being sold for household use
Format
- odor
- Odor score 
- temp
- 
Temperature coded as -1, 0 and 1 
- gas
- Gas/Liquid ratio coded as -1, 0 and 1 
- pack
- Packing height coded as -1, 0 and 1 
Source
"Statistical Design and Analysis of Experiments" by P. John, Macmillan, 1971
Ohio Children Wheeze Status
Description
The ohio data frame has 2148 rows and 4 columns. The dataset is a
subset of the six-city study, a longitudinal study of the health effects of
air pollution.
Format
This data frame contains the following columns:
- resp
- an indicator of wheeze status (1=yes, 0=no) 
- id
- a numeric vector for subject id 
- age
- a numeric vector of age, 0 is 9 years old 
- smoke
- an indicator of maternal smoking at the first year of the study 
References
Fitzmaurice, G.M. and Laird, N.M. (1993) A likelihood-based method for analyzing longitudinal binary responses, Biometrika 80: 141–151.
Spache Shuttle Challenger O-rings
Description
The 1986 crash of the space shuttle Challenger was linked to failure of O-ring seals in the rocket engines. Data was collected on the 23 previous shuttle missions. The launch temperature on the day of the crash was 31F.
Format
A data frame with 23 observations on the following 2 variables.
- temp
- temperature at launch in degrees F 
- damage
- number of damage incidents out of 6 possible 
Source
Presidential Commission on the Space Shuttle Challenger Accident, Vol. 1, 1986: 129-131.
References
S. Dalal, E. Fowlkes and B. Hoadley (1989) "Risk Analysis of the Space Shuttle: Pre-Challenger Prediction of Failure." Journal of the American Statistical Association. 84: 945-957.
Ozone in LA in 1976
Description
A study the relationship between atmospheric ozone concentration and meteorology in the Los Angeles Basin in 1976. A number of cases with missing variables have been removed for simplicity.
Format
A data frame with 330 observations on the following 10 variables.
- O3
- Ozone conc., ppm, at Sandbug AFB. 
- vh
- a numeric vector 
- wind
- wind speed 
- humidity
- a numeric vector 
- temp
- temperature 
- ibh
- inversion base height 
- dpg
- Daggett pressure gradient 
- ibt
- a numeric vector 
- vis
- visibility 
- doy
- day of the year 
Source
Breiman, L. and J. H. Friedman (1985). Estimating optimal transformations for multiple regression and correlation. Journal of the American Statistical Association 80, 580-598.
Examples
data(ozone)
## maybe str(ozone) ; plot(ozone) ...
Marijuana and parent alcohol and drug use
Description
445 college students were classified according to both frequency of marijuana use and parental use of alcohol and psychoactive drugs.
Format
A data frame with 9 observations on the following 3 variables.
- parent
- Number of parents using drugs or alcohol - a factor with levels - Both- Neither- One
- student
- Student usage of marijuana - a factor with levels - Never- Occasional- Regular
- count
- the number of cases 
Source
Ellis, Godfrey J. and Stone, Lorene H. (1979) Marijuana Use in College: "An Evaluation of a Modeling Explanation" Youth and Society 10, 323-34
Carbon dioxide effects on peanut oil extraction
Description
The peanut data frame has 16 rows and 6 columns. Carbon dioxide
effects on peanut oil extraction
Format
This data frame contains the following columns:
- press
- CO2 pressure - two levels low=0, high=1 
- temp
- CO2 temperature - two levels low=0, high=1 
- moist
- peanut moisture - two levels low=0, high=1 
- flow
- CO2 flow rate - two levels low=0, high=1 
- size
- peanut particle size - two levels low=0, high=1 
- solubility
- the amount of oil that could dissolve in the CO2 
Source
Kilgo, M (1989) "An Application of Fractional Factorial Experimental Designs" Quality Engineering, 1, 45-54
Penicillin yield by block and treatment
Description
The production of penicillin uses a raw material, corn steep liquor, is quite variable and can only be made in blends sufficient for four runs. There are four processes, A, B, C and D, for the production.
Format
A data frame with 20 observations on the following 3 variables.
- treat
- a factor with levels - A- B- C- D
- blend
- a factor with levels - Blend1- Blend2- Blend3- Blend4- Blend5
- yield
- a numeric vector 
Source
Box, G., W. Hunter, and J. Hunter (1978). Statistics for Experimenters. New York: Wiley.
Examples
data(penicillin)
## maybe str(penicillin) ; plot(penicillin) ...
Birth weights in Philadelphia
Description
Data based on a 5
Format
A data frame with 1115 observations on the following 5 variables.
- black
- is the mother Black? 
- educ
- mother's years of education 
- smoke
- does the mother smoke during pregnancy? 
- gestate
- gestational age in weeks 
- grams
- birth weight in grams 
Source
I. T. Elo, G. Rodriguez and H. Lee (2001). Racial and Neighborhood Disparities in Birthweight in Philadelphia. Paper presented at the Annual Meeting of the Population Association of America, Washington, DC 2001.
Examples
data(phbirths)
## maybe str(phbirths) ; plot(phbirths) ...
Diabetes survey on Pima Indians
Description
The National Institute of Diabetes and Digestive and Kidney Diseases conducted a study on 768 adult female Pima Indians living near Phoenix.
Format
The dataset contains the following variables
- pregnant
- Number of times pregnant 
- glucose
- 
Plasma glucose concentration at 2 hours in an oral glucose tolerance test 
- diastolic
- Diastolic blood pressure (mm Hg) 
- triceps
- Triceps skin fold thickness (mm) 
- insulin
- 2-Hour serum insulin (mu U/ml) 
- bmi
- 
Body mass index (weight in kg/(height in metres squared)) 
- diabetes
- Diabetes pedigree function 
- age
- Age (years) 
- test
- test whether the patient shows signs of diabetes (coded 0 if negative, 1 if positive) 
Source
The data may be obtained from UCI Repository of machine learning databases at http://archive.ics.uci.edu/ml/
NIST data on ultrasonic measurements of defects in the Alaska pipeline
Description
Researchers at National Institutes of Standards and Technology (NIST) collected data on ultrasonic measurements of the depths of defects in the Alaska pipeline in the field. The depth of the defects were then remeasured in the laboratory. These measurements were performed in six different batches. The laboratory measurements are more accurate than the in-field measurements, but more time consuming and expensive.
Format
A data frame with 107 observations on the following 3 variables.
- Field
- measurement of depth of defect on site 
- Lab
- measurement of depth of defect in the lab 
- Batch
- the batch of measurements 
Source
Office of the Director of the Institute of Materials Research (now the Materials Science and Engineering Laboratory) of NIST
Pneumonoconiosis in coal miners
Description
The data for this example contains the number of coal miners classified by radiological examination into one of three categories of pneumonoultramicroscopicosilicovolcanoconiosis (known as pneumonoconiosis for short) and by number of years spent working at the coal face divided into eight categories.
Format
A data frame with 24 observations on the following 3 variables.
- Freq
- number of miners 
- status
- pneumoconiosis status - a factor with levels - mild- normal- severe
- year
- number of years service (midpoint of interval) 
Source
M. Aitkin and D. Anderson and B. Francis and J. Hinde (1989) "Statistical Modelling in GLIM" Oxford University Press.
Marijuana usage by youth
Description
The National Youth Survey collected a sample of 11 to 17 year olds - 117 boys and 120 girls - asking questions about marijuana usage.
Format
A data frame with 486 observations on the following 7 variables.
- sex
- 1=Male, 2=Female 
- year.76
- 1=never used, 2=used no more than once a month, 3=used more than once a month in 1976 
- year.77
- 1=never used, 2=used no more than once a month, 3=used more than once a month in 1977 
- year.78
- 1=never used, 2=used no more than once a month, 3=used more than once a month in 1978 
- year.79
- 1=never used, 2=used no more than once a month, 3=used more than once a month in 1979 
- year.80
- 1=never used, 2=used no more than once a month, 3=used more than once a month in 1980 
- count
- Number of cases in this category 
Source
ICPSR, University of Michigan
References
Lang J., McDonald, J and Smith P. (1999) "Association-Marginal Modeling of Mutlivariate Categorical Responses: A Maximum Likelihood Approach" JASA 94, 1161-
Prostate cancer surgery
Description
The prostate data frame has 97 rows and 9 columns. A study on 97 men
with prostate cancer who were due to receive a radical prostatectomy.
Format
This data frame contains the following columns:
- lcavol
- log(cancer volume) 
- lweight
- 
log(prostate weight) 
- age
- age 
- lbph
- log(benign prostatic hyperplasia amount) 
- svi
- seminal vesicle invasion 
- lcp
- log(capsular penetration) 
- gleason
- 
Gleason score 
- pgg45
- percentage Gleason scores 4 or 5 
- lpsa
- log(prostate specific antigen) 
Source
Andrews DF and Herzberg AM (1985): Data. New York: Springer-Verlag
Partial Residual Plot
Description
Makes a Partial Residual plot
Usage
prplot(g, i)
Arguments
| g | An object returned from lm() | 
| i | index of predictor | 
Value
none
Author(s)
Julian Faraway
Examples
data(stackloss)
g <- lm(stack.loss ~ .,stackloss)
prplot(g,1)
Panel Study of Income Dynamics subset
Description
The Panel Study of Income Dynamics (PSID), begun in 1968, is a longitudinal study of a representative sample of U.S. individuals. The study is conducted at the Survey Research Center, Institute for Social Research, University of Michigan and is still continuing. The data represents a small subset of the total data.
Format
A data frame with 1661 observations on the following 6 variables.
- age
- age in 1968 
- educ
- years of education 
- sex
- sex of individual, - For- M
- income
- annual income in dollars 
- year
- calendar year 
- person
- ID number for individual 
Source
Martha S. Hill, The Panel Study of Income Dynamics: A User's Guide, Sage Publications, 1992,Newbury Park, CA.
Brightness of paper pulp depending on shift operator
Description
The pulp data frame has 20 rows and 2 columns. Data comes from an
experiment to test the paper brightness depending on a shift operator.
Format
This data frame contains the following columns:
- bright
- Brightness of the pulp as measured by a reflectance meter 
- operator
- Shift operator a-d 
Source
"Statistical techniques applied to production situations" F. Sheldon (1960) Industrial and Engineering Chemistry, 52, 507-509
Leg strength and punting
Description
Investigators studied physical characteristics and ability in 13 (American) football punters. Each volunteer punted a football ten times. The investigators recorded the average distance for the ten punts, in feet.
Format
A data frame with 13 observations on the following 7 variables.
- Distance
- average distance over 10 punts 
- Hang
- hang time 
- RStr
- right leg strength in pounds 
- LStr
- left leg strength in pounds 
- RFlex
- right hamstring muscle flexibility in degrees 
- LFlex
- left hamstring muscle flexibility in degrees 
- OStr
- overall leg strength in foot pounds 
Source
Unknown
Examples
data(punting)
## maybe str(punting) ; plot(punting) ...
Production of PVC by operator and resin railcar
Description
Data from an experiment to study factors affecting the production of the plastic PVC, 3 operators used 8 different devices called resin railcars to produce PVC. For each of the 24 combinations, two samples were produced.
Format
Dataset contains the following variables
- psize
- Particle size 
- operator
- Operator number 1, 2 or 3 
- resin
- Resin railcar 1-8 
Source
R. Morris and E. Watson (1998) "A comparison of the techniques used to evaluate the measurement process" Quality Engineering, 11, 213-219
Activity in pyrimidines
Description
Structural information on 74 2,4-diamino- 5-(substituted benzyl) pyrimidines used as inhibitors of DHFR in E. coli. There are 3 positions where chemical activity occurs and 9 attributes per position leading to 27 total predictors. One predictor had no variability and was removed from the data set. 26 chemical properties of 74 compounds and an activity level
Format
A data frame with 74 observations on the following 27 variables.
- p1.polar
- measured on a [0,1] scale 
- p1.size
- measured on a [0,1] scale 
- p1.flex
- measured on a [0,1] scale 
- p1.h.doner
- measured on a [0,1] scale 
- p1.h.acceptor
- measured on a [0,1] scale 
- p1.pi.doner
- measured on a [0,1] scale 
- p1.pi.acceptor
- measured on a [0,1] scale 
- p1.polarisable
- measured on a [0,1] scale 
- p1.sigma
- measured on a [0,1] scale 
- p2.polar
- measured on a [0,1] scale 
- p2.size
- measured on a [0,1] scale 
- p2.flex
- measured on a [0,1] scale 
- p2.h.doner
- measured on a [0,1] scale 
- p2.h.acceptor
- measured on a [0,1] scale 
- p2.pi.doner
- measured on a [0,1] scale 
- p2.pi.acceptor
- measured on a [0,1] scale 
- p2.polarisable
- measured on a [0,1] scale 
- p2.sigma
- measured on a [0,1] scale 
- p3.polar
- measured on a [0,1] scale 
- p3.size
- measured on a [0,1] scale 
- p3.flex
- measured on a [0,1] scale 
- p3.h.doner
- measured on a [0,1] scale 
- p3.h.acceptor
- measured on a [0,1] scale 
- p3.pi.doner
- measured on a [0,1] scale 
- p3.polarisable
- measured on a [0,1] scale 
- p3.sigma
- measured on a [0,1] scale 
- activity
- log 1/Ki, where Ki is the inhibition constant as experimentally assayed, scaled to [0,1] 
Source
Jonathan D. Hirst, Ross D. King, Michael J. E. Sternberg (1994) Quantitative structure-activity relationships by neural networks and inductive logic programming. I. The inhibition of dihydrofolate reductase by pyrimidines doi:10.1007/BF00125375
Examples
data(pyrimidines)
## maybe str(pyrimidines) ; plot(pyrimidines) ...
Labeled QQ plot
Description
Makes a labeled QQ plot
Usage
qqnorml(
  y,
  main = "Normal Q-Q Plot",
  xlab = "Theoretical Quantiles",
  ylab = "Sample Quantiles",
  ...
)
Arguments
| y | A numeric vector | 
| main | main label | 
| xlab | x-axis label | 
| ylab | y-axis label | 
| ... | arguments passed to plot() | 
Value
none
Author(s)
Julian Faraway
See Also
qqnorm
Examples
qqnorml(rnorm(16))
Rabbit weight gain by diet and litter
Description
A nutritionist studied the effects of six diets, on weight gain of domestic rabbits. From past experience with sizes of litters, it was felt that only 3 uniform rabbits could be selected from each available litter. There were ten litters available forming blocks of size three.
Format
The variables in the dataset were
- treat
- 
Diet a through f 
- gain
- Weight gain 
- block
- Block (10 litters) 
Source
"Experimental Design and Analysis" by M. Lentner and T. Bishop, Valley Book Company, 1986
Rat growth weights affected by additives
Description
The data consist of 5 weekly measurements of body weight for 27 rats. The first 10 rats are on a control treatment while 7 rats have thyroxine added to their drinking water. 10 Rats have thiouracil added to their water.
Format
A data frame with 135 observations on the following 4 variables.
- wt
- Weight of the rat 
- weeks
- Week of the study from 0 to 4 
- subject
- the rat code number 
- treat
- treatment applied to the rat drinking water - a factor with levels - control- thiouracil- thyroxine
Source
Unknown
Effect of toxic agents on rats
Description
An experiment was conducted as part of an investigation to combat the effects of certain toxic agents.
Format
A data frame with 48 observations on the following 3 variables.
- time
- survival time in tens of hours 
- poison
- the poison type - a factor with levels - I- II- III
- treat
- the treatment - a factor with levels - A- B- C- D
Source
Box G and Cox D (1964) "An analysis of transformations" J. Roy. Stat. Soc. Series B. 26 211.
Shape and plate effects on current noise in resistors
Description
The resceram data frame has 12 rows and 3 columns. Shape and plate
effects on current noise in resistors
Format
This data frame contains the following columns:
- noise
- current noise 
- shape
- the geometrical shape of the resistor, A, B, C or D 
- plate
- the ceramic plate on which the resistor was mounted. Only three resistors will fit on one plate. 
Source
Natrella, M (1963) "Experimental Statistics" National Bureau of Standards Handbook 91, Gaithersburg MD.
Salmonella reverse mutagenicity assay
Description
The data was collected in a salmonella reverse mutagenicity assay where the numbers of revertant colonies of TA98 Salmonella observed on each of three replicate plates for different doses of quinoline
Format
A data frame with 18 observations on the following 2 variables.
- colonies
- numbers of revertant colonies of TA98 Salmonella 
- dose
- dose level of quinoline 
Source
Breslow N.E. (1984), Extra-Poisson Variation in Log-linear Models, ApplStat, pp. 38-44.
School expenditure and test scores from USA in 1994-95
Description
The sat data frame has 50 rows and 7 columns.  Data were collected to
study the relationship between expenditures on public education and test
results.
Format
This data frame contains the following columns:
- expend
- Current expenditure per pupil in average daily attendance in public elementary and secondary schools, 1994-95 (in thousands of dollars) 
- ratio
- Average pupil/teacher ratio in public elementary and secondary schools, Fall 1994 
- salary
- 
Estimated average annual salary of teachers in public elementary and secondary schools, 1994-95 (in thousands of dollars) 
- takers
- Percentage of all eligible students taking the SAT, 1994-95 
- verbal
- Average verbal SAT score, 1994-95 
- math
- Average math SAT score, 1994-95 
- total
- 
Average total score on the SAT, 1994-95 
Source
"Getting What You Pay For: The Debate Over Equity in Public School Expenditures" D. Guber, Journal of Statistics Education, 1999
Savings rates in 50 countries
Description
The savings data frame has 50 rows and 5 columns.  The data is
averaged over the period 1960-1970.
Format
This data frame contains the following columns:
- sr
- savings rate - personal saving divided by disposable income 
- pop15
- percent population under age of 15 
- pop75
- percent population over age of 75 
- dpi
- per-capita disposable income in dollars 
- ddpi
- percent growth rate of dpi 
Details
Now also appears as LifeCycleSavings in the datasets package
Source
Belsley, D., Kuh. E. and Welsch, R. (1980) "Regression Diagnostics" Wiley.
See Also
LifeCycleSavings
Car seat position depending driver size
Description
Car drivers like to adjust the seat position for their own comfort. Car designers would find it helpful to know where different drivers will position the seat depending on their size and age. Researchers at the HuMoSim laboratory at the University of Michigan collected data on 38 drivers.
Format
The dataset contains the following variables
- Age
- Age in years 
- Weight
- Weight in lbs 
- HtShoes
- Height in shoes in cm 
- Ht
- Height bare foot in cm 
- Seated
- Seated height in cm 
- Arm
- lower arm length in cm 
- Thigh
- Thigh length in cm 
- Leg
- Lower leg length in cm 
- hipcenter
- horizontal distance of the midpoint of the hips from a fixed location in the car in mm 
Source
"Linear Models in R" by Julian Faraway, CRC Press, 2004
Germination of seeds depending on moisture and covering
Description
A Biologist analyzed an experiment to determine the effect of moisture content on seed germination. Eight boxes of 100 seeds each were treated with the same moisture level. 4 boxes were covered and 4 left uncovered. The process was repeated at 6 different moisture levels (nonlinear scale).
Format
A data frame with 48 observations on the following 3 variables.
- germ
- percentage germinated 
- moisture
- moisture level 
- covered
- a factor with levels - no- yes
Source
Chatfield C. (1995) Problem Solving: A Statistician's Guide, 2ed Chapman Hall.
Examples
data(seeds)
## maybe str(seeds) ; plot(seeds) ...
Semiconductor split-plot experiment
Description
The semicond data frame has 48 rows and 5 columns.
Format
This data frame contains the following columns:
- resistance
- a numeric vector 
- ET
- a factor with levels - 1to- 4representing etch time.
- Wafer
- a factor with levels - 1to- 3
- position
- a factor with levels - 1to- 4
- Grp
- an ordered factor with levels - 1/1<- 1/2<- 1/3<- 2/1<- 2/2<- 2/3<- 3/1<- 3/2<- 3/3<- 4/1<- 4/2<- 4/3
Details
Also found in the SASmixed package
Source
Littel, R. C., Milliken, G. A., Stroup, W. W., and Wolfinger, R. D. (1996), SAS System for Mixed Models, SAS Institute (Data Set 2.2(b)).
Post traumatic stress disorder in abused adult females
Description
The data for this example come from a study of the effects of childhood sexual abuse on adult females. 45 women being treated at a clinic, who reported childhood sexual abuse, were measured for post traumatic stress disorder and childhood physical abuse both on standardized scales. 31 women also being treated at the same clinic, who did not report childhood sexual abuse were also measured. The full study was more complex than reported here and so readers interested in the subject matter should refer to the original article.
Format
The variables in the dataset are
- cpa
- Childhood physical abuse on standard scale 
- ptsd
- Post-traumatic stress disorder on standard scale 
- csa
- Childhood sexual abuse - abused or not abused 
Source
N. Rodriguez and S. Ryan and H. Vande Kemp and D. Foy (1997) "Postraumatic stress disorder in adult female survivors of childhood sexual abuse: A comparison study", Journal of Consulting and Clinical Pyschology, 65, 53-59
Marital sex ratings
Description
Data from a questionaire from 91 couples in the Tucson, Arizona area. Subjects answered the question "Sex is fun for me and my partner". The possible answers were "never or occasionally","fairly often","very often" and "almost always"
Format
A data frame with 16 observations on the following 3 variables.
- y
- the count 
- husband
- a factor with levels - never- fairly- very- always
- wife
- a factor with levels - never- fairly- very- always
Source
Hout, M., Duncan, O. and Sobel M. (1987) Association and heterogeneity: Structural models of similarities and differences. Sociological Methods. 17, 145-184.
Snail production
Description
A study was conducted to optimize snail production for consumption. The percentage water content of the tissues of snails grown under three different levels of relative humidity and two different temperatures was recorded. For each combination, 4 snails were observed.
Format
A data frame with 24 observations on the following 3 variables.
- water
- percentage water content 
- temp
- temperature in C 
- humid
- relative humidity 
Source
Unknown
Examples
data(snail)
## maybe str(snail) ; plot(snail) ...
Solder skips in printing circuit boards
Description
ATT ran an experiment varying five factors relevant to a wave-soldering procedure for mounting components on printed circuit boards. The response variable, skips, is a count of how many solder skips appeared to a visual inspection.
Format
A data frame with 900 observations on the following 6 variables.
- Opening
- a factor with levels - L- M- S
- Solder
- a factor with levels - Thick- Thin
- Mask
- a factor with levels - A1.5- A3- A6- B3- B6
- PadType
- a factor with levels - D4- D6- D7- L4- L6- L7- L8- L9- W4- W9
- Panel
- a numeric vector 
- skips
- count of how many solder skips appeared to a visual inspection 
Source
Comizzoli, R. B., J. M. Landwehr, and J. D. Sinclair (1990). Robust materials and processes: Key to reliability. AT&T Technical Journal 69(6), 113-128.
Examples
data(solder)
## maybe str(solder) ; plot(solder) ...
Block design task testing child ability
Description
Behavioural scientists at Macquarie University conducted an experiment to test the time taken to perform a block design task with 24 fifth grade children (12 boys and 12 girls).
Format
A data frame with 24 observations and 3 variables.
- group
- Solution attempted first by row(r) or corner(c) 
- time
- Time taken to complete the task in seconds 
- eft
- Score on the embedded figures test which is a measure of difficulty in abstracting logical structure of a problem from its context. 
Source
Statistical Modelling in GLIM (1989) M. Aitkin and D. Anderson and B. Francis and J. Hinde Oxford University Press
Sonoluminescence
Description
The sono data frame has 16 rows and 8 columns.  Sonoluminescence is
the process of turning sound energy into light.  An experiment was conducted
to study factors affecting this process.
Format
This data frame contains the following columns:
- Intensity
- Sonoluminescent light intensity 
- Molarity
- Amount of Solute. The coding is "low" for 0.10 mol and "high" for 0.33 mol. 
- Solute
- Solute type. The coding is "low" for sugar and "high" for glycerol. 
- pH
- The coding is "low" for 3 and "high" for 11. 
- Gas
- Gas type in water. The coding is "low" for helium and "high" for air. 
- Water
- Water depth. The coding is "low" for half and "high" for full. 
- Horn
- Horn depth. The coding is "low" for 5 mm and "high" for 10 mm. 
- Flask
- Flask clamping. The coding is "low" for unclamped and "high" for clamped. 
Source
Eva Wilcox and Ken Inn of the NIST Physics Laboratory conducted this experiment during 1999 and published in NIST/SEMATECH e-Handbook of Statistical Methods, http://www.itl.nist.gov/div898/handbook/
Germination failures for soybean seeds
Description
An experiment was conducted to compare the germination rates of the five varieties of soybean. Five plots were available.
Format
A data frame with 25 observations on the following 3 variables.
- variety
- the variety - a factor with levels - arasan- check- fermate- semesan- spergon
- replicate
- the plot - a factor with levels - 1- 2- 3- 4- 5
- failure
- the number of failures out of 100 planted seeds 
Source
Snedecor G. and Cochran W. (1967) Statistical Methods (6th Ed) Iowa State University Press
Teaching methods in Economics
Description
A study to determine the effectiveness of a new teaching method in Economics
Format
A data frame with 32 observations on the following 4 variables.
- grade
- 1 = exam grades improved, 0 = not improved 
- psi
- 1 = student exposed to PSI (a new teach method), 0 = not exposed 
- tuce
- a measure of ability when entering the class 
- gpa
- grade point average 
Source
Spector, L. and Mazzeo, M. (1980), "Probit Analysis and Economic Education", Journal of Economic Education, 11, 37 - 44.
Speedometer cable shrinkage
Description
Speedometer cables can be noisy because of shrinkage in the plastic casing material. An experiment was conducted to find out what caused shrinkage by screening a large number of factors. The engineers started with 15 different factors.
Format
The dataset contains the following variables: (variables a-o are 2 level factors, coded "+" and "-" where "+" indicates a higher value where appropriate)
- a
- liner outer diameter 
- b
- liner die 
- c
- liner material 
- d
- liner line speed 
- e
- wire braid type 
- f
- braiding tension 
- g
- wire diameter 
- h
- liner tension 
- i
- liner temperature 
- j
- coating material 
- k
- coating die type 
- l
- melt temperature 
- m
- screen pack 
- n
- cooling method 
- o
- line speed 
- y
- percentage shrinkage per specimen 
Source
G. P. Box and S. Bisgaard and C. Fung (1988) "An explanation and critque of Taguchi's contributions to quality engineering", Quality and reliability engineering international, 4, 123-131
Star temperatures and light intensites
Description
Data on the log of the surface temperature and the log of the light intensity of 47 stars in the star cluster CYG OB1, which is in the direction of Cygnus,
Format
A data frame with 47 observations on the following 3 variables.
- index
- a numeric vector 
- temp
- temperature 
- light
- light intensity 
Source
Rousseeuw, P. and A. Leroy (1987). Robust Regression and Outlier Detection. New York: Wiley.
Examples
data(star)
## maybe str(star) ; plot(star) ...
Marks in a statistics class
Description
Marks from Statistics 500 one year at the University of Michigan
Format
A data frame with 55 observations on the following 4 variables.
- midterm
- a numeric vector 
- final
- a numeric vector 
- hw
- a numeric vector 
- total
- a numeric vector 
Source
Julian Faraway
Examples
data(stat500)
## maybe str(stat500) ; plot(stat500) ...
Stepping and effect on heart rate
Description
An experiment was conducted to explore the nature of the relationship between a person's heart rate and the frequency at which that person stepped up and down on steps of various heights.
Format
A data frame with 30 observations on the following 6 variables.
- Order
- running order within the experiment 
- Block
- Experimenter used 
- Height
- 0 if step at the low (5.75in) height, 1 if at the high (11.5in) height 
- Frequency
- the rate of stepping. 0 if slow (14 steps/min), 1 if medium (21 steps/min), 2 if high (28 steps/min) 
- RestHR
- the resting heart rate of the subject before a trial, in beats per minute 
- HR
- the final heart rate of the subject after a trial, in beats per minute 
Source
Unknown
Examples
data(stepping)
## maybe str(stepping) ; plot(stepping) ...
Strong interaction experiment data
Description
Example Dataset from "Practical Regression and Anova"
Format
Dataframe with 10 cases
- momentum
- energy
- inverse total energy 
- crossx
- Scattering cross-section/sec 
- sd
- standard deviation 
Source
Weisberg, H., Beier, H., Brody, H., Patton, R., Raychaudhari, K., Takeda, H., Thern, R. and Van Berg, R. (1978). s-dependence of proton fragmentation by hadrons. II. Incident laboratory momenta, 30–250 GeV/c. Physics Review D, 17, 2875–2887.
References
Weisberg, S. (2014). Applied Linear Regression, 4th edition. Hoboken NJ: Wiley.
Suicide method data from the UK
Description
One year of suicide data from the United Kingdom crossclassified by sex, age and method.
Format
A data frame with 36 observations on the following 4 variables.
- y
- number of people 
- cause
- method used - a factor with levels - drug(suicide by solid or liquid matter),- gas,- gun(guns, knives or explosives)- hang(hanging, strangling, suffocating or drowning,- jump- other
- age
- a factor with levels - m(middle-aged)- o(old)- y(young)
- sex
- a factor with levels - f- m
Source
Everitt B. & Dunn G. (1991) "Applied Multivariate Data Analysis" Edward Arnold
Abbreviated Regression Summary
Description
Generic summaries for lm, glm and mer objects
Usage
sumary(object, ...)
Arguments
| object | An lm, glm or mer object returned from lm(), glm() or lmer() respectively | 
| ... | further arguments passed to or from other methods. | 
Details
This generic function provides an abbreviated regression output containing
the more useful information. Users wanting to see more are advised to use
summary()
Value
returns the same as summary()
Author(s)
Julian Faraway
References
This function is adapted from the display() function in
the arm package
See Also
Examples
data(stackloss)
object <- lm(stack.loss ~ .,stackloss)
sumary(object)
Study of teenage gambling in Britain
Description
The teengamb data frame has 47 rows and 5 columns. A survey was
conducted to study teenage gambling in Britain.
Format
This data frame contains the following columns:
- sex
- 0=male, 1=female 
- status
- Socioeconomic status score based on parents' occupation 
- income
- in pounds per week 
- verbal
- verbal score in words out of 12 correctly defined 
- gamble
- expenditure on gambling in pounds per year 
Source
Ide-Smith & Lea, 1988, Journal of Gambling Behavior, 4, 110-118
Toenail infection treatment study
Description
The data come from a Multicenter study comparing two oral treatments for toenail infection. Patients were evaluated for the degree of separation of the nail. Patients were randomized into two treatments and were followed over seven visits - four in the first year and yearly thereafter. The patients have not been treated prior to the first visit so this should be regarded as the baseline.
Format
A data frame with 1908 observations on the following 5 variables.
- ID
- ID of patient 
- outcome
- 0=none or mild seperation, 1=moderate or severe 
- treatment
- the treatment A=0 or B=1 
- month
- time of the visit (not exactly monthly intervals hence not round numbers) 
- visit
- the number of the visit 
Source
De Backer, M., De Vroey, C., Lesaffre, E., Scheys, I., and De Keyser, P. (1998). Twelve weeks of continuous oral therapy for toenail onychomycosis caused by dermatophytes: A double-blind comparative trial of terbinafine 250 mg/day versus itraconazole 200 mg/day. Journal of the American Academy of Dermatology, 38, 57-63.
References
Lesaffre, E. and Spiessens, B. (2001). On the effect of the number of quadrature points in a logistic random-effects model: An example. Journal of the Royal Statistical Society, Series C, 50, 325-335.
G. Fitzmaurice, N. Laird and J. Ware (2004) Applied Longitudinal Analysis, Wiley
Survival of trout eggs depending on time and location
Description
Boxes of trout eggs were buried at five different stream locations and retrieved at 4 different times. The number of surviving eggs was recorded. The box was not returned to the stream.
Format
A data frame with 20 observations on the following 4 variables.
- survive
- the number of surviving eggs 
- total
- the number of eggs in the box 
- location
- the location in the stream with levels - 1- 2- 3- 4- 5
- period
- the number of weeks after placement that the box was withdrawn levels - 4- 7- 8- 11
Source
Manly B. (1978) Regression models for proportions with extraneous variance. Biometrie-Praximetrie, 18, 1-18.
References
Hinde J. and Demetrio C. (1988) Overdispersion: Models and estimation. Computational Statistics and Data Analysis. 27, 151-170.
Truck leaf spring experiment
Description
Data on an experiment concerning the production of leaf springs for trucks.
A 2^{5-1} fractional factorial experiment with 3 replicates was
carried out with objective of recommending production settings to achieve a
free height as close as possible to 8 inches.
Format
A data frame with 48 observations on the following 6 variables.
- B
- furnace temperature - a factor with levels - +- -
- C
- heating time - a factor with levels - +- -
- D
- transfer time - a factor with levels - +- -
- E
- hold-down time - a factor with levels - +- -
- O
- quench oil temperature - a factor with levels - +- -
- height
- leaf spring free height in inches 
Source
J. J. Pignatiello and J. S. Ramberg (1985) Contribution to discussion of offline quality control, parameter design and the Taguchi method, Journal of Quality Technology, 17 198-206.
References
P. McCullagh and J. Nelder (1989) "Generalized Linear Models" Chapman and Hall, 2nd ed.
Incubation temperature and the sex of turtles
Description
Incubation temperature can affect the sex of turtles. There are 3 independent replicates for each temperature.
Format
A data frame with 15 observations on the following 3 variables.
- temp
- temperature in degrees centigrade 
- male
- number of male turtles hatched 
- female
- number of female turtles hatched 
Source
Beyond Traditional Statistical Methods Copyright 2000 D. Cook, P. Dixon, W. M. Duckworth, M. S. Kaiser, K. Koehler, W. Q. Meeker and W. R. Stephenson. Developed as part of NSF/ILI grant DUE9751644.
Examples
data(turtle)
Life, TVs and Doctors
Description
Life expectancy, doctors and televisions collected on 38 countries in 1993
Format
A data frame with 38 observations on the following 3 variables.
- life
- Life expectancy in years 
- tv
- Number of people per television set 
- doctor
- Number of people per doctor 
Source
Unknown, data for illustration purposes only
Examples
data(tvdoctor)
## maybe str(tvdoctor) ; plot(tvdoctor) ...
Twin IQs from Burt
Description
Study of IQ in twins reared apart
Format
A dataframe with the following variables:
- Foster
- IQ of the fostered child 
- Biological
- IQ of the biological child 
- Social
- social class of natural parents 
Source
Burt, C. (1966). The genetic estimation of differences in intelligence: A study of monozygotic twins reared together and apart. Br. J. Psych., 57, 147-153.
References
Weisberg, S. (2014). Applied Linear Regression, 4th edition. Hoboken NJ: Wiley.
UNC student opinions about the Vietnam War
Description
A student newspaper conducted a survey of student opinions about the Vietnam War in May 1967. Responses were classified by sex, year in the program and one of four opinions. The survey was voluntary.
Format
A data frame with 40 observations on the following 4 variables.
- y
- the count 
- policy
- a factor with levels - A(defeat power of North Vietnam by widespread bombing and land invasion)- B(follow the present policy)- C(withdraw troops to strong points and open negotiations on elections involving the Viet Cong)- D(immediate withdrawal of all U.S. troops)
- sex
- a factor with levels - Female- Male
- year
- a factor with levels - Fresh- Grad- Junior- Senior- Soph
Source
M. Aitkin and D. Anderson and B. Francis and J. Hinde (1989) "Statistical Modelling in GLIM" Oxford University Press.
Weekly wages of US male workers in 1988
Description
The uswages data frame has 2000 rows and 10 columns. Weekly Wages for
US male workers sampled from the Current Population Survey in 1988.
Format
This data frame contains the following columns:
- wage
- Real weekly wages in dollars (deflated by personal consumption expenditures - 1992 base year) 
- educ
- Years of education 
- exper
- Years of experience 
- race
- 1 if Black, 0 if White (other races not in sample) 
- smsa
- 1 if living in Standard Metropolitan Statistical Area, 0 if not 
- ne
- 1 if living in the North East 
- mw
- 1 if living in the Midwest 
- we
- 1 if living in the West 
- so
- 1 if living in the South 
- pt
- 1 if working part time, 0 if not 
Source
Bierens, H.J., and D. Ginther (2001): "Integrated Conditional Moment Testing of Quantile Regression Models", Empirical Economics 26, 307-324
vif
Description
vif
Usage
vif(object)
## Default S3 method:
vif(object)
## S3 method for class 'lm'
vif(object)
Arguments
| object | a data matrix (design matrix without intercept) or a model object | 
Details
Computes the variance inflation factors
Value
variance inflation factors
Author(s)
Julian Faraway
Examples
data(stackloss)
vif(stackloss[,-4])
#  Air.Flow Water.Temp Acid.Conc.
#    2.9065     2.5726     1.3336
Acuity of vision in response to light flash
Description
The acuity of vision for seven subjects was tested. The response is the lag in milliseconds between a light flash and a response in the cortex of the eye. Each eye is tested at four different powers of lens. An object at the distance of the second number appears to be at distance of the first number.
Format
A data frame with 56 observations on the following 4 variables.
- acuity
- a numeric vector 
- power
- a factor with levels - 6/6- 6/18- 6/36- 6/60
- eye
- a factor with levels - left- right
- subject
- a factor with levels - 1- 2- 3- 4- 5- 6- 7
Source
Crowder, M. J. and D. J. Hand (1990). Analysis of Repeated Measures. London: Chapman & Hall.
Examples
data(vision)
## maybe str(vision) ; plot(vision) ...
resitivity of wafer in semiconductor experiment
Description
A full factorial experiment with four two-level predictors.
Format
A data frame with 16 observations on the following 5 variables.
- x1
- a factor with levels - -- +
- x2
- a factor with levels - -- +
- x3
- a factor with levels - -- +
- x4
- a factor with levels - -- +
- resist
- Resistivity of the wafer 
Source
Myers, R. and Montgomery D. (1997) A tutorial on generalized linear models, Journal of Quality Technology, 29, 274-291.
Defects in a wave soldering process
Description
Components are attached to an electronic circuit card assembly by a
wave-soldering process. The soldering process involves baking and preheating
the circuit card and then passing it through a solder wave by conveyor.
Defect arise during the process. Design is 2^{7-3} with 3 replicates.
Format
A data frame with 16 observations on the following 10 variables.
- y1
- Number of defects in the first replicate 
- y2
- Number of defects in the second replicate 
- y3
- Number of defects in the third replicate 
- prebake
- prebake condition - a factor with levels - 1- 2
- flux
- flux density - a factor with levels - 1- 2
- speed
- conveyor speed - a factor with levels - 1- 2
- preheat
- preheat condition - a factor with levels - 1- 2
- cooling
- cooling time - a factor with levels - 1- 2
- agitator
- ultrasonic solder agitator - a factor with levels - 1- 2
- temp
- solder temperature - facctor with levels - 1- 2
Source
L. Condra (1993) Reliability improvement with design of experiments. Marcel Dekker, NY.
References
M. Hamada and J. Nelder (1997) Generalized linear models for quality improvement experiments, Journal of Quality Technology, 29, 292-304
Wisconsin breast cancer database
Description
Data come from a study of breast cancer in Wisconsin. There are 681 cases of potentially cancerous tumors of which 238 are actually malignant. Determining whether a tumor is really malignant is traditionally determined by an invasive surgical procedure. The purpose of this study was to determine whether a new procedure called fine needle aspiration which draws only a small sample of tissue could be effective in determining tumor status.
Format
A data frame with 681 observations on the following 10 variables.
- Class
- 0 if malignant, 1 if benign 
- Adhes
- marginal adhesion 
- BNucl
- bare nuclei 
- Chrom
- bland chromatin 
- Epith
- epithelial cell size 
- Mitos
- mitoses 
- NNucl
- normal nucleoli 
- Thick
- clump thickness 
- UShap
- cell shape uniformity 
- USize
- cell size uniformity 
Details
The predictor values are determined by a doctor observing the cells and rating them on a scale from 1 (normal) to 10 (most abnormal) with respect to the particular characteristic.
Source
Bennett, K.,P., and Mangasarian, O.L., Neural network training via linear programming. In P. M. Pardalos, editor, Advances in Optimization and Parallel Computing, pages 56-57. Elsevier Science, 1992
Western Collaborative Group Study
Description
3154 healthy young men aged 39-59 from the San Francisco area were assessed for their personality type. All were free from coronary heart disease at the start of the research. Eight and a half years later change in this situation was recorded.
Format
A data frame with 3154 observations on the following 13 variables.
- age
- age in years 
- height
- height in inches 
- weight
- weight in pounds 
- sdp
- systolic blood pressure in mm Hg 
- dbp
- diastolic blood pressure in mm Hg 
- chol
- Fasting serum cholesterol in mm % 
- behave
- behavior type which is a factor with levels - A1- A2- B3- B4
- cigs
- number of cigarettes smoked per day 
- dibep
- behavior type a factor with levels - A(Agressive)- B(Passive)
- chd
- coronary heat disease developed is a factor with levels - no- yes
- typechd
- type of coronary heart disease is a factor with levels - angina- infdeath- none- silent
- timechd
- Time of CHD event or end of follow-up 
- arcus
- arcus senilis is a factor with levels - absent- present
Details
The WCGS began in 1960 with 3,524 male volunteers who were employed by 11 California companies. Subjects were 39 to 59 years old and free of heart disease as determined by electrocardiogram. After the initial screening, the study population dropped to 3,154 and the number of companies to 10 because of various exclusions. The cohort comprised both blue- and white-collar employees. At baseline the following information was collected: socio-demographic including age, education, marital status, income, occupation; physical and physiological including height, weight, blood pressure, electrocardiogram, and corneal arcus; biochemical including cholesterol and lipoprotein fractions; medical and family history and use of medications; behavioral data including Type A interview, smoking, exercise, and alcohol use. Later surveys added data on anthropometry, triglycerides, Jenkins Activity Survey, and caffeine use. Average follow-up continued for 8.5 years with repeat examinations
Source
Statistics for Epidemiology by N. Jewell (2004)
References
Coronary Heart Disease in the Western Collaborative Group Study Final Follow-up Experience of 8 1/2 Years Ray H. Rosenman, MD; Richard J. Brand, PhD; C. David Jenkins, PhD; Meyer Friedman, MD; Reuben Straus, MD; Moses Wurm, MD JAMA. 1975;233(8):872-877. doi:10.1001/jama.1975.03260080034016.
Examples
data(wcgs)
## maybe str(wcgs) ; plot(wcgs) ...
welding strength DOE
Description
An experiment to investigate factors affecting welding strength.
Format
A data frame with 16 observations on the following 10 variables.
- Rods
- a 0-1 predictor 
- Drying
- a 0-1 predictor 
- Material
- a 0-1 predictor 
- Thickness
- a 0-1 predictor 
- Angle
- a 0-1 predictor 
- Opening
- a 0-1 predictor 
- Current
- a 0-1 predictor 
- Method
- a 0-1 predictor 
- Preheating
- a 0-1 predictor 
- Strength
- The welding strength 
Source
G. Box and R. Meyer (1986) Dispersion effects from fractional designs, Technometrics, 28, 19-27
Percentage of Body Fat and Body Measurements in Women
Description
Age, weight, height, and 10 body circumference measurements are recorded for 184 women. Each woman's percentage of body fat was accurately estimated by an underwater weighing technique.
Format
A data frame with 184 observations on the following 19 variables.
- siri
- Percent body fat using Siri's equation 
- weight
- Weight (lbs) 
- height
- Height (inches) 
- bmi
- Body Mass Index 
- age
- Age (yrs) 
- neck
- Neck circumference (cm) 
- chest
- Chest circumference (cm) 
- calf
- Calf circumference (cm) 
- biceps
- Extended biceps circumference (cm) 
- hip
- Hip circumference (cm) 
- abdom
- Horizontal minimal measurement, at the end of a normal expiration (cm) 
- forearm
- Forearm circumference (cm) 
- thigh
- (Proximal Thigh) Horizontal measurement immediately distal to the gluteal furrow (cm) 
- mthigh
- (Middle Thigh) Measurement midway between the midpoint of the inguinal crease and the proximal border of the patella (cm) 
- dthigh
- (Distal Thigh) Measurement proximal to the femoral epicondyles (cm) 
- wrist
- Wrist circumference (cm) distal to the styloid processes 
- knee
- Knee circumference (cm) 
- elbow
- A minimal circumference measurement with the elbow extended (cm) 
- ankle
- Ankle circumference (cm) 
Source
Roger W. Johnson (2021): Fitting Percentage of Body Fat to Simple Body Measurements: College Women, Journal of Statistics and Data Science Education, DOI: 10.1080/26939169.2021.1971585 (Note that I have changed some of the variable names to correspond with the older fat data for men)
See Also
Insect damage to wheat by variety
Description
Insect damage to wheat by variety
Format
A data frame with 13 observations on the following 2 variables.
- damage
- a numeric vector 
- variety
- a factor with levels - A- B- C- D
Source
Unknown
Examples
data(wheat)
## maybe str(wheat) ; plot(wheat) ...
Data on players from the 2010 World Cup
Description
Data on players from the 2010 World Cup
Format
A data frame with 595 observations on the following 7 variables.
- Team
- Country 
- Position
- a factor with levels - Defender- Forward- Goalkeeper- Midfielder
- Time
- Time played in minutes 
- Shots
- Number of shots attempted 
- Passes
- Number of passes made 
- Tackles
- Number of tackles made 
- Saves
- Number of saves made 
Details
None
Source
Lost
Examples
data(worldcup)
## maybe str(worldcup) ; plot(worldcup) ...