考虑以下示例:


df = pd.read_csv('myFile.txt',delim_whitespace=True,header=None)
df.columns=['vary','vax1','varx2']
y,X = ps.dmatrices('vary ~ varx1 + varx2',data=df, return_type='dataframe')
model = sm.OLS(y,X) # Describe Model
results = model.fit() # Fit model
print results.summary()

hypotheses = 'varx1 = 0.0'
t_test = results.t_test(hypotheses)
print(t_test)



这是针对替代假设“ H1:varX1的系数不同于零”测试“ H0:varX1的系数为零”

我想使“ H1:VarX1的系数小于-c,c> 0”(单侧替代)。

可能吗?

最佳答案

Statsmodels似乎不允许单尾测试。您可以在拟合模型后调用scipy.stats

from scipy import stats
c = 0.5 #lower bound
print(stats.t.cdf(var_x1, df, loc=c, scale=1))


这将为您提供c右边的t-dist下的区域。

关于python - Python和Statsmodels:如何包含替代test_t假设?,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/35456624/

10-12 18:51