R: Converts measurements to standard deviation scores (SDS)
y2z
R Documentation
Converts measurements to standard deviation scores (SDS)
Description
Converts measurements into age- and sex-conditional standard deviation score
(SDS) using an external reference.
Usage
y2z(y = c(75, 80, 85), x = 1, sex = "M", sub = "N",
ref = get("nl4.hgt"), dist = "LMS", dec = 3, sex.fallback = "M",
sub.fallback = "N", tail.adjust = FALSE)
Arguments
y
A numerical vector containing the outcome measurements. The length
length(y) determines the size of the output vector.
x
A vector containing the values of the numerical covariate (typically
decimal age or height) at which conversion is desired. Values are replicated
to match length(y).
sex
A character vector indicating whether the male ("M") of
female ("F")reference should be used. Values are replicated to match
length(y).
sub
A character vector indicating the level of the sub field of
the reference standard defined in ref
ref
A data frame containing a factor sex, a numerical variable
age containing the tabulated decimal point ages, and two or more
numerical variables with reference values. See details.
dist
A string identifying the type of distribution. Values values are:
"NO", "BCCG", "LMS", "BCPE" and "BCT".
The default is "LMS".
dec
A scalar value indicating the number of decimals used to round the
value.
sex.fallback
The level of the sex field used when no match is
found. The default is "M" for males. Specify sex.fallback="NA"
if unmatched entries should receive a NA value.
sub.fallback
The level of the sub field used when no match is
found. The default is "N" for normal. Specify
sub.fallback="NA" if unmatched entries should receive a NA
value.
tail.adjust
Logical. If TRUE then the WHO method for
tail adjustment is applied. The default is FALSE.
Details
Functions z2y() and y2z() are the inverse of each other.
The argument dist determines the statistical distribution. The
possibilities are as follows:
list(""NO"")
ref
should contain columns mean and sd, containing the mean and the
standard deviation in the external reference population.
list(""LMS"")
ref should contain columns L, S
and M containing the LMS parameters.
list(""BCCG"")
ref should contain columns mu,
sigma and nu containing the Box-Cox Cole-Green parameters.
list(""BCPE"")
ref should contain columns mu,
sigma, nu and tau containing the Box-Cox Power
Exponential parameters.
list(""BCT"")
ref should contain
columns mu, sigma, nu and tau containing the
Box-Cox T distribution parameters.
Value
For y2z(): A vector with length(y) elements containing
the standard deviation score. For z2y(): A vector with
length(z) elements containing quantiles.
Author(s)
Stef van Buuren, 2010
See Also
z2y
Examples
boys <- boys7482
# SDS of height 115 cm at age 5 years,
# relative to Dutch boys reference
y2z(y=115, x=5)
# same relative to Dutch girls
y2z(y=115, x=5, sex="F")
# SDS of IOTF BMI cut-off value for overweight (boys 2-18)
# relative to Dutch boys reference
cutoff <- c(
18.41, 18.15, 17.89, 17.72, 17.55, 17.49, 17.42, 17.49, 17.55, 17.74,
17.92, 18.18, 18.44, 18.77, 19.10, 19.47, 19.84, 20.20, 20.55, 20.89,
21.22, 21.57, 21.91, 22.27, 22.62, 22.96, 23.29, 23.60, 23.90, 24.18,
24.46, 24.73, 25.00)
age <- seq(2, 18, by=0.5)
(z <- y2z(y=cutoff, x=age, sex="M", ref=nl4.bmi))
# apply inverse transformation to check calculations
round(z2y(z, age, ref=nl4.bmi), 2)
cutoff
# calculate percentiles of weight 12 kg at 2 years (boys, girls)
100*round(pnorm(y2z(y=c(12,12), x=2, sex=c("M","F"), ref=nl4.wgt)),2)
# # percentage of children lighter than 15kg at ages 2-5
e <- expand.grid(age=2:5, sex=c("M","F"))
z <- y2z(y=rep(15,nrow(e)), x=e$age, sex=e$sex, ref=nl4.wgt)
w <- matrix(100*round(pnorm(z),2), nrow=2, byrow=TRUE)
dimnames(w) <- list(c("boys","girls"),2:5)
w
# analysis in Z scale
hgt.z <- y2z(y=boys$hgt, x=boys$age, sex="M", ref=nl4.hgt)
wgt.z <- y2z(y=boys$wgt, x=boys$age, sex="M", ref=nl4.wgt)
plot(hgt.z, wgt.z, col="blue")
# z2y
# quantile at SD=0 of age 2 years,
# height Dutch boys
z2y(z=0, x=2)
# same for Dutch girls
z2y(z=0, x=2, sex="F")
# quantile at SD=c(-1,0,1) of age 2 years, BMI Dutch boys
z2y(z=c(-1,0,+1), x=2, ref=nl4.bmi)
# 0SD line (P50) in kg of weight for age in 5-10 year, Dutch boys
z2y(z=rep(0,6), x=5:10, ref=nl4.wgt)
# 95th percentile (P95), age 10 years, wfa, Dutch boys
z2y(z=qnorm(0.95), x=10, ref=nl4.wgt)
# table of P3, P10, P50, P90, P97 of weight for 5-10 year old dutch boys
# age per year
age <- 5:10
p <- c(0.03,0.1,0.5,0.9,0.97)
z <- rep(qnorm(p), length(age))
x <- rep(age, each=length(p))
w <- matrix(z2y(z, x=x, sex="M", ref=nl4.wgt), ncol=length(p),
byrow=TRUE)
dimnames(w) <- list(age, p)
round(w,1)
# standard set of Z-scores of weight for all tabulated ages, boys & girls
# and three etnicities
sds <- c(-2.5, -2, -1, 0, 1, 2, 2.5)
age <- nl4.wgt$x
z <- rep(sds, times=length(age))
x <- rep(age, each=length(sds))
sex <- rep(c("M","F"), each=length(z)/2)
w <- z2y(z=z, x=x, sex=sex, ref=nl4.wgt)
w <- matrix(w, ncol=length(sds), byrow=TRUE)
dimnames(w) <- list(age, sds)
data.frame(sub=nl4.wgt$sub,sex=nl4.wgt$sex,round(w,2), row.names=NULL)
# P85 of BMI in 5-8 year old Dutch boys and girls
e <- expand.grid(age=5:8, sex=c("M","F"))
w <- z2y(z=rep(qnorm(0.85),nrow(e)), x=e$age, sex=e$sex, ref=nl4.bmi)
w <- matrix(w, nrow=2, byrow=TRUE)
dimnames(w) <- list(c("boys","girls"),5:8)
w
# data transformation of height z-scores to cm-scale
z <- c(-1.83, 0.09, 2.33, 0.81, -1.20)
x <- c(8.33, 0.23, 19.2, 24.3, 10)
sex <- c("M", "M", "F", "M", "F")
round(z2y(z=z, x=x, sex=sex, ref=nl4.hgt), 1)
# interpolate published height standard
# to daily values, days 0-31, boys
# on centiles -2SD, 0SD and +2SD
days <- 0:31
sds <- c(-2, 0, +2)
z <- rep(sds, length(days))
x <- rep(round(days/365.25,4), each=length(sds))
w <- z2y(z, x, sex="M", ref=nl4.hgt)
w <- matrix(w, ncol=length(sds), byrow=TRUE)
dimnames(w) <- list(days, sds)
w
Results
R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> library(AGD)
AGD 0.35 2015-05-27
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/AGD/y2z.Rd_%03d_medium.png", width=480, height=480)
> ### Name: y2z
> ### Title: Converts measurements to standard deviation scores (SDS)
> ### Aliases: y2z
> ### Keywords: distribution
>
> ### ** Examples
>
> boys <- boys7482
>
> # SDS of height 115 cm at age 5 years,
> # relative to Dutch boys reference
> y2z(y=115, x=5)
[1] 0.424
>
> # same relative to Dutch girls
> y2z(y=115, x=5, sex="F")
[1] 0.706
>
> # SDS of IOTF BMI cut-off value for overweight (boys 2-18)
> # relative to Dutch boys reference
> cutoff <- c(
+ 18.41, 18.15, 17.89, 17.72, 17.55, 17.49, 17.42, 17.49, 17.55, 17.74,
+ 17.92, 18.18, 18.44, 18.77, 19.10, 19.47, 19.84, 20.20, 20.55, 20.89,
+ 21.22, 21.57, 21.91, 22.27, 22.62, 22.96, 23.29, 23.60, 23.90, 24.18,
+ 24.46, 24.73, 25.00)
> age <- seq(2, 18, by=0.5)
> (z <- y2z(y=cutoff, x=age, sex="M", ref=nl4.bmi))
[1] 1.448 1.459 1.411 1.371 1.299 1.255 1.189 1.190 1.173 1.209 1.223 1.255
[13] 1.273 1.308 1.342 1.381 1.416 1.446 1.459 1.470 1.470 1.471 1.464 1.462
[25] 1.452 1.442 1.426 1.412 1.395 1.382 1.375 1.368 1.372
>
> # apply inverse transformation to check calculations
> round(z2y(z, age, ref=nl4.bmi), 2)
[1] 18.41 18.15 17.89 17.72 17.55 17.49 17.42 17.49 17.55 17.74 17.92 18.18
[13] 18.44 18.77 19.10 19.47 19.84 20.20 20.55 20.89 21.22 21.57 21.91 22.27
[25] 22.62 22.96 23.29 23.60 23.90 24.18 24.46 24.73 25.00
> cutoff
[1] 18.41 18.15 17.89 17.72 17.55 17.49 17.42 17.49 17.55 17.74 17.92 18.18
[13] 18.44 18.77 19.10 19.47 19.84 20.20 20.55 20.89 21.22 21.57 21.91 22.27
[25] 22.62 22.96 23.29 23.60 23.90 24.18 24.46 24.73 25.00
>
> # calculate percentiles of weight 12 kg at 2 years (boys, girls)
> 100*round(pnorm(y2z(y=c(12,12), x=2, sex=c("M","F"), ref=nl4.wgt)),2)
[1] 24 41
>
> # # percentage of children lighter than 15kg at ages 2-5
> e <- expand.grid(age=2:5, sex=c("M","F"))
> z <- y2z(y=rep(15,nrow(e)), x=e$age, sex=e$sex, ref=nl4.wgt)
> w <- matrix(100*round(pnorm(z),2), nrow=2, byrow=TRUE)
> dimnames(w) <- list(c("boys","girls"),2:5)
> w
2 3 4 5
boys 89 46 11 1
girls 96 57 16 2
>
> # analysis in Z scale
> hgt.z <- y2z(y=boys$hgt, x=boys$age, sex="M", ref=nl4.hgt)
> wgt.z <- y2z(y=boys$wgt, x=boys$age, sex="M", ref=nl4.wgt)
> plot(hgt.z, wgt.z, col="blue")
>
>
> # z2y
>
> # quantile at SD=0 of age 2 years,
> # height Dutch boys
> z2y(z=0, x=2)
[1] 88.85
>
> # same for Dutch girls
> z2y(z=0, x=2, sex="F")
[1] 87.49
>
> # quantile at SD=c(-1,0,1) of age 2 years, BMI Dutch boys
> z2y(z=c(-1,0,+1), x=2, ref=nl4.bmi)
[1] 15.173 16.420 17.770
>
> # 0SD line (P50) in kg of weight for age in 5-10 year, Dutch boys
> z2y(z=rep(0,6), x=5:10, ref=nl4.wgt)
[1] 19.82 22.37 25.03 27.86 30.76 33.79
>
> # 95th percentile (P95), age 10 years, wfa, Dutch boys
> z2y(z=qnorm(0.95), x=10, ref=nl4.wgt)
[1] 45.278
>
> # table of P3, P10, P50, P90, P97 of weight for 5-10 year old dutch boys
> # age per year
> age <- 5:10
> p <- c(0.03,0.1,0.5,0.9,0.97)
> z <- rep(qnorm(p), length(age))
> x <- rep(age, each=length(p))
> w <- matrix(z2y(z, x=x, sex="M", ref=nl4.wgt), ncol=length(p),
+ byrow=TRUE)
> dimnames(w) <- list(age, p)
> round(w,1)
0.03 0.1 0.5 0.9 0.97
5 15.8 16.9 19.8 23.4 25.4
6 17.6 19.0 22.4 26.8 29.3
7 19.5 21.1 25.0 30.3 33.4
8 21.6 23.3 27.9 34.1 37.8
9 23.6 25.5 30.8 38.1 42.5
10 25.6 27.8 33.8 42.2 47.5
>
> # standard set of Z-scores of weight for all tabulated ages, boys & girls
> # and three etnicities
> sds <- c(-2.5, -2, -1, 0, 1, 2, 2.5)
> age <- nl4.wgt$x
> z <- rep(sds, times=length(age))
> x <- rep(age, each=length(sds))
> sex <- rep(c("M","F"), each=length(z)/2)
> w <- z2y(z=z, x=x, sex=sex, ref=nl4.wgt)
> w <- matrix(w, ncol=length(sds), byrow=TRUE)
> dimnames(w) <- list(age, sds)
> data.frame(sub=nl4.wgt$sub,sex=nl4.wgt$sex,round(w,2), row.names=NULL)
sub sex X.2.5 X.2 X.1 X0 X1 X2 X2.5
1 N M 2.58 2.77 3.15 3.55 3.97 4.40 4.62
2 N M 2.73 2.93 3.33 3.75 4.19 4.65 4.89
3 N M 2.87 3.08 3.50 3.94 4.40 4.89 5.14
4 N M 3.01 3.23 3.67 4.13 4.62 5.12 5.39
5 N M 3.16 3.38 3.84 4.33 4.84 5.37 5.65
6 N M 3.31 3.54 4.02 4.53 5.06 5.62 5.91
7 N M 3.46 3.70 4.20 4.72 5.27 5.85 6.15
8 N M 3.61 3.86 4.37 4.92 5.50 6.10 6.41
9 N M 3.75 4.01 4.54 5.11 5.71 6.34 6.66
10 N M 3.89 4.16 4.71 5.30 5.92 6.57 6.91
11 N M 4.03 4.30 4.87 5.48 6.12 6.80 7.15
12 N M 4.30 4.59 5.19 5.84 6.53 7.25 7.63
13 N M 4.56 4.86 5.51 6.19 6.92 7.69 8.09
14 N M 4.81 5.13 5.80 6.51 7.27 8.07 8.49
15 N M 5.05 5.38 6.07 6.82 7.62 8.46 8.90
16 N M 5.27 5.61 6.33 7.11 7.94 8.83 9.29
17 N M 5.48 5.83 6.58 7.39 8.25 9.18 9.66
18 N M 5.68 6.04 6.82 7.65 8.54 9.50 10.01
19 N M 5.86 6.23 7.03 7.89 8.81 9.80 10.33
20 N M 6.03 6.42 7.24 8.12 9.07 10.10 10.63
21 N M 6.38 6.78 7.63 8.55 9.54 10.62 11.18
22 N M 6.68 7.10 7.99 8.95 9.99 11.12 11.72
23 N M 6.95 7.38 8.30 9.30 10.39 11.56 12.19
24 N M 7.21 7.65 8.60 9.63 10.76 11.98 12.63
25 N M 7.45 7.90 8.88 9.94 11.10 12.37 13.04
26 N M 7.67 8.14 9.13 10.23 11.43 12.74 13.44
27 N M 7.86 8.34 9.37 10.50 11.74 13.11 13.83
28 N M 8.05 8.54 9.59 10.75 12.02 13.42 14.17
29 N M 8.24 8.74 9.82 11.00 12.30 13.74 14.51
30 N M 8.85 9.38 10.52 11.78 13.18 14.73 15.57
31 N M 9.79 10.36 11.62 13.02 14.59 16.35 17.30
32 N M 10.59 11.21 12.58 14.12 15.86 17.83 18.91
33 N M 11.39 12.06 13.52 15.19 17.09 19.26 20.45
34 N M 12.18 12.89 14.47 16.28 18.36 20.76 22.10
35 N M 12.98 13.74 15.45 17.42 19.72 22.40 23.91
36 N M 13.83 14.64 16.47 18.60 21.11 24.07 25.75
37 N M 14.67 15.54 17.50 19.82 22.57 25.87 27.77
38 N M 15.54 16.47 18.57 21.08 24.09 27.76 29.89
39 N M 16.39 17.38 19.64 22.37 25.69 29.79 32.21
40 N M 17.28 18.33 20.74 23.68 27.30 31.82 34.52
41 N M 18.14 19.26 21.85 25.03 28.99 34.02 37.07
42 N M 19.08 20.26 23.02 26.43 30.73 36.25 39.64
43 N M 19.97 21.23 24.18 27.86 32.54 38.65 42.44
44 N M 20.86 22.20 25.35 29.30 34.38 41.11 45.32
45 N M 21.78 23.20 26.54 30.76 36.22 43.51 48.12
46 N M 22.66 24.16 27.73 32.25 38.15 46.09 51.15
47 N M 23.59 25.18 28.97 33.79 40.10 48.62 54.07
48 N M 24.55 26.25 30.28 35.43 42.17 51.29 57.12
49 N M 25.53 27.36 31.69 37.23 44.48 54.28 60.53
50 N M 26.65 28.62 33.29 39.24 47.00 57.39 63.97
51 N M 27.85 29.99 35.06 41.48 49.79 60.77 67.64
52 N M 29.23 31.56 37.07 43.99 52.84 64.35 71.44
53 N M 30.74 33.29 39.29 46.77 56.22 68.31 75.64
54 N M 32.59 35.35 41.80 49.77 59.69 72.17 79.61
55 N M 34.62 37.58 44.46 52.87 63.22 76.06 83.62
56 N M 36.88 40.00 47.20 55.94 66.59 79.64 87.25
57 N M 39.20 42.43 49.87 58.82 69.65 82.80 90.42
58 N M 41.54 44.84 52.41 61.45 72.32 85.43 92.98
59 N M 43.65 47.01 54.66 63.77 74.67 87.77 95.30
60 N M 45.58 48.96 56.64 65.75 76.61 89.61 97.06
61 N M 47.22 50.62 58.32 67.43 78.26 91.18 98.58
62 N M 48.71 52.11 59.79 68.85 79.59 92.38 99.68
63 N M 49.95 53.35 61.02 70.06 80.75 93.45 100.70
64 N M 51.07 54.47 62.12 71.10 81.70 94.27 101.43
65 N M 52.01 55.41 63.06 72.04 82.62 95.15 102.28
66 N M 52.90 56.30 63.94 72.89 83.43 95.89 102.98
67 N M 53.88 57.26 64.84 73.70 84.10 96.37 103.33
68 N M 54.73 58.10 65.67 74.49 84.83 97.02 103.92
69 N M 55.58 58.95 66.50 75.28 85.56 97.65 104.50
70 N M 55.58 58.95 66.50 75.28 85.56 97.65 104.50
71 N F 2.58 2.77 3.15 3.55 3.97 4.40 4.62
72 N F 2.73 2.93 3.33 3.75 4.19 4.65 4.89
73 N F 2.87 3.08 3.50 3.94 4.40 4.89 5.14
74 N F 3.01 3.23 3.67 4.13 4.62 5.12 5.39
75 N F 3.16 3.38 3.84 4.33 4.84 5.37 5.65
76 N F 3.31 3.54 4.02 4.53 5.06 5.62 5.91
77 N F 3.46 3.70 4.20 4.72 5.27 5.85 6.15
78 N F 3.61 3.86 4.37 4.92 5.50 6.10 6.41
79 N F 3.75 4.01 4.54 5.11 5.71 6.34 6.66
80 N F 3.89 4.16 4.71 5.30 5.92 6.57 6.91
81 N F 4.03 4.30 4.87 5.48 6.12 6.80 7.15
82 N F 4.30 4.59 5.19 5.84 6.53 7.25 7.63
83 N F 4.56 4.86 5.51 6.19 6.92 7.69 8.09
84 N F 4.81 5.13 5.80 6.51 7.27 8.07 8.49
85 N F 5.05 5.38 6.07 6.82 7.62 8.46 8.90
86 N F 5.27 5.61 6.33 7.11 7.94 8.83 9.29
87 N F 5.48 5.83 6.58 7.39 8.25 9.18 9.66
88 N F 5.68 6.04 6.82 7.65 8.54 9.50 10.01
89 N F 5.86 6.23 7.03 7.89 8.81 9.80 10.33
90 N F 6.03 6.42 7.24 8.12 9.07 10.10 10.63
91 N F 6.38 6.78 7.63 8.55 9.54 10.62 11.18
92 N F 6.68 7.10 7.99 8.95 9.99 11.12 11.72
93 N F 6.95 7.38 8.30 9.30 10.39 11.56 12.19
94 N F 7.21 7.65 8.60 9.63 10.76 11.98 12.63
95 N F 7.45 7.90 8.88 9.94 11.10 12.37 13.04
96 N F 7.67 8.14 9.13 10.23 11.43 12.74 13.44
97 N F 7.86 8.34 9.37 10.50 11.74 13.11 13.83
98 N F 8.05 8.54 9.59 10.75 12.02 13.42 14.17
99 N F 8.24 8.74 9.82 11.00 12.30 13.74 14.51
100 N F 8.85 9.38 10.52 11.78 13.18 14.73 15.57
101 N F 9.79 10.36 11.62 13.02 14.59 16.35 17.30
102 N F 10.59 11.21 12.58 14.12 15.86 17.83 18.91
103 N F 11.39 12.06 13.52 15.19 17.09 19.26 20.45
104 N F 12.18 12.89 14.47 16.28 18.36 20.76 22.10
105 N F 12.98 13.74 15.45 17.42 19.72 22.40 23.91
106 N F 13.83 14.64 16.47 18.60 21.11 24.07 25.75
107 N F 14.67 15.54 17.50 19.82 22.57 25.87 27.77
108 N F 15.54 16.47 18.57 21.08 24.09 27.76 29.89
109 N F 16.39 17.38 19.64 22.37 25.69 29.79 32.21
110 N F 17.28 18.33 20.74 23.68 27.30 31.82 34.52
111 N F 18.14 19.26 21.85 25.03 28.99 34.02 37.07
112 N F 19.08 20.26 23.02 26.43 30.73 36.25 39.64
113 N F 19.97 21.23 24.18 27.86 32.54 38.65 42.44
114 N F 20.86 22.20 25.35 29.30 34.38 41.11 45.32
115 N F 21.78 23.20 26.54 30.76 36.22 43.51 48.12
116 N F 22.66 24.16 27.73 32.25 38.15 46.09 51.15
117 N F 23.59 25.18 28.97 33.79 40.10 48.62 54.07
118 N F 24.55 26.25 30.28 35.43 42.17 51.29 57.12
119 N F 25.53 27.36 31.69 37.23 44.48 54.28 60.53
120 N F 26.65 28.62 33.29 39.24 47.00 57.39 63.97
121 N F 27.85 29.99 35.06 41.48 49.79 60.77 67.64
122 N F 29.23 31.56 37.07 43.99 52.84 64.35 71.44
123 N F 30.74 33.29 39.29 46.77 56.22 68.31 75.64
124 N F 32.59 35.35 41.80 49.77 59.69 72.17 79.61
125 N F 34.62 37.58 44.46 52.87 63.22 76.06 83.62
126 N F 36.88 40.00 47.20 55.94 66.59 79.64 87.25
127 N F 39.20 42.43 49.87 58.82 69.65 82.80 90.42
128 N F 41.54 44.84 52.41 61.45 72.32 85.43 92.98
129 N F 43.65 47.01 54.66 63.77 74.67 87.77 95.30
130 N F 45.58 48.96 56.64 65.75 76.61 89.61 97.06
131 N F 47.22 50.62 58.32 67.43 78.26 91.18 98.58
132 N F 48.71 52.11 59.79 68.85 79.59 92.38 99.68
133 N F 49.95 53.35 61.02 70.06 80.75 93.45 100.70
134 N F 51.07 54.47 62.12 71.10 81.70 94.27 101.43
135 N F 52.01 55.41 63.06 72.04 82.62 95.15 102.28
136 N F 52.90 56.30 63.94 72.89 83.43 95.89 102.98
137 N F 53.88 57.26 64.84 73.70 84.10 96.37 103.33
138 N F 54.73 58.10 65.67 74.49 84.83 97.02 103.92
139 N F 55.58 58.95 66.50 75.28 85.56 97.65 104.50
140 N F 55.58 58.95 66.50 75.28 85.56 97.65 104.50
141 T M 2.58 2.77 3.15 3.55 3.97 4.40 4.62
142 T M 2.73 2.93 3.33 3.75 4.19 4.65 4.89
143 T M 2.87 3.08 3.50 3.94 4.40 4.89 5.14
144 T M 3.01 3.23 3.67 4.13 4.62 5.12 5.39
145 T M 3.16 3.38 3.84 4.33 4.84 5.37 5.65
146 T M 3.31 3.54 4.02 4.53 5.06 5.62 5.91
147 T M 3.46 3.70 4.20 4.72 5.27 5.85 6.15
148 T M 3.61 3.86 4.37 4.92 5.50 6.10 6.41
149 T M 3.75 4.01 4.54 5.11 5.71 6.34 6.66
150 T M 3.89 4.16 4.71 5.30 5.92 6.57 6.91
151 T M 4.03 4.30 4.87 5.48 6.12 6.80 7.15
152 T M 4.30 4.59 5.19 5.84 6.53 7.25 7.63
153 T M 4.56 4.86 5.51 6.19 6.92 7.69 8.09
154 T M 4.81 5.13 5.80 6.51 7.27 8.07 8.49
155 T M 5.05 5.38 6.07 6.82 7.62 8.46 8.90
156 T M 5.27 5.61 6.33 7.11 7.94 8.83 9.29
157 T M 5.48 5.83 6.58 7.39 8.25 9.18 9.66
158 T M 5.68 6.04 6.82 7.65 8.54 9.50 10.01
159 T M 5.86 6.23 7.03 7.89 8.81 9.80 10.33
160 T M 6.03 6.42 7.24 8.12 9.07 10.10 10.63
161 T M 6.38 6.78 7.63 8.55 9.54 10.62 11.18
162 T M 6.68 7.10 7.99 8.95 9.99 11.12 11.72
163 T M 6.95 7.38 8.30 9.30 10.39 11.56 12.19
164 T M 7.21 7.65 8.60 9.63 10.76 11.98 12.63
165 T M 7.45 7.90 8.88 9.94 11.10 12.37 13.04
166 T M 7.67 8.14 9.13 10.23 11.43 12.74 13.44
167 T M 7.86 8.34 9.37 10.50 11.74 13.11 13.83
168 T M 8.05 8.54 9.59 10.75 12.02 13.42 14.17
169 T M 8.24 8.74 9.82 11.00 12.30 13.74 14.51
170 T M 8.85 9.38 10.52 11.78 13.18 14.73 15.57
171 T M 9.79 10.36 11.62 13.02 14.59 16.35 17.30
172 T M 10.59 11.21 12.58 14.12 15.86 17.83 18.91
173 T M 11.39 12.06 13.52 15.19 17.09 19.26 20.45
174 T M 12.18 12.89 14.47 16.28 18.36 20.76 22.10
175 T M 12.98 13.74 15.45 17.42 19.72 22.40 23.91
176 T M 13.83 14.64 16.47 18.60 21.11 24.07 25.75
177 T M 14.67 15.54 17.50 19.82 22.57 25.87 27.77
178 T M 15.54 16.47 18.57 21.08 24.09 27.76 29.89
179 T M 16.39 17.38 19.64 22.37 25.69 29.79 32.21
180 T M 17.28 18.33 20.74 23.68 27.30 31.82 34.52
181 T M 18.14 19.26 21.85 25.03 28.99 34.02 37.07
182 T M 19.08 20.26 23.02 26.43 30.73 36.25 39.64
183 T M 19.97 21.23 24.18 27.86 32.54 38.65 42.44
184 T M 20.86 22.20 25.35 29.30 34.38 41.11 45.32
185 T M 21.78 23.20 26.54 30.76 36.22 43.51 48.12
186 T M 22.66 24.16 27.73 32.25 38.15 46.09 51.15
187 T M 23.59 25.18 28.97 33.79 40.10 48.62 54.07
188 T M 24.55 26.25 30.28 35.43 42.17 51.29 57.12
189 T M 25.53 27.36 31.69 37.23 44.48 54.28 60.53
190 T M 26.65 28.62 33.29 39.24 47.00 57.39 63.97
191 T M 27.85 29.99 35.06 41.48 49.79 60.77 67.64
192 T M 29.23 31.56 37.07 43.99 52.84 64.35 71.44
193 T M 30.74 33.29 39.29 46.77 56.22 68.31 75.64
194 T M 32.59 35.35 41.80 49.77 59.69 72.17 79.61
195 T M 34.62 37.58 44.46 52.87 63.22 76.06 83.62
196 T M 36.88 40.00 47.20 55.94 66.59 79.64 87.25
197 T M 39.20 42.43 49.87 58.82 69.65 82.80 90.42
198 T M 41.54 44.84 52.41 61.45 72.32 85.43 92.98
199 T M 43.65 47.01 54.66 63.77 74.67 87.77 95.30
200 T M 45.58 48.96 56.64 65.75 76.61 89.61 97.06
201 T M 47.22 50.62 58.32 67.43 78.26 91.18 98.58
202 T M 48.71 52.11 59.79 68.85 79.59 92.38 99.68
203 T M 49.95 53.35 61.02 70.06 80.75 93.45 100.70
204 T M 51.07 54.47 62.12 71.10 81.70 94.27 101.43
205 T M 52.01 55.41 63.06 72.04 82.62 95.15 102.28
206 T M 52.90 56.30 63.94 72.89 83.43 95.89 102.98
207 T M 53.88 57.26 64.84 73.70 84.10 96.37 103.33
208 T M 48.33 50.76 56.51 63.85 73.58 87.13 96.14
209 T F 2.46 2.63 2.97 3.34 3.74 4.17 4.39
210 T F 2.59 2.76 3.12 3.51 3.93 4.38 4.62
211 T F 2.73 2.90 3.28 3.69 4.13 4.60 4.85
212 T F 2.85 3.04 3.43 3.86 4.32 4.82 5.08
213 T F 2.98 3.17 3.58 4.03 4.51 5.03 5.30
214 T F 3.11 3.31 3.75 4.21 4.71 5.26 5.54
215 T F 3.25 3.46 3.90 4.38 4.90 5.46 5.76
216 T F 3.38 3.60 4.06 4.56 5.10 5.68 5.99
217 T F 3.51 3.74 4.21 4.73 5.29 5.90 6.22
218 T F 3.65 3.88 4.37 4.90 5.48 6.10 6.43
219 T F 3.78 4.01 4.52 5.07 5.67 6.31 6.66
220 T F 4.03 4.28 4.82 5.40 6.04 6.73 7.09
221 T F 4.28 4.54 5.11 5.72 6.39 7.12 7.50
222 T F 4.52 4.79 5.38 6.03 6.74 7.50 7.91
223 T F 4.74 5.03 5.64 6.32 7.06 7.87 8.30
224 T F 4.96 5.25 5.89 6.59 7.36 8.19 8.64
225 T F 5.16 5.46 6.12 6.85 7.65 8.52 8.98
226 T F 5.35 5.66 6.35 7.10 7.93 8.83 9.32
227 T F 5.54 5.86 6.56 7.33 8.18 9.10 9.60
228 T F 5.71 6.04 6.76 7.55 8.42 9.38 9.89
229 T F 6.02 6.37 7.12 7.95 8.87 9.88 10.43
230 T F 6.30 6.66 7.45 8.32 9.28 10.35 10.92
231 T F 6.56 6.94 7.75 8.66 9.66 10.78 11.38
232 T F 6.82 7.21 8.04 8.97 10.00 11.14 11.76
233 T F 7.04 7.44 8.30 9.27 10.35 11.55 12.20
234 T F 7.25 7.66 8.55 9.55 10.66 11.90 12.57
235 T F 7.46 7.88 8.79 9.81 10.95 12.23 12.92
236 T F 7.65 8.08 9.01 10.06 11.23 12.54 13.26
237 T F 7.84 8.27 9.23 10.30 11.50 12.85 13.58
238 T F 8.41 8.88 9.90 11.06 12.37 13.84 14.65
239 T F 9.35 9.87 11.02 12.32 13.81 15.51 16.45
240 T F 10.23 10.80 12.07 13.52 15.19 17.12 18.20
241 T F 11.06 11.68 13.07 14.68 16.55 18.74 19.97
242 T F 11.86 12.54 14.05 15.82 17.90 20.36 21.77
243 T F 12.64 13.37 15.01 16.94 19.24 21.99 23.57
244 T F 13.38 14.16 15.94 18.05 20.59 23.68 25.48
245 T F 14.15 14.99 16.90 19.21 22.02 25.48 27.53
246 T F 14.96 15.86 17.93 20.44 23.54 27.42 29.75
247 T F 15.80 16.77 19.01 21.77 25.22 29.62 32.29
248 T F 16.72 17.76 20.18 23.19 27.00 31.92 34.96
249 T F 17.61 18.73 21.37 24.68 28.93 34.55 38.07
250 T F 18.55 19.76 22.61 26.22 30.92 37.24 41.26
251 T F 19.50 20.79 23.87 27.80 32.98 40.03 44.59
252 T F 20.39 21.78 25.11 29.40 35.12 43.02 48.19
253 T F 21.31 22.80 26.39 31.04 37.28 46.00 51.76
254 T F 22.24 23.85 27.71 32.75 39.55 49.10 55.43
255 T F 23.21 24.94 29.09 34.52 41.85 52.15 58.97
256 T F 24.16 26.03 30.52 36.40 44.34 55.48 62.85
257 T F 25.28 27.29 32.13 38.46 46.96 58.77 66.50
258 T F 26.52 28.70 33.93 40.73 49.78 62.19 70.19
259 T F 28.05 30.39 36.00 43.21 52.69 65.44 73.53
260 T F 29.93 32.41 38.30 45.80 55.49 68.27 76.23
261 T F 32.00 34.58 40.67 48.32 58.09 70.76 78.54
262 T F 34.15 36.78 42.95 50.64 60.37 72.87 80.48
263 T F 36.26 38.89 45.05 52.70 62.33 74.66 82.15
264 T F 38.19 40.81 46.91 54.48 64.00 76.21 83.65
265 T F 39.83 42.43 48.48 56.00 65.49 77.74 85.24
266 T F 41.24 43.81 49.82 57.29 66.77 79.09 86.71
267 T F 42.42 44.98 50.93 58.37 67.86 80.31 88.08
268 T F 43.48 46.00 51.90 59.28 68.74 81.24 89.10
269 T F 44.23 46.75 52.64 60.04 69.59 82.33 90.42
270 T F 44.99 47.48 53.33 60.70 70.24 83.04 91.22
271 T F 45.59 48.07 53.90 61.27 70.85 83.81 92.14
272 T F 46.14 48.61 54.43 61.79 71.41 84.48 92.95
273 T F 46.63 49.10 54.90 62.25 71.89 85.06 93.64
274 T F 47.11 49.56 55.34 62.68 72.33 85.57 94.24
275 T F 47.55 49.99 55.75 63.08 72.73 86.03 94.77
276 T F 47.90 50.34 56.11 63.47 73.20 86.71 95.66
277 T F 48.33 50.76 56.51 63.85 73.58 87.13 96.14
278 M M 2.46 2.63 2.97 3.34 3.74 4.17 4.39
279 M M 2.59 2.76 3.12 3.51 3.93 4.38 4.62
280 M M 2.73 2.90 3.28 3.69 4.13 4.60 4.85
281 M M 2.85 3.04 3.43 3.86 4.32 4.82 5.08
282 M M 2.98 3.17 3.58 4.03 4.51 5.03 5.30
283 M M 3.11 3.31 3.75 4.21 4.71 5.26 5.54
284 M M 3.25 3.46 3.90 4.38 4.90 5.46 5.76
285 M M 3.38 3.60 4.06 4.56 5.10 5.68 5.99
286 M M 3.51 3.74 4.21 4.73 5.29 5.90 6.22
287 M M 3.65 3.88 4.37 4.90 5.48 6.10 6.43
288 M M 3.78 4.01 4.52 5.07 5.67 6.31 6.66
289 M M 4.03 4.28 4.82 5.40 6.04 6.73 7.09
290 M M 4.28 4.54 5.11 5.72 6.39 7.12 7.50
291 M M 4.52 4.79 5.38 6.03 6.74 7.50 7.91
292 M M 4.74 5.03 5.64 6.32 7.06 7.87 8.30
293 M M 4.96 5.25 5.89 6.59 7.36 8.19 8.64
294 M M 5.16 5.46 6.12 6.85 7.65 8.52 8.98
295 M M 5.35 5.66 6.35 7.10 7.93 8.83 9.32
296 M M 5.54 5.86 6.56 7.33 8.18 9.10 9.60
297 M M 5.71 6.04 6.76 7.55 8.42 9.38 9.89
298 M M 6.02 6.37 7.12 7.95 8.87 9.88 10.43
299 M M 6.30 6.66 7.45 8.32 9.28 10.35 10.92
300 M M 6.56 6.94 7.75 8.66 9.66 10.78 11.38
301 M M 6.82 7.21 8.04 8.97 10.00 11.14 11.76
302 M M 7.04 7.44 8.30 9.27 10.35 11.55 12.20
303 M M 7.25 7.66 8.55 9.55 10.66 11.90 12.57
304 M M 7.46 7.88 8.79 9.81 10.95 12.23 12.92
305 M M 7.65 8.08 9.01 10.06 11.23 12.54 13.26
306 M M 7.84 8.27 9.23 10.30 11.50 12.85 13.58
307 M M 8.41 8.88 9.90 11.06 12.37 13.84 14.65
308 M M 9.35 9.87 11.02 12.32 13.81 15.51 16.45
309 M M 10.23 10.80 12.07 13.52 15.19 17.12 18.20
310 M M 11.06 11.68 13.07 14.68 16.55 18.74 19.97
311 M M 11.86 12.54 14.05 15.82 17.90 20.36 21.77
312 M M 12.64 13.37 15.01 16.94 19.24 21.99 23.57
313 M M 13.38 14.16 15.94 18.05 20.59 23.68 25.48
314 M M 14.15 14.99 16.90 19.21 22.02 25.48 27.53
315 M M 14.96 15.86 17.93 20.44 23.54 27.42 29.75
316 M M 15.80 16.77 19.01 21.77 25.22 29.62 32.29
317 M M 16.72 17.76 20.18 23.19 27.00 31.92 34.96
318 M M 17.61 18.73 21.37 24.68 28.93 34.55 38.07
319 M M 18.55 19.76 22.61 26.22 30.92 37.24 41.26
320 M M 19.50 20.79 23.87 27.80 32.98 40.03 44.59
321 M M 20.39 21.78 25.11 29.40 35.12 43.02 48.19
322 M M 21.31 22.80 26.39 31.04 37.28 46.00 51.76
323 M M 22.24 23.85 27.71 32.75 39.55 49.10 55.43
324 M M 23.21 24.94 29.09 34.52 41.85 52.15 58.97
325 M M 24.16 26.03 30.52 36.40 44.34 55.48 62.85
326 M M 25.28 27.29 32.13 38.46 46.96 58.77 66.50
327 M M 26.52 28.70 33.93 40.73 49.78 62.19 70.19
328 M M 28.05 30.39 36.00 43.21 52.69 65.44 73.53
329 M M 29.93 32.41 38.30 45.80 55.49 68.27 76.23
330 M M 32.00 34.58 40.67 48.32 58.09 70.76 78.54
331 M M 34.15 36.78 42.95 50.64 60.37 72.87 80.48
332 M M 36.26 38.89 45.05 52.70 62.33 74.66 82.15
333 M M 38.19 40.81 46.91 54.48 64.00 76.21 83.65
334 M M 39.83 42.43 48.48 56.00 65.49 77.74 85.24
335 M M 41.24 43.81 49.82 57.29 66.77 79.09 86.71
336 M M 42.42 44.98 50.93 58.37 67.86 80.31 88.08
337 M M 43.48 46.00 51.90 59.28 68.74 81.24 89.10
338 M M 44.23 46.75 52.64 60.04 69.59 82.33 90.42
339 M M 44.99 47.48 53.33 60.70 70.24 83.04 91.22
340 M M 45.59 48.07 53.90 61.27 70.85 83.81 92.14
341 M M 46.14 48.61 54.43 61.79 71.41 84.48 92.95
342 M M 46.63 49.10 54.90 62.25 71.89 85.06 93.64
343 M M 47.11 49.56 55.34 62.68 72.33 85.57 94.24
344 M M 47.55 49.99 55.75 63.08 72.73 86.03 94.77
345 M M 48.33 50.76 56.51 63.85 73.58 87.13 96.14
346 M F 2.46 2.63 2.97 3.34 3.74 4.17 4.39
347 M F 2.59 2.76 3.12 3.51 3.93 4.38 4.62
348 M F 2.73 2.90 3.28 3.69 4.13 4.60 4.85
349 M F 2.85 3.04 3.43 3.86 4.32 4.82 5.08
350 M F 2.98 3.17 3.58 4.03 4.51 5.03 5.30
351 M F 3.11 3.31 3.75 4.21 4.71 5.26 5.54
352 M F 3.25 3.46 3.90 4.38 4.90 5.46 5.76
353 M F 3.38 3.60 4.06 4.56 5.10 5.68 5.99
354 M F 3.51 3.74 4.21 4.73 5.29 5.90 6.22
355 M F 3.65 3.88 4.37 4.90 5.48 6.10 6.43
356 M F 3.78 4.01 4.52 5.07 5.67 6.31 6.66
357 M F 4.03 4.28 4.82 5.40 6.04 6.73 7.09
358 M F 4.28 4.54 5.11 5.72 6.39 7.12 7.50
359 M F 4.52 4.79 5.38 6.03 6.74 7.50 7.91
360 M F 4.74 5.03 5.64 6.32 7.06 7.87 8.30
361 M F 4.96 5.25 5.89 6.59 7.36 8.19 8.64
362 M F 5.16 5.46 6.12 6.85 7.65 8.52 8.98
363 M F 5.35 5.66 6.35 7.10 7.93 8.83 9.32
364 M F 5.54 5.86 6.56 7.33 8.18 9.10 9.60
365 M F 5.71 6.04 6.76 7.55 8.42 9.38 9.89
366 M F 6.02 6.37 7.12 7.95 8.87 9.88 10.43
367 M F 6.30 6.66 7.45 8.32 9.28 10.35 10.92
368 M F 6.56 6.94 7.75 8.66 9.66 10.78 11.38
369 M F 6.82 7.21 8.04 8.97 10.00 11.14 11.76
370 M F 7.04 7.44 8.30 9.27 10.35 11.55 12.20
371 M F 7.25 7.66 8.55 9.55 10.66 11.90 12.57
372 M F 7.46 7.88 8.79 9.81 10.95 12.23 12.92
373 M F 7.65 8.08 9.01 10.06 11.23 12.54 13.26
374 M F 7.84 8.27 9.23 10.30 11.50 12.85 13.58
375 M F 8.41 8.88 9.90 11.06 12.37 13.84 14.65
376 M F 9.35 9.87 11.02 12.32 13.81 15.51 16.45
377 M F 10.23 10.80 12.07 13.52 15.19 17.12 18.20
378 M F 11.06 11.68 13.07 14.68 16.55 18.74 19.97
379 M F 11.86 12.54 14.05 15.82 17.90 20.36 21.77
380 M F 12.64 13.37 15.01 16.94 19.24 21.99 23.57
381 M F 13.38 14.16 15.94 18.05 20.59 23.68 25.48
382 M F 14.15 14.99 16.90 19.21 22.02 25.48 27.53
383 M F 14.96 15.86 17.93 20.44 23.54 27.42 29.75
384 M F 15.80 16.77 19.01 21.77 25.22 29.62 32.29
385 M F 16.72 17.76 20.18 23.19 27.00 31.92 34.96
386 M F 17.61 18.73 21.37 24.68 28.93 34.55 38.07
387 M F 18.55 19.76 22.61 26.22 30.92 37.24 41.26
388 M F 19.50 20.79 23.87 27.80 32.98 40.03 44.59
389 M F 20.39 21.78 25.11 29.40 35.12 43.02 48.19
390 M F 21.31 22.80 26.39 31.04 37.28 46.00 51.76
391 M F 22.24 23.85 27.71 32.75 39.55 49.10 55.43
392 M F 23.21 24.94 29.09 34.52 41.85 52.15 58.97
393 M F 24.16 26.03 30.52 36.40 44.34 55.48 62.85
394 M F 25.28 27.29 32.13 38.46 46.96 58.77 66.50
395 M F 26.52 28.70 33.93 40.73 49.78 62.19 70.19
396 M F 28.05 30.39 36.00 43.21 52.69 65.44 73.53
397 M F 29.93 32.41 38.30 45.80 55.49 68.27 76.23
398 M F 32.00 34.58 40.67 48.32 58.09 70.76 78.54
399 M F 34.15 36.78 42.95 50.64 60.37 72.87 80.48
400 M F 36.26 38.89 45.05 52.70 62.33 74.66 82.15
401 M F 38.19 40.81 46.91 54.48 64.00 76.21 83.65
402 M F 39.83 42.43 48.48 56.00 65.49 77.74 85.24
403 M F 41.24 43.81 49.82 57.29 66.77 79.09 86.71
404 M F 42.42 44.98 50.93 58.37 67.86 80.31 88.08
405 M F 43.48 46.00 51.90 59.28 68.74 81.24 89.10
406 M F 44.23 46.75 52.64 60.04 69.59 82.33 90.42
407 M F 44.99 47.48 53.33 60.70 70.24 83.04 91.22
408 M F 45.59 48.07 53.90 61.27 70.85 83.81 92.14
409 M F 46.14 48.61 54.43 61.79 71.41 84.48 92.95
410 M F 46.63 49.10 54.90 62.25 71.89 85.06 93.64
411 M F 47.11 49.56 55.34 62.68 72.33 85.57 94.24
412 M F 47.55 49.99 55.75 63.08 72.73 86.03 94.77
413 M F 47.90 50.34 56.11 63.47 73.20 86.71 95.66
414 M F 48.33 50.76 56.51 63.85 73.58 87.13 96.14
>
> # P85 of BMI in 5-8 year old Dutch boys and girls
> e <- expand.grid(age=5:8, sex=c("M","F"))
> w <- z2y(z=rep(qnorm(0.85),nrow(e)), x=e$age, sex=e$sex, ref=nl4.bmi)
> w <- matrix(w, nrow=2, byrow=TRUE)
> dimnames(w) <- list(c("boys","girls"),5:8)
> w
5 6 7 8
boys 17.152 17.282 17.508 17.857
girls 17.090 17.386 17.837 18.327
>
> # data transformation of height z-scores to cm-scale
> z <- c(-1.83, 0.09, 2.33, 0.81, -1.20)
> x <- c(8.33, 0.23, 19.2, 24.3, 10)
> sex <- c("M", "M", "F", "M", "F")
> round(z2y(z=z, x=x, sex=sex, ref=nl4.hgt), 1)
[1] 123.8 60.7 185.3 189.7 135.6
>
> # interpolate published height standard
> # to daily values, days 0-31, boys
> # on centiles -2SD, 0SD and +2SD
> days <- 0:31
> sds <- c(-2, 0, +2)
> z <- rep(sds, length(days))
> x <- rep(round(days/365.25,4), each=length(sds))
> w <- z2y(z, x, sex="M", ref=nl4.hgt)
> w <- matrix(w, ncol=length(sds), byrow=TRUE)
> dimnames(w) <- list(days, sds)
> w
-2 0 2
0 47.163 51.320 55.477
1 47.265 51.430 55.594
2 47.371 51.543 55.716
3 47.474 51.653 55.833
4 47.580 51.767 55.954
5 47.682 51.877 56.071
6 47.784 51.986 56.188
7 47.890 52.100 56.310
8 47.993 52.210 56.427
9 48.096 52.321 56.545
10 48.203 52.435 56.667
11 48.305 52.545 56.785
12 48.412 52.659 56.907
13 48.515 52.770 57.024
14 48.618 52.880 57.142
15 48.727 52.995 57.264
16 48.832 53.106 57.381
17 48.937 53.217 57.498
18 49.046 53.333 57.619
19 49.151 53.444 57.736
20 49.260 53.559 57.857
21 49.366 53.670 57.974
22 49.469 53.781 58.093
23 49.577 53.896 58.216
24 49.681 54.007 58.334
25 49.784 54.118 58.453
26 49.892 54.234 58.576
27 49.996 54.345 58.694
28 50.103 54.460 58.817
29 50.213 54.573 58.933
30 50.324 54.686 59.049
31 50.438 54.803 59.169
>
>
>
>
>
> dev.off()
null device
1
>