University of Michigan Center for Statistical 
Genetics
Search
 
 

 
 

PSEUDO -- Family Variance Weights

Basic summaries

Once the main simulation runs, PSEUDO also produces a number of summaries of family variance weights . By default, the program will output two tables for each hypothesis you've tested. The first table will be a listing of maximum and minimum variance weights for each chromosome.


MAXIMUM AND MINIMUM WEIGHTS BY CHROMOSOME
=========================================

Hypothesis  :  Hypothesis_5
Affections  :  trait [QTL]
Target Lod  :  1.81
Target Rank :  2

  CHR     MAX WT     FAMILY  |     MIN WT     FAMILY |    AVERAGE
----------------------------------------------------------------
    1    0.04655        148  |    0.04288        113 |    0.04406
    2    0.04601          9  |    0.04284        107 |    0.04406
    3    0.04536        194  |    0.04282        164 |    0.04402
    4    0.04614         91  |    0.04270         28 |    0.04394
    5    0.04587         67  |    0.04295        184 |    0.04394
    6    0.04562         26  |    0.04282         59 |    0.04399
    7    0.04557         65  |    0.04278         94 |    0.04406
    8    0.04567         94  |    0.04288         45 |    0.04397
    9    0.04637         35  |    0.04291        125 |    0.04403
   10    0.04591         79  |    0.04289        195 |    0.04399
   11    0.04568         12  |    0.04259          3 |    0.04403
   12    0.04676        188  |    0.04262        151 |    0.04407
   13    0.04588        125  |    0.04277         41 |    0.04400
   14    0.04566         97  |    0.04265         10 |    0.04398
   15    0.04606        189  |    0.04284         72 |    0.04406
   16    0.04565         53  |    0.04287        102 |    0.04403
   17    0.04561        132  |    0.04288        193 |    0.04401
   18    0.04565        185  |    0.04276         72 |    0.04399
   19    0.04592         95  |    0.04280        103 |    0.04402
   20    0.04622         49  |    0.04260         85 |    0.04401
   21    0.04587        101  |    0.04264        186 |    0.04405
   22    0.04587        161  |    0.04288         38 |    0.04404
The second presents the same information, on a per-family basis.
MAXIMUM AND MINIMUM WEIGHTS BY FAMILY
=====================================

Hypothesis  :  Hypothesis_0
Affections  :  affection [PAIRS], affection [ALL]
Target Lod  :  3.60
Target Rank :  1

    FAMILY |     AVERAGE  |      MAX WT        CHR  |      MIN WT        CHR
----------------------------------------------------------------------------
         1 |     0.00647  |     0.00669          5  |     0.00635         17 
         2 |     0.00649  |     0.00664         18  |     0.00639         12 
         3 |     0.00648  |     0.00663         15  |     0.00641          5 
         4 |     0.00650  |     0.00669          8  |     0.00639         12 
         5 |     0.00649  |     0.00658          5  |     0.00638          7 
         6 |     0.00647  |     0.00660         20  |     0.00636          4 
         7 |     0.00649  |     0.00663         16  |     0.00639         19 
         8 |     0.00648  |     0.00656          5  |     0.00638         10 
         9 |     0.00646  |     0.00652          4  |     0.00637         16 
        10 |     0.00647  |     0.00660         11  |     0.00637         13 
        11 |     0.00645  |     0.00653         10  |     0.00638         11 
        12 |     0.00648  |     0.00660          9  |     0.00639          4 
        13 |     0.00648  |     0.00661          7  |     0.00640          4 
        14 |     0.00647  |     0.00664         22  |     0.00638         11 
        15 |     0.00649  |     0.00659          4  |     0.00641         21 
        16 |     0.00650  |     0.00662          3  |     0.00638         17 
        17 |     0.00650  |     0.00663          4  |     0.00636         12 
        18 |     0.00647  |     0.00661          5  |     0.00635          2 

Optional weight summaries

If you'd like, you can also tell PSEUDO to produce output files that provide more information on family variance weights. If you're interested only in obtaining the weights for each family and chromosome, try the --weights option. When this option is selected PSEUDO will write the estimated weight for each family and chromosome to the file pseudo.weights. For the example data, PSEUDO produced the following weights output:


VARIANCE WEIGHTS BY FAMILY AND CHROMOSOME FOR Hypothesis_0
==========================================================
        
CHROMOSOME     FAMILY     WEIGHT
---------------------------------
         1          1    0.00645
         1          2    0.00655
         1          3    0.00650
         1          4    0.00643
         1          5    0.00657
         1          6    0.00644
         1          7    0.00651
         1          8    0.00644
         1          9    0.00645
         1         10    0.00648
         1         11    0.00643
         1         12    0.00643
         1         13    0.00648
         1         14    0.00645
         1         15    0.00642
         1         16    0.00648
         1         17    0.00639
         1         18    0.00650
         .
         .
         .

If you're interested in more detailed information, the --reps option will be useful. For this option, PSEUDO output all proportions used to calculate each family variance weight. The table below shows a small portion of the output produced for our example simulation using the pseudo_scan.ped data set. If you recall, we generated a 20 gene-dropping replicates of the original data to be used for z-score sampling. If you look at the look at the column under the "SELECTIONS" heading, the number of times each of the 20 z-score replicates in the pool for family 1 on chromosome 1 was chosen is listed; here replicate 1 was chosen 4952 times. Of the simulations where replicate 1 was chosen, 30 were positive for Hypothesis_0, with overall success proportion of 0.00606.


REPLICATE SUMMARY RESULTS FOR Hypothesis_0
==========================================

[CHROMOSOME 1 FAMILY 5001]


    REPLICATE   SELECTIONS     POSITIVE     PROP POS
----------------------------------------------------

            1         4952           30      0.00606
            2         5041           25      0.00496
            3         4962           37      0.00746
            4         5040           26      0.00516
            5         4930           30      0.00609
            6         5078           35      0.00689
            7         5143           39      0.00758
            8         5021           36      0.00717
            9         4982           33      0.00662
           10         5018           34      0.00678
           11         5000           27      0.00540
           12         5093           28      0.00550
           13         4918           20      0.00407 
           14         4987           39      0.00782
           15         5072           27      0.00532
           16         4866           35      0.00719
           17         4948           35      0.00707
           18         4951           35      0.00707
           19         4995           28      0.00561
           20         5003           30      0.00600
                                                                       

 
 

University of Michigan | School of Public Health | Abecasis Lab