1.基于改进函数因果似然度的糖尿病危险因素因果发现方法,其特征在于,包括:获取糖尿病危险因素变量子集的联合分布,具体包括:
设{X1,X2,···,XN}为糖尿病危险因素变量集合,其中N为糖尿病危险因素变量个数;
令G为子集X={X1,X2,···,Xn}上的因果结构,n为具有因果关系的糖尿病危险因素变量个数;P(Xi=x)为变量Xi=x的概率,x为对应的观测值; 为变量Xi在其所有因变量条件下的概率,其中1≤i≤n≤N,i为具有因果关系的糖尿病危险因素变量下标;设因果结构G满足因果马尔可夫假设、因果忠诚性假设,则其联合分布P(X)表示为:其中 为变量Xi的所有因变量;
根据所述联合分布与因果结构求观察数据的对数似然度,具体包括:给定 危险 因素 观察 数 据 其中 为 n维 向量 ,即令 表示为包含 的观测值,观察数据的对数似然度
为:
将所述观察数据的对数似然度转化成观察数据噪声的对数似然度,进而建立糖尿病危险因素FCL模型,具体包括:设Ei表示Xi对应的随机噪声变量,采用加性噪声模型 作为因果关系生成机制,其中Fi是Xi的因果函数,随机噪声变量Ei独立于因变量 因此,有如下等式:其中, 表示为包含 的观测值, 为Xi的所有因变量;
由式(2)和(3)可知,观察数据的对数似然度等于观察数据噪声的对数似然度;又设S=
迭代后的糖尿病危险因素因果结构对数似然度表示为:
其中,α′是迭代后的总阈值;由(7)式可知,初始糖尿病危险因素IFCL模型的总阈值为nα,视为每个因果节点对数似然度都具有一个阈值α,即每次使用式(6)进行迭代后,更新的节
点对数似然度都会有一个阈值α减少,这会致使迭代后的总阈值小于初始糖尿病危险因素IFCL模型的总阈值,即α′<nα;因此,迭代过程中必须搜索到对数似然度更大的因果节点,才能达到迭代终止条件 这是基于IFCL模型的糖尿病危险因素因果发现方法能输出更优化因果结构的根本原因。