python关于q检验

关于statsmodels模块里的q_stat

python里面自带Q检验，
氮素！
小伙伴们注意一下，python里的q检验最好用acf的接口进去
什么意思呢？
就是说，不要直接用statsmodel.tsa.stattools.q_stat
因为这个函数里面的x序列应该是你的数据的自相关系数序列，
不是数据序列！
对，我就是这样做出来的东西奇奇怪怪。。后来才发现
所以不要和我一样踩到坑里。
如果需要q检验的话，
可以调用statsmodel.tsa.stattools.acf(x,qstat = True),
x就是数据序列
这样就会自动输出qstat的值以及相应的p值

具体一点的例子就是

from scipy import stats 
import numpy as np
def q_stat(x, nobs, g=0 ,type="ljungbox"):
    #g是在做AR模型总体检验时，对自由度的修正，是AR参数的个数
    
    x = np.asarray(x)
    mean = np.mean(x)
    var