Package 'cNORM'

Description

Selects polynomial degrees for the two components of a beta-binomial model ($\alpha$/$\beta$ when mode = 2 or $\mu$/$\sigma$ when mode = 1) by minimizing BIC.

Usage

autoselect.betabinomial(
  age,
  score,
  n = NULL,
  weights = NULL,
  mode = 2,
  max_alpha = 4,
  max_beta = 4,
  min_alpha = 1,
  min_beta = 1,
  control = NULL,
  scale = "T",
  parallel = TRUE,
  n_cores = NULL,
  plot = TRUE,
  verbose = TRUE
)

Arguments

Details

Parallel execution is attempted by default. If the workers cannot access the cNORM namespace, the function transparently falls back to sequential execution.

Value

The selected fitted model

Description

Selects the polynomial degrees for the four components of a cnorm.shash model ($\mu$, $\sigma$, $\epsilon$, $\delta$) by minimizing BIC over a full grid.

Usage

autoselect.shash(
  age,
  score,
  weights = NULL,
  max_mu = 4,
  max_sigma = 3,
  max_epsilon = 3,
  max_delta = 1,
  min_mu = 1,
  min_sigma = 1,
  min_epsilon = 1,
  min_delta = 0,
  delta = 1,
  control = NULL,
  scale = "T",
  parallel = TRUE,
  n_cores = NULL,
  plot = TRUE,
  verbose = TRUE
)

Arguments

Details

A degree of 0 for $\delta$ is interpreted as a **fixed** $\delta$ (passed to cnorm.shash as delta_degree = NULL) with value taken from the delta argument. This allows the grid search to choose between a fixed and a varying tail-weight model.

Parallel execution is attempted by default. If the workers cannot access the cNORM namespace, the function transparently falls back to sequential execution.

Value

The selected fitted cnormShash model with an additional element $selection containing:

• evaluated: data frame of every combination tried, sorted by BIC.

• selected: list with the chosen degrees and BIC.

See Also

cnorm.shash, autoselect.betabinomial

Examples

## Not run: 
m <- autoselect.shash(elfe$group, elfe$raw)
m$selection$evaluated
summary(m, age = elfe$group, score = elfe$raw)

## End(Not run)

Description

Computes Taylor polynomial regression models by evaluating a series of models with increasing predictors. It aims to find a consistent model that effectively captures the variance in the data. It draws on the regsubsets function from the leaps package and builds up to 20 models for each number of predictors, evaluates these models regarding model consistency and selects consistent model with the highest R^2. This automatic model selection should usually be accompanied with visual inspection of the percentile plots and assessment of fit statistics. Set R^2 or number of terms manually to retrieve a more parsimonious model, if desired.

Usage

bestModel(
  data,
  raw = NULL,
  R2 = NULL,
  k = NULL,
  t = NULL,
  predictors = NULL,
  terms = 0,
  weights = NULL,
  force.in = NULL,
  plot = TRUE,
  extensive = TRUE,
  subsampling = FALSE
)

Arguments

Details

The functions rankBySlidingWindow, rankByGroup, bestModel, computePowers and prepareData are usually not called directly, but accessed through other functions like cnorm.

Additional functions like plotSubset(model) and cnorm.cv can aid in model evaluation.

Value

The model. Further exploration can be done using plotSubset(model) and plotPercentiles(data, model).

See Also

plotSubset, plotPercentiles, plotPercentileSeries, checkConsistency

Other model: checkConsistency(), cnorm.cv(), derive(), modelSummary(), print.cnorm(), printSubset(), rangeCheck(), regressionFunction(), summary.cnorm()

Examples

# Example with sample data
## Not run: 
# It is not recommended to directly use this function. Rather use 'cnorm' instead.
normData <- prepareData(elfe)
model <- bestModel(normData)
plotSubset(model)
plotPercentiles(buildCnormObject(normData, model))

# Specifying variables explicitly
preselectedModel <- bestModel(normData, predictors = c("L1", "L3", "L1A3", "A2", "A3"))
print(regressionFunction(preselectedModel))

## End(Not run)

Description

This function calculates the $\alpha$ (a) and $\beta$ (b) parameters of a beta binomial distribution, along with the mean (m), variance (var) based on the input vector 'x' and the maximum number 'n'.

Usage

betaCoefficients(x, n = NULL)

Arguments

Details

The beta-binomial distribution is a discrete probability distribution that models the number of successes in a fixed number of trials, where the probability of success varies from trial to trial. This variability in success probability is modeled by a beta distribution. Such a calculation is particularly relevant in scenarios where there is heterogeneity in success probabilities across trials, which is common in real-world situations, as for example the number of correct solutions in a psychometric test, where the test has a fixed number of items.

Value

A numeric vector containing the calculated parameters in the following order: alpha (a), beta (b), mean (m), standard deviation (sd), and the maximum number (n).

Description

Helper function to build a cnorm object from a data object and a model object from the bestModel function for compatibility reasons.

Usage

buildCnormObject(data, model)

Arguments

Value

A cnorm object

Examples

## Not run: 
  data <- prepareData(elfe)
  model <- bestModel(data, k = 4)
  model.cnorm <- buildCnormObject(data, model)

## End(Not run)

Description

Build regression function for bestModel

Usage

buildFunction(raw, k, t, age)

Arguments

Value

regression function

Description

The function is an inline for searching zeros in the inverse regression function. It collapses the regression function at a specific age and simplifies the coefficients.

Usage

calcPolyInL(raw, age, model)

Arguments

Value

The coefficients

Description

The function is an inline for searching zeros in the inverse regression function. It collapses the regression function at a specific age and simplifies the coefficients. Optimized version of the prior 'calcPolyInLBase'

Usage

calcPolyInLBase2(raw, age, coeff, k)

Arguments

Value

The coefficients

Description

By the courtesy of the Center of Disease Control (CDC), cNORM includes human growth data for children and adolescents age 2 to 25 that can be used to model trajectories of the body mass index and to estimate percentiles for clinical definitions of under- and overweight. The data stems from the NHANES surveys in the US and was published in 2012 as public domain. The data was cleaned by removing missing values and it includes the following variables from or based on the original dataset.

Usage

CDC

Format

A data frame with 45053 rows and 7 variables:

age

continuous age in years, based on the month variable

group

age group; chronological age in years at the time of examination

month

chronological age in month at the time of examination

sex

sex of the participant, 1 = male, 2 = female

height

height of the participants in cm

weight

weight of the participants in kg

bmi

the body mass index, computed by (weight in kg)/(height in m)^2

A data frame with 45035 rows and 7 columns

References

Center for Disease Control and Prevention (2012). National Health and Nutrition Examination Survey: Questionnaires, datasets and related documentation. U.S. Department of Health and Human Services (original source not available anymore).

Description

While abilities increase and decline over age, within one age group, the norm scores always have to show a monotonic increase or decrease with increasing raw scores. Violations of this assumption are an indication for problems in modeling the relationship between raw and norm scores. There are several reasons, why this might occur:

1. Vertical extrapolation: Choosing extreme norm scores, e. g. values -3 <= x and x >= 3 In order to model these extreme values, a large sample dataset is necessary.

2. Horizontal extrapolation: Taylor polynomials converge in a certain radius. Using the model values outside the original dataset may lead to inconsistent results.

3. The data cannot be modeled with Taylor polynomials, or you need another power parameter (k) or R2 for the model.

Usage

checkConsistency(
  model,
  minAge = NULL,
  maxAge = NULL,
  minNorm = NULL,
  maxNorm = NULL,
  minRaw = NULL,
  maxRaw = NULL,
  stepAge = NULL,
  stepNorm = 1,
  warn = FALSE,
  silent = FALSE
)

Arguments

Details

In general, extrapolation (point 1 and 2) can carefully be done to a certain degree outside the original sample, but it should in general be handled with caution. Please note that at extreme values, the models most likely become independent and it is thus recommended to restrict the norm score range to the relevant range of abilities, e.g. +/- 2.5 SD via the minNorm and maxNorm parameter.

Value

Boolean, indicating model violations (TRUE) or no problems (FALSE)

See Also

Other model: bestModel(), cnorm.cv(), derive(), modelSummary(), print.cnorm(), printSubset(), rangeCheck(), regressionFunction(), summary.cnorm()

Examples

## Not run: 
  model <- cnorm(raw = elfe$raw, group = elfe$group, plot = FALSE)
  modelViolations <- checkConsistency(model, minNorm = 25, maxNorm = 75)
  plotDerivative(model, minNorm = 25, maxNorm = 75)

## End(Not run)

Description

Check, if NA or values <= 0 occur and issue warning

Usage

checkWeights(weights)

Arguments

Description

Conducts continuous norming in one step and returns an object including ranked raw data and the continuous norming model. Please consult the function description of 'rankByGroup', 'rankBySlidingWindow' and 'bestModel' for specifics of the steps in the data preparation and modeling process. In addition to the raw scores, either provide

• a numeric vector for the grouping information (group)

• a numeric age vector and the width of the sliding window (age, width)

for the ranking of the raw scores. You can adjust the grade of smoothing of the regression model by setting the k and terms parameter. In general, increasing k to more than 4 and the number of terms lead to a higher fit, while lower values lead to more smoothing. The power parameter for the age trajectory can be specified independently by 't'. If both parameters are missing, cnorm uses k = 5 and t = 3 by default.

Usage

cnorm(
  raw = NULL,
  group = NULL,
  age = NULL,
  width = NA,
  weights = NULL,
  scale = "T",
  method = 4,
  descend = FALSE,
  k = NULL,
  t = NULL,
  terms = 0,
  R2 = NULL,
  plot = TRUE,
  extensive = TRUE,
  subsampling = FALSE
)

