bnlearn (version 3.1)

coronary: Coronary Heart Disease data set

Description

Probable risk factors for coronary trombosis, comprising data from 1841 men.

Usage

data(coronary)

format

The coronary data set contains the following 6 variables:

• Smoking(smoking): a two-level factor with levelsnoandyes.
• M. Work(strenuous mental work): a two-level factor with levelsnoandyes.
• P. Work(strenuous physical work): a two-level factor with levelsnoandyes.
• Pressure(systolic blood pressure): a two-level factor with levels<140< code="">and>140.
• Proteins(ratio of beta and alpha lipoproteins): a two-level factor with levels<3< code="">and>3.
• Family(family anamnesis of coronary heart disease): a two-level factor with levelsnegandpos.

source

Edwards DI (2000). Introduction to Graphical Modelling. Springer, 2nd edition.

Reinis Z, Pokorny J, Basika V, Tiserova J, Gorican K, Horakova D, Stuchlikova E, Havranek T, Hrabovsky F (1981). "Prognostic Significance of the Risk Profile in the Prevention of Coronary Heart Disease". Bratisl Lek Listy, 76, 137-150. Published on Bratislava Medical Journal, in Czech.

Whittaker J (1990). Graphical Models in Applied Multivariate Statistics. Wiley.

Examples

# This is the undirected graphical model from Whittaker (1990).
data(coronary)
ug = empty.graph(names(coronary))
arcs(ug, ignore.cycles = TRUE) = matrix(
  c("Family", "M. Work", "M. Work", "Family",
    "M. Work", "P. Work", "P. Work", "M. Work",
    "M. Work", "Proteins", "Proteins", "M. Work",
    "M. Work", "Smoking", "Smoking", "M. Work",
    "P. Work", "Smoking", "Smoking", "P. Work",
    "P. Work", "Proteins", "Proteins", "P. Work",
    "Smoking", "Proteins", "Proteins", "Smoking",
    "Smoking", "Pressure", "Pressure", "Smoking",
    "Pressure", "Proteins", "Proteins", "Pressure"),
  ncol = 2, byrow = TRUE,
  dimnames = list(c(), c("from", "to")))