Sib-pair Command: test locus


ClassAnalysis and data manipulation command
Nametest locus
Arguments [<marker 1>] [...<marker N>]]

Carry out the usual Sib-pair Mendelian error screen on the specified locus and only on active pedigrees. Useful if checking was turned off when the pedigree was first read in.

Example:

>> set checking off
>> set imputation -1
>> set loc a marker
>> read ped inline
1 1 x x m 3/3
1 2 x x f x/x
1 3 1 2 f 2/4
;;;;
>> run
>> test locus a

NOTE:  inconsistency due child 1-3 at locus a {  2/4  }

Locus "a"
------------------
Sibship: 1-1 x 1-2

Inconsistency between parent and child genotypes.

                        [1]                 (2)    
                        3/3                 x/x  
                         |                   |
                         +=========+=========+
                                   |
                                   |
                                  (3)    
                                  2/4  


Pedigree    Loci with Mendel errors
----------  -----------------------
1                  1   100.0%

Locus       Peds with Mendel errors
----------  -----------------------
a                  1   100.0%                                         

>> wri
Writing 1 pedigrees:

!
!       S        
!       e        
!       x    a   
!
1 1 x x m   3/3  
1 2 x x f   x/x  
1 3 1 2 f   2/4  

>> set error_dropping on
>> test locus a
>> wri
Writing 1 pedigrees:

!
!       S        
!       e        
!       x    a   
!
1 1 x x m   x/x  
1 2 x x f   x/x  
1 3 1 2 f   x/x  


<< (test haploid)Up to index>> (test map)