Self-reported and measured BMI
Dataset containing height and weight data (measured, self-reported) from two studies.
A data frame with 2060 rows and 15 variables:
Study, either krul or mgg (factor)
Person identification number
Population, all NL (factor)
Age of respondent in years
Sex of respondent (factor)
Height measured (cm)
Weight measured (kg)
Height reported (cm)
Weight reported (kg)
Pregnancy (factor), all Not pregnant
Educational level (factor)
Ethnicity (factor)
Obtained through web survey (factor)
BMI measured (kg/m2)
BMI reported (kg/m2)
This dataset combines two datasets: krul data (Krul, 2010) (1257
persons) and the mgg data (Van Keulen 2011; Van der Klauw 2011) (803
persons). The krul dataset contains height and weight (both measures
and self-reported) from 1257 Dutch adults, whereas the mgg dataset
contains self-reported height and weight for 803 Dutch adults. Section 7.3 in
Van Buuren (2012) shows how the missing measured data can be imputed in the
mgg data, so corrected prevalence estimates can be calculated.
Krul, A., Daanen, H. A. M., Choi, H. (2010). Self-reported and measured weight, height and body mass index (BMI) in Italy, The Netherlands and North America. European Journal of Public Health, 21(4), 414-419.
Van Keulen, H.M.,, Chorus, A.M.J., Verheijden, M.W. (2011). Monitor Convenant Gezond Gewicht Nulmeting (determinanten van) beweeg- en eetgedrag van kinderen (4-11 jaar), jongeren (12-17 jaar) en volwassenen (18+ jaar). TNO/LS 2011.016. Leiden: TNO.
Van der Klauw, M., Van Keulen, H.M., Verheijden, M.W. (2011). Monitor Convenant Gezond Gewicht Beweeg- en eetgedrag van kinderen (4-11 jaar), jongeren (12-17 jaar) en volwassenen (18+ jaar) in 2010 en 2011. TNO/LS 2011.055. Leiden: TNO. (in Dutch)
Van Buuren, S. (2018). Flexible Imputation of Missing Data. Second Edition. Chapman & Hall/CRC. Boca Raton, FL.
md.pattern(selfreport[, c("age", "sex", "hm", "hr", "wm", "wr")])
### FIMD Section 7.3.5 Application
bmi <- function(h, w) {
return(w / (h / 100)^2)
}
init <- mice(selfreport, maxit = 0)
meth <- init$meth
meth["bm"] <- "~bmi(hm,wm)"
pred <- init$pred
pred[, c("src", "id", "web", "bm", "br")] <- 0
imp <- mice(selfreport, pred = pred, meth = meth, seed = 66573, maxit = 2, m = 1)
## imp <- mice(selfreport, pred=pred, meth=meth, seed=66573, maxit=20, m=10)
### Like FIMD Figure 7.6
cd <- complete(imp, 1)
xy <- xy.coords(cd$bm, cd$br - cd$bm)
plot(xy,
col = mdc(2), xlab = "Measured BMI", ylab = "Reported - Measured BMI",
xlim = c(17, 45), ylim = c(-5, 5), type = "n", lwd = 0.7
)
polygon(x = c(30, 20, 30), y = c(0, 10, 10), col = "grey95", border = NA)
polygon(x = c(30, 40, 30), y = c(0, -10, -10), col = "grey95", border = NA)
abline(0, 0, lty = 2, lwd = 0.7)
idx <- cd$src == "krul"
xyc <- xy
xyc$x <- xy$x[idx]
xyc$y <- xy$y[idx]
xys <- xy
xys$x <- xy$x[!idx]
xys$y <- xy$y[!idx]
points(xyc, col = mdc(1), cex = 0.7)
points(xys, col = mdc(2), cex = 0.7)
lines(lowess(xyc), col = mdc(4), lwd = 2)
lines(lowess(xys), col = mdc(5), lwd = 2)
text(1:4, x = c(40, 28, 20, 32), y = c(4, 4, -4, -4), cex = 3)
box(lwd = 1)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.