Running 'test-bmi_for_age.R' [208s]
  Running 'test-bp_distributions.R' [5s]
  Running 'test-bp_parameters.R' [0s]
  Running 'test-est_norm.R' [0s]
  Running 'test-gemelli1990.R' [156s]
  Running 'test-head_circ_for_age.R' [6s]
  Running 'test-length_for_age_inf.R' [173s]
  Running 'test-lo2013.R' [0s]
  Running 'test-nhlbi_bp_norms.R' [0s]
  Running 'test-stature_for_age.R' [202s]
  Running 'test-weight_for_age.R' [31s]
  Running 'test-weight_for_age_inf.R' [6s]
  Running 'test-weight_for_length_inf.R' [175s]
  Running 'test-weight_for_stature.R' [8s]
Running the tests in 'tests/test-bmi_for_age.R' failed.
Complete output:
  > library(pedbp)
  > 
  > d <- pedbp:::cdc_lms_data
  > 
  > 
  > # Testing p_bmi_for_age {{{
  > test_p03 <-
  +   p_bmi_for_age(
  +            q    = d[d$set == "bmi_for_age", "p03"]
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > test_p05 <-
  +   p_bmi_for_age(
  +            q    = d[d$set == "bmi_for_age", "p05"]
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > test_p10 <-
  +   p_bmi_for_age(
  +            q    = d[d$set == "bmi_for_age", "p10"]
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > test_p25 <-
  +   p_bmi_for_age(
  +            q    = d[d$set == "bmi_for_age", "p25"]
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > test_p50 <-
  +   p_bmi_for_age(
  +            q    = d[d$set == "bmi_for_age", "p50"]
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > test_p75 <-
  +   p_bmi_for_age(
  +            q    = d[d$set == "bmi_for_age", "p75"]
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > test_p85 <-
  +   p_bmi_for_age(
  +            q    = d[d$set == "bmi_for_age", "p85"]
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > test_p90 <-
  +   p_bmi_for_age(
  +            q    = d[d$set == "bmi_for_age", "p90"]
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > test_p95 <-
  +   p_bmi_for_age(
  +            q    = d[d$set == "bmi_for_age", "p95"]
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > test_p97 <-
  +   p_bmi_for_age(
  +            q    = d[d$set == "bmi_for_age", "p97"]
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > 
  > stopifnot(all(round(test_p03, 2) == 0.03))
  > stopifnot(all(round(test_p05, 2) == 0.05))
  > stopifnot(all(round(test_p10, 2) == 0.10))
  > stopifnot(all(round(test_p25, 2) == 0.25))
  > stopifnot(all(round(test_p50, 2) == 0.50))
  > stopifnot(all(round(test_p75, 2) == 0.75))
  > stopifnot(all(round(test_p85, 2) == 0.85))
  > stopifnot(all(round(test_p90, 2) == 0.90))
  > stopifnot(all(round(test_p95, 2) == 0.95))
  > stopifnot(all(round(test_p97, 2) == 0.97))
  > 
  > # }}}
  > 
  > # Tesing q_bmi_for_age {{{
  > test_q03 <-
  +   q_bmi_for_age(
  +            p    = c(0.03)
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > test_q05 <-
  +   q_bmi_for_age(
  +            p    = c(0.05)
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > test_q10 <-
  +   q_bmi_for_age(
  +            p    = c(0.10)
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > test_q25 <-
  +   q_bmi_for_age(
  +            p    = c(0.25)
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > test_q50 <-
  +   q_bmi_for_age(
  +            p    = c(0.50)
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > test_q75 <-
  +   q_bmi_for_age(
  +            p    = c(0.75)
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > test_q85 <-
  +   q_bmi_for_age(
  +            p    = c(0.85)
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > test_q90 <-
  +   q_bmi_for_age(
  +            p    = c(0.90)
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > test_q95 <-
  +   q_bmi_for_age(
  +            p    = c(0.95)
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > test_q97 <-
  +   q_bmi_for_age(
  +            p    = c(0.97)
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > stopifnot(all.equal(test_q03, d[d$set == "bmi_for_age", "p03"]))
  > stopifnot(all.equal(test_q05, d[d$set == "bmi_for_age", "p05"]))
  > stopifnot(all.equal(test_q10, d[d$set == "bmi_for_age", "p10"]))
  > stopifnot(all.equal(test_q25, d[d$set == "bmi_for_age", "p25"]))
  > stopifnot(all.equal(test_q50, d[d$set == "bmi_for_age", "p50"]))
  > stopifnot(all.equal(test_q75, d[d$set == "bmi_for_age", "p75"]))
  > stopifnot(all.equal(test_q85, d[d$set == "bmi_for_age", "p85"]))
  > stopifnot(all.equal(test_q90, d[d$set == "bmi_for_age", "p90"]))
  > stopifnot(all.equal(test_q95, d[d$set == "bmi_for_age", "p95"]))
  > stopifnot(all.equal(test_q97, d[d$set == "bmi_for_age", "p97"]))
  > 
  > # end testing q_bmi_for_age }}}
  > 
  > # Testing z_bmi_for_age {{{
  > test_z03 <-
  +   z_bmi_for_age(
  +            q    = d[d$set == "bmi_for_age", "p03"]
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > test_z05 <-
  +   z_bmi_for_age(
  +            q    = d[d$set == "bmi_for_age", "p05"]
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > test_z10 <-
  +   z_bmi_for_age(
  +            q    = d[d$set == "bmi_for_age", "p10"]
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > test_z25 <-
  +   z_bmi_for_age(
  +            q    = d[d$set == "bmi_for_age", "p25"]
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > test_z50 <-
  +   z_bmi_for_age(
  +            q    = d[d$set == "bmi_for_age", "p50"]
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > test_z75 <-
  +   z_bmi_for_age(
  +            q    = d[d$set == "bmi_for_age", "p75"]
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > test_z85 <-
  +   z_bmi_for_age(
  +            q    = d[d$set == "bmi_for_age", "p85"]
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > test_z90 <-
  +   z_bmi_for_age(
  +            q    = d[d$set == "bmi_for_age", "p90"]
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > test_z95 <-
  +   z_bmi_for_age(
  +            q    = d[d$set == "bmi_for_age", "p95"]
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > test_z97 <-
  +   z_bmi_for_age(
  +            q    = d[d$set == "bmi_for_age", "p97"]
  +          , age  = d[d$set == "bmi_for_age", "age"]
  +          , male = d[d$set == "bmi_for_age", "male"]
  +          )
  > 
  > 
  > stopifnot(all.equal(test_z03, rep(qnorm(0.03), length(test_z03))))
  > stopifnot(all.equal(test_z05, rep(qnorm(0.05), length(test_z05))))
  > stopifnot(all.equal(test_z10, rep(qnorm(0.10), length(test_z10))))
  > stopifnot(all.equal(test_z25, rep(qnorm(0.25), length(test_z25))))
  > stopifnot(all.equal(test_z50, rep(qnorm(0.50), length(test_z50))))
  > stopifnot(all.equal(test_z75, rep(qnorm(0.75), length(test_z75))))
  > stopifnot(all.equal(test_z85, rep(qnorm(0.85), length(test_z85))))
  > stopifnot(all.equal(test_z90, rep(qnorm(0.90), length(test_z90))))
  > stopifnot(all.equal(test_z95, rep(qnorm(0.95), length(test_z95))))
  > stopifnot(all.equal(test_z97, rep(qnorm(0.97), length(test_z97))))
  > 
  > # }}}
  > 
  > 
  > proc.time()
     user  system elapsed 
    35.65    2.87   38.53 
Running the tests in 'tests/test-gemelli1990.R' failed.
Complete output:
  > library(pedbp)
  > 
  > stopifnot(class(gemelli1990) == "data.frame")
  > 
  > # expected names:
  > stopifnot(all(c("age", "male", "sbp_mean", "dbp_mean", "sbp_sd", "dbp_sd") %in% names(gemelli1990)))
  > 
  > # testing charateristics of the data set
  > stopifnot( dim(gemelli1990) == c(8, 6))
  > 
  > stopifnot(all(gemelli1990$male %in% c(0, 1)))
  > 
  > proc.time()
     user  system elapsed 
     0.17    0.09    0.25 
Running the tests in 'tests/test-length_for_age_inf.R' failed.
Complete output:
  > library(pedbp)
  > 
  > d <- pedbp:::cdc_lms_data
  > 
  > 
  > # Testing p_length_for_age_inf {{{
  > test_p03 <-
  +   p_length_for_age_inf(
  +            q    = d[d$set == "length_for_age_inf", "p03"]
  +          , age  = d[d$set == "length_for_age_inf", "age"]
  +          , male = d[d$set == "length_for_age_inf", "male"]
  +          )
  > 
  > test_p05 <-
  +   p_length_for_age_inf(
  +            q    = d[d$set == "length_for_age_inf", "p05"]
  +          , age  = d[d$set == "length_for_age_inf", "age"]
  +          , male = d[d$set == "length_for_age_inf", "male"]
  +          )
  > test_p10 <-
  +   p_length_for_age_inf(
  +            q    = d[d$set == "length_for_age_inf", "p10"]
  +          , age  = d[d$set == "length_for_age_inf", "age"]
  +          , male = d[d$set == "length_for_age_inf", "male"]
  +          )
  > 
  > test_p25 <-
  +   p_length_for_age_inf(
  +          q    = d[d$set == "length_for_age_inf", "p25"]
  +        , age  = d[d$set == "length_for_age_inf", "age"]
  +        , male = d[d$set == "length_for_age_inf", "male"]
  +        )
  > 
  > test_p50 <-
  +   p_length_for_age_inf(
  +            q    = d[d$set == "length_for_age_inf", "p50"]
  +          , age  = d[d$set == "length_for_age_inf", "age"]
  +          , male = d[d$set == "length_for_age_inf", "male"]
  +          )
  > 
  > test_p75 <-
  +   p_length_for_age_inf(
  +            q    = d[d$set == "length_for_age_inf", "p75"]
  +          , age  = d[d$set == "length_for_age_inf", "age"]
  +          , male = d[d$set == "length_for_age_inf", "male"]
  +          )
  > 
  > test_p85 <-
  +   p_length_for_age_inf(
  +            q    = d[d$set == "length_for_age_inf", "p85"]
  +          , age  = d[d$set == "length_for_age_inf", "age"]
  +          , male = d[d$set == "length_for_age_inf", "male"]
  +          )
  > 
  > test_p90 <-
  +   p_length_for_age_inf(
  +            q    = d[d$set == "length_for_age_inf", "p90"]
  +          , age  = d[d$set == "length_for_age_inf", "age"]
  +          , male = d[d$set == "length_for_age_inf", "male"]
  +          )
  > 
  > test_p95 <-
  +   p_length_for_age_inf(
  +            q    = d[d$set == "length_for_age_inf", "p95"]
  +          , age  = d[d$set == "length_for_age_inf", "age"]
  +          , male = d[d$set == "length_for_age_inf", "male"]
  +          )
  > 
  > test_p97 <-
  +   p_length_for_age_inf(
  +            q    = d[d$set == "length_for_age_inf", "p97"]
  +          , age  = d[d$set == "length_for_age_inf", "age"]
  +          , male = d[d$set == "length_for_age_inf", "male"]
  +          )
  > 
  > 
  > stopifnot(all(round(test_p03, 2) == 0.03))
  > stopifnot(all(round(test_p05, 2) == 0.05))
  > stopifnot(all(round(test_p10, 2) == 0.10))
  > stopifnot(all(round(test_p25, 2) == 0.25))
  > stopifnot(all(round(test_p50, 2) == 0.50))
  > stopifnot(all(round(test_p75, 2) == 0.75))
  > # stopifnot(all(round(test_p85, 2) == 0.85))
  > stopifnot(all(is.na(test_p85)))
  > stopifnot(all(round(test_p90, 2) == 0.90))
  > stopifnot(all(round(test_p95, 2) == 0.95))
  > stopifnot(all(round(test_p97, 2) == 0.97))
  > 
  > # }}}
  > 
  > # Tesing q_length_for_age_inf {{{
  > test_q03 <-
  +   q_length_for_age_inf(
  +            p    = c(0.03)
  +          , age  = d[d$set == "length_for_age_inf", "age"]
  +          , male = d[d$set == "length_for_age_inf", "male"]
  +          )
  > 
  > test_q05 <-
  +   q_length_for_age_inf(
  +            p    = c(0.05)
  +          , age  = d[d$set == "length_for_age_inf", "age"]
  +          , male = d[d$set == "length_for_age_inf", "male"]
  +          )
  > 
  > test_q10 <-
  +   q_length_for_age_inf(
  +            p    = c(0.10)
  +          , age  = d[d$set == "length_for_age_inf", "age"]
  +          , male = d[d$set == "length_for_age_inf", "male"]
  +          )
  > 
  > test_q25 <-
  +   q_length_for_age_inf(
  +            p    = c(0.25)
  +          , age  = d[d$set == "length_for_age_inf", "age"]
  +          , male = d[d$set == "length_for_age_inf", "male"]
  +          )
  > 
  > test_q50 <-
  +   q_length_for_age_inf(
  +            p    = c(0.50)
  +          , age  = d[d$set == "length_for_age_inf", "age"]
  +          , male = d[d$set == "length_for_age_inf", "male"]
  +          )
  > 
  > test_q75 <-
  +   q_length_for_age_inf(
  +            p    = c(0.75)
  +          , age  = d[d$set == "length_for_age_inf", "age"]
  +          , male = d[d$set == "length_for_age_inf", "male"]
  +          )
  > 
  > # no need to test 85th as it is not in the provided data set
  > stopifnot(all(is.na(d[d$set == "length_for_age_inf", "p85"])))
  > # test_q85 <-
  > #   q_length_for_age_inf(
  > #            p    = c(0.85)
  > #          , age  = d[d$set == "length_for_age_inf", "age"]
  > #          , male = d[d$set == "length_for_age_inf", "male"]
  > #          )
  > 
  > test_q90 <-
  +   q_length_for_age_inf(
  +            p    = c(0.90)
  +          , age  = d[d$set == "length_for_age_inf", "age"]
  +          , male = d[d$set == "length_for_age_inf", "male"]
  +          )
  > 
  > test_q95 <-
  +   q_length_for_age_inf(
  +            p    = c(0.95)
  +          , age  = d[d$set == "length_for_age_inf", "age"]
  +          , male = d[d$set == "length_for_age_inf", "male"]
  +          )
  > 
  > test_q97 <-
  +   q_length_for_age_inf(
  +            p    = c(0.97)
  +          , age  = d[d$set == "length_for_age_inf", "age"]
  +          , male = d[d$set == "length_for_age_inf", "male"]
  +          )
  > 
  > stopifnot(all.equal(test_q03, d[d$set == "length_for_age_inf", "p03"], tol = 1e-6))
  > stopifnot(all.equal(test_q05, d[d$set == "length_for_age_inf", "p05"], tol = 1e-6))
  > stopifnot(all.equal(test_q10, d[d$set == "length_for_age_inf", "p10"], tol = 1e-6))
  > stopifnot(all.equal(test_q25, d[d$set == "length_for_age_inf", "p25"], tol = 1e-6))
  > stopifnot(all.equal(test_q50, d[d$set == "length_for_age_inf", "p50"], tol = 1e-6))
  > stopifnot(all.equal(test_q75, d[d$set == "length_for_age_inf", "p75"], tol = 1e-6))
  > # stopifnot(all.equal(test_q85, d[d$set == "length_for_age_inf", "p85"], tol = 1e-6))
  > stopifnot(all.equal(test_q90, d[d$set == "length_for_age_inf", "p90"], tol = 1e-6))
  > stopifnot(all.equal(test_q95, d[d$set == "length_for_age_inf", "p95"], tol = 1e-6))
  > stopifnot(all.equal(test_q97, d[d$set == "length_for_age_inf", "p97"], tol = 1e-6))
  > 
  > # end testing q_length_for_age_inf }}}
  > 
  > # Testing z_length_for_age_inf {{{
  > test_z03 <-
  +   z_length_for_age_inf(
  +            q    = d[d$set == "length_for_age_inf", "p03"]
  +          , age  = d[d$set == "length_for_age_inf", "age"]
  +          , male = d[d$set == "length_for_age_inf", "male"]
  +          )
  > 
  > test_z05 <-
  +   z_length_for_age_inf(
  +            q    = d[d$set == "length_for_age_inf", "p05"]
  +          , age  = d[d$set == "length_for_age_inf", "age"]
  +          , male = d[d$set == "length_for_age_inf", "male"]
  +          )
  > 
  > test_z10 <-
  +   z_length_for_age_inf(
  +            q    = d[d$set == "length_for_age_inf", "p10"]
  +          , age  = d[d$set == "length_for_age_inf", "age"]
  +          , male = d[d$set == "length_for_age_inf", "male"]
  +          )
  > 
  > test_z25 <-
  +   z_length_for_age_inf(
  +            q    = d[d$set == "length_for_age_inf", "p25"]
  +          , age  = d[d$set == "length_for_age_inf", "age"]
  +          , male = d[d$set == "length_for_age_inf", "male"]
  +          )
  > 
  > test_z50 <-
  +   z_length_for_age_inf(
  +            q    = d[d$set == "length_for_age_inf", "p50"]
  +          , age  = d[d$set == "length_for_age_inf", "age"]
  +          , male = d[d$set == "length_for_age_inf", "male"]
  +          )
  > 
  > test_z75 <-
  +   z_length_for_age_inf(
  +            q    = d[d$set == "length_for_age_inf", "p75"]
  +          , age  = d[d$set == "length_for_age_inf", "age"]
  +          , male = d[d$set == "length_for_age_inf", "male"]
  +          )
  > 
  > stopifnot(all(is.na(d[d$set == "length_for_age_inf", "p85"])))
  > # test_z85 <-
  > #   z_length_for_age_inf(
  > #            q    = d[d$set == "length_for_age_inf", "p85"]
  > #          , age  = d[d$set == "length_for_age_inf", "age"]
  > #          , male = d[d$set == "length_for_age_inf", "male"]
  > #          )
  > 
  > test_z90 <-
  +   z_length_for_age_inf(
  +            q    = d[d$set == "length_for_age_inf", "p90"]
  +          , age  = d[d$set == "length_for_age_inf", "age"]
  +          , male = d[d$set == "length_for_age_inf", "male"]
  +          )
  > 
  > test_z95 <-
  +   z_length_for_age_inf(
  +            q    = d[d$set == "length_for_age_inf", "p95"]
  +          , age  = d[d$set == "length_for_age_inf", "age"]
  +          , male = d[d$set == "length_for_age_inf", "male"]
  +          )
  > 
  > test_z97 <-
  +   z_length_for_age_inf(
  +            q    = d[d$set == "length_for_age_inf", "p97"]
  +          , age  = d[d$set == "length_for_age_inf", "age"]
  +          , male = d[d$set == "length_for_age_inf", "male"]
  +          )
  > 
  > 
  > stopifnot(all.equal(test_z03, rep(qnorm(0.03), length(test_z03)), tol = 2.2e-6))
  > stopifnot(all.equal(test_z05, rep(qnorm(0.05), length(test_z05)), tol = 1e-6))
  > stopifnot(all.equal(test_z10, rep(qnorm(0.10), length(test_z10)), tol = 1e-6))
  > stopifnot(all.equal(test_z25, rep(qnorm(0.25), length(test_z25)), tol = 1e-6))
  > stopifnot(all.equal(test_z50, rep(qnorm(0.50), length(test_z50)), tol = 1e-6))
  > stopifnot(all.equal(test_z75, rep(qnorm(0.75), length(test_z75)), tol = 1e-6))
  > # stopifnot(all.equal(test_z85, rep(qnorm(0.85), length(test_z85)), tol = 1e-6))
  > stopifnot(all.equal(test_z90, rep(qnorm(0.90), length(test_z90)), tol = 1e-6))
  > stopifnot(all.equal(test_z95, rep(qnorm(0.95), length(test_z95)), tol = 1e-6))
  > stopifnot(all.equal(test_z97, rep(qnorm(0.97), length(test_z97)), tol = 2.2e-6))
  > 
  > # }}}
  > 
  > 
  > proc.time()
     user  system elapsed 
     5.62    0.12    5.70 
Running the tests in 'tests/test-stature_for_age.R' failed.
Complete output:
  > library(pedbp)
  > 
  > d <- pedbp:::cdc_lms_data
  > 
  > 
  > # Testing p_stature_for_age {{{
  > test_p03 <-
  +   p_stature_for_age(
  +            q    = d[d$set == "stature_for_age", "p03"]
  +          , age  = d[d$set == "stature_for_age", "age"]
  +          , male = d[d$set == "stature_for_age", "male"]
  +          )
  > 
  > test_p05 <-
  +   p_stature_for_age(
  +            q    = d[d$set == "stature_for_age", "p05"]
  +          , age  = d[d$set == "stature_for_age", "age"]
  +          , male = d[d$set == "stature_for_age", "male"]
  +          )
  > test_p10 <-
  +   p_stature_for_age(
  +            q    = d[d$set == "stature_for_age", "p10"]
  +          , age  = d[d$set == "stature_for_age", "age"]
  +          , male = d[d$set == "stature_for_age", "male"]
  +          )
  > 
  > test_p25 <-
  +   p_stature_for_age(
  +            q    = d[d$set == "stature_for_age", "p25"]
  +          , age  = d[d$set == "stature_for_age", "age"]
  +          , male = d[d$set == "stature_for_age", "male"]
  +          )
  > 
  > test_p50 <-
  +   p_stature_for_age(
  +            q    = d[d$set == "stature_for_age", "p50"]
  +          , age  = d[d$set == "stature_for_age", "age"]
  +          , male = d[d$set == "stature_for_age", "male"]
  +          )
  > 
  > test_p75 <-
  +   p_stature_for_age(
  +            q    = d[d$set == "stature_for_age", "p75"]
  +          , age  = d[d$set == "stature_for_age", "age"]
  +          , male = d[d$set == "stature_for_age", "male"]
  +          )
  > 
  > stopifnot(all(is.na( d[d$set == "stature_for_age", "p85"])))
  > # test_p85 <-
  > #   p_stature_for_age(
  > #            q    = d[d$set == "stature_for_age", "p85"]
  > #          , age  = d[d$set == "stature_for_age", "age"]
  > #          , male = d[d$set == "stature_for_age", "male"]
  > #          )
  > 
  > test_p90 <-
  +   p_stature_for_age(
  +            q    = d[d$set == "stature_for_age", "p90"]
  +          , age  = d[d$set == "stature_for_age", "age"]
  +          , male = d[d$set == "stature_for_age", "male"]
  +          )
  > 
  > test_p95 <-
  +   p_stature_for_age(
  +            q    = d[d$set == "stature_for_age", "p95"]
  +          , age  = d[d$set == "stature_for_age", "age"]
  +          , male = d[d$set == "stature_for_age", "male"]
  +          )
  > 
  > test_p97 <-
  +   p_stature_for_age(
  +            q    = d[d$set == "stature_for_age", "p97"]
  +          , age  = d[d$set == "stature_for_age", "age"]
  +          , male = d[d$set == "stature_for_age", "male"]
  +          )
  > 
  > 
  > stopifnot(all(round(test_p03, 2) == 0.03))
  > stopifnot(all(round(test_p05, 2) == 0.05))
  > stopifnot(all(round(test_p10, 2) == 0.10))
  > stopifnot(all(round(test_p25, 2) == 0.25))
  > stopifnot(all(round(test_p50, 2) == 0.50))
  > stopifnot(all(round(test_p75, 2) == 0.75))
  > # stopifnot(all(round(test_p85, 2) == 0.85))
  > stopifnot(all(round(test_p90, 2) == 0.90))
  > stopifnot(all(round(test_p95, 2) == 0.95))
  > stopifnot(all(round(test_p97, 2) == 0.97))
  > 
  > # }}}
  > 
  > # Tesing q_stature_for_age {{{
  > test_q03 <-
  +   q_stature_for_age(
  +            p    = c(0.03)
  +          , age  = d[d$set == "stature_for_age", "age"]
  +          , male = d[d$set == "stature_for_age", "male"]
  +          )
  > 
  > test_q05 <-
  +   q_stature_for_age(
  +            p    = c(0.05)
  +          , age  = d[d$set == "stature_for_age", "age"]
  +          , male = d[d$set == "stature_for_age", "male"]
  +          )
  > 
  > test_q10 <-
  +   q_stature_for_age(
  +            p    = c(0.10)
  +          , age  = d[d$set == "stature_for_age", "age"]
  +          , male = d[d$set == "stature_for_age", "male"]
  +          )
  > 
  > test_q25 <-
  +   q_stature_for_age(
  +            p    = c(0.25)
  +          , age  = d[d$set == "stature_for_age", "age"]
  +          , male = d[d$set == "stature_for_age", "male"]
  +          )
  > 
  > test_q50 <-
  +   q_stature_for_age(
  +            p    = c(0.50)
  +          , age  = d[d$set == "stature_for_age", "age"]
  +          , male = d[d$set == "stature_for_age", "male"]
  +          )
  > 
  > test_q75 <-
  +   q_stature_for_age(
  +            p    = c(0.75)
  +          , age  = d[d$set == "stature_for_age", "age"]
  +          , male = d[d$set == "stature_for_age", "male"]
  +          )
  > 
  > # not need to test 85th
  > stopifnot(all(is.na( d[d$set == "stature_for_age", "p85"])))
  > # test_q85 <-
  > #   q_stature_for_age(
  > #            p    = c(0.85)
  > #          , age  = d[d$set == "stature_for_age", "age"]
  > #          , male = d[d$set == "stature_for_age", "male"]
  > #          )
  > 
  > test_q90 <-
  +   q_stature_for_age(
  +            p    = c(0.90)
  +          , age  = d[d$set == "stature_for_age", "age"]
  +          , male = d[d$set == "stature_for_age", "male"]
  +          )
  > 
  > test_q95 <-
  +   q_stature_for_age(
  +            p    = c(0.95)
  +          , age  = d[d$set == "stature_for_age", "age"]
  +          , male = d[d$set == "stature_for_age", "male"]
  +          )
  > 
  > test_q97 <-
  +   q_stature_for_age(
  +            p    = c(0.97)
  +          , age  = d[d$set == "stature_for_age", "age"]
  +          , male = d[d$set == "stature_for_age", "male"]
  +          )
  > 
  > stopifnot(all.equal(test_q03, d[d$set == "stature_for_age", "p03"]))
  > stopifnot(all.equal(test_q05, d[d$set == "stature_for_age", "p05"]))
  > stopifnot(all.equal(test_q10, d[d$set == "stature_for_age", "p10"]))
  > stopifnot(all.equal(test_q25, d[d$set == "stature_for_age", "p25"]))
  > stopifnot(all.equal(test_q50, d[d$set == "stature_for_age", "p50"]))
  > stopifnot(all.equal(test_q75, d[d$set == "stature_for_age", "p75"]))
  > # stopifnot(all.equal(test_q85, d[d$set == "stature_for_age", "p85"]))
  > stopifnot(all.equal(test_q90, d[d$set == "stature_for_age", "p90"]))
  > stopifnot(all.equal(test_q95, d[d$set == "stature_for_age", "p95"]))
  > stopifnot(all.equal(test_q97, d[d$set == "stature_for_age", "p97"]))
  > 
  > # end testing q_stature_for_age }}}
  > 
  > # Testing z_stature_for_age {{{
  > test_z03 <-
  +   z_stature_for_age(
  +            q    = d[d$set == "stature_for_age", "p03"]
  +          , age  = d[d$set == "stature_for_age", "age"]
  +          , male = d[d$set == "stature_for_age", "male"]
  +          )
  > 
  > test_z05 <-
  +   z_stature_for_age(
  +            q    = d[d$set == "stature_for_age", "p05"]
  +          , age  = d[d$set == "stature_for_age", "age"]
  +          , male = d[d$set == "stature_for_age", "male"]
  +          )
  > 
  > test_z10 <-
  +   z_stature_for_age(
  +            q    = d[d$set == "stature_for_age", "p10"]
  +          , age  = d[d$set == "stature_for_age", "age"]
  +          , male = d[d$set == "stature_for_age", "male"]
  +          )
  > 
  > test_z25 <-
  +   z_stature_for_age(
  +            q    = d[d$set == "stature_for_age", "p25"]
  +          , age  = d[d$set == "stature_for_age", "age"]
  +          , male = d[d$set == "stature_for_age", "male"]
  +          )
  > 
  > test_z50 <-
  +   z_stature_for_age(
  +            q    = d[d$set == "stature_for_age", "p50"]
  +          , age  = d[d$set == "stature_for_age", "age"]
  +          , male = d[d$set == "stature_for_age", "male"]
  +          )
  > 
  > test_z75 <-
  +   z_stature_for_age(
  +            q    = d[d$set == "stature_for_age", "p75"]
  +          , age  = d[d$set == "stature_for_age", "age"]
  +          , male = d[d$set == "stature_for_age", "male"]
  +          )
  > 
  > stopifnot(all(is.na(d[d$set == "stature_for_age", "p85"])))
  > # test_z85 <-
  > #   z_stature_for_age(
  > #            q    = d[d$set == "stature_for_age", "p85"]
  > #          , age  = d[d$set == "stature_for_age", "age"]
  > #          , male = d[d$set == "stature_for_age", "male"]
  > #          )
  > 
  > test_z90 <-
  +   z_stature_for_age(
  +            q    = d[d$set == "stature_for_age", "p90"]
  +          , age  = d[d$set == "stature_for_age", "age"]
  +          , male = d[d$set == "stature_for_age", "male"]
  +          )
  > 
  > test_z95 <-
  +   z_stature_for_age(
  +            q    = d[d$set == "stature_for_age", "p95"]
  +          , age  = d[d$set == "stature_for_age", "age"]
  +          , male = d[d$set == "stature_for_age", "male"]
  +          )
  > 
  > test_z97 <-
  +   z_stature_for_age(
  +            q    = d[d$set == "stature_for_age", "p97"]
  +          , age  = d[d$set == "stature_for_age", "age"]
  +          , male = d[d$set == "stature_for_age", "male"]
  +          )
  > 
  > 
  > stopifnot(all.equal(test_z03, rep(qnorm(0.03), length(test_z03))))
  > stopifnot(all.equal(test_z05, rep(qnorm(0.05), length(test_z05))))
  > stopifnot(all.equal(test_z10, rep(qnorm(0.10), length(test_z10))))
  > stopifnot(all.equal(test_z25, rep(qnorm(0.25), length(test_z25))))
  > stopifnot(all.equal(test_z50, rep(qnorm(0.50), length(test_z50))))
  > stopifnot(all.equal(test_z75, rep(qnorm(0.75), length(test_z75))))
  > # stopifnot(all.equal(test_z85, rep(qnorm(0.85), length(test_z85))))
  > stopifnot(all.equal(test_z90, rep(qnorm(0.90), length(test_z90))))
  > stopifnot(all.equal(test_z95, rep(qnorm(0.95), length(test_z95))))
  > stopifnot(all.equal(test_z97, rep(qnorm(0.97), length(test_z97))))
  > 
  > # }}}
  > 
  > 
  > proc.time()
     user  system elapsed 
    31.25    2.48   33.76 
Running the tests in 'tests/test-weight_for_length_inf.R' failed.
Complete output:
  > library(pedbp)
  > 
  > d <- pedbp:::cdc_lms_data
  > 
  > # Testing p_weight_for_length_inf {{{
  > test_p03 <-
  +   p_weight_for_length_inf(
  +            q    = d[d$set == "weight_for_length_inf", "p03"]
  +          , length  = d[d$set == "weight_for_length_inf", "length"]
  +          , male = d[d$set == "weight_for_length_inf", "male"]
  +          )
  > 
  > test_p05 <-
  +   p_weight_for_length_inf(
  +            q    = d[d$set == "weight_for_length_inf", "p05"]
  +          , length  = d[d$set == "weight_for_length_inf", "length"]
  +          , male = d[d$set == "weight_for_length_inf", "male"]
  +          )
  > test_p10 <-
  +   p_weight_for_length_inf(
  +            q    = d[d$set == "weight_for_length_inf", "p10"]
  +          , length  = d[d$set == "weight_for_length_inf", "length"]
  +          , male = d[d$set == "weight_for_length_inf", "male"]
  +          )
  > 
  > test_p25 <-
  +   p_weight_for_length_inf(
  +            q    = d[d$set == "weight_for_length_inf", "p25"]
  +          , length  = d[d$set == "weight_for_length_inf", "length"]
  +          , male = d[d$set == "weight_for_length_inf", "male"]
  +          )
  > 
  > test_p50 <-
  +   p_weight_for_length_inf(
  +            q    = d[d$set == "weight_for_length_inf", "p50"]
  +          , length  = d[d$set == "weight_for_length_inf", "length"]
  +          , male = d[d$set == "weight_for_length_inf", "male"]
  +          )
  > 
  > test_p75 <-
  +   p_weight_for_length_inf(
  +            q    = d[d$set == "weight_for_length_inf", "p75"]
  +          , length  = d[d$set == "weight_for_length_inf", "length"]
  +          , male = d[d$set == "weight_for_length_inf", "male"]
  +          )
  > 
  > stopifnot(all(is.na(d[d$set == "weight_for_length_inf", "p85"])))
  > # test_p85 <-
  > #   p_weight_for_length_inf(
  > #            q    = d[d$set == "weight_for_length_inf", "p85"]
  > #          , length  = d[d$set == "weight_for_length_inf", "length"]
  > #          , male = d[d$set == "weight_for_length_inf", "male"]
  > #          )
  > 
  > test_p90 <-
  +   p_weight_for_length_inf(
  +            q    = d[d$set == "weight_for_length_inf", "p90"]
  +          , length  = d[d$set == "weight_for_length_inf", "length"]
  +          , male = d[d$set == "weight_for_length_inf", "male"]
  +          )
  > 
  > test_p95 <-
  +   p_weight_for_length_inf(
  +            q    = d[d$set == "weight_for_length_inf", "p95"]
  +          , length  = d[d$set == "weight_for_length_inf", "length"]
  +          , male = d[d$set == "weight_for_length_inf", "male"]
  +          )
  > 
  > test_p97 <-
  +   p_weight_for_length_inf(
  +            q    = d[d$set == "weight_for_length_inf", "p97"]
  +          , length  = d[d$set == "weight_for_length_inf", "length"]
  +          , male = d[d$set == "weight_for_length_inf", "male"]
  +          )
  > 
  > 
  > stopifnot(all(round(test_p03, 2) == 0.03))
  > stopifnot(all(round(test_p05, 2) == 0.05))
  > stopifnot(all(round(test_p10, 2) == 0.10))
  > stopifnot(all(round(test_p25, 2) == 0.25))
  > stopifnot(all(round(test_p50, 2) == 0.50))
  > stopifnot(all(round(test_p75, 2) == 0.75))
  > # stopifnot(all(round(test_p85, 2) == 0.85))
  > stopifnot(all(round(test_p90, 2) == 0.90))
  > stopifnot(all(round(test_p95, 2) == 0.95))
  > stopifnot(all(round(test_p97, 2) == 0.97))
  > 
  > # }}}
  > 
  > # Tesing q_weight_for_length_inf {{{
  > test_q03 <-
  +   q_weight_for_length_inf(
  +            p    = c(0.03)
  +          , length  = d[d$set == "weight_for_length_inf", "length"]
  +          , male = d[d$set == "weight_for_length_inf", "male"]
  +          )
  > 
  > test_q05 <-
  +   q_weight_for_length_inf(
  +            p    = c(0.05)
  +          , length  = d[d$set == "weight_for_length_inf", "length"]
  +          , male = d[d$set == "weight_for_length_inf", "male"]
  +          )
  > 
  > test_q10 <-
  +   q_weight_for_length_inf(
  +            p    = c(0.10)
  +          , length  = d[d$set == "weight_for_length_inf", "length"]
  +          , male = d[d$set == "weight_for_length_inf", "male"]
  +          )
  > 
  > test_q25 <-
  +   q_weight_for_length_inf(
  +            p    = c(0.25)
  +          , length  = d[d$set == "weight_for_length_inf", "length"]
  +          , male = d[d$set == "weight_for_length_inf", "male"]
  +          )
  > 
  > test_q50 <-
  +   q_weight_for_length_inf(
  +            p    = c(0.50)
  +          , length  = d[d$set == "weight_for_length_inf", "length"]
  +          , male = d[d$set == "weight_for_length_inf", "male"]
  +          )
  > 
  > test_q75 <-
  +   q_weight_for_length_inf(
  +            p    = c(0.75)
  +          , length  = d[d$set == "weight_for_length_inf", "length"]
  +          , male = d[d$set == "weight_for_length_inf", "male"]
  +          )
  > 
  > stopifnot(all(is.na(d[d$set == "weight_for_length_inf", "p85"])))
  > # test_q85 <-
  > #   q_weight_for_length_inf(
  > #            p    = c(0.85)
  > #          , length  = d[d$set == "weight_for_length_inf", "length"]
  > #          , male = d[d$set == "weight_for_length_inf", "male"]
  > #          )
  > 
  > test_q90 <-
  +   q_weight_for_length_inf(
  +            p    = c(0.90)
  +          , length  = d[d$set == "weight_for_length_inf", "length"]
  +          , male = d[d$set == "weight_for_length_inf", "male"]
  +          )
  > 
  > test_q95 <-
  +   q_weight_for_length_inf(
  +            p    = c(0.95)
  +          , length  = d[d$set == "weight_for_length_inf", "length"]
  +          , male = d[d$set == "weight_for_length_inf", "male"]
  +          )
  > 
  > test_q97 <-
  +   q_weight_for_length_inf(
  +            p    = c(0.97)
  +          , length  = d[d$set == "weight_for_length_inf", "length"]
  +          , male = d[d$set == "weight_for_length_inf", "male"]
  +          )
  > 
  > stopifnot(all.equal(test_q03, d[d$set == "weight_for_length_inf", "p03"]))
  > stopifnot(all.equal(test_q05, d[d$set == "weight_for_length_inf", "p05"]))
  > stopifnot(all.equal(test_q10, d[d$set == "weight_for_length_inf", "p10"]))
  > stopifnot(all.equal(test_q25, d[d$set == "weight_for_length_inf", "p25"]))
  > stopifnot(all.equal(test_q50, d[d$set == "weight_for_length_inf", "p50"]))
  > stopifnot(all.equal(test_q75, d[d$set == "weight_for_length_inf", "p75"]))
  > # stopifnot(all.equal(test_q85, d[d$set == "weight_for_length_inf", "p85"]))
  > stopifnot(all.equal(test_q90, d[d$set == "weight_for_length_inf", "p90"]))
  > stopifnot(all.equal(test_q95, d[d$set == "weight_for_length_inf", "p95"]))
  > stopifnot(all.equal(test_q97, d[d$set == "weight_for_length_inf", "p97"]))
  > 
  > # end testing q_weight_for_length_inf }}}
  > 
  > # Testing z_weight_for_length_inf {{{
  > test_z03 <-
  +   z_weight_for_length_inf(
  +            q    = d[d$set == "weight_for_length_inf", "p03"]
  +          , length  = d[d$set == "weight_for_length_inf", "length"]
  +          , male = d[d$set == "weight_for_length_inf", "male"]
  +          )
  > 
  > test_z05 <-
  +   z_weight_for_length_inf(
  +            q    = d[d$set == "weight_for_length_inf", "p05"]
  +          , length  = d[d$set == "weight_for_length_inf", "length"]
  +          , male = d[d$set == "weight_for_length_inf", "male"]
  +          )
  > 
  > test_z10 <-
  +   z_weight_for_length_inf(
  +            q    = d[d$set == "weight_for_length_inf", "p10"]
  +          , length  = d[d$set == "weight_for_length_inf", "length"]
  +          , male = d[d$set == "weight_for_length_inf", "male"]
  +          )
  > 
  > test_z25 <-
  +   z_weight_for_length_inf(
  +            q    = d[d$set == "weight_for_length_inf", "p25"]
  +          , length  = d[d$set == "weight_for_length_inf", "length"]
  +          , male = d[d$set == "weight_for_length_inf", "male"]
  +          )
  > 
  > test_z50 <-
  +   z_weight_for_length_inf(
  +            q    = d[d$set == "weight_for_length_inf", "p50"]
  +          , length  = d[d$set == "weight_for_length_inf", "length"]
  +          , male = d[d$set == "weight_for_length_inf", "male"]
  +          )
  > 
  > test_z75 <-
  +   z_weight_for_length_inf(
  +            q    = d[d$set == "weight_for_length_inf", "p75"]
  +          , length  = d[d$set == "weight_for_length_inf", "length"]
  +          , male = d[d$set == "weight_for_length_inf", "male"]
  +          )
  > 
  > stopifnot(all(is.na(d[d$set == "weight_for_length_inf", "p85"])))
  > # test_z85 <-
  > #   z_weight_for_length_inf(
  > #            q    = d[d$set == "weight_for_length_inf", "p85"]
  > #          , length  = d[d$set == "weight_for_length_inf", "length"]
  > #          , male = d[d$set == "weight_for_length_inf", "male"]
  > #          )
  > 
  > test_z90 <-
  +   z_weight_for_length_inf(
  +            q    = d[d$set == "weight_for_length_inf", "p90"]
  +          , length  = d[d$set == "weight_for_length_inf", "length"]
  +          , male = d[d$set == "weight_for_length_inf", "male"]
  +          )
  > 
  > test_z95 <-
  +   z_weight_for_length_inf(
  +            q    = d[d$set == "weight_for_length_inf", "p95"]
  +          , length  = d[d$set == "weight_for_length_inf", "length"]
  +          , male = d[d$set == "weight_for_length_inf", "male"]
  +          )
  > 
  > test_z97 <-
  +   z_weight_for_length_inf(
  +            q    = d[d$set == "weight_for_length_inf", "p97"]
  +          , length  = d[d$set == "weight_for_length_inf", "length"]
  +          , male = d[d$set == "weight_for_length_inf", "male"]
  +          )
  > 
  > 
  > stopifnot(all.equal(test_z03, rep(qnorm(0.03), length(test_z03))))
  > stopifnot(all.equal(test_z05, rep(qnorm(0.05), length(test_z05))))
  > stopifnot(all.equal(test_z10, rep(qnorm(0.10), length(test_z10))))
  > stopifnot(all.equal(test_z25, rep(qnorm(0.25), length(test_z25))))
  > stopifnot(all.equal(test_z50, rep(qnorm(0.50), length(test_z50))))
  > stopifnot(all.equal(test_z75, rep(qnorm(0.75), length(test_z75))))
  > # stopifnot(all.equal(test_z85, rep(qnorm(0.85), length(test_z85))))
  > stopifnot(all.equal(test_z90, rep(qnorm(0.90), length(test_z90))))
  > stopifnot(all.equal(test_z95, rep(qnorm(0.95), length(test_z95))))
  > stopifnot(all.equal(test_z97, rep(qnorm(0.97), length(test_z97))))
  > 
  > # }}}
  > 
  > 
  > proc.time()
     user  system elapsed 
     8.29    0.31    8.59