Skip to main content

Table 3 Folded Skellam mixture parameter estimates and results of LRTs for equal Poisson mean values

From: Analyzing allele specific RNA expression using mixture models

Parameter Mix1 Mix2 Mix3 Mix4 Mix5 Mix6
π i 0.54 (0.54, 0.55) 0.1 (0.10, 0.11) 0.0065 (0.0064, 0.0066) 0.037 (0.036, 0.038) 0.0003 (0.0003, 0.00035) 0.3 (0.3, 0.31)
λ i,1 65.7 (65.4, 66.5) 83.8 (82.6, 84.2) 268 (263.3, 269.4) 92.7 (91.4, 93.1) 214.8 (212.2, 216.3) 4.81 (4.75, 4.84)
λ i,2 69.2 (69.2, 70.2) 106 (105, 107) 80.3 (79.9, 81.5) 166.0 (165.9, 169.1) 78.1 (77.0, 78.5) 5.39 (5.29, 5.40)
L 0 −17852 −2074 NA −650 NA −7860
L 1 −17864 −1967 −522 −8233
p -value 1 <0.00001 <0.00001 1
No. of SNPs 5459 482 3 130 2 4626
No. of Genes 471 165 3 72 2 407
  1. Only SNPs on 3’ UTR and classified into Poisson mixture component 1 were used for fitting the folded Skellam mixture (overall log-likelihood = -34979; BIC = 70117; sample-size = 10702; (λ i,1, λ i,2) is estimate of the ordered pair (λ i,P , λ i,M ). NAs indicate insufficient sample sizes for LRTs