Sib-pair Command: homoz

Class	Analysis and data manipulation command
Name	homoz
Arguments	[<binary trait> [|<quantitative trait> >|>=|<|<=|==|^= <threshold>]]

Prints the asymptotic Z statistic and a one-sided MC P-value for whether homozygosity at each marker locus is increased in probands, either affected if the trait is binary, or above or below the given threshold if the trait is quantitative.

Example:

#
# Increased homozygosity in cases consistent with a recessive
# disease gene
#
>> include linclex.in
>> homoz lincl

--------------------------------------------------
Marker homozygosity in cases of trait "lincl"
--------------------------------------------------

Marker        N    Obs    Exp    Fis     Z  Asy P  Emp P  Iters
-------- ------ ------ ------ ------ ------ ------ ------ -----
D13S162      16 0.8750 0.4262 0.7822    3.7 0.0002 0.0100   201 HOM *  
D13S160      16 0.8125 0.3262 0.7217    4.3 2e-5   0.0004   201 HOM ** 
D13S170      12 0.7500 0.1885 0.6919    5.2 2e-7   0.0008   201 HOM ** 
--------                      ------
Multilocus                    0.7268