Arguments

Value

cnorm object including the ranked raw data and the regression model.

References

1. Gary, S. & Lenhard, W. (2021). In norming we trust. Diagnostica.

2. Gary, S., Lenhard, W. & Lenhard, A. (2021). Modelling Norm Scores with the cNORM Package in R. Psych, 3(3), 501-521. https://doi.org/10.3390/psych3030033

3. Gary, S., Lenhard, W., Lenhard, A., & Herzberg, D. (2023). A tutorial on automatic post-stratification and weighting in conventional and regression-based norming of psychometric tests. Behavior Research Methods. https://doi.org/10.3758/s13428-023-02207-0

4. Gary, S., Lenhard, A., Lenhard, W., & Herzberg, D. S. (2023). Reducing the Bias of Norm Scores in Non-Representative Samples: Weighting as an Adjunct to Continuous Norming Methods. Assessment, 30(8), 2491–2509. https://doi.org/10.1177/10731911231153832

5. Lenhard, A., Lenhard, W., Suggate, S. & Segerer, R. (2016). A continuous solution to the norming problem. Assessment, Online first, 1-14. doi:10.1177/1073191116656437

6. Lenhard, A., Lenhard, W., Gary, S. (2018). Continuous Norming (cNORM). The Comprehensive R Network, Package cNORM, available: https://CRAN.R-project.org/package=cNORM

7. Lenhard, A., Lenhard, W., Gary, S. (2019). Continuous norming of psychometric tests: A simulation study of parametric and semi-parametric approaches. PLoS ONE, 14(9), e0222279. doi:10.1371/journal.pone.0222279

8. Lenhard, W., & Lenhard, A. (2020). Improvement of Norm Score Quality via Regression-Based Continuous Norming. Educational and Psychological Measurement(Online First), 1-33. https://doi.org/10.1177/0013164420928457

See Also

rankByGroup, rankBySlidingWindow, computePowers, bestModel

Examples

## Not run: 
# Conventional norming
cnorm(raw = elfe$raw)

# Continuous norming via group
m1 <- cnorm(raw = elfe$raw, group = elfe$group)

# Continuous norming via continuous age + sliding window
m2 <- cnorm(raw = ppvt$raw, age = ppvt$age, width = 1)

# Norm tables with confidence intervals
normTable(c(2.125, 2.375, 2.625), m1, CI = .9, reliability = .95)
rawTable (c(2.125, 2.375, 2.625), m1, CI = .9, reliability = .95)

## End(Not run)

Description

This function fits a beta-binomial regression model where both the $\alpha$ and $\beta$ parameters of the beta-binomial distribution are modeled as polynomial functions of the predictor variable (typically age). Setting mode to 1 fits a beta-binomial model on the basis of $\mu$ and $\sigma$, setting it to 2 (default) fits a beta-binomial model directly on the basis of $\alpha$ and $\beta$.

Usage

cnorm.betabinomial(
  age,
  score,
  n = NULL,
  weights = NULL,
  mode = 2,
  alpha = 3,
  beta = 3,
  control = NULL,
  scale = "T",
  plot = T
)

Arguments

Details

The function standardizes the input variables, fits polynomial models for both the alpha and beta parameters, and uses maximum likelihood estimation to find the optimal parameters. The optimization is performed using the L-BFGS-B method.

Value

A list of class "cnormBetaBinomial" or "cnormBetaBinomial2". In case of mode 2 containing:

Examples

## Not run: 
# Fit a beta-binomial regression model to the PPVT data
model <- cnorm.betabinomial(ppvt$age, ppvt$raw, n = 228)
summary(model)

# Use weights for post-stratification
marginals <- data.frame(var = c("sex", "sex", "migration", "migration"),
                        level = c(1,2,0,1),
                        prop = c(0.51, 0.49, 0.65, 0.35))
weights <- computeWeights(ppvt, marginals)
model <- cnorm.betabinomial(ppvt$age, ppvt$raw, n = 228, weights = weights)

## End(Not run)

Description

This function fits a beta-binomial regression model where both the alpha and beta parameters of the beta-binomial distribution are modeled as polynomial functions of the predictor variable (typically age). While 'cnorm-betabinomial' fits a beta-binomial model on the basis of $\mu$ and $\sigma$, this function fits a beta-binomial model directly on the basis of $\gamma$ and $\beta$.

Usage

cnorm.betabinomial2(
  age,
  score,
  n = NULL,
  weights = NULL,
  alpha_degree = 3,
  beta_degree = 3,
  control = NULL,
  scale = "T",
  plot = TRUE
)

Arguments

Details

The function standardizes the input variables, fits polynomial models for both the alpha and beta parameters, and uses maximum likelihood estimation to find the optimal parameters. The optimization is performed using the L-BFGS-B method.

Value

A list of class "cnormBetaBinomial2" containing:

Description

Assists in determining the optimal number of terms for the regression model using repeated Monte Carlo cross-validation. It leverages an 80-20 split between training and validation data, with stratification by norm group or random sample in case of using sliding window ranking.

Usage

cnorm.cv(
  data,
  formula = NULL,
  repetitions = 5,
  norms = TRUE,
  min = 1,
  max = 12,
  cv = "full",
  pCutoff = NULL,
  width = NA,
  raw = NULL,
  group = NULL,
  age = NULL,
  weights = NULL
)

Arguments

Details

Successive models, with an increasing number of terms, are evaluated, and the RMSE for raw scores plotted. This encompasses the training, validation, and entire dataset. If 'norms' is set to TRUE (default), the function will also calculate the mean norm score reliability and crossfit measures. Note that due to the computational requirements of norm score calculations, execution can be slow, especially with numerous repetitions or terms.

When 'cv' is set to "full" (default), both test and validation datasets are ranked separately, providing comprehensive cross-validation. For a more streamlined validation process focused only on modeling, a pre-ranked dataset can be used. The output comprises RMSE for raw score models, norm score R^2, delta R^2, crossfit, and the norm score SE according to Oosterhuis, van der Ark, & Sijtsma (2016).

This function is not yet prepared for the 'extensive' search strategy, introduced in version 3.3, but instead relies on the first model per number of terms, without consistency check.

For assessing overfitting:

$CROSSFIT = R(Training; Model)^2 / R(Validation; Model)^2$

A CROSSFIT > 1 suggests overfitting, < 1 suggests potential underfitting, and values around 1 are optimal, given a low raw score RMSE and high norm score validation R^2.

Suggestions for ideal model selection:

• Visual inspection of percentiles with 'plotPercentiles' or 'plotPercentileSeries'.

• Pair visual inspection with repeated cross-validation (e.g., 10 repetitions).

• Aim for low raw score RMSE and high norm score R^2, avoiding terms with significant overfit (e.g., crossfit > 1.1).

Value

Table with results per term number: RMSE for raw scores, R^2 for norm scores, and crossfit measure.

References

Oosterhuis, H. E. M., van der Ark, L. A., & Sijtsma, K. (2016). Sample Size Requirements for Traditional and Regression-Based Norms. Assessment, 23(2), 191–202. https://doi.org/10.1177/1073191115580638

See Also

Other model: bestModel(), checkConsistency(), derive(), modelSummary(), print.cnorm(), printSubset(), rangeCheck(), regressionFunction(), summary.cnorm()

Examples

## Not run: 
# Example: Plot cross-validation RMSE by number of terms (up to 9) with three repetitions.
result <- cnorm(raw = elfe$raw, group = elfe$group)
cnorm.cv(result$data, min = 2, max = 9, repetitions = 3)

# Using a cnorm object examines the predefined formula.
cnorm.cv(result, repetitions = 1)

## End(Not run)

Description

Launcher for the graphical user interface of cNORM for distribution free continuous norming

Usage

cNORM.GUI(launch.browser = TRUE)

Arguments

Examples

## Not run: 
# Launch graphical user interface
cNORM.GUI()

## End(Not run)

Description

Launch the cNORM Parametric Modeling Shiny Application

Usage

cNORM.GUI2(launch.browser = TRUE)

Arguments

Examples

## Not run: 
# Launch graphical user interface
cNORM.GUI2()

## End(Not run)

Description

This function fits a Sinh-Arcsinh (shash; Jones & Pewsey, 2009) regression model for continuous norm score modeling, where the distribution parameters vary smoothly as polynomial functions of age or other predictors. The shash distribution is well-suited for psychometric data as it can flexibly model skewness and tail weight independently, making it ideal for handling floor effects, ceiling effects, and varying degrees of individual differences across age groups. In a simulation study (Lenhard et al, 2019), the shash model demonstrated superior performance compared to other parametric approaches from the Box Cox family of functions. In contrast to Box Cox, Sinh-Arcsinh can model distributions including zero and negativ values.

Usage

cnorm.shash(
  age,
  score,
  weights = NULL,
  mu_degree = 3,
  sigma_degree = 2,
  epsilon_degree = 2,
  delta_degree = 1,
  delta = 1,
  control = NULL,
  scale = "T",
  plot = TRUE
)

Arguments

Details

This implementation uses the Jones & Pewsey (2009) parameterization of the Sinh-Arcsinh distribution. Parameters are estimated using maximum likelihood via the L-BFGS-B algorithm. In case, optimization fails, try reducing model complexity by reducing polynomial degrees or fixing the delta parameter.

The Sinh-Arcsinh Distribution

The shash distribution is defined by the transformation:

