symptoms {randomLCA} | R Documentation |
Allergy and respiratory symptoms for infants 0 to 2 years in six month periods. Outcome is presence or absence of symptom in the six months. Original data was collected at Visits 1-7 over the 2 year period which were summarised to six month periods.
This is data relating to a paper which will hopefully be published shortly. I will be updating this entry then.
Note that these models can be slow to fit, with the "symptoms.lca2random2" model taking about 3 hours. I have a couple more possibilities to speed things up. These are to code the algorithms completely in C (they are currently a mix of C and R), use a Bayesian methodology to prevent the estimates moving towards the boundaries (this can also be described as penalized likelihood) and to look at other optimisation functions.
Thanks to the investigators of the CAPS study for making the data available.
data(symptoms)
A data frame with 444 observations on the following 17 variables.
Nightcough.13
Wheeze.13
Itchyrash.13
FlexDerma.13
Nightcough.45
Wheeze.45
Itchyrash.45
FlexDerma.45
Nightcough.6
Wheeze.6
Itchyrash.6
FlexDerma.6
Nightcough.7
Wheeze.7
Itchyrash.7
FlexDerma.7
Freq
Mihrshai et al (2001)
Mihrshahi, S., Peat, J.K., Webb, K., Tovey, R.E., Marks, G.B., Mellis, C.M. and Leeder S.R. (2001) The Childhood Asthma Prevention Study (CAPS): Design and research protocol of a randomized trial for the primary prevention of asthma. Control led Clinical Trials, 22:333–354.
data(symptoms) ## Not run: symptoms.lca2 <- randomLCA(symptoms[,1:16],freq=symptoms$Freq,nclass=2) symptoms.lca2random <- randomLCA(symptoms[,1:16],freq=symptoms$Freq,initmodel=symptoms.lca2,random=TRUE,nclass=2,blocksize=4) symptoms.lca2random2 <- randomLCA(symptoms[,1:16],freq=symptoms$Freq,initmodel=symptoms.lca2random,random=TRUE,level2=TRUE,nclass=2,blocksize=4) ## End(Not run)