1.一种miRNA因果调控网络识别方法,其特征在于,包括:获取给定匹配样本的miRNA单细胞转录组矩阵和靶基因单细胞转录组矩阵,其中,所述匹配样本中包括多个单细胞匹配样本,每个所述单细胞匹配样本中包括多个miRNA和靶基因;
获取所述miRNA单细胞转录组矩阵和所述靶基因单细胞转录组矩阵的miRNA因果调控网络,其中,所述因果调控网络用于表示所述miRNA和所述靶基因之间的因果调控关系;
将所述miRNA单细胞转录组矩阵和所述靶基因单细胞转录组矩阵中的第i个单细胞剔除,获取剔除第i个单细胞后的miRNA因果调控网络,包括:将所述miRNA单细胞转录组矩阵和所述靶基因单细胞转录组矩阵中的第i个单细胞剔除,得到剔除第i个单细胞后的miRNA和靶基因单细胞转录组矩阵;
根据剔除第i个单细胞后的miRNA和靶基因单细胞转录组矩阵,获取所述剔除第i个单细胞后的miRNA和靶基因单细胞转录组矩阵组成的完全部分有向无环CPDAG图;
计算所述CPDAG图中每条边的miRNA‑靶基因的因果效应,得到所述剔除第i个单细胞后的miRNA和靶基因单细胞转录组矩阵的miRNA因果调控网络,计算方法是:先计算剔除单细胞前的CPDAG图中一条边的因果效应,再计算剔除第i个单细胞后的CPDAG图中这条边的因果效应,计算这条边剔除第i个单细胞前的因果效应的显著性p值和剔除第i个单细胞后的因果效应的显著性p值;
采用上述方法计算CPDAG图中每条边的miRNA‑靶基因的因果效应,得到剔除第i个单细胞前的因果效应的显著性p值矩阵和剔除第i个单细胞后的因果效应的显著性p值矩阵;
给定p阈值,将每条边剔除i个单细胞前的因果效应的显著性p值和剔除第i个单细胞后的因果效应的显著性p值分别与p阈值进行比较,比较结果组成miRNA因果调控关系矩阵,从剔除第i个单细胞前的因果效应的显著性p值矩阵得到未剔除第i个单细胞前的miRNA因果调控关系矩阵,从剔除第i个单细胞后的因果效应的显著性p值矩阵得到剔除第i个单细胞后的miRNA因果调控关系矩阵;
其中,i为大于1的整数;
根据所述miRNA因果调控网络以及所述剔除第i个单细胞后的miRNA因果调控网络,确定第i个单细胞的miRNA因果调控网络。
2.如权利要求1所述的方法,其特征在于,获取所述miRNA单细胞转录组矩阵和所述靶基因单细胞转录组矩阵的miRNA因果调控网络,包括:根据所述miRNA单细胞转录组矩阵和所述靶基因单细胞转录组矩阵,获取所述miRNA单细胞转录组矩阵和所述靶基因单细胞转录组矩阵组成的完全部分有向无环CPDAG图;
计算所述CPDAG图中每条边的miRNA‑靶基因的因果效应,得到所述miRNA单细胞转录组矩阵和所述靶基因单细胞转录组矩阵的miRNA因果调控网络。
3.如权利要求1所述的方法,其特征在于,根据所述miRNA因果调控网络以及所述剔除第i个单细胞后的miRNA因果调控网络,确定第i个单细胞的miRNA因果调控网络,包括:根据所述miRNA单细胞转录组矩阵和所述靶基因单细胞转录组矩阵的miRNA因果调控网络、以及所述剔除第i个单细胞后的miRNA和靶基因单细胞转录组矩阵的miRNA因果调控网络,获取第i个单细胞的差异性miRNA因果调控关系;
融合所有第i个单细胞的所述差异性miRNA因果调控关系,确定所述第i个单细胞的miRNA因果调控网络。
4.一种miRNA因果调控网络识别装置,其特征在于,包括:获取模块,用于获取给定匹配样本的miRNA单细胞转录组矩阵和靶基因单细胞转录组矩阵,其中,所述匹配样本中包括多个单细胞匹配样本,每个所述单细胞匹配样本中包括多个miRNA和靶基因;
所述获取模块,还用于获取所述miRNA单细胞转录组矩阵和所述靶基因单细胞转录组矩阵的miRNA因果调控网络,其中,所述因果调控网络用于表示所述miRNA和所述靶基因之间的因果调控关系;
所述获取模块,还用于将所述miRNA单细胞转录组矩阵和所述靶基因单细胞转录组矩阵中的第i个单细胞剔除,获取剔除第i个单细胞后的miRNA因果调控网络,其中,i为大于1的整数;
确定模块,用于根据所述miRNA因果调控网络以及所述剔除第i个单细胞后的miRNA因果调控网络,确定第i个单细胞的miRNA因果调控网络。
5.如权利要求4所述的装置,其特征在于,所述获取模块,具体用于根据所述miRNA单细胞转录组矩阵和所述靶基因单细胞转录组矩阵,获取所述miRNA单细胞转录组矩阵和所述靶基因单细胞转录组矩阵组成的完全部分有向无环图CPDAG;
计算所述CPDAG中每条边的miRNA‑靶基因的因果效应,得到所述miRNA单细胞转录组矩阵和所述靶基因单细胞转录组矩阵的miRNA因果调控网络。
6.如权利要求5所述的装置,其特征在于,所述获取模块,具体用于将所述miRNA单细胞转录组矩阵和所述靶基因单细胞转录组矩阵中的第i个单细胞剔除,得到剔除第i个单细胞后的miRNA和靶基因单细胞转录组矩阵;
根据剔除第i个单细胞后的miRNA和靶基因单细胞转录组矩阵,获取所述剔除第i个单细胞后的miRNA和靶基因单细胞转录组矩阵组成的完全部分有向无环图CPDAG;
计算所述CPDAG中每条边的miRNA‑靶基因的因果效应,得到所述剔除第i个单细胞后的miRNA和靶基因单细胞转录组矩阵的miRNA因果调控网络。
7.如权利要求6所述的装置,其特征在于,所述确定模块,具体用于根据所述miRNA单细胞转录组矩阵和所述靶基因单细胞转录组矩阵的miRNA因果调控网络、以及所述剔除第i个单细胞后的miRNA和靶基因单细胞转录组矩阵的miRNA因果调控网络,获取第i个单细胞的差异性miRNA因果调控关系;
融合所有第i个单细胞的所述差异性miRNA因果调控关系,确定所述第i个单细胞的miRNA因果调控网络。
8.一种电子设备,其特征在于,包括处理器、存储介质和总线,所述存储介质存储有所述处理器可执行的机器可读指令,当所述电子设备运行时,所述处理器与所述存储介质之间通过总线通信,所述处理器执行所述机器可读指令,以执行权利要求1‑3任一项提供的方法。
9.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器运行时执行权利要求1‑3任一项提供的方法。