$X = \mu + \sigma \cdot \sinh\left(\frac{\text{arcsinh}(Y) - \epsilon}{\delta}\right)$

where Y is a standard normal variable, Y ~ N(0,1).

This transformation provides:

• mu: Location parameter (similar to mean)

• sigma: Scale parameter (similar to standard deviation)

• epsilon: Skewness parameter (epsilon = 0 for symmetry)

• delta: Tail weight parameter (delta = 1 for normal-like tails)

Model Selection

Choose polynomial degrees based on:

• Sample size (higher degrees need more data)

• Theoretical expectations about developmental trajectories

• Model comparison criteria (AIC, BIC)

• Visual inspection of fitted curves

For most applications, mu_degree = 3, sigma_degree = 2, epsilon_degree = 2, delta_degree = 1 provides a good balance of flexibility and parsimony.

Value

An object of class "cnormShash" containing the fitted model results. This is a list with components:

Note

• The function requires the input data to have sufficient variability. Very small datasets or datasets with little age spread may cause convergence problems.

• Polynomial models can exhibit edge effects at the boundaries of the age range. Predictions outside the observed age range should be made cautiously.

• If convergence fails, try: (1) reducing polynomial degrees, (2) adjusting the delta parameter, (3) providing custom control parameters, or (4) checking for data quality issues.

• The tail weight parameter delta is fixed across ages by default. For applications where tail behavior changes substantially with age, consider setting the delta_degree paramerer to 1 or 2.

Author(s)

Wolfgang Lenhard

References

Jones, M. C., & Pewsey, A. (2009). Sinh-arcsinh distributions. *Biometrika*, 96(4), 761-780.

Lenhard, A., Lenhard, W., Gary, S. (2019). Continuous norming of psychometric tests: A simulation study of parametric and semi-parametric approaches. *PLoS ONE*, 14(9), e0222279. https://doi.org/10.1371/journal.pone.0222279

See Also

plot for plotting fitted models, predict for generating predictions, cnorm.betabinomial2 for discrete beta-binomial alternative

Examples

## Not run: 
# Basic usage with default settings
model <- cnorm.shash(age = children$age, score = children$raw_score)

# Custom polynomial degrees for complex developmental pattern
model_complex <- cnorm.shash(
  age = adolescents$age,
  score = adolescents$vocabulary_score,
  mu_degree = 4,         # Complex mean trajectory
  sigma_degree = 3,      # Changing variability pattern
  epsilon_degree = 2,    # Skewness shifts
  delta_degree = NULL, # set to NULL to activate fixed delta
  delta = 1.3            # Slightly heavy tails
)

# With sampling weights
model_weighted <- cnorm.shash(
  age = survey$age,
  score = survey$score,
  weights = survey$sample_weight
)

# Custom optimization control for difficult convergence
model_robust <- cnorm.shash(
  age = mixed$age,
  score = mixed$score,
  control = list(factr = 1e-6, maxit = 2000),
  delta = 1.5
)

# Compare model fit
compare(model, model_complex)

## End(Not run)

Description

This function creates a visualization comparing two norm models by displaying their percentile curves. The first model is shown with solid lines, the second with dashed lines. If age and score vectors are provided, manifest percentiles are displayed as dots. The function works with regular cnorm models, beta-binomial models, and shash models, allowing comparison between different model types.

Usage

compare(
  model1,
  model2,
  percentiles = c(0.025, 0.1, 0.25, 0.5, 0.75, 0.9, 0.975),
  age = NULL,
  score = NULL,
  weights = NULL,
  title = NULL,
  subtitle = NULL
)

Arguments

Value

A ggplot object showing the comparison of both models

See Also

Other plot: plot.cnorm(), plot.cnormBetaBinomial(), plot.cnormBetaBinomial2(), plotDensity(), plotDerivative(), plotNorm(), plotNormCurves(), plotPercentileSeries(), plotPercentiles(), plotRaw(), plotSubset()

Examples

## Not run: 
# Compare different types of models
model1 <- cnorm(group = elfe$group, raw = elfe$raw)
model2 <- cnorm.betabinomial(elfe$group, elfe$raw)
model3 <- cnorm.shash(elfe$group, elfe$raw)

# Compare traditional cnorm with shash
compare(model1, model3, age = elfe$group, score = elfe$raw)

# Compare beta-binomial with shash
compare(model2, model3, age = elfe$group, score = elfe$raw)

## End(Not run)

Description

The function computes powers of the norm variable e. g. T scores (location, L), an explanatory variable, e. g. age or grade of a data frame (age, A) and the interactions of both (L X A). The k variable indicates the degree up to which powers and interactions are build. These predictors can be used later on in the bestModel function to model the norm sample. Higher values of k allow for modeling the norm sample closer, but might lead to over-fit. In general k = 3 or k = 4 (default) is sufficient to model human performance data. For example, k = 2 results in the variables L1, L2, A1, A2, and their interactions L1A1, L2A1, L1A2 and L2A2 (but k = 2 is usually not sufficient for the modeling). Please note, that you do not need to use a normal rank transformed scale like T or IQ; you can use the percentiles for the 'normValue' as well.

Usage

computePowers(data, k = 5, norm = NULL, age = NULL, t = 3, silent = FALSE)

Arguments

Details

The functions rankBySlidingWindow, rankByGroup, bestModel, computePowers and prepareData are usually not called directly, but accessed through other functions like cnorm.

Value

data.frame with the powers and interactions of location and explanatory variable / age

See Also

bestModel

Other prepare: prepareData(), rankByGroup(), rankBySlidingWindow()

Examples

## Not run: 
# Dataset with grade levels as grouping
data.elfe <- rankByGroup(elfe)
data.elfe <- computePowers(data.elfe)

# Dataset with continuous age variable and k = 5
data.ppvt <- rankByGroup(ppvt)
data.ppvt <- computePowers(data.ppvt, age = "age", k = 5)

## End(Not run)

Description

Computes and standardizes weights via raking to compensate for non-stratified samples. It is based on the implementation in the survey R package. It reduces data collection #' biases in the norm data by the means of post stratification, thus reducing the effect of unbalanced data in percentile estimation and norm data modeling.

Usage

computeWeights(data, population.margins, standardized = TRUE)

Arguments

Details

This function computes standardized raking weights to overcome biases in norm samples. It generates weights, by drawing on the information of population shares (e. g. for sex, ethnic group, region ...) and subsequently reduces the influence of over-represented groups or increases underrepresented cases. The returned weights are either raw or standardized and scaled to be larger than 0.

Raking in general has a number of advantages over post stratification and it additionally allows cNORM to draw on larger datasets, since less cases have to be removed during stratification. To use this function, additionally to the data, a data frame with stratification variables has to be specified. The data frame should include a row with (a) the variable name, (b) the level of the variable and (c) the according population proportion.

Value

a vector with the standardized weights

Examples

## Not run: 
# cNORM features a dataset on vocabulary development (ppvt)
# that includes variables like sex or migration. In order
# to weight the data, we have to specify the population shares.
# According to census, the population includes 52% boys
# (factor level 1 in the ppvt dataset) and 70% / 30% of persons
# without / with a a history of migration (= 0 / 1 in the dataset).
# First we set up the popolation margins with all shares of the
# different levels:

margins <- data.frame(variables = c("sex", "sex",
                                    "migration", "migration"),
                      levels = c(1, 2, 0, 1),
                      share = c(.52, .48, .7, .3))
head(margins)

# Now we use the population margins to generate weights
# through raking

weights <- computeWeights(ppvt, margins)


# There are as many different weights as combinations of
# factor levels, thus only four in this specific case

unique(weights)


# To include the weights in the cNORM modelling, we have
# to pass them as weights. They are then used to set up
# weighted quantiles and as weights in the regession.

model <- cnorm(raw = ppvt$raw,
               group=ppvt$group,
               weights = weights)

## End(Not run)

Description

In order to check model assumptions, a table of the first order derivative of the model coefficients is created.

Usage

derivationTable(A, model, minNorm = NULL, maxNorm = NULL, step = 0.1)

Arguments

Value

data.frame with norm scores and the predicted scores based on the derived regression function

See Also

plotDerivative, derive

Other predict: getNormCurve(), normTable(), predict.cnormBetaBinomial(), predict.cnormBetaBinomial2(), predict.cnormShash(), predictNorm(), predictRaw(), rawTable()

Examples

## Not run: 
# Generate cnorm object from example data
cnorm.elfe <- cnorm(raw = elfe$raw, group = elfe$group)

# retrieve function for time point 6
d <- derivationTable(6, cnorm.elfe, step = 0.5)

## End(Not run)

Description

Calculates the derivative of the location / norm value from the regression model with the first derivative as the default. This is useful for finding violations of model assumptions and problematic distribution features as f. e. bottom and ceiling effects, non-progressive norm scores within an age group or in general #' intersecting percentile curves.

Usage

derive(model, order = 1)

Arguments

Value

The derived coefficients

See Also

Other model: bestModel(), checkConsistency(), cnorm.cv(), modelSummary(), print.cnorm(), printSubset(), rangeCheck(), regressionFunction(), summary.cnorm()

Examples

## Not run: 
  m <- cnorm(raw = elfe$raw, group = elfe$group)
  derivedCoefficients <- derive(m, order=1)

## End(Not run)

Description

This function provides diagnostic information for a fitted beta-binomial model from the cnorm.betabinomial function. It returns various metrics related to model convergence, fit, and complexity. In case, age and raw scores are provided, the function as well computes R2, rmse and bias for the norm scores based on the manifest and predicted norm scores.

Usage

