早上好。
我正在使用CausalImpact python3软件包(https://nbviewer.jupyter.org/github/dafiti/causalimpact/blob/master/examples/getting_started.ipynb),并且希望将摘要统计信息显示在数据框中。这是我的代码。
from causalimpact import CausalImpact
import boto3
import numpy as np
import pandas as pd
from statsmodels.tsa.arima_process import ArmaProcess
import matplotlib
import os
from datetime import datetime
#setting up now time from datetime
now = datetime.now()
#generate our data with random values to check for lift
np.random.seed(12345)
ar = np.r_[1, 0.9]
ma = np.array([1])
arma_process = ArmaProcess(ar, ma)
X = 100 + arma_process.generate_sample(nsample=100)
y = 1.2 * X + np.random.normal(size=100)
y[70:] += 5
#generate data frame to be used for CausalImpact
data = pd.DataFrame({'X': X, 'y': y}, columns=['y', 'X'])
data.plot()
pre_period = [0, 69]
post_period = [70, 99]
ci = CausalImpact(data, pre_period, post_period)
#ci.plot(figsize=(12, 6))
print(ci.summary())
问题是我希望结果在数据框中
这是上面代码的结果
Average Cumulative
实际125.2 3756.9 预测(标准)120.4(0.3)3611.8(8.0) 95%CI [119.9,120.9] [3596.3,3627.8]
绝对效果(标准差)4.8(0.3)145.1(8.0) 95%CI [4.3,5.4] [129.1,160.6]
相对效果(标准)4.0%(0.2%)4.0%(0.2%) 95%CI [3.6%,4.4%] [3.6%,4.4%]
后尾区概率p:0.0 后验概率。因果关系:100.00%
有关更多详细信息,请运行命令:print(impact.summary('report')) 分析报告{CausalImpact}
如何获取结果并将其放入数据框? (附加图片,因为结果表示在这里看起来不正确)
非常感谢你 哈里斯 enter image description here