Sib-pair Command: alla

Class	Algebraic operators and function
Name	alla | allb
Arguments	<marker>

Return the first or second allele for each individual at the given marker. If you want to test the allele for a given individual is a character value (eg "A"), the comparison has to be to the homozygous genotype ("A/A"), else the parser will complain there is no locus "A".

Example:

>> count alla D14S53 == 151 or allb D14S53 == 151

Count where "alla D14S53 = = 151 or allb D14S53 = = 151":

Pedigree       Con=T   Num  ASPs Trios    4+
-------------- ----- ----- ----- ----- -----
         Total    19    77     1     3     0

>>  tab D14S53

------------------------------------
Tabulation of "D14S53              "
------------------------------------
  D14S53      Count  Percent
-----------------------------
  144/144        1      5.0
  144/151        6     30.0
  148/151        1      5.0
  149/151        1      5.0
  151/151        6     30.0
  151/155        4     20.0
  151/157        1      5.0
-----------------------------
  Total         20 (and  57  missing)

>> macro mkdummy
>> set loc %1_%2 qua
>> %1_%2 = 0
>> if (alla %1 == %2) then %1_%2 = %1_%2 + 1
>> if (allb %1 == %2) then %1_%2 = %1_%2 + 1
>> if (untyp %1) then %1_%2 = x
>> ;;;;
>> mkdummy D14S53 { 144 148 149 151 155 157 }

>>  count alla rs11260964 == "A/A"
Number of matched persons   = 3067 out of 4699 (65.3%)
Number of matched pedigrees = 1755 out of 2050 (85.6%)
>>  tab rs11260964
rs11260964           x: 94285     A/A:  2027     A/G:  1040     G/G:   480
>> 2027+1040
=> 3067.0000000000000