Sib-pair Command: asp


ClassAnalysis and data manipulation command
Nameasp
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)