Use of a random coefficient regression (RCR) model to estimate growth parameters

We used a random coefficient regression (RCR) model to estimate growth parameters for the time series of observed serum glucose levels in the Replicate 1 of the Genetic Analysis Workshop 13 simulated data. For comparison, a two time-point interval was also selected and the slope between these two observations was calculated. This process yielded four phenotypes: the RCR growth phenotype, a two time-point slope phenotype, and Time 1 and Time 2 serum glucose level phenotypes. These four phenotypes were used for linkage analyses on simulated chromosomes 5, 7, 9, and 21, those chromosomes that contained loci affecting the growth course for serum glucose levels. The linkage analysis of the RCR-derived phenotype showed overwhelming evidence for linkage at one locus (LOD 65.78 on chromosome 5), while showing elevated but nonsignificant LOD scores for two other loci (LOD 1.25 on chromosome 7, LOD 1.10 on chromosome 9), and no evidence of linkage for the final locus. The two time-point slope phenotype showed evidence for linkage at one locus (LOD 4.16 on chromosome 5) but no evidence for linkage at any of the other loci. A parallel cross-sectional approach, using as input phenotypes the endpoints of the two-point slope phenotype, gave strong linkage results for the major locus on chromosome 5 (maximal LOD scores of 17.90 and 27.24 for Time 1 and Time 2, respectively) while showing elevated but nonsignificant linkage results on chromosome 7 (maximal LOD scores of 1.71 and 1.48) and no evidence for linkage at the two remaining loci. The RCR growth parameter showed more power to detect linkage to the major locus than either the cross-sectional or two-point slope approach, but the cross-sectional approach gave a higher maximal LOD score for one of the minor loci.