identifyComponentModel {BGmisc}

identifyComponentModel Determine if a variance components model is identified

Description

identifyComponentModel Determine if a variance components model is identified

Usage

identifyComponentModel(..., verbose = TRUE)

Arguments

...

Comma-separated relatedness component matrices representing the variance components of the model.

verbose

logical. If FALSE, suppresses messages about identification; TRUE by default.

Details

This function checks the identification status of a given variance components model by examining the rank of the concatenated matrices of the components. If any components are not identified, their names are returned in the output.

Value

A list of length 2 containing:

• identified: TRUE if the model is identified, FALSE otherwise.

• nidp: A vector of non-identified parameters, specifying the names of components that are not simultaneously identified.

Examples

identifyComponentModel(A = list(matrix(1, 2, 2)), C = list(matrix(1, 2, 2)), E = diag(1, 2))