@@ -128,10 +128,10 @@ class PrePostFit(ExperimentalDesign):
128
128
>>> result = cp.pymc_experiments.PrePostFit(
129
129
... sc,
130
130
... treatment_time,
131
- ... formula="actual ~ 0 + a + b + c + d + e + f + g",
131
+ ... formula="actual ~ 0 + a + g",
132
132
... model=cp.pymc_models.WeightedSumFitter(
133
133
... sample_kwargs={
134
- ... "draws": 2000 ,
134
+ ... "draws": 400 ,
135
135
... "target_accept": 0.95,
136
136
... "random_seed": seed,
137
137
... "progressbar": False
@@ -140,16 +140,11 @@ class PrePostFit(ExperimentalDesign):
140
140
... )
141
141
>>> result.summary() # doctest: +NUMBER
142
142
==================================Pre-Post Fit==================================
143
- Formula: actual ~ 0 + a + b + c + d + e + f + g
143
+ Formula: actual ~ 0 + a + g
144
144
Model coefficients:
145
- a 0.3, 94% HDI [0.3, 0.3]
146
- b 0.0, 94% HDI [0.0, 0.0]
147
- c 0.3, 94% HDI [0.2, 0.3]
148
- d 0.0, 94% HDI [0.0, 0.1]
149
- e 0.0, 94% HDI [0.0, 0.0]
150
- f 0.1, 94% HDI [0.1, 0.2]
151
- g 0.0, 94% HDI [0.0, 0.0]
152
- sigma 0.2, 94% HDI [0.2, 0.3]
145
+ a 0.6, 94% HDI [0.6, 0.6]
146
+ g 0.3, 94% HDI [0.3, 0.3]
147
+ sigma 0.7, 94% HDI [0.6, 0.9]
153
148
"""
154
149
155
150
def __init__ (
@@ -777,28 +772,14 @@ class RegressionDiscontinuity(ExperimentalDesign):
777
772
... formula="y ~ 1 + x + treated + x:treated",
778
773
... model=cp.pymc_models.LinearRegression(
779
774
... sample_kwargs={
780
- ... "draws": 2000 ,
775
+ ... "draws": 100 ,
781
776
... "target_accept": 0.95,
782
777
... "random_seed": seed,
783
778
... "progressbar": False,
784
779
... },
785
780
... ),
786
781
... treatment_threshold=0.5,
787
782
... )
788
- >>> result.summary() # doctest: +NUMBER
789
- ============================Regression Discontinuity============================
790
- Formula: y ~ 1 + x + treated + x:treated
791
- Running variable: x
792
- Threshold on running variable: 0.5
793
- <BLANKLINE>
794
- Results:
795
- Discontinuity at threshold = 0.91
796
- Model coefficients:
797
- Intercept 0.0, 94% HDI [0.0, 0.1]
798
- treated[T.True] 2.4, 94% HDI [1.6, 3.2]
799
- x 1.3, 94% HDI [1.1, 1.5]
800
- x:treated[T.True] -3.0, 94% HDI [-4.1, -2.0]
801
- sigma 0.3, 94% HDI [0.3, 0.4]
802
783
"""
803
784
804
785
def __init__ (
0 commit comments