Python StatTestMC.cdfの例

コード例 #1

ファイル: try_mctools.py プロジェクト: kasunsp/pinalpha_mvp

mc1 = StatTestMC(normalnoisesim, lb)
mc1.run(5000, statindices=lrange(4))

print(mc1.summary_quantiles([1,2,3], stats.chi2([2,3,4]).ppf,
                            varnames=['lag 1', 'lag 2', 'lag 3'],
                            title='acorr_ljungbox'))
print('\n\n')

frac = [0.01, 0.025, 0.05, 0.1, 0.975]
crit = stats.chi2([2,3,4]).ppf(np.atleast_2d(frac).T)
print(mc1.summary_cdf([1,2,3], frac, crit,
                      varnames=['lag 1', 'lag 2', 'lag 3'],
                      title='acorr_ljungbox'))
print(mc1.cdf(crit, [1,2,3])[1])

#----------------------

def randwalksim(nobs=500, drift=0.0):
        return (drift+np.random.randn(nobs)).cumsum()

def adf20(x):
    return unitroot_adf(x, 2, trendorder=0, autolag=None)

print(adf20(np.random.randn(100)))

mc2 = StatTestMC(randwalksim, adf20)
mc2.run(10000, statindices=[0,1])
frac = [0.01, 0.05, 0.1]
#bug

コード例 #2

ファイル: try_mctools.py プロジェクト: yarikoptic/pystatsmodels

mc1.run(5000, statindices=range(4))

print mc1.summary_quantiles([1, 2, 3],
                            stats.chi2([2, 3, 4]).ppf,
                            varnames=['lag 1', 'lag 2', 'lag 3'],
                            title='acorr_ljungbox')
print '\n\n'

frac = [0.01, 0.025, 0.05, 0.1, 0.975]
crit = stats.chi2([2, 3, 4]).ppf(np.atleast_2d(frac).T)
print mc1.summary_cdf([1, 2, 3],
                      frac,
                      crit,
                      varnames=['lag 1', 'lag 2', 'lag 3'],
                      title='acorr_ljungbox')
print mc1.cdf(crit, [1, 2, 3])[1]

#----------------------


def randwalksim(nobs=500, drift=0.0):
    return (drift + np.random.randn(nobs)).cumsum()


def adf20(x):
    return unitroot_adf(x, 2, trendorder=0, autolag=None)


print adf20(np.random.randn(100))

mc2 = StatTestMC(randwalksim, adf20)

コード例 #3

ファイル: try_mctools.py プロジェクト: AnaMP/statsmodels

mc1 = StatTestMC(normalnoisesim, lb)
mc1.run(5000, statindices=range(4))

print mc1.summary_quantiles([1,2,3], stats.chi2([2,3,4]).ppf,
                            varnames=['lag 1', 'lag 2', 'lag 3'],
                            title='acorr_ljungbox')
print '\n\n'

frac = [0.01, 0.025, 0.05, 0.1, 0.975]
crit = stats.chi2([2,3,4]).ppf(np.atleast_2d(frac).T)
print mc1.summary_cdf([1,2,3], frac, crit,
                      varnames=['lag 1', 'lag 2', 'lag 3'],
                      title='acorr_ljungbox')
print mc1.cdf(crit, [1,2,3])[1]

#----------------------

def randwalksim(nobs=500, drift=0.0):
        return (drift+np.random.randn(nobs)).cumsum()

def adf20(x):
    return unitroot_adf(x, 2, trendorder=0, autolag=None)

print adf20(np.random.randn(100))

mc2 = StatTestMC(randwalksim, adf20)
mc2.run(10000, statindices=[0,1])
frac = [0.01, 0.05, 0.1]
#bug

コード例 #4

ファイル: try_mctools.py プロジェクト: alfonsodiecko/PYTHON_DIST

mc1 = StatTestMC(normalnoisesim, lb)
mc1.run(5000, statindices=list(range(4)))

print(
    mc1.summary_quantiles(
        [1, 2, 3], stats.chi2([2, 3, 4]).ppf, varnames=["lag 1", "lag 2", "lag 3"], title="acorr_ljungbox"
    )
)
print("\n\n")

frac = [0.01, 0.025, 0.05, 0.1, 0.975]
crit = stats.chi2([2, 3, 4]).ppf(np.atleast_2d(frac).T)
print(mc1.summary_cdf([1, 2, 3], frac, crit, varnames=["lag 1", "lag 2", "lag 3"], title="acorr_ljungbox"))
print(mc1.cdf(crit, [1, 2, 3])[1])

# ----------------------


def randwalksim(nobs=500, drift=0.0):
    return (drift + np.random.randn(nobs)).cumsum()


def adf20(x):
    return unitroot_adf(x, 2, trendorder=0, autolag=None)


print(adf20(np.random.randn(100)))

mc2 = StatTestMC(randwalksim, adf20)