1.一种细胞特异性miRNA因果调控网络推理方法,其特征在于,包括:
获取目标单细胞的miRNA和靶基因的双重表达谱数据,所述目标单细胞的miRNA和靶基因的双重表达谱数据中包括背景数据和扰动数据;
根据所述目标单细胞的miRNA和靶基因的双重表达谱数据,获取CPDAG图,所述CPDAG图中包括多条边;
计算所述CPDAG图中每条边在背景数据和扰动数据之间的差异因果效应;
根据所述差异因果效应符合预设条件的边,得到所述目标单细胞的细胞特异性miRNA因果调控网络;
所述根据所述目标单细胞的miRNA和靶基因的双重表达谱数据,获取CPDAG图,包括:根据所述目标单细胞的miRNA和靶基因的双重表达谱数据获取所述miRNA靶基因的先验网络图;根据所述先验网络图中每条边的因果关系,获得多个编码相同条件独立关系的多个DAG等价类;通过所述CPDAG图表示所述多个编码相同条件独立关系的多个DAG等价类;
所述计算所述CPDAG图中每条边在背景数据和扰动数据之间的差异因果效应,包括:通过修正潜在混杂因子,计算所述CPDAG图中每条边在所述背景数据中的第一因果效应;通过修正潜在混杂因子,计算所述CPDAG图中每条边在所述扰动数据中的第二因果效应;根据所述第一因果效应和所述第二因果效应,获取所述差异因果效应。
2.根据权利要求1所述的方法,其特征在于,所述获取目标单细胞的miRNA和靶基因的双重表达谱数据,包括:获取所述目标单细胞的转录组数据,根据预设的基因注释信息,获取所述目标单细胞的miRNA和靶基因的双重表达谱数据;
通过统计扰动策略,根据所述目标单细胞的miRNA和靶基因的双重表达谱数据获取所述扰动数据,所述统计扰动策略包括每次统计时敲除一个单细胞;
将获取到的所述目标单细胞的miRNA和靶基因的双重表达谱数据作为所述背景数据。
3.根据权利要求1或2所述的方法,其特征在于,在得到所述目标单细胞的细胞特异性miRNA因果调控网络之后,所述方法还包括:通过网络相似性分析、聚类分析以及验证分析中的至少一种分析方式,对所述目标单细胞的细胞特异性miRNA因果调控网络验证。
4.一种细胞特异性miRNA因果调控网络推理装置,其特征在于,包括:
获取模块,用于获取目标单细胞的miRNA和靶基因的双重表达谱数据,所述目标单细胞的miRNA和靶基因的双重表达谱数据中包括背景数据和扰动数据;
所述获取模块,还用于根据所述目标单细胞的miRNA和靶基因的双重表达谱数据,获取CPDAG图,所述CPDAG图中包括多条边;
计算模块,用于计算所述CPDAG图中每条边在背景数据和扰动数据之间的差异因果效应;
所述获取模块,还用于根据所述差异因果效应符合预设条件的边,得到所述目标单细胞的细胞特异性miRNA因果调控网络;
所述获取模块,具体用于根据所述目标单细胞的miRNA和靶基因的双重表达谱数据获取所述miRNA靶基因的先验网络图;根据所述先验网络图中每条边的因果关系,获得多个编码相同条件独立关系的多个DAG等价类;通过所述CPDAG图表示所述多个编码相同条件独立关系的多个DAG等价类;
所述计算模块,具体用于通过修正潜在混杂因子,计算所述CPDAG图中每条边在所述背景数据中的第一因果效应;通过修正潜在混杂因子,计算所述CPDAG图中每条边在所述扰动数据中的第二因果效应;根据所述第一因果效应和所述第二因果效应,获取所述差异因果效应。
5.根据权利要求4所述的装置,其特征在于,所述获取模块,具体用于获取所述目标单细胞的转录组数据,根据预设的基因注释信息,获取所述目标单细胞的miRNA和靶基因的双重表达谱数据;
通过统计扰动策略,根据所述目标单细胞的miRNA和靶基因的双重表达谱数据获取所述扰动数据,所述统计扰动策略包括每次统计时敲除一个单细胞;
将获取到的所述目标单细胞的miRNA和靶基因的双重表达谱数据作为所述背景数据。
6.一种电子设备,其特征在于,包括处理器、存储介质和总线,所述存储介质存储有所述处理器可执行的机器可读指令,当所述电子设备运行时,所述处理器与所述存储介质之间通过总线通信,所述处理器执行所述机器可读指令,以执行权利要求1‑3任一项所述的方法。
7.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器运行时执行权利要求1‑3任一项所述的方法。