Dear All,
Thank you for a great software package! I'm currently working on a project
where MI is not a viable option. I'm hoping to use pattern mixture models
to estimate propensity scores for cases and controls having missing data
for covariates. Is there then any way to then use MatchIt for matching and
computing balance statistics? Hoping that there might be a way to use
stratification or missing data indicators to do this (also considering
general location models-- just in case these would work better here). I'll
be tremendously appreciative of any thoughts or advice you may have on
this.
Also want to let you know that I can't open the documents that you show
links to in your responses (e.g
http://gking.harvard.edu/node/4355/rbuild_documentation/How_Do_I3.html.<http://gking.harvard.edu/node/4355/rbuild_documentation/How_Do_I3.html>Figured
out the syntax here but would be helpful to be able to see the docs
you reference. Is there anything I can do to access them?
Thank you!
Best wishes,
Sharon Simonton
University of Michigan
Ann Arbor, MI
Hi,
Hoping someone might be able to help me with the above warning message.
I have a dataframe with 7 covariates. 6 out of the 7 run without any issues, but any matching method run involving the 7th returns the following message;
"Warning message: glm.fit: fitted probabilities numerically 0 or 1 occurred"
I've googled the above and despite some fairly detailed explanations I'm none the wiser as to how to progress.
The variable is numerical. Sample as follows:
"[1] 42.108184 1.925875 42.108184 3.540687 2.331532 2.331532 2.781410 2.781410 1.701947 1.701947
[11] 1.701947 1.701947 1.701947 1.966572 1.966572 10.994199 0.000000 0.748349 3.666423 3.666423
[21] 3.666423 3.666423 0.746897 2.046640 1.019349 0.371999 1.019349 0.208058 0.194368 0.194368
[31] 0.952280 0.952280 0.952280 1.202849 1.202849 1.361906 1.361906 1.694844 0.618513 1.361906
[41] 1.361906 1.361906 0.304016 0.304016 0.641256 0.641256 0.641256 0.234019 0.641256 0.641256"
Any suggestions would be very much appreciated.
Cheers,
Anthony
Anthony Dancer
Imperial College London
Silwood Park Campus
Buckhurst Road
Ascot SL5 7PY
t: +44(0) 7968 836 451
e: anthony.dancer09(a)imperial.ac.uk
s: anthony.dancer1
Hi,
I have a question regarding how to stratify a matching problem in MatchIt. I've been using the package optmatch directly, to identify optimal full matches for a stratified dataset (i.e., I have treatment and control groups for two sites, and I'd like to match site-by-site). Extracting the matched data from optmatch is complex, so I've switched to calling optmatch through MatchIt, using the following code:
m1.out<-matchit(Tr~ HH_SIZE+DOM_ETH+HH_GENDER+ED_LEVEL+ED_LEVEL_MISSING+RESIDENT+HH_AGE+MARKET,data=covariates, method="full", distance="linear.logit", discard="hull.both", min.controls=1,max.controls=Inf)
Is there a way to incorporate the strata denoted by "MPAID" when running method='full' in MatchIt? I'd like to be able to view one set of balance statistics, but still enforce matching within sites.
Any advice would be greatly appreciated!
Thank you!
Louise