compute the difference between two estimates of the LIGE
Diff.lige(cfaobj1, cfaobj2, se = T, h)
the first CFA object
the second CFA object
boolean whether or not to compute standard errors
a bandwidth
a CFASE object
# NOT RUN {
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))
Diff.lige(out$cfa1, out$cfa2, h=0.5)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab