def test_intercept1(self):
    prod = (self.terms[0] * self.terms[2])
    f = (self.formula + formula.I)
    icol = f.names().index('intercept')
    f.namespace = self.namespace
    assert_almost_equal(f()[icol], np.ones((40,)))
def test_intercept1(self):
    prod = (self.terms[0] * self.terms[2])
    f = (self.formula + formula.I)
    icol = f.names().index('intercept')
    f.namespace = self.namespace
    assert_almost_equal(f()[icol], np.ones((40,)))