library(corpcor)#compute partial correlations

library(GeneNet)#compute the p_value of corpcor

m.cov = cbind(bai[,i], LOS[,i], T[,i], P[,i]) #數(shù)組

    c = cov(m.cov)#創(chuàng)建協(xié)方差矩陣

    m.pcor = cor2pcor(c) #compute partial correlations

    kappan = dim(T)[1]-4+1#自由度

    p.pcor = cor0.test(m.pcor, kappa = kappan)#計(jì)算p_value