Example #1
0
 def __init__(self,
              kinship,
              design,
              kinship_cm=None,
              kinship_cross=None,
              jitter=1e-4):
     if kinship_cm is None: kinship_cm = kinship
     if kinship_cross is None: kinship_cross = kinship
     ff_dim = 2
     self.covff = FreeFormCov(ff_dim, jitter=1e-4)
     self._K = kinship
     self._ZK = sp.dot(design, kinship_cross.T)
     self._KZ = sp.dot(kinship_cross, design.T)
     self._ZKZ = sp.dot(design, sp.dot(kinship_cm, design.T))
     Covariance.__init__(self, kinship.shape[0])