|
245 | 245 | " T_data,\n", |
246 | 246 | " Y_data,\n", |
247 | 247 | " coords,\n", |
248 | | - " priors={\n", |
249 | | - " \"beta_\": [0, 1],\n", |
250 | | - " \"beta_trt\": [0, 1],\n", |
251 | | - " \"alpha_trt\": [0, 1],\n", |
252 | | - " \"alpha_outcome\": [0, 1],\n", |
253 | | - " \"sigma\": 1,\n", |
254 | | - " \"beta_ps\": [0, 1],\n", |
255 | | - " },\n", |
| 248 | + " priors=None,\n", |
256 | 249 | " noncentred=True,\n", |
257 | 250 | " normal_outcome=True,\n", |
258 | 251 | "):\n", |
| 252 | + " if priors is None:\n", |
| 253 | + " priors = {\n", |
| 254 | + " \"beta_\": [0, 1],\n", |
| 255 | + " \"beta_trt\": [0, 1],\n", |
| 256 | + " \"alpha_trt\": [0, 1],\n", |
| 257 | + " \"alpha_outcome\": [0, 1],\n", |
| 258 | + " \"sigma\": 1,\n", |
| 259 | + " \"beta_ps\": [0, 1],\n", |
| 260 | + " }\n", |
259 | 261 | " with pm.Model(coords=coords) as model:\n", |
260 | 262 | " X_data_trt = pm.Data(\"X\", X_trt, dims=(\"obs\", \"beta_trt\"))\n", |
261 | 263 | " X_data_outcome = pm.Data(\"X_outcome\", X_outcome, dims=(\"obs\", \"betas\"))\n", |
|
606 | 608 | " X_trt,\n", |
607 | 609 | " T_data,\n", |
608 | 610 | " coords,\n", |
609 | | - " priors={\n", |
610 | | - " \"beta_\": [0, 1],\n", |
611 | | - " \"beta_trt\": [0, 1],\n", |
612 | | - " \"alpha_trt\": [0, 1],\n", |
613 | | - " \"alpha_outcome\": [0, 1],\n", |
614 | | - " },\n", |
| 611 | + " priors=None,\n", |
615 | 612 | " noncentred=True,\n", |
616 | 613 | "):\n", |
| 614 | + " if priors is None:\n", |
| 615 | + " priors = {\n", |
| 616 | + " \"beta_\": [0, 1],\n", |
| 617 | + " \"beta_trt\": [0, 1],\n", |
| 618 | + " \"alpha_trt\": [0, 1],\n", |
| 619 | + " \"alpha_outcome\": [0, 1],\n", |
| 620 | + " }\n", |
617 | 621 | " with pm.Model(coords=coords) as model_trt:\n", |
618 | 622 | " X_data_trt = pm.Data(\"X\", X_trt, dims=(\"obs\", \"betas_trt\"))\n", |
619 | 623 | " T_data_ = pm.Data(\"T\", T_data, dims=\"obs\")\n", |
|
649 | 653 | " X_outcome,\n", |
650 | 654 | " Y_data,\n", |
651 | 655 | " coords,\n", |
652 | | - " priors={\n", |
653 | | - " \"beta_\": [0, 1],\n", |
654 | | - " \"beta_trt\": [0, 1],\n", |
655 | | - " \"alpha_trt\": [0, 1],\n", |
656 | | - " \"alpha_outcome\": [0, 1],\n", |
657 | | - " \"sigma\": 1,\n", |
658 | | - " \"beta_ps\": [0, 1],\n", |
659 | | - " },\n", |
| 656 | + " priors=None,\n", |
660 | 657 | " noncentred=True,\n", |
661 | 658 | " spline_component=False,\n", |
662 | 659 | " propensity_score_idata=None,\n", |
663 | 660 | " normal_outcome=True,\n", |
664 | 661 | " winsorize_boundary=0.0,\n", |
665 | 662 | "):\n", |
| 663 | + " if priors is None:\n", |
| 664 | + " priors = {\n", |
| 665 | + " \"beta_\": [0, 1],\n", |
| 666 | + " \"beta_trt\": [0, 1],\n", |
| 667 | + " \"alpha_trt\": [0, 1],\n", |
| 668 | + " \"alpha_outcome\": [0, 1],\n", |
| 669 | + " \"sigma\": 1,\n", |
| 670 | + " \"beta_ps\": [0, 1],\n", |
| 671 | + " }\n", |
666 | 672 | " propensity_scores = az.extract(propensity_score_idata)[\"p\"]\n", |
667 | 673 | " with pm.Model(coords=coords) as model_outcome:\n", |
668 | 674 | " X_data_outcome = pm.Data(\"X_outcome\", X_outcome, dims=(\"obs\", \"betas\"))\n", |
|
908 | 914 | " X_outcome,\n", |
909 | 915 | " Y_data,\n", |
910 | 916 | " coords,\n", |
911 | | - " priors={\n", |
912 | | - " \"beta_\": [0, 1],\n", |
913 | | - " \"alpha_outcome\": [0, 1],\n", |
914 | | - " \"sigma\": 1,\n", |
915 | | - " },\n", |
| 917 | + " priors=None,\n", |
916 | 918 | " noncentred=True,\n", |
917 | 919 | "):\n", |
| 920 | + " if priors is None:\n", |
| 921 | + " priors = {\"beta_\": [0, 1], \"alpha_outcome\": [0, 1], \"sigma\": 1}\n", |
918 | 922 | " with pm.Model(coords=coords) as reg_model:\n", |
919 | 923 | " X_data_outcome = pm.Data(\"X_outcome\", X_outcome, dims=(\"obs\", \"betas\"))\n", |
920 | 924 | " Y_data_ = pm.Data(\"Y\", Y_data, dims=\"obs\")\n", |
|
0 commit comments