.-
help for ^simuped2^ and ^simuped3^                                   (STB-58: dm82)
.-

Simulating two- and three-generation families
---------------------------------------------

^simuped2^ #Age1 #Std1 #Age2 #Std2 [^,^ ^r^eps^(^#^)^ ^sav^ing^(^filename^)^ ^a^lle^(^#^)^ ^s^ib^(^#^)^]

^simuped3^ #Age1 #Std1 #Age2 #Std2 #Age3 #Std3 [^,^ ^r^eps^(^#^)^
    ^sav^ing^(^filename^)^ ^a^lle^(^#^)^ ^s^ib^(^#^)^ ^si3(^#^)^]


Description
-----------

^simuped2^ and ^simuped3^ are immediate commands used for generating two- and
three-generation family data, respectively. For each person in a family, the
sex is generated by a probability of 0.5. The age of a person is generated
according to a normal distribution, with means #Age1, #Age2 and #Age3 for the
first, second and third generation. The standard deviation of the age is given
by #Std1, #Std2 and #Std3 for the three generation, respectively.

The number of siblings in a generation is a random number, distributed
according to a Poisson distribution. The mean size of the siblings in
the second- and third-generation is given by the argument ^sib(^#^)^
and ^si3(^#^)^, respectively.

The Hardy-Weinberg equilibrium is assumed for the genotypic
distribution of people in the first generation (see
e.g. Elandt-Johnson, 1971). The allele frequency of a biallelic locus
A is given by the argument ^alle(^#^)^. The genotype of a person in
the second- and third-generation is generated according the Mendelian
inheritance, that is, a person inherits the allele A from father (or
mother) with a half probability.

The simulated family data are saved in a file specified by the argument
^saving(^filename^)^, and the number of replications is specified by the
argument ^reps(^#^)^.


Options
-------

^reps(^#^)^ specifies the number of simulated families. The default value is
    100.

^saving(^filename^)^ specifies the file into which the simulated data is
    saved.  The default file name is temp.dta.

^alle(^#^)^ specifies the allele frequency of an biallelic locus A. The
    default value is 0.1.

^sib(^#^)^ spcifies the number of siblings in the second generation. The
    default value is 3.

^si3(^#^)^ spcifies the number of siblings in the third generation. The
    default value is 3.



Examples
--------

 . ^simugen2 70 10 40 10, reps(1000) sav(output) alle(0.05) sib(5)^

 . ^simugen3 80 10 50 10 20 10, reps(2000) alle(0.1) sib(4) si3(3.5)^



Author
------

       Jisheng Cui
       University of Melbourne
       j.cui@@gpph.unimelb.edu.au