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