Estimate parameters
Arguments
- psdesign
A psdesign object, it must have risk model and integration model components
- estimate
An estimate object created by ps_estimate
Examples
test <- psdesign(generate_example_data(n = 100), Z = Z, Y = Y.obs, S = S.obs, BIP = BIP)
test + integrate_parametric(S.1 ~ BIP) + risk_binary(D = 50) + ps_estimate(method = "BFGS")
#> Augmented data frame: 100 obs. by 6 variables.
#> Z Y S.1 S.0 cdfweights BIP
#> 1 0 0 NA 0.9050 1 0.8621
#> 2 1 0 0.706 NA 1 -0.2432
#> 3 1 0 0.675 NA 1 -0.2061
#> 4 0 0 NA -0.0366 1 0.0192
#> 5 0 0 NA 0.1348 1 0.0296
#> 6 0 0 NA 0.6176 1 0.5498
#>
#> Empirical TE: 0.775
#>
#> Mapped variables:
#> Z -> Z
#> Y -> Y.obs
#> S -> S.obs
#> BIP -> BIP
#>
#> Integration models:
#> integration model for S.1 :
#> integrate_parametric(formula = S.1 ~ BIP )
#>
#> Risk models:
#> risk_binary(D = 50 )
#>
#> Estimated parameters:
#> (Intercept) S.1 Z S.1:Z
#> -0.510 -0.298 -1.272 -0.941
#> Convergence: TRUE
#>
#> No bootstraps present, see ?ps_bootstrap.