diagnostics.betabinomial(model, age = NULL, score = NULL, weights = NULL)

Arguments

Details

The AIC and BIC are calculated as: AIC = 2k - 2ln(L) BIC = ln(n)k - 2ln(L) where k is the number of parameters, L is the maximum likelihood, and n is the number of observations.

Value

A list containing the following diagnostic information:

• converged: Logical indicating whether the optimization algorithm converged.

• n_evaluations: Number of function evaluations performed during optimization.

• n_gradient: Number of gradient evaluations performed during optimization.

• final_value: Final value of the objective function (negative log-likelihood).

• message: Any message returned by the optimization algorithm.

• AIC: Akaike Information Criterion.

• BIC: Bayesian Information Criterion.

• max_gradient: Maximum absolute gradient at the solution (if available).

Examples

## Not run: 
# Fit a beta-binomial model
model <- cnorm.betabinomial(ppvt$age, ppvt$raw)

# Get diagnostic information
diag_info <- diagnostics.betabinomial(model)

# Print the diagnostic information
print(diag_info)

# Summary the diagnostic information
summary(diag_info)

# Check if the model converged
if(diag_info$converged) {
  cat("Model converged successfully.\n")
} else {
  cat("Warning: Model did not converge.\n")
}

# Compare AIC and BIC
cat("AIC:", diag_info$AIC, "\n")
cat("BIC:", diag_info$BIC, "\n")

## End(Not run)

Description

A dataset containing the raw data of 1400 students from grade 2 to 5 in the sentence comprehension test from ELFE 1-6 (Lenhard & Schneider, 2006). In this test, students are presented lists of sentences with one gap. The student has to fill in the correct solution by selecting from a list of 5 alternatives per sentence. The alternatives include verbs, adjectives, nouns, pronouns and conjunctives. Each item stems from the same word type. The text is speeded, with a time cutoff of 180 seconds. The variables are as follows:

Usage

elfe

Format

A data frame with 1400 rows and 3 variables:

personID

ID of the student

group

grade level, with x.5 indicating the end of the school year and x.0 indicating the middle of the school year

raw

the raw score of the student, spanning values from 0 to 28

A data frame with 1400 rows and 3 columns

Source

https://www.psychometrica.de/elfe2.html

References

Lenhard, W. & Schneider, W.(2006). Ein Leseverstaendnistest fuer Erst- bis Sechstklaesser. Goettingen/Germany: Hogrefe.

Examples

## Not run: 
  # prepare data, retrieve model and plot percentiles
  model <- cnorm(group = elfe$group, raw = elfe$raw)

## End(Not run)

Description

Splits a continuous variable into n groups (either equidistant or with equal number of observations) and returns the group mean for every observation.

Usage

getGroups(x, n = NULL, equidistant = FALSE)

Arguments

Value

Numeric vector of the same length as x with the group mean assigned to each observation.

Examples

## Not run: 
x <- rnorm(1000, mean = 50, sd = 10)
m <- getGroups(x, n = 10)

## End(Not run)

Description

As with this continuous norming regression approach, raw scores are modeled as a function of age and norm score (location), getNormCurve is a straightforward approach to show the raw score development over age, while keeping the norm value constant. This way, e. g. academic performance or intelligence development of a specific ability is shown.

Usage

getNormCurve(
  norm,
  model,
  minAge = NULL,
  maxAge = NULL,
  step = 0.1,
  minRaw = NULL,
  maxRaw = NULL
)

Arguments

Value

data.frame of the variables raw, age and norm

See Also

Other predict: derivationTable(), normTable(), predict.cnormBetaBinomial(), predict.cnormBetaBinomial2(), predict.cnormShash(), predictNorm(), predictRaw(), rawTable()

Examples

## Not run: 
# Generate cnorm object from example data
cnorm.elfe <- cnorm(raw = elfe$raw, group = elfe$group)
getNormCurve(35, cnorm.elfe)

## End(Not run)

Description

Calculates the standard error (SE) or root mean square error (RMSE) of the norm scores In case of large datasets, both results should be almost identical

Usage

getNormScoreSE(model, type = 2)

Arguments

Value

The standard error (SE) of the norm scores sensu Oosterhuis et al. (2016) or the RMSE

References

Oosterhuis, H. E. M., van der Ark, L. A., & Sijtsma, K. (2016). Sample Size Requirements for Traditional and Regression-Based Norms. Assessment, 23(2), 191–202. https://doi.org/10.1177/1073191115580638

Description

Prints the results and regression function of a cnorm model

Usage

modelSummary(object, ...)

Arguments

Value

A report on the regression function, weights, R2 and RMSE

See Also

Other model: bestModel(), checkConsistency(), cnorm.cv(), derive(), print.cnorm(), printSubset(), rangeCheck(), regressionFunction(), summary.cnorm()

Description

This function generates a norm table for a specific age based on the regression model by assigning raw scores to norm scores. Please specify the range of norm scores, you want to cover. A T value of 25 corresponds to a percentile of .6. As a consequence, specifying a range of T = 25 to T = 75 would cover 98.4 the population. Please be careful when extrapolating vertically (at the lower and upper end of the age specific distribution). Depending on the size of your standardization sample, extreme values with T < 20 or T > 80 might lead to inconsistent results. In case a confidence coefficient (CI, default .9) and the reliability is specified, confidence intervals are computed for the true score estimates, including a correction for regression to the mean (Eid & Schmidt, 2012, p. 272).

Usage

normTable(
  A,
  model,
  minNorm = NULL,
  maxNorm = NULL,
  minRaw = NULL,
  maxRaw = NULL,
  step = NULL,
  monotonuous = TRUE,
  CI = 0.9,
  reliability = NULL,
  pretty = TRUE
)

Arguments

Value

either data.frame with norm scores, predicted raw scores and percentiles in case of simple A value or a list of norm tables if vector of A values was provided

References

Eid, M. & Schmidt, K. (2012). Testtheorie und Testkonstruktion. Hogrefe.

See Also

rawTable

Other predict: derivationTable(), getNormCurve(), predict.cnormBetaBinomial(), predict.cnormBetaBinomial2(), predict.cnormShash(), predictNorm(), predictRaw(), rawTable()

Examples

## Not run: 
# Generate cnorm object from example data
cnorm.elfe <- cnorm(raw = elfe$raw, group = elfe$group)

# create single norm table
norms <- normTable(3.5, cnorm.elfe, minNorm = 25, maxNorm = 75, step = 0.5)

# create list of norm tables
norms <- normTable(c(2.5, 3.5, 4.5), cnorm.elfe,
  minNorm = 25, maxNorm = 75,
  step = 1, minRaw = 0, maxRaw = 26
)

# conventional norming, set age to arbitrary value
model <- cnorm(raw = elfe$raw)
normTable(0, model)

## End(Not run)

Description

This function generates a norm table for a specific ages based on the beta binomial regression model. In case a confidence coefficient (CI, default .9) and the reliability is specified, confidence intervals are computed for the true score estimates, including a correction for regression to the mean (Eid & Schmidt, 2012, p. 272).

Usage

normTable.betabinomial(
  model,
  ages,
  n = NULL,
  m = NULL,
  range = 3,
  CI = 0.9,
  reliability = NULL
)

Arguments

Value

A list of data frames with columns: x, Px, Pcum, Percentile, z, norm score and possibly confidence interval

Description

Generates norm tables for specific ages based on a fitted SinH-ArcSinH (shash) regression model. Computes probabilities, percentiles, z-scores, and norm scores for a specified range of raw scores. Optionally includes confidence intervals when reliability is provided.

Usage

normTable.shash(
  model,
  ages,
  start = NULL,
  end = NULL,
  step = 1,
  CI = 0.9,
  reliability = NULL
)

Arguments

Details

For continuous shash distributions, probability densities are computed and converted to cumulative probabilities and percentiles. When reliability is specified, confidence intervals include correction for regression to the mean.

Value

List of data frames (one per age) containing:

Examples

## Not run: 
# Basic norm table
model <- cnorm.shash(age, score)
tables <- normTable.shash(model, ages = c(7, 8, 9), start = 0, end = 50)

# With confidence intervals and finer granularity
tables_ci <- normTable.shash(model, ages = c(8, 9), start = 10, end = 40,
                             step = 0.5, CI = 0.95, reliability = 0.85)

## End(Not run)

Description

S3 function for plotting cnorm objects

Usage

## S3 method for class 'cnorm'
plot(x, y, ...)

Arguments

See Also

Other plot: compare(), plot.cnormBetaBinomial(), plot.cnormBetaBinomial2(), plotDensity(), plotDerivative(), plotNorm(), plotNormCurves(), plotPercentileSeries(), plotPercentiles(), plotRaw(), plotSubset()

Description

This function creates a visualization of a fitted cnormBetaBinomial model, including the original data points manifest percentiles and specified percentile lines. Note that the beta-binomial model aims at discrete raw scores. We decided to display continuous percentile lines nonetheless on order to maintain visual comparability with other modelling techniques. If you prefer discretization, set the "discrete" parameter to TRUE.

Usage

## S3 method for class 'cnormBetaBinomial'
plot(x, ...)

Arguments

Value

A ggplot object.

See Also

Other plot: compare(), plot.cnorm(), plot.cnormBetaBinomial2(), plotDensity(), plotDerivative(), plotNorm(), plotNormCurves(), plotPercentileSeries(), plotPercentiles(), plotRaw(), plotSubset()

Description

This function creates a visualization of a fitted cnormBetaBinomial model, including the original data points manifest percentiles and specified percentile lines.

Usage

