getResDiff.CFA

Get the difference between two CFA objects

getResDiff.CFA(cfaobj1, cfaobj2, fun, se = T, ...)

Arguments

cfaobj1	the first CFA object
cfaobj2	the second CFA object
fun	a function to apply for every value of the treatment in the cfaobj
se	whether or not to compute standard errors
...	can pass additional arguments to fun using this argument

Value

CFASE object

Examples

if (FALSE) {
data(igm)
tvals <- seq(10,12,length.out=8)
yvals <- seq(quantile(igm$lcfincome, .05), quantile(igm$lcfincome, .95), length.out=50)

## obtain counterfactual results
out <- cfa2(lcfincome ~ lfincome, tvals, yvals, igm, method1="qr",
xformla2=~HEDUC, method2="qr", iters=10, tau1=seq(.05,.95,.05),
tau2=seq(.05,.95,.05))

## get the difference between the average that adjusts for covariates and
## the one that does not
getResDiff.CFA(out$cfa1, out$cfa2, E)
}