在Python中查找三向方差分析的P值

Question

我完全被这件事缠住了。有没有办法从多次模拟运行中提取p值？我在这里尝试过这个方法https://pythonfordatascience.org/anova-2-way-n-way/

pValues = np.zeros(100)

for k in range (99) :
       y = simCorn(overallEffect=10, fertilizerEffect=[1,2,3],rowEffect=[1,0,1], colEffect=[0,1,1]) 
       pValues[k] = ols('Yield ~ C(Fertilizer)+C(Row)+C(Column)', y).fit()

但我甚至不能让基本ols函数工作，更不用说分离和提取p值了

所以，我尝试了statsmodels方法

pValues = np.zeros(100)

for k in range (99) :
       y = simCorn(overallEffect=10, fertilizerEffect=[1,2,3],rowEffect=[1,0,1], colEffect=[0,1,1]) 
       pValues[k] = sm.stats.anova_lm(data=y).Pr(F)

……但似乎什么都没用

这里的任何帮助都将不胜感激