anova（一更）

最新推荐文章于 2024-08-28 16:14:39 发布

qdu5er

最新推荐文章于 2024-08-28 16:14:39 发布

阅读量628

点赞数

分类专栏：参数假设检

本文链接：https://blog.csdn.net/qdu5er/article/details/81368703

版权

参数假设检专栏收录该内容

2 篇文章 0 订阅

订阅专栏

import tushare as ts

from anova import Anova
gd=ts.get_hist_data(start='2018-01-01',code='600795')#国电电力
hd=ts.get_hist_data(start='2018-01-01',code='600726')#华电能源
#gd=gd.loc[:,'close']
###############################
print(Anova(gd))

def Anova(factor):
    import numpy as np
    import pandas as pd
    m=len(factor.columns)
    n=len(factor.index)
    fmean=factor.mean()    
    totalmean=fmean.mean()/m
    fstd=factor.std()
    sse=(fstd.sum())*n
    ssb=np.array(fmean)
    ssb=np.power(ssb-totalmean,2)*n
    ssb=pd.Series(ssb)
    ssb=ssb.sum()
    F=ssb*(n-m)/(sse*(m-1))
    return F

作者：隐