Description Usage Arguments Details Value Author(s) Examples
Perform a likelihood ratio test for two asreml-models that have an equal fixed part and a nested random part.
1 | LR.test.asreml(model0, model1)
|
model0 |
An asreml-object, see |
model1 |
An asreml-object with the same fixed part as |
The likelihood ratio test is performed by comparing two times the difference in log-likelihood with the quantiles if the χ^2-distribution with as degrees of freedom the difference in number of variance components estimated in both models.
If the degrees of freedom equal 1, the obtained p-value is divided by 2.
A list with class "htest
" containing the following
components:
The value of the test statistic.
The degrees of freedom for the test statistic.
The p-value for the test.
NULL
NULL
NULL
Extended model outperforms the 'null model'.
Likelihood ratio test (Chi-square)
Names of the two models that were tested.
Ruud Derijcker, Katrien Baert
1 2 3 4 5 6 7 | library(asreml)
data(oats)
model0 <- asreml(fixed = yield ~ Variety*Nitrogen, random = ~ Blocks:Wplots,
data = oats)
model1 <- asreml(fixed = yield ~ Variety*Nitrogen, random = ~ Blocks/Wplots,
data = oats)
LR.test.asreml(model0, model1)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.