Menu

SeparateChromosomes2 output

2018-05-02
2018-05-02
  • Roy Francis

    Roy Francis - 2018-05-02

    I ran SeparateChromosomes2 on my dataset like so:

    java -cp ~/ngs/programs/lep-map3/bin/ SeparateChromosomes2 data="snp-2.txt" lodLimit=10 theta=0.03 numThreads=8 subsample=0.25 > "snp-3-ss.txt"
    

    With the subsample argument, my idea was to use 25% of the total number of SNPs. But my output file has the same number of lines are the original number of SNPs. The output is mostly zeros. Only 4% of the rows have an integer other than zero.

    0
    0
    0
    0
    0
    1513
    0
    0
    0
    ...
    

    What does this mean?

     
  • Pasi Rastas

    Pasi Rastas - 2018-05-02

    Dear Roy,

    Yes, this is expected. In your case, on average at least 75% of rows are zeros. This way the map is compatible with the data and you could add the remaining markers to the map if you later want that.

    If you want to permanently reduce the number of markers, you could use simple scripting on the input data.

    Cheers,
    Pasi

     

Log in to post a comment.