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 : (Q, P) array-like
D: (N, P) array-like :
|
---|---|
Returns: | tf : 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