livertests: Dataset: livertests

In reflimR: Reference Limit Estimation Using Routine Laboratory Data

Dataset: livertests

Description

Example data showing eight different biomarkers (laboratory tests), which are frequently measured in healthy controls and patients with liver diseases.

Usage

livertests

Format

A data frame with 612 rows and 11 columns:

Category

healthy reference individual or patient

Age

age in years

Sex

sex f = female or m = male

ALB

albumin, g/L

ALT

alanine aminotransferase, U/L

AST

aspartate aminotransferase, U/L

BIL

bilirubin, µmol/l

CHE

choline esterase, kU/L

CREA

creatinine, µmol/L

GGT

gamma-glutamyl transferase, U/L

PROT

total protein, mg/L

Source

<https://archive.ics.uci.edu/ml/datasets/HCV+data>

Examples

summary(livertests)
pie(table(livertests$Category), labels = c("patients", "controls"))
plot(livertests$Age, livertests$ALB, xlab = "Age [yr]", ylab = "ALB [g/L]")
grid()
abline(lm(livertests$ALB ~ livertests$Age))

che <- livertests$CHE
ref <- livertests$CHE[livertests$Category == "reference"]
pat <- livertests$CHE[livertests$Category == "patient"]

hist(che, breaks = 1 : 20, col = "white", main = "cholinesterase", xlab = "kU/L")
hist(ref, breaks = 1 : 20, col = rgb(0, 1, 0, 0.5), add = TRUE)
hist(pat, breaks = 1 : 20, col = rgb(1, 0, 0, 0.5), add = TRUE)
legend("topright", fill = c(rgb(1,1,1,1), rgb(0,1,0,0.5), rgb(1,0,0,0.5)),
    legend = c("all", "controls", "patients"))

t.test(ref, pat)
var.test(ref, pat)

che.f <- livertests$CHE[livertests$Sex == "f"]
che.m <- livertests$CHE[livertests$Sex == "m"]
plot(density(che.f), xlim = c(0, 20), col = "red",
   main = "cholinesterase", xlab = "kU/L")
lines(density(che.m), col = "blue")
legend("topright", lty = 1, col = c("red", "blue"), legend = c("females", "males"))

reflim(che.m, main = "CHE (m)", xlab = "kU/L")
reflim(livertests$AST[livertests$Sex == "m"], main = "AST (m)", xlab = "U/L")

reflimR documentation built on May 29, 2024, 12:18 p.m.