Class | Analysis and data manipulation command |
Name | asp |
Arguments | <binary trait>|(<quantitative trait> >|>=|<|<=|==|^= <threshold>) |
Prints IBS-based affected (full and half) sib-pair statistics for all marker loci versus the trait. It also prints the mean IBD sharing for full sibs, along with the exact (binomial) two-tailed P-value for the "mean" test. All possible sib-pairs are used, and are treated as independent.
Example:
>> asp AD
------------------------------------------------
IBS ASP for "AD" v. all marker loci
------------------------------------------------
Marker NPairs mIBS ExpIBS Asy P mIBD Asy P
---------- ------ ------ ------ ------ ------ ------
D14S52 10 0.7000 0.6716 0.4145 0.4994 1.0000 ASP .
D14S43 10 1.0000 0.6376 0.0001 0.6343 0.3284 ASP ***
D14S53 7 0.9286 0.7439 0.1680 0.6002 0.6312 ASP .
>> set ple 1
NOTE: Print level 1
>> asp AD
------------------------------------------------
IBS ASP for "AD" v. all marker loci
------------------------------------------------
----------- ASP analysis for "D14S52" --------------
No. of IBS Sharing Mean IBS sharing
Pairs 2/2 1/2 0/2 Obs Exp
Full-sibs 10 4 6 0 0.7000 0.6716
Half-sibs 0 0 0 0 0.0000 0.5227
Expectd IBS Sharing
2/2 1/2 0/2
Full-sibs 0.421 0.500 0.078
Half-sibs 0.202 0.642 0.156
Full-Sib Chi-square (2 df) = 1.8 (P=0.4145)
Mean full-sib IBD sharing =0.4994 (P=1.0000)
[...]
<< (schaid) | Up to index | >> (pen) |