mklmm.py 文件源码

def getInitialHyps(self, X, C, y):
        self.logdetXX  = np.linalg.slogdet(C.T.dot(C))[1]

        hyp0_sig2e = [0.5*np.log(0.5*y.var())]
        Linreg = sklearn.linear_model.LinearRegression(fit_intercept=False, normalize=False, copy_X=False)
        Linreg.fit(C, y)
        hyp0_fixedEffects = Linreg.coef_        
        return hyp0_sig2e, hyp0_fixedEffects