Performs statistical likelihood-ratio tests that evaluate the goodness-of-fit of a nested model against a more general one.
testMod(ModE,RestMod=ModE@ModelConfig[2]:length(ModE@ModelConfig),FullMod="Next")
An object of class IdtE representing the estimates of a model fitted to a data set of interval-value variables
Indices of the restricted models beeing evaluated in the NULL hypothesis
Either indices of the general models beeing evaluated in the alternative hypothesis or the strings "Next" (default) or "All". In the former case a Restricted model is always compared against the most parsimonious alternative that encompasses it, and in latter all possible comparisons are performed
An object of type ConfTests with the results of the tests performed
# Create an Interval-Data object containing the intervals of temperatures by quarter
# for 899 Chinese meteorological stations.
ChinaT <- IData(ChinaTemp[1:8])
# Estimate by maximum likelihood the parameters of Gaussian models
# for the Winter (1st and 4th) quarter intervals
ChinaWTE <- mle(ChinaT[,c(1,4)])
cat("China maximum likelhiood estimation results for Winter quarters:\n")
print(ChinaWTE)
# Perform Likelihood-Ratio tests comparing models with consecutive nested Configuration
testMod(ChinaWTE)
# Perform Likelihood-Ratio tests comparing all possible models
testMod(ChinaWTE,FullMod="All")
# Compare model with covariance Configuration case 3 (MidPoints independent of LogRanges)
# against model with covariance Configuration 1 (unrestricted covariance)
testMod(ChinaWTE,RestMod=3,FullMod=1)
Run the code above in your browser using DataLab