Method of statistical genomic analysis

A method for statistical genomic analysis, comprising:

• designing a study of a microarray chips of genomic material that utilizes multi-stage testing;
• normalizing a measurement of data received from the study;
• estimating an experimental error of the measurement of data with an error estimator (θi) for the ith gene is {hacek over (θ)}i=(1−Bi){circumflex over (θ)}i+Bi{tilde over (θ)}i, where Bi=((k−r−2)/ (k−r))Si/(Si+Â), where k is the number of genes, where r is the number of predictor variables, and where Si is the variance of the ith gene;
• and
• inferring conclusions about the expression of the microarray chips of genomic material from the data received from the study based on the calculation of the probability (Xi) that a difference between groups during an experiment happened by chance, where the probability density function for Xi is
• f
• X
• i
• ⁡
• (
• x
• i
• ;
• λ
• _
• ,
• r
• _
• ,
• s
• _
• )
• =
• λ
• 0
• ⁢
• I
• (
• 0
• ,
• 1
• )
• ⁡
• (
• x
• i
• )
• +
• ∑
• j
• =
• 1
• v
• ⁢
• λ
• j
• ⁢
• I
• (
• 0
• ,
• 1
• )
• ⁡
• (
• x
• i
• )
• ⁢
• x
• i
• r
• j
• -
• 1
• ⁡
• (
• 1
• -
• x
• i
• )
• s
• j
• -
• 1
• B
• ⁡
• (
• r
• j
• ,
• s
• j
• )
• ,
• where v is a mixture of beta distributions with parameters rj and sj, j=1 . . . v, and
• where B(rj, sj)=, ∫01urj−1(1−u)sj−1du.