## S3 method for class 'cnormBetaBinomial2'
plot(x, ...)

Arguments

Value

A ggplot object.

See Also

Other plot: compare(), plot.cnorm(), plot.cnormBetaBinomial(), plotDensity(), plotDerivative(), plotNorm(), plotNormCurves(), plotPercentileSeries(), plotPercentiles(), plotRaw(), plotSubset()

Description

Plot SinH-ArcSinH Model with Data and Percentile Lines

Usage

## S3 method for class 'cnormShash'
plot(x, ...)

Arguments

Value

A ggplot object

Description

General convenience plotting function

Usage

plotCnorm(x, y, ...)

Arguments

Description

This function plots density curves based on the regression model against the raw scores. It supports both traditional continuous norming models and beta-binomial models. The function allows for customization of the plot range and groups to be displayed.

Usage

plotDensity(
  model,
  minRaw = NULL,
  maxRaw = NULL,
  minNorm = NULL,
  maxNorm = NULL,
  group = NULL
)

Arguments

Details

The function generates density curves for specified age groups, allowing for easy comparison of score distributions across different ages.

For beta-binomial models, the density is based on the probability mass function, while for traditional models, it uses a normal distribution based on the norm scores.

Value

A ggplot object representing the density functions.

Note

Please check for inconsistent curves, especially those showing implausible shapes such as violations of biuniqueness in the cnorm models.

See Also

plotNormCurves, plotPercentiles

Other plot: compare(), plot.cnorm(), plot.cnormBetaBinomial(), plot.cnormBetaBinomial2(), plotDerivative(), plotNorm(), plotNormCurves(), plotPercentileSeries(), plotPercentiles(), plotRaw(), plotSubset()

Examples

## Not run: 
# For traditional continuous norming model
result <- cnorm(raw = elfe$raw, group = elfe$group)
plotDensity(result, group = c(2, 4, 6))

# For beta-binomial model
bb_model <- cnorm.betabinomial(age = ppvt$age, score = ppvt$raw, n = 228)
plotDensity(bb_model)

## End(Not run)

Description

This function plots the scores obtained via the first order derivative of the regression model in dependence of the norm score.

Usage

plotDerivative(
  model,
  minAge = NULL,
  maxAge = NULL,
  minNorm = NULL,
  maxNorm = NULL,
  stepAge = NULL,
  stepNorm = NULL,
  order = 1
)

Arguments

Details

The results indicate the progression of the norm scores within each age group. The regression-based modeling approach relies on the assumption of a linear progression of the norm scores. Negative scores in the first order derivative indicate a violation of this assumption. Scores near zero are typical for bottom and ceiling effects in the raw data.

The regression models usually converge within the range of the original values. In case of vertical and horizontal extrapolation, with increasing distance to the original data, the risk of assumption violation increases as well.

Value

A ggplot object representing the derivative of the regression function.

Note

This function is currently incompatible with reversed raw score scales ('descent' option).

See Also

checkConsistency, bestModel, derive

Other plot: compare(), plot.cnorm(), plot.cnormBetaBinomial(), plot.cnormBetaBinomial2(), plotDensity(), plotNorm(), plotNormCurves(), plotPercentileSeries(), plotPercentiles(), plotRaw(), plotSubset()

Examples

## Not run: 
# For traditional continuous norming model
result <- cnorm(raw = elfe$raw, group = elfe$group)
plotDerivative(result, minAge=2, maxAge=5, stepAge=.2, minNorm=25, maxNorm=75, stepNorm=1)

## End(Not run)

Description

This function plots the manifest norm score against the fitted norm score from the inverse regression model per group. This helps to inspect the precision of the modeling process. The scores should not deviate too far from the regression line. Applicable for Taylor polynomial, beta-binomial, and shash models.

Usage

plotNorm(
  model,
  age = NULL,
  score = NULL,
  width = NULL,
  weights = NULL,
  group = FALSE,
  minNorm = NULL,
  maxNorm = NULL,
  type = 0
)

Arguments

Value

A ggplot object representing the norm scores plot.

See Also

Other plot: compare(), plot.cnorm(), plot.cnormBetaBinomial(), plot.cnormBetaBinomial2(), plotDensity(), plotDerivative(), plotNormCurves(), plotPercentileSeries(), plotPercentiles(), plotRaw(), plotSubset()

Examples

## Not run: 
# Load example data set, compute model and plot results

# Taylor polynomial model
model <- cnorm(raw = elfe$raw, group = elfe$group)
plot(model, "norm")

# Beta binomial models; maximum number of items in elfe is n = 28
model.bb <- cnorm.betabinomial(elfe$group, elfe$raw, n = 28)
plotNorm(model.bb, age = elfe$group, score = elfe$raw)

## End(Not run)

Description

This function plots the norm curves based on the regression model. It supports both Taylor polynomial models, beta-binomial models, and shash models.

Usage

plotNormCurves(
  model,
  normList = NULL,
  minAge = NULL,
  maxAge = NULL,
  step = 0.1,
  minRaw = NULL,
  maxRaw = NULL
)

Arguments

Details

Please check the function for inconsistent curves: The different curves should not intersect. Violations of this assumption are a strong indication of violations of model assumptions in modeling the relationship between raw and norm scores.

Common reasons for inconsistencies include: 1. Vertical extrapolation: Choosing extreme norm scores (e.g., scores <= -3 or >= 3). 2. Horizontal extrapolation: Using the model scores outside the original dataset. 3. The data cannot be modeled with the current approach, or you need another power parameter (k) or R2 for the model.

Value

A ggplot object representing the norm curves.

See Also

checkConsistency, plotDerivative, plotPercentiles

Other plot: compare(), plot.cnorm(), plot.cnormBetaBinomial(), plot.cnormBetaBinomial2(), plotDensity(), plotDerivative(), plotNorm(), plotPercentileSeries(), plotPercentiles(), plotRaw(), plotSubset()

Examples

## Not run: 
# For Taylor continuous norming model
m <- cnorm(raw = ppvt$raw, group = ppvt$group)
plotNormCurves(m, minAge=2, maxAge=5)

# For beta-binomial model
bb_model <- cnorm.betabinomial(age = ppvt$age, score = ppvt$raw, n = 228)
plotNormCurves(bb_model)

## End(Not run)

Description

The function plots the norm curves based on the regression model against the actual percentiles from the raw data. As in 'plotNormCurves', please check for inconsistent curves, especially intersections. Violations of this assumption are a strong indication for problems in modeling the relationship between raw and norm scores. In general, extrapolation (point 1 and 2) can carefully be done to a certain degree outside the original sample, but it should in general be handled with caution. The original percentiles are displayed as distinct points in the according color, the model based projection of percentiles are drawn as lines. Please note, that the estimation of the percentiles of the raw data is done with the quantile function with the default settings. In case, you get 'jagged' or disorganized percentile curve, try to reduce the 'k' and/or 't' parameter in modeling.

Usage

plotPercentiles(
  model,
  minRaw = NULL,
  maxRaw = NULL,
  minAge = NULL,
  maxAge = NULL,
  raw = NULL,
  group = NULL,
  percentiles = c(0.025, 0.1, 0.25, 0.5, 0.75, 0.9, 0.975),
  scale = NULL,
  title = NULL,
  subtitle = NULL,
  points = FALSE
)

Arguments

See Also

plotNormCurves, plotPercentileSeries

Other plot: compare(), plot.cnorm(), plot.cnormBetaBinomial(), plot.cnormBetaBinomial2(), plotDensity(), plotDerivative(), plotNorm(), plotNormCurves(), plotPercentileSeries(), plotRaw(), plotSubset()

Examples

## Not run: 
  # Load example data set, compute model and plot results
  result <- cnorm(raw = elfe$raw, group = elfe$group)
  plotPercentiles(result)

## End(Not run)

Description

This functions makes use of 'plotPercentiles' to generate a series of plots with different number of predictors. It draws on the information provided by the model object to determine the bounds of the modeling (age and standard score range). It can be used as an additional model check to determine the best fitting model. Please have a look at the ' plotPercentiles' function for further information.

Usage

plotPercentileSeries(
  model,
  start = 1,
  end = NULL,
  group = NULL,
  percentiles = c(0.025, 0.1, 0.25, 0.5, 0.75, 0.9, 0.975),
  filename = NULL
)

Arguments

Value

the complete list of plots

See Also

plotPercentiles

Other plot: compare(), plot.cnorm(), plot.cnormBetaBinomial(), plot.cnormBetaBinomial2(), plotDensity(), plotDerivative(), plotNorm(), plotNormCurves(), plotPercentiles(), plotRaw(), plotSubset()

Examples

## Not run: 
  # Load example data set, compute model and plot results
  result <- cnorm(raw = elfe$raw, group = elfe$group)
  plotPercentileSeries(result, start=4, end=6)

## End(Not run)

Description

The function plots the raw data against the fitted scores from the regression model per group. This helps to inspect the precision of the modeling process. The scores should not deviate too far from regression line.

Usage

plotRaw(model, group = FALSE, type = 0)

Arguments

See Also

Other plot: compare(), plot.cnorm(), plot.cnormBetaBinomial(), plot.cnormBetaBinomial2(), plotDensity(), plotDerivative(), plotNorm(), plotNormCurves(), plotPercentileSeries(), plotPercentiles(), plotSubset()

Examples

## Not run: 
  # Compute model with example dataset and plot results
  result <- cnorm(raw = elfe$raw, group = elfe$group)
  plotRaw(result)

## End(Not run)

Description

