limix.stats.linear_kinship

limix.stats.linear_kinship(G, out=None, verbose=True)

Estimate Kinship matrix via linear kernel.

Examples

>>> from numpy.random import RandomState
>>> from numpy import array_str
>>> from limix.stats import linear_kinship
>>>
>>> random = RandomState(1)
>>> X = random.randn(4, 100)
>>> K = linear_kinship(X, verbose=False)
>>> print(array_str(K, precision=4))
[[ 0.9131 -0.1928 -0.3413 -0.379 ]
 [-0.1928  0.8989 -0.2356 -0.4704]
 [-0.3413 -0.2356  0.9578 -0.3808]
 [-0.379  -0.4704 -0.3808  1.2302]]