[发明专利]一种用于计算基因的组织特异表达的鲁棒z-score打分方法

权利要求书

1.一种用于计算基因的组织特异表达的鲁棒z-score打分方法，其特征在于：包括以下步骤：

S1.取某个基因或转录本在N个组织的表达值按从小到达的顺序进行排序；

S2.为每个表达值构建线性拟合模型：

y_(i)＝β₀+β₁*i+e_(i) (1)

其中β₀和β₁是回归系数，其初始值利用最小二乘法生成，i∈[1…N]表示表达值排序的序号，y_(i)表示排序后的第i个表达值，e_(i)表示第i个表达值的残差；

S3.计算每个表达值的权值：

其中，w_(i)是权值，k是加权阈值；

S4.根据计算得到的权值利用加权最小二乘法对β₀、β₁、e_(i)进行更新；并令其中是前次加权的样本标准差；

S5.重复执行步骤S2～S4，直至β₀和β₁恒定，将β₀和β₁代入式(1)中，得到表达值y_(i)的估计值利用y_(i)求取残差，然后利用公式(2)求取权值；

S6.对w_(i)进行标准化：

利用权值产生加权均值和加权标准差

S7.计算鲁棒z-score：

S8.利用鲁棒z-score计算组织特异表达值的统计显著性p_i值：

其中，Φ(x)是标准正态分布的分布函数；

S9.设定一个阈值，若p_i小于设定的阈值，则拒绝零假设，认为该基因或细胞型在序号为i的组织上有特异表达；如果p_i大于或等于设定的阈值，则接受零假设，认为该基因或细胞型在序号为i的组织上没有特异表达。

2.根据权利要求1所述的用于计算基因的组织特异表达的鲁棒z-score打分方法，其特征在于：所述S9步骤中设定的阈值为0.05。