We consider the problem of estimating the expected causal effect \(E[Y|do(X)]\) for a target variable \(Y\) when treatment \(X\) is set by intervention, focusing on continuous random variables. In settings without selection bias or confounding, \(E[Y|do(X)] = E[Y|X]\) , which can be estimated using standard regression methods. However, regression fails when systematic missingness induced by selection bias, or confounding distorts the data. Proxy variables unaffected by the selection process can, under certain constraints, be used to correct for selection bias to recover \(E[Y|X]\) , and hence \(E[Y|do(X)]\) , reliably. When data is additionally affected by confounding, recovering the causal effect from selection-biased data is more challenging and requires access to proxies to both correct for confounding and for the selection mechanism. Assuming access to such proxies from external unbiased observational data, we derive theoretical conditions ensuring identifiability and recoverability of causal effects. We further introduce a linear two-step regression estimator (TSR), which can be extended through adding non-linear basis functions, capable of exploiting proxy variables to adjust for selection bias while accounting for confounding. We show that TSR is consistent with previous estimators when confounding is absent, but achieves a lower variance. Extensive simulation studies validate TSR’s correctness for scenarios that include both selection bias and confounding with proxy variables.