MACS3.Signal.VariantStat module

Module for SAPPER BAMParser class

This code is free software; you can redistribute it and/or modify it under the terms of the BSD License (see the file COPYING included with the distribution).

@status: experimental @version: $Revision$ @author: Tao Liu @contact: tliu4@buffalo.edu

MACS3.Signal.VariantStat.CalModel_Heter_AS(top1_bq_T, top1_bq_C, top2_bq_T, top2_bq_C, max_allowed_ar=0.99)

Return (lnL, BIC)

kc ki AS_alleleratio

MACS3.Signal.VariantStat.CalModel_Heter_noAS(top1_bq_T, top1_bq_C, top2_bq_T, top2_bq_C)

Return (lnL, BIC)

k_T k_C

MACS3.Signal.VariantStat.CalModel_Homo(top1_bq_T, top1_bq_C, top2_bq_T, top2_bq_C): Return (lnL, BIC).

MACS3.Signal.VariantStat.calculate_GQ(lnL1, lnL2, lnL3): GQ1 = -10*log_{10}((L2+L3)/(L1+L2+L3))

MACS3.Signal.VariantStat.calculate_GQ_heterASsig(lnL1, lnL2)