statsmodels.tools.tools.isestimable¶

statsmodels.tools.tools.isestimable(c, d)[source]¶

True if (Q, P) contrast c is estimable for (N, P) design d.

From an Q x P contrast matrix C and an N x P design matrix D, checks if the contrast C is estimable by looking at the rank of vstack([C,D]) and verifying it is the same as the rank of D.

Parameters:¶

c : array_like¶

A contrast matrix with shape (Q, P). If 1 dimensional assume shape is (1, P).

d : array_like¶

The design matrix, (N, P).

Returns:¶

True if the contrast c is estimable on design d.

Return type:¶

bool

Examples

>>> d = np.array([[1, 1, 1, 0, 0, 0],
...               [0, 0, 0, 1, 1, 1],
...               [1, 1, 1, 1, 1, 1]]).T
>>> isestimable([1, 0, 0], d)
False
>>> isestimable([1, -1, 0], d)
True