# Control of variance component testing

### Description

`varCompTest.control`

provide a list of arguments that control the details of the testing method for testing the nullity of variance components in `varComp.test`

.

### Usage

1 2 3 4 5 6 7 8 9 10 | ```
varCompTest.control(test = "LinScore", LinScore.wt = "InvSTD", LinScore.acc = 1e-08,
LinScore.lim = 1e+06L, LinScore.method = c("AS155", "SSAS155"),
VM03.method = c("SSChiBarSq", "pboot"), VM03.nsim = 10000L,
SS95.method = c("pboot", "pboot"), SS95.nsim = 10000L,
RLRT.method = c("exact", "pboot"), RLRT.nsim = 10000L,
information = 'EI'
)
varCompTests
varCompScoreTests
LinScoreWeightingMethods
``` |

### Arguments

`test` |
A character vector, specifying the tests to be performed. Available choices are given in |

`LinScore.wt` |
A character vector, specifying the weighting method for the linear score test. Available choices are given in |

`LinScore.acc` |
A single small positive number, specifying the accuracy when computing p-values from the linear score test. This is the same as the |

`LinScore.lim` |
A large positive integer, specifying the maximum number of integration terms when computing the |

`LinScore.method` |
A character vector with two elements. The first component specifies the method for obtaining null distributions when the null hypothesis contains no variance components other than the error variance; the second component specifies the method for obtaining null distributions when the null hypothesis contains additional variance parameters other than the error variance. Available choices are |

`VM03.method` |
A character vector with two components for the VM03 test, each component being a character string. The first component specifies the method for obtaining null distributions when the null hypothesis contains no variance components other than the error variance; the second component specifies the method for obtaining null distributions when the null hypothesis contains additional variance parameters other than the error variance. Currently, the 2nd component is discarded. Available choices are |

`VM03.nsim` |
A large positive integer, specifying the number of simulations used to approximate the null, when |

`SS95.method` |
A vector of two components for the SS95 test, each component being a character string. The first component specifies the method for obtaining null distributions when the null hypothesis contains no variance components other than the error variance; the second component specifies the method for obtaining null distributions when the null hypothesis contains additional variance parameters other than the error variance. Currently available chioces are |

`SS95.nsim` |
A large positive integer, specifying the number of simulations used to approximate the null, when |

`RLRT.method` |
A vector of two components for the RLRT test, each component being a character string. The first component specifies the method for obtaining null distributions when the null hypothesis contains no variance components other than the error variance; the second component specifies the method for obtaining null distributions when the null hypothesis contains additional variance parameters other than the error variance. Currently available choices are "exact" (for either the exact null of Crainiceanu and Ruppert [2003] or the pseudo-likelihood heuristic method of Greven et al [2008]) and "pboot" (for the 1st component, this is an alias of "exact"; for the 2nd component, this is not yet implemented). See also |

`RLRT.nsim` |
A large positive integer, specifying the number of simulations used to approximate the null, when |

`information` |
A character scalar, specifying the information matrix to be used in score-based tests. Available choices are given in |

### Value

A list of class `varCompTest.control`

with length equal to the length `test`

argument. This can be passed as the `control`

or `fit.control`

arguments to `varComp.test`

### Author(s)

Long Qu

### References

Qu L, Guennel T, Marshall SL. (2013) Linear Score Tests for Variance Components in Linear Mixed Models and Applications to Genetic Association Studies. *Biometrics*, Volume 69, Issue 4, pages 883–892.

Daniel B. Hall and Jens T. Praestgaard (2001) Order-restricted tests for homogeneity in generalised linear and nonlinear mixed models. *Biometrika* 88, 739–751.

Mervyn J. Silvapulle and Paramsothy Silvapulle (1995) A score test against one-sided alternatives. *Journal of the American Statistical Association* 90, 342–349.

Geert Verbeke and Geert Molenberghs (2003) The use of score tests for inference on variance components. *Biometrics* 59, 254–262.

Geert Verbeke and Geert Molenberghs (2007) What Can Go Wrong With the Score Test?. *The American Statistician* 61, 289–290.

Ciprian M. Crainiceanu and David Ruppert (2003) Likelihood ratio tests in linear mixed models with one variance component. *Journal of the Royal Statistical Society: Series B (Statistical Methodology)* 66, 165–185.

Greven S, Crainiceanu CM Kuchenhoff H and Peters A (2008) Restricted Likelihood Ratio Testing for Zero Variance Components in Linear Mixed Models. *Journal of Computational and Graphical Statistics* 17, 870–891.

Robert B. Davies (1980) The Distribution of a Linear Combination of *χ^2* Random variables. *Journal of the Royal Statistical Society. Series C (Applied Statistics)* 29, 323–333.

### See Also

`varComp.test`

, `CompQuadForm::davies`

, `RLRsim::RLRTSim`