Fitting Haschka's copula-based Instrumental Variable Model (IMA) using Gaussian copulas to address endogeneity without external instruments, while allowing dependence between endogenous and exogenous regressors.
The copula-based IMA approach corrects endogeneity by constructing control functions from the marginal distributions of the endogenous regressors and incorporating them into the estimation equation. Inference is based on bootstrapping.
The key contribution of IMA over the original Park and Gupta (2012) copula
correction (copulaCorrection) is that it does not require the
endogenous and exogenous regressors to be independent. When this independence
assumption is violated, which is common in practice, the original copula
correction yields severely biased estimates of all model coefficients.
IMA corrects for this by exploiting the content of the exogenous
variables through a first-stage auxiliary regression, in a similar way to
instrumental variable identification.
The IMA method requires a normally distributed structural error. At least one endogenous regressor \(P_{i,k}\) is continuously and nonnormally distributed. It does not require that the exogenous regressors \(X_i\) be independent of the endogenous regressors \(P_{i,k}\) or to satisfy any distributional requirement, i.e., \(X_i\) may be binary or non-continuously distributed. This method supports only continuous endogenous regressors.
copulaIMA(
formula,
data,
cdf = c("adj.ecdf", "resc.ecdf", "ecdf", "kde"),
num.boots = 1000,
verbose = TRUE
)A symbolic description of the model to be fitted. See the "Details" section for the exact notation.
A data.frame containing the data of all parts specified in the formula parameter.
Character string specifying the method used to estimate the
marginal distribution functions of the endogenous regressors.
One of "adj.ecdf", "resc.ecdf", "ecdf", or "kde".
"adj.ecdf"Adjusted empirical CDF with midrank correction (Liengaard et al., 2024). Keeps values strictly inside (0,1).
"resc.ecdf"Rescaled empirical CDF via
copula::pobs (Qian et al., 2024).
"ecdf"Empirical CDF with boundary replacement (Becker et al., 2022).
"kde"Integral of a density estimator via ks::kcde
used in (Park and Gupta, 2012).
Positive integer giving the number of bootstrap replications used for standard error estimation. A minimum of 1000 is recommended.
Show details about the running of the function.
An object of class rendo.copula.ima which is a list that contains:
formulaThe formula given to specify the fitted model.
modelThe model.frame used for model fitting.
termsThe terms object used for model fitting.
coefficientsA named vector of all coefficients resulting from model fitting.
names.main.coefsA vector specifying which coefficients are from the model. For internal usage.
fitted.valuesFitted values of the structural model.
residualsThe structural residuals.
boots.paramsThe bootstrapped coefficients.
n.boots.attemptedThe number of bootstrap iterations that were attempted.
n.boots.failedThe number of bootstrap iterations resulting in a failed fit.
cdfThe used cdf function.
names.endo.regsThe names of the continuous endogenous regressors.
res.lm.augmentedThe fitted augmented regression model, including the control function terms.
Consider the structural regression model with \(K\) endogenous regressors:
$$Y_i = \mu + \sum_{k=1}^{K} P_{i,k} \alpha_k + X_i' \beta + \varepsilon_i$$
where \(i = 1, \ldots, n\) is the number of observations, \(Y_i\) is the dependent variable, \(P_{i,k}\) are continuous endogenous regressors that may be correlated with both the structural error \(\varepsilon_i\) and the exogenous regressors \(X_i\), \(X_i\) is a vector of exogenous regressors uncorrelated with \(\varepsilon_i\), and \(\mu, \alpha_k, \beta\) are the structural model parameters.
IMA is an augmented OLS estimator. The estimation proceeds in three steps:
For every explanatory variable \(W \in \{X_1, \ldots, X_J, P_1, \ldots, P_K\}\), compute the normal score \(W^* = \Phi^{-1}(\hat{F}_W(W))\), where \(\hat{F}_W\) is the estimated marginal CDF of \(W\) and \(\Phi^{-1}\) is the standard normal quantile function.
For each endogenous regressor \(P_{i,k}^*\), regress it on the normal scores of the exogenous regressors \(X_1^*, \ldots, X_J^*\) without intercept and retain the residuals \(\hat{\varepsilon}_k\) as the copula correction terms. This first-stage regression exploits the dependence between endogenous and exogenous regressors in a way analogous to instrumental variable identification.
Augment the structural model with \(\hat{\varepsilon}_1, \ldots, \hat{\varepsilon}_{K}\) as additional regressors and estimate by OLS.
Haschka (2025) reported that simulation results showed that copula-based IMA estimator may exhibit a slightly larger bias than alternative two-stage approaches when the sample size is very small (e.g., n = 100) and with intercept. However, the bias decreases rapidly as sample size increases and become negligible for moderate sample sizes (around n >= 600). The estimator appears asymptotically unbiased.
Bootstrap inference is performed by resampling the data with replacement. Degenerate bootstrap samples (e.g. singular design matrices or failed model estimation) are discarded and resampled until the requested number of valid bootstrap replications is obtained. The percentage of discarded samples is reported as a warning. Confidence intervals are computed using only successful bootstrap replications.
The formula argument follows a two-part notation separated by
|. The first part specifies the structural model. The second part
identifies the continuous endogenous regressors using continuous():
y ~ X + P | continuous(P) # typical use: with intercepty ~ X + P1 + P2 | continuous(P1) + continuous(P2) # two endogenous regressorsy ~ X + P - 1 | continuous(P) # no intercept (in simulation settings)In typical applied settings the model includes an intercept. The
no-intercept specification (-1) is primarily used in simulation
studies following the design of (Haschka (2025) Section 4.1)
and is unlikely to be appropriate for real datasets.
Haschka, R. E. (2025). Robustness of copula-correction models in causal analysis: Exploiting between-regressor correlation. IMA Journal of Management Mathematics, 36, 161-180. doi:10.1093/imaman/dpae018
Park, S. and Gupta, S., (2012), "Handling Endogenous Regressors by Joint Estimation Using Copulas", Marketing Science, 31(4), 567-86. doi:10.1287/mksc.1120.0718
Becker, JM., Proksch, D., & Ringle, C.M. Revisiting Gaussian copulas to handle endogenous regressors. Journal of the Academy of Marketing Science, 50, 46–66 (2022). doi:10.1007/s11747-021-00805-y
Liengaard, B. D., Becker, J.-M., Bennedsen, M., Heiler, P., Taylor, L. N., & Ringle, C. M. (2025). Dealing with regression models' endogeneity by means of an adjusted estimator for the Gaussian copula approach. Journal of the Academy of Marketing Science, 53(1), 279-299 doi:10.1007/s11747-024-01055-4
Qian, Y., Koschmann, A., and Xie, H. (2024). A practical guide to endogeneity correction using copulas. NBER Working Paper No. w32231. doi:10.2139/ssrn.4754776
Other copula-based methods:
copulaCorrection()
#------------------------------------------------------------------------
# Example 1: Single endogenous regressor with continuous
# and normal exogenous regressor (Haschka 2025, Section 4.1 Scenario 1)
# True values: alpha = 1 (P), beta = 1 (X), no intercept
#------------------------------------------------------------------------
data(dataCopIMAContExo)
res <- copulaIMA(
y ~ X + P - 1 | continuous(P),
data = dataCopIMAContExo,
cdf = "adj.ecdf",
num.boots = 1000
)
#> Fitting Haschka's copula-based IMA model for 1 continuous endogenous regressor(s).
#> Running 1000 bootstraps.
#>
|
| | 0%
|
| | 1%
|
|= | 1%
|
|= | 2%
|
|== | 2%
|
|== | 3%
|
|== | 4%
|
|=== | 4%
|
|=== | 5%
|
|==== | 5%
|
|==== | 6%
|
|===== | 6%
|
|===== | 7%
|
|===== | 8%
|
|====== | 8%
|
|====== | 9%
|
|======= | 9%
|
|======= | 10%
|
|======= | 11%
|
|======== | 11%
|
|======== | 12%
|
|========= | 12%
|
|========= | 13%
|
|========= | 14%
|
|========== | 14%
|
|========== | 15%
|
|=========== | 15%
|
|=========== | 16%
|
|============ | 16%
|
|============ | 17%
|
|============ | 18%
|
|============= | 18%
|
|============= | 19%
|
|============== | 19%
|
|============== | 20%
|
|============== | 21%
|
|=============== | 21%
|
|=============== | 22%
|
|================ | 22%
|
|================ | 23%
|
|================ | 24%
|
|================= | 24%
|
|================= | 25%
|
|================== | 25%
|
|================== | 26%
|
|=================== | 26%
|
|=================== | 27%
|
|=================== | 28%
|
|==================== | 28%
|
|==================== | 29%
|
|===================== | 29%
|
|===================== | 30%
|
|===================== | 31%
|
|====================== | 31%
|
|====================== | 32%
|
|======================= | 32%
|
|======================= | 33%
|
|======================= | 34%
|
|======================== | 34%
|
|======================== | 35%
|
|========================= | 35%
|
|========================= | 36%
|
|========================== | 36%
|
|========================== | 37%
|
|========================== | 38%
|
|=========================== | 38%
|
|=========================== | 39%
|
|============================ | 39%
|
|============================ | 40%
|
|============================ | 41%
|
|============================= | 41%
|
|============================= | 42%
|
|============================== | 42%
|
|============================== | 43%
|
|============================== | 44%
|
|=============================== | 44%
|
|=============================== | 45%
|
|================================ | 45%
|
|================================ | 46%
|
|================================= | 46%
|
|================================= | 47%
|
|================================= | 48%
|
|================================== | 48%
|
|================================== | 49%
|
|=================================== | 49%
|
|=================================== | 50%
|
|=================================== | 51%
|
|==================================== | 51%
|
|==================================== | 52%
|
|===================================== | 52%
|
|===================================== | 53%
|
|===================================== | 54%
|
|====================================== | 54%
|
|====================================== | 55%
|
|======================================= | 55%
|
|======================================= | 56%
|
|======================================== | 56%
|
|======================================== | 57%
|
|======================================== | 58%
|
|========================================= | 58%
|
|========================================= | 59%
|
|========================================== | 59%
|
|========================================== | 60%
|
|========================================== | 61%
|
|=========================================== | 61%
|
|=========================================== | 62%
|
|============================================ | 62%
|
|============================================ | 63%
|
|============================================ | 64%
|
|============================================= | 64%
|
|============================================= | 65%
|
|============================================== | 65%
|
|============================================== | 66%
|
|=============================================== | 66%
|
|=============================================== | 67%
|
|=============================================== | 68%
|
|================================================ | 68%
|
|================================================ | 69%
|
|================================================= | 69%
|
|================================================= | 70%
|
|================================================= | 71%
|
|================================================== | 71%
|
|================================================== | 72%
|
|=================================================== | 72%
|
|=================================================== | 73%
|
|=================================================== | 74%
|
|==================================================== | 74%
|
|==================================================== | 75%
|
|===================================================== | 75%
|
|===================================================== | 76%
|
|====================================================== | 76%
|
|====================================================== | 77%
|
|====================================================== | 78%
|
|======================================================= | 78%
|
|======================================================= | 79%
|
|======================================================== | 79%
|
|======================================================== | 80%
|
|======================================================== | 81%
|
|========================================================= | 81%
|
|========================================================= | 82%
|
|========================================================== | 82%
|
|========================================================== | 83%
|
|========================================================== | 84%
|
|=========================================================== | 84%
|
|=========================================================== | 85%
|
|============================================================ | 85%
|
|============================================================ | 86%
|
|============================================================= | 86%
|
|============================================================= | 87%
|
|============================================================= | 88%
|
|============================================================== | 88%
|
|============================================================== | 89%
|
|=============================================================== | 89%
|
|=============================================================== | 90%
|
|=============================================================== | 91%
|
|================================================================ | 91%
|
|================================================================ | 92%
|
|================================================================= | 92%
|
|================================================================= | 93%
|
|================================================================= | 94%
|
|================================================================== | 94%
|
|================================================================== | 95%
|
|=================================================================== | 95%
|
|=================================================================== | 96%
|
|==================================================================== | 96%
|
|==================================================================== | 97%
|
|==================================================================== | 98%
|
|===================================================================== | 98%
|
|===================================================================== | 99%
|
|======================================================================| 99%
|
|======================================================================| 100%
summary(res)
#>
#> Call:
#> copulaIMA(formula = y ~ X + P - 1 | continuous(P), data = dataCopIMAContExo,
#> cdf = "adj.ecdf", num.boots = 1000)
#>
#> Coefficients:
#> Point Estimate Boots SE Lower Boots CI (95%) Upper Boots CI (95%)
#> X 0.96636 0.03188 0.90603 1.03104
#> P 1.05741 0.04575 0.96998 1.14269
#> P_cop 0.61766 0.03563 0.54597 0.68712
#> Number of bootstraps: 1000
#> Num bootstraps attempted: 1000
#> Num bootstraps failed and discarded: 0 (0%)
#>
#> Continuous endogenous variables: P
#> Used cdf: adj.ecdf
#------------------------------------------------------------------------
# Example 2: Two endogenous regressors with intercept and no exogenous regressor
# True values: mu=10, alpha1 = 1 (P1), alpha2 = 1 (P2)
# Extension of the first example
#------------------------------------------------------------------------
data("dataCopIMAMultiEndo")
res2 <- copulaIMA(
# Alternative: y ~ P1 + P2 | continuous(P1, P2)
y ~ P1 + P2 | continuous(P1) + continuous(P2),
data = dataCopIMAMultiEndo,
cdf = "adj.ecdf",
num.boots = 1000
)
#> Fitting Haschka's copula-based IMA model for 2 continuous endogenous regressor(s).
#> Running 1000 bootstraps.
#>
|
| | 0%
|
| | 1%
|
|= | 1%
|
|= | 2%
|
|== | 2%
|
|== | 3%
|
|== | 4%
|
|=== | 4%
|
|=== | 5%
|
|==== | 5%
|
|==== | 6%
|
|===== | 6%
|
|===== | 7%
|
|===== | 8%
|
|====== | 8%
|
|====== | 9%
|
|======= | 9%
|
|======= | 10%
|
|======= | 11%
|
|======== | 11%
|
|======== | 12%
|
|========= | 12%
|
|========= | 13%
|
|========= | 14%
|
|========== | 14%
|
|========== | 15%
|
|=========== | 15%
|
|=========== | 16%
|
|============ | 16%
|
|============ | 17%
|
|============ | 18%
|
|============= | 18%
|
|============= | 19%
|
|============== | 19%
|
|============== | 20%
|
|============== | 21%
|
|=============== | 21%
|
|=============== | 22%
|
|================ | 22%
|
|================ | 23%
|
|================ | 24%
|
|================= | 24%
|
|================= | 25%
|
|================== | 25%
|
|================== | 26%
|
|=================== | 26%
|
|=================== | 27%
|
|=================== | 28%
|
|==================== | 28%
|
|==================== | 29%
|
|===================== | 29%
|
|===================== | 30%
|
|===================== | 31%
|
|====================== | 31%
|
|====================== | 32%
|
|======================= | 32%
|
|======================= | 33%
|
|======================= | 34%
|
|======================== | 34%
|
|======================== | 35%
|
|========================= | 35%
|
|========================= | 36%
|
|========================== | 36%
|
|========================== | 37%
|
|========================== | 38%
|
|=========================== | 38%
|
|=========================== | 39%
|
|============================ | 39%
|
|============================ | 40%
|
|============================ | 41%
|
|============================= | 41%
|
|============================= | 42%
|
|============================== | 42%
|
|============================== | 43%
|
|============================== | 44%
|
|=============================== | 44%
|
|=============================== | 45%
|
|================================ | 45%
|
|================================ | 46%
|
|================================= | 46%
|
|================================= | 47%
|
|================================= | 48%
|
|================================== | 48%
|
|================================== | 49%
|
|=================================== | 49%
|
|=================================== | 50%
|
|=================================== | 51%
|
|==================================== | 51%
|
|==================================== | 52%
|
|===================================== | 52%
|
|===================================== | 53%
|
|===================================== | 54%
|
|====================================== | 54%
|
|====================================== | 55%
|
|======================================= | 55%
|
|======================================= | 56%
|
|======================================== | 56%
|
|======================================== | 57%
|
|======================================== | 58%
|
|========================================= | 58%
|
|========================================= | 59%
|
|========================================== | 59%
|
|========================================== | 60%
|
|========================================== | 61%
|
|=========================================== | 61%
|
|=========================================== | 62%
|
|============================================ | 62%
|
|============================================ | 63%
|
|============================================ | 64%
|
|============================================= | 64%
|
|============================================= | 65%
|
|============================================== | 65%
|
|============================================== | 66%
|
|=============================================== | 66%
|
|=============================================== | 67%
|
|=============================================== | 68%
|
|================================================ | 68%
|
|================================================ | 69%
|
|================================================= | 69%
|
|================================================= | 70%
|
|================================================= | 71%
|
|================================================== | 71%
|
|================================================== | 72%
|
|=================================================== | 72%
|
|=================================================== | 73%
|
|=================================================== | 74%
|
|==================================================== | 74%
|
|==================================================== | 75%
|
|===================================================== | 75%
|
|===================================================== | 76%
|
|====================================================== | 76%
|
|====================================================== | 77%
|
|====================================================== | 78%
|
|======================================================= | 78%
|
|======================================================= | 79%
|
|======================================================== | 79%
|
|======================================================== | 80%
|
|======================================================== | 81%
|
|========================================================= | 81%
|
|========================================================= | 82%
|
|========================================================== | 82%
|
|========================================================== | 83%
|
|========================================================== | 84%
|
|=========================================================== | 84%
|
|=========================================================== | 85%
|
|============================================================ | 85%
|
|============================================================ | 86%
|
|============================================================= | 86%
|
|============================================================= | 87%
|
|============================================================= | 88%
|
|============================================================== | 88%
|
|============================================================== | 89%
|
|=============================================================== | 89%
|
|=============================================================== | 90%
|
|=============================================================== | 91%
|
|================================================================ | 91%
|
|================================================================ | 92%
|
|================================================================= | 92%
|
|================================================================= | 93%
|
|================================================================= | 94%
|
|================================================================== | 94%
|
|================================================================== | 95%
|
|=================================================================== | 95%
|
|=================================================================== | 96%
|
|==================================================================== | 96%
|
|==================================================================== | 97%
|
|==================================================================== | 98%
|
|===================================================================== | 98%
|
|===================================================================== | 99%
|
|======================================================================| 99%
|
|======================================================================| 100%
summary(res2)
#>
#> Call:
#> copulaIMA(formula = y ~ P1 + P2 | continuous(P1) + continuous(P2),
#> data = dataCopIMAMultiEndo, cdf = "adj.ecdf", num.boots = 1000)
#>
#> Coefficients:
#> Point Estimate Boots SE Lower Boots CI (95%) Upper Boots CI (95%)
#> (Intercept) 9.835738 0.600372 8.680275 11.090617
#> P1 1.417652 0.424032 0.573655 2.270830
#> P2 0.715764 0.385053 0.007068 1.518044
#> P1_cop 0.258177 0.117059 0.019112 0.491788
#> P2_cop 0.497003 0.110308 0.263097 0.700726
#> Number of bootstraps: 1000
#> Num bootstraps attempted: 1000
#> Num bootstraps failed and discarded: 0 (0%)
#>
#> Continuous endogenous variables: P1, P2
#> Used cdf: adj.ecdf
#------------------------------------------------------------------------
# Example 3: Single endogenous regressor with binary exogenous regressor
# (Haschka 2025, section 4.1 scenario 2)
# This example shows one of the key example of IMA method, i.e., the exogenous
# regressor does not need to be continuous or normally distributed.
# X the exogenous regressor is binary (0 or 1).
# True values: alpha = 1 (P), beta = 1 (X), no intercept.
#------------------------------------------------------------------------
data("dataCopIMABinExo")
res3 <- copulaIMA(
y ~ X + P - 1 | continuous(P),
data = dataCopIMABinExo,
cdf = "adj.ecdf",
num.boots = 1000
)
#> Fitting Haschka's copula-based IMA model for 1 continuous endogenous regressor(s).
#> Running 1000 bootstraps.
#>
|
| | 0%
|
| | 1%
|
|= | 1%
|
|= | 2%
|
|== | 2%
|
|== | 3%
|
|== | 4%
|
|=== | 4%
|
|=== | 5%
|
|==== | 5%
|
|==== | 6%
|
|===== | 6%
|
|===== | 7%
|
|===== | 8%
|
|====== | 8%
|
|====== | 9%
|
|======= | 9%
|
|======= | 10%
|
|======= | 11%
|
|======== | 11%
|
|======== | 12%
|
|========= | 12%
|
|========= | 13%
|
|========= | 14%
|
|========== | 14%
|
|========== | 15%
|
|=========== | 15%
|
|=========== | 16%
|
|============ | 16%
|
|============ | 17%
|
|============ | 18%
|
|============= | 18%
|
|============= | 19%
|
|============== | 19%
|
|============== | 20%
|
|============== | 21%
|
|=============== | 21%
|
|=============== | 22%
|
|================ | 22%
|
|================ | 23%
|
|================ | 24%
|
|================= | 24%
|
|================= | 25%
|
|================== | 25%
|
|================== | 26%
|
|=================== | 26%
|
|=================== | 27%
|
|=================== | 28%
|
|==================== | 28%
|
|==================== | 29%
|
|===================== | 29%
|
|===================== | 30%
|
|===================== | 31%
|
|====================== | 31%
|
|====================== | 32%
|
|======================= | 32%
|
|======================= | 33%
|
|======================= | 34%
|
|======================== | 34%
|
|======================== | 35%
|
|========================= | 35%
|
|========================= | 36%
|
|========================== | 36%
|
|========================== | 37%
|
|========================== | 38%
|
|=========================== | 38%
|
|=========================== | 39%
|
|============================ | 39%
|
|============================ | 40%
|
|============================ | 41%
|
|============================= | 41%
|
|============================= | 42%
|
|============================== | 42%
|
|============================== | 43%
|
|============================== | 44%
|
|=============================== | 44%
|
|=============================== | 45%
|
|================================ | 45%
|
|================================ | 46%
|
|================================= | 46%
|
|================================= | 47%
|
|================================= | 48%
|
|================================== | 48%
|
|================================== | 49%
|
|=================================== | 49%
|
|=================================== | 50%
|
|=================================== | 51%
|
|==================================== | 51%
|
|==================================== | 52%
|
|===================================== | 52%
|
|===================================== | 53%
|
|===================================== | 54%
|
|====================================== | 54%
|
|====================================== | 55%
|
|======================================= | 55%
|
|======================================= | 56%
|
|======================================== | 56%
|
|======================================== | 57%
|
|======================================== | 58%
|
|========================================= | 58%
|
|========================================= | 59%
|
|========================================== | 59%
|
|========================================== | 60%
|
|========================================== | 61%
|
|=========================================== | 61%
|
|=========================================== | 62%
|
|============================================ | 62%
|
|============================================ | 63%
|
|============================================ | 64%
|
|============================================= | 64%
|
|============================================= | 65%
|
|============================================== | 65%
|
|============================================== | 66%
|
|=============================================== | 66%
|
|=============================================== | 67%
|
|=============================================== | 68%
|
|================================================ | 68%
|
|================================================ | 69%
|
|================================================= | 69%
|
|================================================= | 70%
|
|================================================= | 71%
|
|================================================== | 71%
|
|================================================== | 72%
|
|=================================================== | 72%
|
|=================================================== | 73%
|
|=================================================== | 74%
|
|==================================================== | 74%
|
|==================================================== | 75%
|
|===================================================== | 75%
|
|===================================================== | 76%
|
|====================================================== | 76%
|
|====================================================== | 77%
|
|====================================================== | 78%
|
|======================================================= | 78%
|
|======================================================= | 79%
|
|======================================================== | 79%
|
|======================================================== | 80%
|
|======================================================== | 81%
|
|========================================================= | 81%
|
|========================================================= | 82%
|
|========================================================== | 82%
|
|========================================================== | 83%
|
|========================================================== | 84%
|
|=========================================================== | 84%
|
|=========================================================== | 85%
|
|============================================================ | 85%
|
|============================================================ | 86%
|
|============================================================= | 86%
|
|============================================================= | 87%
|
|============================================================= | 88%
|
|============================================================== | 88%
|
|============================================================== | 89%
|
|=============================================================== | 89%
|
|=============================================================== | 90%
|
|=============================================================== | 91%
|
|================================================================ | 91%
|
|================================================================ | 92%
|
|================================================================= | 92%
|
|================================================================= | 93%
|
|================================================================= | 94%
|
|================================================================== | 94%
|
|================================================================== | 95%
|
|=================================================================== | 95%
|
|=================================================================== | 96%
|
|==================================================================== | 96%
|
|==================================================================== | 97%
|
|==================================================================== | 98%
|
|===================================================================== | 98%
|
|===================================================================== | 99%
|
|======================================================================| 99%
|
|======================================================================| 100%
summary(res3)
#>
#> Call:
#> copulaIMA(formula = y ~ X + P - 1 | continuous(P), data = dataCopIMABinExo,
#> cdf = "adj.ecdf", num.boots = 1000)
#>
#> Coefficients:
#> Point Estimate Boots SE Lower Boots CI (95%) Upper Boots CI (95%)
#> X 0.89119 0.06468 0.76524 1.01695
#> P 1.07922 0.04596 0.98358 1.16870
#> P_cop 0.52237 0.03281 0.45512 0.58199
#> Number of bootstraps: 1000
#> Num bootstraps attempted: 1000
#> Num bootstraps failed and discarded: 0 (0%)
#>
#> Continuous endogenous variables: P
#> Used cdf: adj.ecdf