# 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
# with configurations C1, C3, C4 and C5 for the Winter (1st and 4th) quarter intervals
#ChinaWTE <- mle(ChinaT[,c(1,4)],Config=c(1,3,4,5))
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 4
# against model with covariance Configuration 1 (unrestricted covariance)
testMod(ChinaWTE,RestMod=4,FullMod=1)
Run the code above in your browser using DataLab