From b6d5d31974dd3d2c96e392fed1ac37c8df51ac98 Mon Sep 17 00:00:00 2001
From: Stef van Buuren The columns The script below creates a figure with -2SD, 0SD and +2SD centiles
@@ -996,7 +996,7 @@ Van Buuren S, Eekhout I, Huizing A (2024). D-score for Child Development. The dscore R package, version 1.9.0. van Buuren S, Eekhout I (2023).
+Child Development with the D-score.
+CRC Press, Boca Raton, FL.
+doi:10.1201/9781003216315.
+ The D-score is a numerical score that measures generic development in children. You may use the D-score to analyze and predict development of children similar to measures like height and weight. The D-score is a numerical score that measures generic development in children. Use the D-score to analyze and predict early development of children similar to measures like height and weight. The You may estimate the D-score and the D-score age-adjusted Z-score (DAZ) from child data on developmental milestones. Four steps are needed: You may estimate the D-score and the Development-for-Age Z-score (DAZ) from child data on developmental milestones. Four steps are needed: This study was supported by the Bill & Melinda Gates Foundation. The contents are the sole responsibility of the authors and may not necessarily represent the official views of the Bill & Melinda Gates Foundation or other agencies that may have supported the primary data studies used in the present study. The authors wish to recognize the principal investigators and their study team members for their generous contribution of the data that made this tool possible and the members of the Ki team who directly or indirectly contributed to the study: Amina Abubakar, Claudia R. Lindgren Alves, Orazio Attanasio, Maureen M. Black, Maria Caridad Araujo, Susan M. Chang-Lopez, Gary L. Darmstadt, Bernice M. Doove, Wafaie Fawzi, Lia C.H. Fernald, Günther Fink, Emanuela Galasso, Melissa Gladstone, Sally M. Grantham-McGregor, Cristina Gutierrez de Pineres, Pamela Jervis, Jena Derakhshani Hamadani, Charlotte Hanlon, Simone M. Karam, Gillian Lancaster, Betzy Lozoff, Gareth McCray, Jeffrey R Measelle, Girmay Medhin, Ana M. B. Menezes, Lauren Pisani, Helen Pitchik, Muneera Rasheed, Lisy Ratsifandrihamanana, Sarah Reynolds, Linda Richter, Marta Rubio-Codina, Norbert Schady, Limbika Sengani, Chris Sudfeld, Marcus Waldman, Susan P. Walker, Ann M. Weber and Aisha K. Yousafzai. This study was supported by the Bill & Melinda Gates Foundation. The contents are the sole responsibility of the authors and may not necessarily represent the official views of the Bill & Melinda Gates Foundation or other agencies that may have supported the primary data studies used in the present study. Updates the Use Simplifies the spelling of the term “D-score” to improve consistency and readability Removes dependency on Replaces Make style more consistent with This is a major update of the BREAKING CHANGE: On May 31, 2024 we detected a long-time error in the calculation of the D-score resulting from an incorrect scale factor that led us to believe that item characteristic curves are steeper than the actually are. The impact of the error on the result is as follows: 1) There is no effect on the difficulty estimates of the Rasch models, 2) The D-score estimates are slightly altered but changes are small, 3) The references are largely unaffected and need not to be redone, 4) The estimates of the SEMs can differ substantially, so inferences based on the SEMs should be re-evaluated, 5) When there were changes in the analyses, the results in the newer method look smoother and are preferred. The correction appeared in the development version BREAKING CHANGE: NEW DEFAULT KEY: Adds a new reference 20240609 SvB: Added A A The Here are more details on the reference population:
+The Returns the age-interpolated median of the D-score of the default
+reference for a given key. Decimal age, numeric vector Character, key of the reference population A vector of length Do not use this function if you want the median D-score for a specific
+reference. A Summarise D-score and DAZlibrary(dplyr)
r <- builtin_references |>
- filter(name == "dutch") |>
+ filter(population == "dutch") |>
select(age, SDM2, SD0, SDP2)
ggplot(md, aes(x = a, y = d, group = id, color = sex)) +
diff --git a/docs/articles/index.html b/docs/articles/index.html
index 397d05b..ea94690 100644
--- a/docs/articles/index.html
+++ b/docs/articles/index.html
@@ -17,7 +17,7 @@
diff --git a/docs/articles/scoring_GSED.html b/docs/articles/scoring_GSED.html
index cab818a..c83db22 100644
--- a/docs/articles/scoring_GSED.html
+++ b/docs/articles/scoring_GSED.html
@@ -33,7 +33,7 @@
@@ -101,10 +101,10 @@
D-score and DAZ
-
R
. The
-app contains online documentation and instructions and will not be
-further discussed here.R
. The app contains online documentation and instructions
+and will not be further discussed here.R
package dscore
. The R
package dscore
at https://CRAN.R-project.org/package=dscore is a flexible
@@ -118,8 +118,8 @@ PreliminariesGetting
-Started with R site;
R
package dscore
on your local machine;data.frame
, a
@@ -976,16 +976,16 @@ Phase 1 references and DAZ
-library(dplyr, warn.conflicts = FALSE, quietly = TRUE)
ref <- builtin_references |>
- filter(name == "phase1") |>
- select(name, age, mu, sigma, nu, tau, SDM2, SD0, SDP2)
+ filter(population == "phase1") |>
+ select(population, age, mu, sigma, nu, tau, SDM2, SD0, SDP2)
head(ref)
+## name age mu sigma nu tau SDM2 SD0 SDP2
-## 1 phase1 0.0383 13.68 0.2456 1.1731 15.422 6.042690 13.68707 20.71620
-## 2 phase1 0.0575 14.36 0.2324 1.2062 15.540 6.699401 14.36568 21.30880
-## 3 phase1 0.0767 15.02 0.2206 1.2375 15.652 7.354450 15.02457 21.88555
-## 4 phase1 0.0958 15.68 0.2100 1.2670 15.758 8.014796 15.68368 22.47481
-## 5 phase1 0.1150 16.35 0.2005 1.2951 15.860 8.680976 16.35299 23.09062
-## 6 phase1 0.1342 17.03 0.1917 1.3218 15.957 9.363085 17.03241 23.72268
## population age mu sigma nu tau SDM2 SD0 SDP2
+## 1 phase1 0.0383 13.68 0.2456 1.1731 15.422 6.042690 13.68707 20.71620
+## 2 phase1 0.0575 14.36 0.2324 1.2062 15.540 6.699401 14.36568 21.30880
+## 3 phase1 0.0767 15.02 0.2206 1.2375 15.652 7.354450 15.02457 21.88555
+## 4 phase1 0.0958 15.68 0.2100 1.2670 15.758 8.014796 15.68368 22.47481
+## 5 phase1 0.1150 16.35 0.2005 1.2951 15.860 8.680976 16.35299 23.09062
+## 6 phase1 0.1342 17.03 0.1917 1.3218 15.957 9.363085 17.03241 23.72268
mu
, sigma
, nu
and
tau
are the age-varying parameters of a Box-Cox \(t\) (BCT) distribution.Phase 1 references and DAZlibrary(patchwork)
r <- builtin_references |>
- filter(name == "phase1" & age <= 3.5) |>
+ filter(population == "phase1" & age <= 3.5) |>
mutate(m = age * 12)
lf2$ins <- "lf"
diff --git a/docs/authors.html b/docs/authors.html
index 472bac5..3f1c2af 100644
--- a/docs/authors.html
+++ b/docs/authors.html
@@ -17,7 +17,7 @@
@@ -82,13 +82,18 @@
Citation
- @Manual{dscore-package,
- title = {D-score for Child Development},
- author = {S. {van Buuren} and I. Eekhout and A. Huizing},
- year = {2024},
- note = {R package version 1.9.0},
- url = {https://github.com/d-score/dscore},
+
@Book{,
+ title = {Child Development with the D-score},
+ author = {Stef {van Buuren} and Iris Eekhout},
+ year = {2023},
+ publisher = {CRC Press},
+ address = {Boca Raton, FL},
+ doi = {https://doi.org/10.1201/9781003216315},
}
diff --git a/docs/index.html b/docs/index.html
index 6a90f82..6f7f09b 100644
--- a/docs/index.html
+++ b/docs/index.html
@@ -12,7 +12,7 @@
-
+
-dscore
package contains tools to
diff --git a/docs/news/index.html b/docs/news/index.html
index af39720..fb7f96c 100644
--- a/docs/news/index.html
+++ b/docs/news/index.html
@@ -17,7 +17,7 @@
@@ -61,60 +61,57 @@ Installation
Overview
-
dscore
package covers your measurement instrument;Keys
dutch
key, 0-2 years: van Buuren (2014)gcdg
key: Weber et al. (2019)gsed
key: GSED team (Maureen Black, Kieran Bromley, Vanessa Cavallera (lead author), Jorge Cuartas, Tarun Dua (corresponding author), Iris Eekhout, Günther Fink, Melissa Gladstone, Katelyn Hepworth, Magdalena Janus, Patricia Kariger, Gillian Lancaster, Dana McCoy, Gareth McCray, Abbie Raikes, Marta Rubio-Codina, Stef van Buuren, Marcus Waldman, Susan Walker and Ann Weber) (2019)gsed
keys: World Health Organization (WHO) (2023)Shiny app
Acknowledgement
-Literature
@@ -160,9 +160,6 @@ Literature
Fernald, L. C. H., E. Prado, P. Kariger, and A. Raikes. 2017. “A Toolkit for Measuring Early Childhood Development in Low and Middle-Income Countries.” https://documents.worldbank.org/en/publication/documents-reports/documentdetail/384681513101293811/a-toolkit-for-measuring-early-childhood-development-in-low-and-middle-income-countries.
Literature
Weber, A. M., M. Rubio-Codina, S. P. Walker, S. van Buuren, I. Eekhout, S. Grantham-McGregor, M. C. Araujo, et al. 2019. “The D-Score: A Metric for Interpreting the Early Development of Infants and Toddlers Across Global Settings.” BMJ Global Health 4: e001724. https://gh.bmj.com/content/bmjgh/4/6/e001724.full.pdf.
Dev status
Changelog
dscore 1.8.9
-phase1_healthy
age-conditional D-score references with three changes:population = "phase1_healthy"
to calculate DAZ relative to this new reference. The new reference is based on healthy children only, and is currently experimental. The current reference "phase1"
remains the default for now, but may be superseeded by "phase1_healthy"
in a future version.tidyselect
magrittr
pipe %>%
by base pipe |>
styler
dscore 1.8.8
-algorithm
argument to dscore()
and dscore_posterior()
to select the algorithm used to calculate the D-score. The default algorithm, called "current"
, correctly uses a scale factor of 1 for calculation in the logistic metric. Previously, a scale factor of 1 was errroneously applied to the D-score metric. The current algorithm results in different D-score estimates and higher SEM estimates, but is preferred to the former method. For backward-compatibility to dscore 1.8.7
and earlier, use algorithm = "1.8.7"
.dscore 1.8.7
-"phase1_healthy"
calculated from selective subsample of the GSED Phase 1 data using the “gsed2212” key. This reference is based on the same data as the "phase1"
reference, but only includes children who were developing well at the time of the assessment. This reference is intended for use in studies where the population of interest is healthy children. Note: This is a temporary name and will be deprecated.dscore 1.8.6
-"dutch_gsed2212"
calculated from Dutch data using the “gsed2212” key. Note: This is a temporary name and will be deprecated.dscore 1.8.5
-dscore 1.9.0
+Overview
+dscore
package featuring: - a new default reference "preliminary_standards"
- a correction of an issue with the scaling factor - a major clean-up of the itembank, references, keys, and R code - improved documentation and examplesMajor issues
+dscore 1.8.8
, and is now incorporated into release dscore 1.9.0
. For backward compatibility to dscore 1.8.7
and earlier, use the argument algorithm = "1.8.7"
in calls to the dscore()
function.dscore_posterior()
now returns a data.frame
with column names that indicate the quadrature point. This was an unnamed matrix
. Code that expects a matrix
as the return of dscore_posterior()
may need to be adapted."preliminary_standards"
calculated from selected subsample of the GSED Phase 1 data, and makes these the default in this release. The reference is a temporary stand-in for a future norm-based standard for normal early child development. This reference replaces the temporary reference "phase1_healthy"
that was introduced in dscore 1.8.7
. Compared to the "phase1"
reference, the "preliminary_standards"
reference has the following differences: 1) D-score estimation uses the new model 20240601 with correct scale factor, 2) Calculates the D-score for SF and LF separately (not combined), 3) Tunes the GAMSLSS model to fit the healthy subsample. The "phase1_healthy"
object is removed.Major changes
+"dutch"
calculated using the "gsed2212"
key."gsed2406"
to accomodate for the changed prior mean because of the adoption of the new reference "preliminary_standards"
as the base population. The key "gsed2406"
is identical to "gsed2212"
, and is the default key in this release.builtin_keys
table that contains proper defaults for the base reference, transformation and quadrature points per keykey
and population
. Previously the index was based on only population
. This change allows for multiple references per key, and for references created for the same population under different keys. The key
field is now mandatory in the reference table.verbose
option to dscore()
, dscore_posterior()
, get_age_equivalent()
, get_reference()
, get_tau()
, daz()
and zad()
to print progress messages to the console on the values of key
, population
, transform
, qp
and algorithm
. This is useful for debugging and for understanding the behavior of the functions.builtin_keys
table. This makes the code more readable and maintainable.sf2206
, lf2206
, 294_0
, gsed2206
, gsed2208
and removes from the builtin_itembank
.dscore()
and dscore_posterior()
can now copy variables from the input data
into the output through the prepend
argument. (#46)dscore_posterior()
now returns a data.frame
with column names that indicate the quadrature point. This was an unnamed matrix
. Code that expects a matrix
as the return of dscore_posterior()
may need to be adapted.dscore 1.8.4
-Specified C++11: please drop specification unless essential
), removes a C++11 specificationdscore 1.8.3
-get_reference()
to phase1
to remain in sync with the default key = "gsed"
+Minor changes
+init_key()
and set_default_xxx()
functions to regulate values for key
, population
, transform
and qp
argumentsdata-raw/data/keys
to more consistent names, adapts data-raw/R/save_builtin_itembank.R
to reflect model history, and rebuilds builtin_itembank
+tibble
and tidyselect
, and replaces the dependency on stringr
by the lighter stringi
+reference
to reference_table
in daz()
and zad()
to avoid confusion with the references
argument in get_reference()
magrittr
pipe %>%
by base pipe |>
+styler
+Specified C++11: please drop specification unless essential
), removes a C++11 specificationpBCT()
daz()
and zad()
using the BCT transformation for positive valuesnames
attribute from the return value of daz()
and zad()
dscore 1.8.2
-pBCT()
when is.na(nu)
is TRUE
+pBCT()
when is.na(nu)
is TRUE
bibtex
packagedata-raw/data/keys
folder and renamed them to improve readabilitydscore 1.8.02023-01-20
Available keys for calculating the D-score
Format
- builtin_keys
is a data.frame
with variables:Name Label key
String. Name of the key indicating the Rasch model default_reference
String. Name of the default reference for the key n_items
Number of items in the key n_instruments
Number of instruments in the key intercept
Intercept to convert logit into D-score slope
Slope to convert logit into D-score from
Starting value of the quadrature points to
Stopping value of the quadrature points by
Increment of the quadrature points retired
Has the key been retired? builtin_keys
is a data.frame
with variables:Name Label key
String. Name of the key indicating the Rasch model base_population
String. Name of the base population for the key n_items
Number of items in the key n_instruments
Number of instruments in the key intercept
Intercept to convert logit into D-score slope
Slope to convert logit into D-score from
Starting value of the quadrature points to
Stopping value of the quadrature points by
Increment of the quadrature points retired
Has the key been retired? Note
builtin_keys
table by
diff --git a/docs/reference/builtin_references.html b/docs/reference/builtin_references.html
index 5a96e49..a441e1f 100644
--- a/docs/reference/builtin_references.html
+++ b/docs/reference/builtin_references.html
@@ -23,7 +23,7 @@
Collection of age-conditional reference distributions
Format
- data.frame
with the following variables:Name Label name
Name of the reference key
D-score key: "dutch"
, "gcdg"
or "gsed2212"
distribution
Distribution family: "LMS"
or "BCT"
age
Decimal age in years mu
M-curve, median D-score, P50 sigma
S-curve, spread expressed as coefficient of variation nu
L-curve, the lambda coefficient of the LMS/BCT model for skewness tau
Kurtosis parameter in the BCT model P3
P3 percentile P10
P10 percentile P25
P25 percentile P50
P50 percentile P75
P75 percentile P90
P90 percentile P97
P97 percentile SDM2
-2SD centile SDM1
-1SD centile SD0
0SD centile, median SDP1
+1SD centile SDP2
+2SD centile data.frame
with the following variables:Name Label population
Name of the reference population key
D-score key, e.g., "dutch"
, "gcdg"
or "gsed"
distribution
Distribution family: "LMS"
or "BCT"
age
Decimal age in years mu
M-curve, median D-score, P50 sigma
S-curve, spread expressed as coefficient of variation nu
L-curve, the lambda coefficient of the LMS/BCT model for skewness tau
Kurtosis parameter in the BCT model P3
P3 percentile P10
P10 percentile P25
P25 percentile P50
P50 percentile P75
P75 percentile P90
P90 percentile P97
P97 percentile SDM2
-2SD centile SDM1
-1SD centile SD0
0SD centile, median SDP1
+1SD centile SDP2
+2SD centile Details
- "dutch"
references were calculated from the SMOCC data, and cover
+ "dutch"
references were calculated from the SMOCC data, and cover
age range 0-2.5 years (van Buuren, 2014).
The "gcdg"
references were calculated from the 15 cohorts of the
GCDG-study, and cover age range 0-5 years (Weber, 2019).
@@ -125,7 +126,7 @@ See also
Examples
# get an overview of available references per key
-table(builtin_references$name, builtin_references$key)
+table(builtin_references$population, builtin_references$key)
#>
#> 293_0 dutch gcdg gsed1912 gsed2212 gsed2406
#> dutch 0 144 0 0 185 185
diff --git a/docs/reference/calculate_posterior.html b/docs/reference/calculate_posterior.html
index 031ba70..612eb15 100644
--- a/docs/reference/calculate_posterior.html
+++ b/docs/reference/calculate_posterior.html
@@ -18,7 +18,7 @@
Median D-score from the default references for the given key
+ Source: R/count_mu.R
+ count_mu.Rd
count_mu(t, key)
Arguments
+ Value
+
+
+length(t)
with the median of the default reference
+population for the key.Details
+ Arguments
data.frame
with the LMS or BCT reference values.
The default NULL
selects the default reference belonging to the key
,
-as specified in the default_reference
field in dscore::builtin_keys
.base_population
field in dscore::builtin_keys
.
String. The name of the reference to calculate DAZ.
-Use unique(builtin_references$name)
to obtain the set of currently
-available references.
String. The name of the reference population to calculate DAZ.
+Use unique(builtin_references$population)
to obtain the set of currently
+available reference populations.
data
with user-supplied values of the prior mean
for each observation. If a numeric scalar, it is used as the prior mean
for all observations. The default (NULL
) will consult the
-default_reference
field in builtin_keys
, and use the corresponding
+base_population
field in builtin_keys
, and use the corresponding
median of that reference as prior mean for the D-score.
@@ -205,11 +203,6 @@ builtin_keys
.
-A string describing the population of the -reference distribution.
A vector of two integers specifying the number of decimals for rounding the D-score and DAZ, respectively. @@ -261,15 +254,12 @@
"gsed2406"
.
The default starting prior is a mean calculated from a so-called
"Count model" that describes mean D-score as a function of age. The
-Count models are stored as functions
-count_mu_preliminary_standards()
-count_mu_phase1()
,
-count_mu_gcdg()
and
-count_mu_dutch()
. The spread of the starting prior
-is 5 D-score points around this mean D-score, which corresponds to
+The Count models are implemented in the function [count_mu()]
.
+By default, the spread of the starting prior
+is 5 D-score points around the mean D-score, which corresponds to
approximately 1.5 to 2 times the normal spread of child of a given age. The
-starting prior is somewhat informative for very short test (say <5 items)
-and uninformative for large number of items (say >10 items).
String. The name of the reference to calculate DAZ.
-Use unique(builtin_references$name)
to obtain the set of currently
-available references.
String. The name of the reference population to calculate DAZ.
+Use unique(builtin_references$population)
to obtain the set of currently
+available reference populations.
transform
and qp
arguments.
-A string describing the population of the -reference distribution.
A string specifying the unit in which age is measured
(either "decimal"
, "days"
or "months"
).
diff --git a/docs/reference/get_itemnames.html b/docs/reference/get_itemnames.html
index 4223325..0316dba 100644
--- a/docs/reference/get_itemnames.html
+++ b/docs/reference/get_itemnames.html
@@ -19,7 +19,7 @@
diff --git a/docs/reference/get_itemtable.html b/docs/reference/get_itemtable.html
index 6c2caed..ea32a51 100644
--- a/docs/reference/get_itemtable.html
+++ b/docs/reference/get_itemtable.html
@@ -21,7 +21,7 @@
diff --git a/docs/reference/get_labels.html b/docs/reference/get_labels.html
index b2b94c8..16d0275 100644
--- a/docs/reference/get_labels.html
+++ b/docs/reference/get_labels.html
@@ -18,7 +18,7 @@
diff --git a/docs/reference/get_reference.html b/docs/reference/get_reference.html
index c430daf..1c1eba8 100644
--- a/docs/reference/get_reference.html
+++ b/docs/reference/get_reference.html
@@ -18,7 +18,7 @@
@@ -69,9 +69,8 @@
get_reference(
- name = NULL,
- key = NULL,
population = NULL,
+ key = NULL,
references = dscore::builtin_references,
verbose = FALSE,
...
@@ -80,10 +79,10 @@ Get D-score reference
Arguments
- - name
-String. The name of the reference to calculate DAZ.
-Use unique(builtin_references$name)
to obtain the set of currently
-available references.
+ - population
+String. The name of the reference population to calculate DAZ.
+Use unique(builtin_references$population)
to obtain the set of currently
+available reference populations.
- key
@@ -92,11 +91,6 @@ Arguments
of the available keys.
-- population
-A string describing the population of the
-reference distribution.
-
-
- references
A data.frame
with the same structure as builtin_references
.
The default is to use builtin_references
.
@@ -150,8 +144,8 @@ See also
Examples
- # see key-name combinations of builtin_references
-table(builtin_references$key, builtin_references$name)
+ # see key-population combinations of builtin_references
+table(builtin_references$key, builtin_references$population)
#>
#> dutch gcdg phase1 preliminary_standards
#> 293_0 0 0 186 0
@@ -164,7 +158,7 @@ Examples
# get the default reference
reftab <- get_reference()
head(reftab, 2)
-#> name key distribution age mu sigma nu
+#> population key distribution age mu sigma nu
#> 944 preliminary_standards gsed2406 BCT 0.0383 13.18 0.2090 1.4212
#> 945 preliminary_standards gsed2406 BCT 0.0575 14.04 0.2016 1.4272
#> tau P3 P10 P25 P50 P75 P90 P97
@@ -176,18 +170,18 @@ Examples
# get the default reference for the key "gsed2212"
reftab <- get_reference(key = "gsed2212", verbose = TRUE)
-#> key: gsed2212
-#> name: phase1
+#> key: gsed2212
+#> population: phase1
# get dutch reference for default key
-reftab <- get_reference(name = "dutch", verbose = TRUE)
-#> key: gsed2406
-#> name: dutch
+reftab <- get_reference(population = "dutch", verbose = TRUE)
+#> key: gsed2406
+#> population: dutch
# loading a non-existing reference yields zero rows
-reftab <- get_reference(name = "france", verbose = TRUE)
-#> key: gsed2406
-#> name: france
+reftab <- get_reference(population = "france", verbose = TRUE)
+#> key: gsed2406
+#> population: france
#> Warning: Reference 'france' for key 'gsed2406' not found.
nrow(reftab)
#> [1] 0
diff --git a/docs/reference/get_tau.html b/docs/reference/get_tau.html
index a33ab68..a78001b 100644
--- a/docs/reference/get_tau.html
+++ b/docs/reference/get_tau.html
@@ -19,7 +19,7 @@
diff --git a/docs/reference/gsample.html b/docs/reference/gsample.html
index 8b64fd3..9b5e083 100644
--- a/docs/reference/gsample.html
+++ b/docs/reference/gsample.html
@@ -18,7 +18,7 @@
diff --git a/docs/reference/index.html b/docs/reference/index.html
index 4ce0c42..4474401 100644
--- a/docs/reference/index.html
+++ b/docs/reference/index.html
@@ -17,7 +17,7 @@
@@ -136,6 +136,10 @@ Internal functions For research purposes only.
+
+
+ Median D-score from the default references for the given key
+
Median of Dutch references
diff --git a/docs/reference/milestones.html b/docs/reference/milestones.html
index 4a55759..5523337 100644
--- a/docs/reference/milestones.html
+++ b/docs/reference/milestones.html
@@ -18,7 +18,7 @@
diff --git a/docs/reference/normalize.html b/docs/reference/normalize.html
index 261be83..c45c8a4 100644
--- a/docs/reference/normalize.html
+++ b/docs/reference/normalize.html
@@ -17,7 +17,7 @@
diff --git a/docs/reference/posterior.html b/docs/reference/posterior.html
index 5f9d61a..5e18697 100644
--- a/docs/reference/posterior.html
+++ b/docs/reference/posterior.html
@@ -17,7 +17,7 @@
diff --git a/docs/reference/rename_gcdg_gsed.html b/docs/reference/rename_gcdg_gsed.html
index 9ada805..6868615 100644
--- a/docs/reference/rename_gcdg_gsed.html
+++ b/docs/reference/rename_gcdg_gsed.html
@@ -18,7 +18,7 @@
diff --git a/docs/reference/sample_hf.html b/docs/reference/sample_hf.html
index ea963f3..0cedb0a 100644
--- a/docs/reference/sample_hf.html
+++ b/docs/reference/sample_hf.html
@@ -18,7 +18,7 @@
diff --git a/docs/reference/sample_lf.html b/docs/reference/sample_lf.html
index 5f92145..d3b194d 100644
--- a/docs/reference/sample_lf.html
+++ b/docs/reference/sample_lf.html
@@ -18,7 +18,7 @@
diff --git a/docs/reference/sample_sf.html b/docs/reference/sample_sf.html
index dce2f10..ac3da23 100644
--- a/docs/reference/sample_sf.html
+++ b/docs/reference/sample_sf.html
@@ -18,7 +18,7 @@
diff --git a/docs/reference/sort_itemnames.html b/docs/reference/sort_itemnames.html
index f6e06ff..5d2ed69 100644
--- a/docs/reference/sort_itemnames.html
+++ b/docs/reference/sort_itemnames.html
@@ -19,7 +19,7 @@
diff --git a/docs/sitemap.xml b/docs/sitemap.xml
index 500acfa..fe7504e 100644
--- a/docs/sitemap.xml
+++ b/docs/sitemap.xml
@@ -39,6 +39,9 @@
https://d-score.github.io/dscore/reference/calculate_posterior.html
+
+ https://d-score.github.io/dscore/reference/count_mu.html
+
https://d-score.github.io/dscore/reference/count_mu_dutch.html