This function plots various information criteria and model fit statistics against the number of predictors or adjusted R-squared, depending on the type of plot selected. It helps in model selection by visualizing different aspects of model performance. Models, which did not pass the initial consistency check are depicted with an empty circle.

Usage

plotSubset(model, type = 0)

Arguments

Details

The function generates different plots to help in model selection:

- For types 1 and 2 (Mallow's Cp and BIC), look for the "elbow" in the curve where the information criterion begins to drop. This often indicates a good balance between model fit and complexity. - For type 0 (Adjusted R2), higher values indicate better fit, but be cautious of overfitting with values approaching 1. - For types 3 and 4 (RMSE and RSS), lower values indicate better fit. - For type 5 (F-test), higher values suggest significant improvement with added predictors. - For type 6 (p-values), values below the significance level (typically 0.05) suggest significant improvement with added predictors.

Value

A ggplot object representing the selected information criterion plot.

Note

It's important to balance statistical measures with practical considerations and to visually inspect the model fit using functions like plotPercentiles.

See Also

bestModel, plotPercentiles, printSubset

Other plot: compare(), plot.cnorm(), plot.cnormBetaBinomial(), plot.cnormBetaBinomial2(), plotDensity(), plotDerivative(), plotNorm(), plotNormCurves(), plotPercentileSeries(), plotPercentiles(), plotRaw()

Examples

## Not run: 
# Compute model with example data and plot information function
cnorm.model <- cnorm(raw = elfe$raw, group = elfe$group)
plotSubset(cnorm.model)

# Plot BIC against adjusted R-squared
plotSubset(cnorm.model, type = 2)

# Plot RMSE against number of predictors
plotSubset(cnorm.model, type = 3)

## End(Not run)

Description

A dataset based on an unstratified sample of PPVT4 data (German adaption). The PPVT4 consists of blocks of items with 12 items each. Each item consists of 4 pictures. The test taker is given a word orally and he or she has to point out the picture matching the oral word. Bottom and ceiling blocks of items are determined according to age and performance. For instance, when a student knows less than 4 word from a block of 12 items, the testing stops. The sample is not identical with the norm sample and includes doublets of cases in order to align the sample size per age group. It is primarily intended for running the cNORM analyses with regard to modeling and stratification.

Usage

ppvt

Format

A data frame with 4542 rows and 6 variables:

age

the chronological age of the child

sex

the sex of the test taker, 1=male, 2=female

migration

migration status of the family, 0=no, 1=yes

region

factor specifying the region, the data were collected; grouped into south, north, east and west

raw

the raw score of the student, spanning values from 0 to 228

group

age group of the child, determined by the getGroups()-function with 12 equidistant age groups

A data frame with 5600 rows and 9 columns

Source

https://www.psychometrica.de/ppvt4.html

References

Lenhard, A., Lenhard, W., Segerer, R. & Suggate, S. (2015). Peabody Picture Vocabulary Test - Revision IV (Deutsche Adaption). Frankfurt a. M./Germany: Pearson Assessment.

Examples

## Not run: 
# Example with continuous age variable, ranked with sliding window
model.ppvt.sliding <- cnorm(age=ppvt$age, raw=ppvt$raw, width=1)

# Example with age groups; you might first want to experiment with
# the granularity of the groups via the 'getGroups()' function
model.ppvt.group <- cnorm(group=ppvt$group, raw=ppvt$raw) # with predefined groups
model.ppvt.group <- cnorm(group=getGroups(ppvt$age, n=15, equidistant = T),
                          raw=ppvt$raw) # groups built 'on the fly'


# plot information function
plot(model.ppvt.group, "subset")

# check model consistency
checkConsistency(model.ppvt.group)

# plot percentiles
plot(model.ppvt.group, "percentiles")

## End(Not run)

Description

This function calculates norm scores based on raw scores, age, and a fitted cnormBetaBinomial model.

Usage

## S3 method for class 'cnormBetaBinomial'
predict(object, ...)

Arguments

Details

The function first predicts the alpha and beta parameters of the beta-binomial distribution for each age using the provided model. It then calculates the cumulative probability for each raw score given these parameters. Finally, it converts these probabilities to the norm scale specified in the model.

Value

A numeric vector of norm scores.

See Also

Other predict: derivationTable(), getNormCurve(), normTable(), predict.cnormBetaBinomial2(), predict.cnormShash(), predictNorm(), predictRaw(), rawTable()

Examples

## Not run: 
# Assuming you have a fitted model named 'bb_model':
model <- cnorm.betabinomial(ppvt$age, ppvt$raw)
raw <- c(100, 121, 97, 180)
ages <- c(7, 8, 9, 10)
norm_scores <- predict(model, ages, raw)

## End(Not run)

Description

This function calculates norm scores based on raw scores, age, and a fitted cnormBetaBinomial model.

Usage

## S3 method for class 'cnormBetaBinomial2'
predict(object, ...)

Arguments

Details

The function first predicts the alpha and beta parameters of the beta-binomial distribution for each age using the provided model. It then calculates the cumulative probability for each raw score given these parameters. Finally, it converts these probabilities to the norm scale specified in the model.

Value

A numeric vector of norm scores.

See Also

Other predict: derivationTable(), getNormCurve(), normTable(), predict.cnormBetaBinomial(), predict.cnormShash(), predictNorm(), predictRaw(), rawTable()

Examples

## Not run: 
# Assuming you have a fitted model named 'bb_model':
model <- cnorm.betabinomial(ppvt$age, ppvt$raw)
raw <- c(100, 121, 97, 180)
ages <- c(7, 8, 9, 10)
norm_scores <- predict(model, ages, raw)

## End(Not run)

Description

This function calculates norm scores based on raw scores, age, and a fitted cnormShash model.

Usage

## S3 method for class 'cnormShash'
predict(object, ...)

Arguments

Details

The function predicts the SinH-ArcSinH (shash) distribution parameters (mu, sigma, epsilon, delta) for each age using the provided model. It then calculates the cumulative probability for each raw score given these parameters using the continuous shash distribution. Finally, it converts these probabilities to the norm scale specified in the model.

Value

A numeric vector of norm scores.

See Also

Other predict: derivationTable(), getNormCurve(), normTable(), predict.cnormBetaBinomial(), predict.cnormBetaBinomial2(), predictNorm(), predictRaw(), rawTable()

Examples

## Not run: 
# Assuming you have a fitted model named 'shash_model':
model <- cnorm.shash(children$age, children$score)
raw_scores <- c(25.5, 30.2, 18.7, 45.3)
ages <- c(7, 8, 9, 10)
norm_scores <- predict(model, ages, raw_scores)

## End(Not run)

Description

This functions numerically determines the norm score for raw scores depending on the level of the explanatory variable A, e. g. norm scores for raw scores at given ages.

Usage

predictNorm(
  raw,
  A,
  model,
  minNorm = NULL,
  maxNorm = NULL,
  force = TRUE,
  silent = FALSE
)

Arguments

Value

The predicted norm score for a raw score, either single value or vector

See Also

Other predict: derivationTable(), getNormCurve(), normTable(), predict.cnormBetaBinomial(), predict.cnormBetaBinomial2(), predict.cnormShash(), predictRaw(), rawTable()

Examples

## Not run: 
# Generate cnorm object from example data
cnorm.elfe <- cnorm(raw = elfe$raw, group = elfe$group)

# return norm value for raw value 21 for grade 2, month 9
specificNormValue <- predictNorm(raw = 21, A = 2.75, cnorm.elfe)

## End(Not run)

Description

Most elementary function to predict raw score based on Location (L, T score), Age (grouping variable) and the coefficients from a regression model.

Usage

predictRaw(norm, age, coefficients, minRaw = -Inf, maxRaw = Inf)

Arguments

Value

the predicted raw score or a data.frame of scores in case, lists of norm scores or age is used

See Also

Other predict: derivationTable(), getNormCurve(), normTable(), predict.cnormBetaBinomial(), predict.cnormBetaBinomial2(), predict.cnormShash(), predictNorm(), rawTable()

Examples

## Not run: 
# Prediction of single scores
model <- cnorm(raw = elfe$raw, group = elfe$group)
predictRaw(35, 3.5, model)

## End(Not run)

Description

This is a convenience method to either load the inbuilt sample dataset, or to provide a data frame with the variables "raw" (for the raw scores) and "group" The function ranks the data within groups, computes norm values, powers of the norm scores and interactions. Afterwards, you can use these preprocessed data to determine the best fitting model.

Usage

prepareData(
  data = NULL,
  group = "group",
  raw = "raw",
  age = "group",
  k = 4,
  t = NULL,
  width = NA,
  weights = NULL,
  scale = "T",
  descend = FALSE,
  silent = FALSE
)

Arguments

Details

The functions rankBySlidingWindow, rankByGroup, bestModel, computePowers and prepareData are usually not called directly, but accessed through other functions like cnorm.

Value

data frame including the norm scores, powers and interactions of the norm score and grouping variable

See Also

Other prepare: computePowers(), rankByGroup(), rankBySlidingWindow()

Examples

## Not run: 
# conducts ranking and computation of powers and interactions with the 'elfe' dataset
data.elfe <- prepareData(elfe)

# use vectors instead of data frame
data.elfe <- prepareData(raw=elfe$raw, group=elfe$group)

# variable names can be specified as well, here with the BMI data included in the package
data.bmi <- prepareData(CDC, group = "group", raw = "bmi", age = "age")

## End(Not run)

# modeling with only one group with the 'elfe' dataset as an example
# this results in conventional norming
data.elfe2 <- prepareData(data = elfe, group = FALSE)
m <- bestModel(data.elfe2)

Description

After conducting the model fitting procedure on the data set, the best fitting model has to be chosen. The print function shows the R2 and other information on the different best fitting models with increasing number of predictors.

Usage

## S3 method for class 'cnorm'
print(x, ...)

Arguments

Value

A table with information criteria

See Also

Other model: bestModel(), checkConsistency(), cnorm.cv(), derive(), modelSummary(), printSubset(), rangeCheck(), regressionFunction(), summary.cnorm()

Description

Print method for SinH-ArcSinH objects

Usage

## S3 method for class 'cnormShash'
print(x, ...)

Arguments

Description

Displays R^2 and other metrics for models with varying predictors, aiding in choosing the best-fitting model after model fitting.

Usage

printSubset(x, ...)

Arguments

Value

Table with model information criteria.

See Also

Other model: bestModel(), checkConsistency(), cnorm.cv(), derive(), modelSummary(), print.cnorm(), rangeCheck(), regressionFunction(), summary.cnorm()

Examples

## Not run: 
  # Using cnorm object from sample data
  result <- cnorm(raw = elfe$raw, group = elfe$group)
  printSubset(result)

## End(Not run)

Description

Regression model only work in a specific range and extrapolation horizontally (outside the original range) or vertically (extreme norm scores) might lead to inconsistent results. The function generates a message, indicating extrapolation and the range of the original data.

Usage

rangeCheck(
  object,
  minAge = NULL,
  maxAge = NULL,
  minNorm = NULL,
  maxNorm = NULL,
  digits = 3,
  ...
)

Arguments

Value

the report

See Also

Other model: bestModel(), checkConsistency(), cnorm.cv(), derive(), modelSummary(), print.cnorm(), printSubset(), regressionFunction(), summary.cnorm()

Examples

## Not run: 
  m <- cnorm(raw = elfe$raw, group = elfe$group)
  rangeCheck(m)

## End(Not run)

Description

This is the initial step, usually done in all kinds of test norming projects, after the scale is constructed and the norm sample is established. First, the data is grouped according to a grouping variable and afterwards, the percentile for each raw value is retrieved. The percentile can be used for the modeling procedure, but in case, the samples to not deviate too much from normality, T, IQ or z scores can be computed via a normal rank procedure based on the inverse cumulative normal distribution. In case of bindings, we use the medium rank and there are different methods for estimating the percentiles (default RankIt).

Usage

rankByGroup(
  data = NULL,
  group = "group",
  raw = "raw",
  weights = NULL,
  method = 4,
  scale = "T",
  descend = FALSE,
  descriptives = TRUE,
  na.rm = TRUE,
  silent = FALSE
)

Arguments

Value

the dataset with the percentiles and norm scales per group

Remarks on using covariates

So far the inclusion of a binary covariate is experimental and far from optimized. The according variable name has to be specified in the ranking procedure and the modeling includes this in the further process. At the moment, during ranking the data are split into the according cells group x covariate, which leads to small sample sizes. Please take care to have enough cases in each combination. Additionally, covariates can lead to unstable modeling solutions. The question, if it is really reasonable to include covariates when norming a test is a decision beyond the pure data modeling. Please use with care or alternatively split the dataset into the two groups beforehand and model them separately.

The functions rankBySlidingWindow, rankByGroup, bestModel, computePowers and prepareData are usually not called directly, but accessed through other functions like cnorm.

See Also

rankBySlidingWindow, computePowers, computeWeights, weighted.rank

Other prepare: computePowers(), prepareData(), rankBySlidingWindow()

Examples

## Not run: 
# Transformation with default parameters: RankIt and converting to T scores
data.elfe <- rankByGroup(elfe, group = "group") # using a data frame with vector names
data.elfe2 <- rankByGroup(raw=elfe$raw, group=elfe$group) # use vectors for raw score and group

# Transformation into Wechsler scores with Yu & Huang (2001) ranking procedure
data.elfe <- rankByGroup(raw = elfe$raw, group = elfe$group, method = 7, scale = c(10, 3))

# cNORM can as well be used for conventional norming, in case no group is given
d <- rankByGroup(raw = elfe$raw)
d <- computePowers(d)
m <- bestModel(d)
rawTable(0, m) # please use an arbitrary value for age when generating the tables

## End(Not run)

Description

The function retrieves all individuals in the predefined age range (x +/- width/2) around each case and ranks that individual based on this individually drawn sample. This function can be directly used with a continuous age variable in order to avoid grouping. When collecting data on the basis of a continuous age variable, cases located far from the mean age of the group receive distorted percentiles when building discrete groups and generating percentiles with the traditional approach. The distortion increases with distance from the group mean and this effect can be avoided by the sliding window. Nonetheless, please ensure, that the optional grouping variable in fact represents the correct mean age of the respective age groups, as this variable is later on used for displaying the manifest data in the percentile plots.

Usage

rankBySlidingWindow(
  data = NULL,
  age = "age",
  raw = "raw",
  weights = NULL,
  width,
  method = 4,
  scale = "T",
  descend = FALSE,
  descriptives = TRUE,
  nGroup = 0,
  group = NA,
  na.rm = TRUE,
  silent = FALSE
)

Arguments

Details

In case of bindings, the function uses the medium rank and applies the algorithms already described in the rankByGroup function. At the upper and lower end of the data sample, the sliding stops and the sample is drawn from the interval min + width and max - width, respectively.

Value

the dataset with the individual percentiles and norm scores

Remarks on using covariates

So far the inclusion of a binary covariate is experimental and far from optimized. The according variable name has to be specified in the ranking procedure and the modeling includes this in the further process. At the moment, during ranking the data are split into the according degrees of the covariate and the ranking is done separately. This may lead to small sample sizes. Please take care to have enough cases in each combination. Additionally, covariates can lead to unstable modeling solutions. The question, if it is really reasonable to include covariates when norming a test is a decision beyond the pure data modeling. Please use with care or alternatively split the dataset into the two groups beforehand and model them separately.

The functions rankBySlidingWindow, rankByGroup, bestModel, computePowers and prepareData are usually not called directly, but accessed through other functions like cnorm.

See Also

rankByGroup, computePowers, computeWeights, weighted.rank, weighted.quantile

Other prepare: computePowers(), prepareData(), rankByGroup()

Examples

## Not run: 
# Transformation using a sliding window
data.elfe2 <- rankBySlidingWindow(elfe, raw = "raw", age = "group", width = 0.5)

# Comparing this to the traditional approach should give us exactly the same
# values, since the sample dataset only has a grouping variable for age
data.elfe <- rankByGroup(elfe, group = "group")
mean(data.elfe$normValue - data.elfe2$normValue)

## End(Not run)

Description

This function is comparable to 'normTable', despite it reverses the assignment: A table with raw scores and the according norm scores for a specific age based on the regression model is generated. This way, the inverse function of the regression model is solved numerically with brute force. Please specify the range of raw values, you want to cover. With higher precision and smaller stepping, this function becomes computational intensive. In case a confidence coefficient (CI, default .9) and the reliability is specified, confidence intervals are computed for the true score estimates, including a correction for regression to the mean (Eid & Schmidt, 2012, p. 272).

Usage

rawTable(
  A,
  model,
  minRaw = NULL,
  maxRaw = NULL,
  minNorm = NULL,
  maxNorm = NULL,
  step = 1,
  monotonuous = TRUE,
  CI = 0.9,
  reliability = NULL,
  pretty = TRUE
)

Arguments

Value

either data.frame with raw scores and the predicted norm scores in case of simple A value or a list of norm tables if vector of A values was provided

References

Eid, M. & Schmidt, K. (2012). Testtheorie und Testkonstruktion. Hogrefe.

See Also

normTable

Other predict: derivationTable(), getNormCurve(), normTable(), predict.cnormBetaBinomial(), predict.cnormBetaBinomial2(), predict.cnormShash(), predictNorm(), predictRaw()

Examples

## Not run: 
# Generate cnorm object from example data
cnorm.elfe <- cnorm(raw = elfe$raw, group = elfe$group)
# generate a norm table for the raw value range from 0 to 28 for the time point month 7 of grade 3
table <- rawTable(3 + 7 / 12, cnorm.elfe, minRaw = 0, maxRaw = 28)

# generate several raw tables
table <- rawTable(c(2.5, 3.5, 4.5), cnorm.elfe, minRaw = 0, maxRaw = 28)

# additionally compute confidence intervals
table <- rawTable(c(2.5, 3.5, 4.5), cnorm.elfe, minRaw = 0, maxRaw = 28, CI = .9, reliability = .94)

# conventional norming, set age to arbitrary value
model <- cnorm(raw = elfe$raw)
rawTable(0, model)

## End(Not run)

Description

The method builds the regression function for the regression model, including the beta weights. It can be used to predict the raw scores based on age and location.

Usage

regressionFunction(model, raw = NULL, digits = NULL)

Arguments

Value

The regression formula as a string

See Also

Other model: bestModel(), checkConsistency(), cnorm.cv(), derive(), modelSummary(), print.cnorm(), printSubset(), rangeCheck(), summary.cnorm()

Examples

## Not run: 
  result <- cnorm(raw = elfe$raw, group = elfe$group)
  regressionFunction(result)

## End(Not run)

Description

Density, distribution function, quantile function and random generation for the Sinh-Arcsinh distribution with location parameter mu, scale parameter sigma, skewness parameter epsilon, and tail weight parameter delta.

Usage

dshash(x, mu = 0, sigma = 1, epsilon = 0, delta = 1, log = FALSE)

pshash(
  q,
  mu = 0,
  sigma = 1,
  epsilon = 0,
  delta = 1,
  lower.tail = TRUE,
  log.p = FALSE
)

qshash(
  p,
  mu = 0,
  sigma = 1,
  epsilon = 0,
  delta = 1,
  lower.tail = TRUE,
  log.p = FALSE
)

rshash(n, mu = 0, sigma = 1, epsilon = 0, delta = 1)

Arguments

Details

The Sinh-Arcsinh distribution (Jones & Pewsey, 2009) is defined by the transformation:

$X = \mu + \sigma \cdot \sinh\left(\frac{\text{asinh}(Z) - \epsilon}{\delta}\right)$

where $Z \sim N(0,1)$ is a standard normal variable.

The four parameters control:

• mu: Location (similar to mean)

• sigma: Scale (similar to standard deviation)

• epsilon: Skewness (epsilon = 0 gives symmetry)

• delta: Tail weight (delta = 1 gives normal-like tails, delta > 1 gives heavier tails, delta < 1 gives lighter tails)

Value

dshash gives the density, pshash gives the distribution function, qshash gives the quantile function, and rshash generates random deviates.

The length of the result is determined by n for rshash, and is the maximum of the lengths of the numerical arguments for the other functions.

References

Jones, M. C., & Pewsey, A. (2009). Sinh-arcsinh distributions. Biometrika, 96(4), 761-780. doi:10.1093/biomet/asp053

See Also

Normal for the normal distribution.

Examples

## Not run: 
# Generate random samples
x <- rshash(1000, mu = 0, sigma = 1, epsilon = 0.5, delta = 1.2)

# Density
plot(density(x))
curve(dshash(x, mu = 0, sigma = 1, epsilon = 0.5, delta = 1.2),
      add = TRUE, col = "red")

# Cumulative probability
pshash(0, mu = 0, sigma = 1, epsilon = 0.5, delta = 1.2)

# Quantiles
qshash(c(0.025, 0.5, 0.975), mu = 0, sigma = 1, epsilon = 0.5, delta = 1.2)

# Compare with normal distribution (epsilon = 0, delta = 1)
par(mfrow = c(2, 2))
x_vals <- seq(-4, 4, length.out = 200)
plot(x_vals, dshash(x_vals), type = "l", main = "Symmetric (like normal)")
plot(x_vals, dshash(x_vals, epsilon = 1), type = "l", main = "Right skewed")
plot(x_vals, dshash(x_vals, delta = 2), type = "l", main = "Heavy tails")
plot(x_vals, dshash(x_vals, delta = 0.5), type = "l", main = "Light tails")

## End(Not run)

Description

Simulate mean per age

Usage

simMean(age)

Arguments

Value

return predicted means

Examples

## Not run: 
x <- simMean(a)

## End(Not run)

Description

Simulate sd per age

Usage

simSD(age)

Arguments

Value

return predicted sd

Examples

## Not run: 
x <- simSD(a)

## End(Not run)

Description

For testing purposes only: The function simulates raw test scores based on a virtual Rasch based test with n results per age group, an evenly distributed age variable, items.n test items with a simulated difficulty and standard deviation. The development trajectories over age group are modeled by a curve linear function of age, with at first fast progression, which slows down over age, and a slightly increasing standard deviation in order to model a scissor effects. The item difficulties can be accessed via $theta and the raw data via $data of the returned object.

Usage

simulateRasch(
  data = NULL,
  n = 100,
  minAge = 1,
  maxAge = 7,
  items.n = 21,
  items.m = 0,
  items.sd = 1,
  Theta = "random",
  width = 1
)

Arguments

Value

a list containing the simulated data and thetas

data

the data.frame with only age, group and raw

sim

the complete simulated data with item level results

theta

the difficulty of the items

Examples

## Not run: 
# simulate data for a rather easy test (m = -1.0)
sim <- simulateRasch(n=150, minAge=1,
                     maxAge=7, items.n = 30, items.m = -1.0,
                     items.sd = 1, Theta = "random", width = 1.0)

# Show item difficulties
mean(sim$theta)
sd(sim$theta)
hist(sim$theta)

# Plot raw scores
boxplot(raw~group, data=sim$data)

# Model data
data <- prepareData(sim$data, age="age")
model <- bestModel(data, k = 4)
printSubset(model)
plotSubset(model, type=0)

## End(Not run)

Description

This function standardizes a numeric vector by subtracting the mean and dividing by the standard deviation. The resulting vector will have a mean of 0 and a standard deviation of 1.

Usage

standardize(x)

Arguments

Value

A numeric vector of the same length as x, containing the standardized values.

Examples

## Not run: 
data <- c(1, 2, 3, 4, 5)
standardized_data <- standardize(data)
print(standardized_data)

## End(Not run)

Description

Function for standardizing raking weights Raking weights get divided by the smallest weight. Thereby, all weights become larger or equal to 1 without changing the ratio of the weights to each other.

Usage

standardizeRakingWeights(weights)

Arguments

Value

the standardized weights

Description

Performs k-fold resampling to estimate averaged coefficients for linear regression. The coefficients are averaged across k different subsets of the data to provide more stable estimates. For small samples (n < 100), returns a standard linear model instead.

Usage

subsample_lm(text, data, weights, k = 10)

Arguments

Details

The function splits the data into k subsets, fits a linear model on k-1 subsets, and stores the coefficients. This process is repeated k times, and the final coefficients are averaged across all iterations to provide more stable estimates.

Value

An object of class 'lm' with averaged coefficients from k-fold resampling. For small samples, returns a standard lm object.

Description

S3 method for printing the results and regression function of a cnorm model

Usage

## S3 method for class 'cnorm'
summary(object, ...)

Arguments

Value

A report on the regression function, weights, R2 and RMSE

See Also

Other model: bestModel(), checkConsistency(), cnorm.cv(), derive(), modelSummary(), print.cnorm(), printSubset(), rangeCheck(), regressionFunction()

Description

This function provides a summary of a fitted beta-binomial continuous norming model, including model fit statistics, convergence information, and parameter estimates.

Usage

## S3 method for class 'cnormBetaBinomial'
summary(object, ...)

Arguments

Details

The summary includes:

• Basic model information (type, number of observations, number of parameters)

• Model fit statistics (log-likelihood, AIC, BIC)

• R-squared, RMSE, and bias (if age and raw scores are provided) in comparison to manifest norm scores

• Convergence information

• Parameter estimates with standard errors, z-values, and p-values

Value

Invisibly returns a list containing detailed diagnostic information about the model. The function primarily produces printed output summarizing the model.

See Also

cnorm.betabinomial, diagnostics.betabinomial

Examples

## Not run: 
model <- cnorm.betabinomial(ppvt$age, ppvt$raw, n = 228)
summary(model)

# Including R-squared, RMSE, and bias in the summary:
summary(model, age = ppvt$age, score = ppvt$raw)

## End(Not run)

Description

This function provides a summary of a fitted beta-binomial continuous norming model, including model fit statistics, convergence information, and parameter estimates.

Usage

## S3 method for class 'cnormBetaBinomial2'
summary(object, ...)

Arguments

Details

The summary includes:

• Basic model information (type, number of observations, number of parameters)

• Model fit statistics (log-likelihood, AIC, BIC)

• R-squared, RMSE, and bias (if age and raw scores are provided) in comparison to manifest norm scores

• Convergence information

• Parameter estimates with standard errors, z-values, and p-values

Value

Invisibly returns a list containing detailed diagnostic information about the model. The function primarily produces printed output summarizing the model.

See Also

cnorm.betabinomial, diagnostics.betabinomial

Examples

## Not run: 
model <- cnorm.betabinomial(ppvt$age, ppvt$raw, n = 228)
summary(model)

# Including R-squared, RMSE, and bias in the summary:
summary(model, age = ppvt$age, raw = ppvt$raw)

## End(Not run)

Description

This function provides a summary of a fitted SinH-ArcSinH (shash) continuous norming model, including model fit statistics, convergence information, and parameter estimates.

Usage

## S3 method for class 'cnormShash'
summary(object, ...)

Arguments

Details

The summary includes:

• Basic model information (polynomial degrees, delta parameter, number of observations)

• Model fit statistics (log-likelihood, AIC, BIC)

• R-squared, RMSE, and bias (if age and raw scores are provided)

• Convergence information

• Parameter estimates with standard errors, z-values, and p-values

Value

Invisibly returns a list containing detailed diagnostic information about the model. The function primarily produces printed output summarizing the model.

See Also

cnorm.shash, diagnostics.shash

Examples

## Not run: 
model <- cnorm.shash(children$age, children$score)
summary(model)

# Including R-squared, RMSE, and bias in the summary:
summary(model, age = children$age, score = children$score)

## End(Not run)

Description

Conducts distribution free continuous norming and aims to find a fitting model. Raw data are modelled as a Taylor polynomial of powers of age and location and their interactions. In addition to the raw scores, either provide a numeric vector for the grouping information (group) for the ranking of the raw scores. You can adjust the grade of smoothing of the regression model by setting the k, t and terms parameter. In general, increasing k and t leads to a higher fit, while lower values lead to more smoothing. If both parameters are missing, taylorSwift uses k = 5 and t = 3 by default.

Usage

taylorSwift(
  raw = NULL,
  group = NULL,
  age = NULL,
  width = NA,
  weights = NULL,
  scale = "T",
  method = 4,
  descend = FALSE,
  k = NULL,
  t = NULL,
  terms = 0,
  R2 = NULL,
  plot = TRUE,
  extensive = TRUE,
  subsampling = FALSE
)