Stata Technical Bulletin
17
|
age3 |
.01487978 |
.01428464 |
.01388494 |
.01356744 |
.01398739 |
|
age4 |
.01041847 |
.01041711 |
.01041616 |
.01041539 |
.01041641 |
|
age5 |
.00831241 |
.00830772 |
.00830446 |
.0083018 |
.0083053 |
|
age6 |
.00621676 |
.00621788 |
.00621866 |
.00621929 |
.00621846 |
. clear
. use Ieecart
(US Death Rates by Year: 1900-1995)
. Ieecart year age m if (age<=100 & year<=1995)
. return list
matrices :
r(mhat) : IOl x 96
r(kt) : 96 x 1
r(bx) : IOl x 1
r(ax) : 101 x 1
Saved Results
leecart saves the following matrices in r():
|
r(ax) |
bæ kt estimated age-specific death rates by year |
Acknowledgment
The U.S. mortality dataset was obtained from the Berkeley Mortality Database (http://demog.berkeley.edu/wilmoth/mortality)
at the University of California at Berkeley.
References
Lee, R. D. and L. Carter. 1992. Modeling and forecasting the time series of U.S. mortality. Journal of the American Statistical Association 87:
659-671.
Good, I. J. 1969. Some applications of the singular value decomposition of a matrix. Technometrics 11: 823-831.
sg150 Hardy-Weinberg equilibrium test in case-control studies
Jisheng Cui, University of Melbourne, Australia, [email protected]
Introduction
Cleves (1999) proposed the Stata command genhwi for testing the Hardy-Weinberg equilibrium (HWE) of one sample of
individuals. However, in case-control studies, two samples of individuals are collected; the cases and the controls. Usually the
genotypic counts of controls are under HWE. We modified genhwi to enable it to be applicable to data in case-control studies.
The new command genhwcci is given for testing whether the genotypic counts of the cases are under HWE, given the controls
are under HWE. This test is more powerful than that given by genhwi because data from controls are utilized.
Syntax
genhwcci #aa1 #Aa1 #aa1 #AA2 #Aa2 #aa2 [. 1 abel (genotypes) binvar ]
Description
genhwcci is an immediate command used for estimating allele frequency, genotype frequencies, disequilibrium coefficients,
and the associated standard error for codominant traits or data of completely known genotypes in case-control studies. For both
genotypic counts of cases and controls, it performs asymptotic HWE tests. It also tests the HWE for genotypic counts of cases,
under the assumption that the genotypic counts of controls are under HWE; where #aa∣, #Aa1, and #aa1 are the counts for the
AA, Aa and aa of the cases; while #AA2, #Aa2, and #aa2 are the genotypic counts of the controls. This command works for
biallelic loci only.