I am using cem to match treatment and control group members from 5 multiply
imputed datasets from Amelia. I've encountered two types of problems.
First, when I run a logit model after calling cem on the multiply imputed
datasets:
> run <- att(new1, Pass ~ SEQ + O_GPA + O_HRS_EARNED + AGE + PELL1 +
LOG_CA_SCORE + LOG_CM_SCORE, data = imputed, model = "logit")
I get the following warning messages, one for each multiply imputed
dataset:
Warning messages:
1: In eval(expr, envir, enclos) :
non-integer #successes in a binomial glm!
2: In eval(expr, envir, enclos) :
non-integer #successes in a binomial glm!
3: In eval(expr, envir, enclos) :
non-integer #successes in a binomial glm!
4: In eval(expr, envir, enclos) :
non-integer #successes in a binomial glm!
5: In eval(expr, envir, enclos) :
non-integer #successes in a binomial glm!
My dependent variable had no missing cases, so the values for it were not
imputed for any observations, and I have checked that the values of it are
either 0 or 1. Any ideas why I am getting this warning? I've noticed from
searching online that this warning arises when weights are included in a
binomial glm model, so I wonder if it might have something to do with how
cem is weighting the control observations.
(2) The second issue relates to displaying the output of the model. When I
type:
> run
I get:
Logistic model on CEM matched data:
SATT point estimate: 1.509872 (p.value=0.001640)
95% conf. interval: [0.570028, 2.449717]
which is fine, but when I type:
> summary(run)
I get:
Treatment effect estimation for data:
NULL
Logistic model estimated on matched data only
Coefficients:
Error in symnum(pv, corr = FALSE, na = FALSE, cutpoints = c(0, 0.001, :
'x' must be between 0 and 1
Any idea why I'm getting this error?
Thanks for any help you can provide,
Bill
Hello,
I am wondering how I can get a composite weight from cem_weights and the sampling weight (the inverse of the probability of a unit to be selected) from my survey data.
Shall I simply use cem_weight*sampling_weight as a new sampling weight? If not, is there any proper way to deal with it?
I am using STATA. cem_weight is supposed to be used as "iweight" while sampling weight is supposed to be used as "pweight". I am worried about the compatibility issue to use cem_weight*sampling_weight as a weight.
Thanks.
Best.
Sifan Zhou
Ph.D. student
Economics Department
University at Albany, SUNY
Hello all,
I am attempting to implement the coarsened exact matching technique
described in Iacus, King, and Porro (2011). The "treatment" in this
case has three categories. It is clear from section 5.1 of the paper
that this is not a problem for the technique generally. The question I
have is how to calculate the weights, or whether the weights are
strictly necessary. If anyone can you point me in the direction of
resources that might illuminate this issue I would be grateful.
Best,
Rose
--
Rose Anne Medeiros, Ph.D.
Quantitative Methodologist
Department of Sociology
Rice University
(713)348-4397
Hi,
I'm get an error message when I run imb in Stata (cem runs fine):
L1meas(): 3301 subscript invalid
imbalance(): - function returned error
<istmt>: - function returned error
r(3301);
My code is something like:
imb new_age borninbeijing female hometown01 BJ_hukou Rural_hukoutype
fatheredu motheredu ifpreschooling schoolingathome,
treatment(public_yn)
Any help would be most appreciated!
Thanks,
Prashant
There is also a 'pair' function in the cem package, which I think produces a paired sample out of a CEM match solution.
However, I realized that pair function pairs up instances that are not necessarily matched.
So I modified the pair.r file in the cem package, and force the pair function to pair up instances that are matched already.
Below is what I diid and hope it helps:
pair <- function(obj, data, method=NULL, mpower=2, verbose=0){
nm <- NULL
strataID <- unique(obj$strata)
reservoir <- NULL
for(i in strataID){
####################################################
###########added by hoktay
########use only the instances that are matched
#### used to be
######### idx <- which(obj$strata==i )
####################################################
idx <- which(obj$strata==i & obj$matched==TRUE)
.....
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Sun, 12 Feb 2012 23:44:05 -0800
> From: "Ashkan Moazzez, MD" <ashkanmoazzez(a)ucla.edu>
> To: cem(a)lists.gking.harvard.edu
> Subject: [cem] CEM
> Message-ID: <20120212234405.79195ckdlsb7deth(a)mail.ucla.edu>
> Content-Type: text/plain; charset=ISO-8859-1; DelSp="Yes";
> format="flowed"
>
> Dear members of the CEM list,
> I was wondering if there is a way to extract the data after you match
> them with CEM? Is there an equivalent command like match.data in
> matchit() for CEM?
>
> Best regards,
> Ashkan
>
>
> Ashkan Moazzez, MD
>
>
> ------------------------------
>
Dear members of the CEM list,
I was wondering if there is a way to extract the data after you match
them with CEM? Is there an equivalent command like match.data in
matchit() for CEM?
Best regards,
Ashkan
Ashkan Moazzez, MD
Hello all,
I am using spacegraph technology for comparing PSM and CEM. I am basically using 2 different implementations of PSM (one in the CEM package, and the other is in the MatchIt package).
I realized that there are many parameters for PSM that can be turned on and off. To be more concrete below is a short list of parameters and I have some preliminary results indicating that these parameters change the performance of the PSM matching (i.e., the L1 imbalance value):
--match with/without replacement
--breaking ties when the distance is same for multiple control instances
--matching ordering
My question is that, given there are many knobs especially for PSM, what are the right settings for these parameters to fairly compare the performance against the performance of CEM?
Thanks.
Huseyin Oktay
----------------------------------------------------
Ph.D. Candidate, Knowledge Discovery Laboratory
Department of Computer Science,
University of Massachusetts Amherst
140 Governors Dr.
Amherst, MA 01003
http://www.cs.umass.edu/~hoktay/http://kdl.cs.umass.edu/
Hi Marco,
One thing to make sure is that you have run each of the installers using the "Run As Administrator" option in Windows (of course, it might be in Italian on your system). In addition, you should also run SPSS as an administrator when you install the extension command. Then you must close SPSS and reopen it.
If you are still having problems, please send me mail and we can try to resolve the problem.
Cheers,
matt.
~~~~~~~~~~~
Matthew Blackwell
PhD Candidate
Institute for Quantitative Social Science
Department of Government
Harvard University
url: http://www.mattblackwell.org
On Monday, October 24, 2011 at 8:24 AM, Matt Blackwell wrote:
> Hi Marco,
>
> One thing to make sure is that you have run each of the installers using the "Run As Administrator" option in Windows (of course, it might be in Italian on your system). In addition, you should also run SPSS as an administrator when you install the extension command. Then you must close SPSS and reopen it.
>
> If you are still having problems, please send me mail and we can try to resolve the problem.
>
> Cheers,
> matt.
>
> ~~~~~~~~~~~
> Matthew Blackwell
> PhD Candidate
> Institute for Quantitative Social Science
> Department of Government
> Harvard University
> url: http://www.mattblackwell.org
>
>
> On Monday, October 24, 2011 at 5:57 AM, Cantalupi Marco wrote:
>
> >
> >
> >
> >
> >
> > Da: Cantalupi Marco
> > Inviato: lunedì 24 ottobre 2011 11.50
> > A: 'cem(a)lists.gking.harvard.edu (mailto:cem@lists.gking.harvard.edu)'
> > Oggetto: cem for spps
> >
> >
> >
> >
> >
> >
> > Hi,
> >
> >
> > the installation of all the plugins goes well till the very end, but apparently the module doesn’t run properly.
> >
> >
> > I hereby paste the spps output after having launched cem via the dialog menu described in “How to use CEM for spss”.
> >
> >
> > Any suggestions?
> >
> >
> > Kind regards
> >
> >
> > M. Cantalupi
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> > GET
> >
> >
> > FILE='C:\Documents and Settings\600481\Desktop\SPPS_CEM\lalonde.sav'.
> >
> >
> >
> >
> >
> > Avviso. Nome comando: GET FILE
> >
> >
> > Il file di sistema PASW Statistics "C:\Documents and Settings\600481\Desktop\SPPS_CEM\lalonde.sav" è scritto in una codifica di caratteri
> >
> >
> > (ISO_8859-1:1987) non compatibile con il valore LOCALE corrente. È possibile che non sia
> >
> >
> > leggibile. Si consiglia di modificare il valore LOCALE o di impostare UNICODE
> >
> >
> > su ON. (DATA 1721)
> >
> >
> > DATASET NAME InsiemeDati1 WINDOW=FRONT.
> >
> >
> > CEM TREATMENT=treated VARIABLES=education black married nodegree re74 re75.
> >
> >
> > Extension command CEM could not be loaded. The module or a module that it requires may be missing, or there may be syntax errors in it.
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> > Agenzia Regionale del Lavoro
> >
> >
> > Marco Cantalupi
> >
> >
> > Ricercatore - Area Mercato del Lavoro
> >
> >
> > via San Francesco, 37 - 34133 Trieste
> >
> >
> > Tel.: 040-377 51 36
> >
> >
> > Cell.: 335 74 34 969
> >
> >
> > Fax: 040-377 51 97
> >
> >
> > e-mail: marco.cantalupi(a)agelav.fvg.it (mailto:amarco.cantalupi@agelav.fvg.it)
> >
> >
> > www.agenzialavoro.regione.fvg.it (http://www.agenzialavoro.regione.fvg.it)
> >
> >
> > skype: marco.cantalupi.agelav
> >
> >
> > P Save a tree - please do not print this email unless you really need to
> >
> > AVVISO DI RISERVATEZZA
> > “Informazioni riservate possono essere contenute nel messaggio o nei suoi allegati. Se non siete i destinatari indicati nel messaggio, o responsabili per la sua consegna alla persona, o se avete ricevuto il messaggio per errore, non dovete trascriverlo, copiarlo o inviarlo a nessuno. In tal caso, dovete cancellare/distruggere il messaggio ed i suoi allegati. Grazie.”
> > CONFIDENTIALITY NOTICE
> > “Confidential information may be contained in this message or in its attachments. If you are not the addressee indicated in this message, or responsible for message delivering to that person, or if you have received this message in error, you may not transcribe, copy or deliver this message to anyone. In that case, you should delete/destroy this message and its attachments. Thank you”.
> >
> >
> >
> >
> >
> >
> >
>
>
________________________________
Da: Cantalupi Marco
Inviato: lunedì 24 ottobre 2011 11.50
A: 'cem(a)lists.gking.harvard.edu'
Oggetto: cem for spps
Hi,
the installation of all the plugins goes well till the very end, but apparently the module doesn't run properly.
I hereby paste the spps output after having launched cem via the dialog menu described in "How to use CEM for spss".
Any suggestions?
Kind regards
M. Cantalupi
GET
FILE='C:\Documents and Settings\600481\Desktop\SPPS_CEM\lalonde.sav'.
Avviso. Nome comando: GET FILE
Il file di sistema PASW Statistics "C:\Documents and Settings\600481\Desktop\SPPS_CEM\lalonde.sav" è scritto in una codifica di caratteri
(ISO_8859-1:1987) non compatibile con il valore LOCALE corrente. È possibile che non sia
leggibile. Si consiglia di modificare il valore LOCALE o di impostare UNICODE
su ON. (DATA 1721)
DATASET NAME InsiemeDati1 WINDOW=FRONT.
CEM TREATMENT=treated VARIABLES=education black married nodegree re74 re75.
Extension command CEM could not be loaded. The module or a module that it requires may be missing, or there may be syntax errors in it.
Agenzia Regionale del Lavoro
Marco Cantalupi
Ricercatore - Area Mercato del Lavoro
via San Francesco, 37 - 34133 Trieste
Tel.: 040-377 51 36
Cell.: 335 74 34 969
Fax: 040-377 51 97
e-mail: marco.cantalupi(a)agelav.fvg.it<mailto:amarco.cantalupi@agelav.fvg.it>
www.agenzialavoro.regione.fvg.it<http://www.agenzialavoro.regione.fvg.it>
skype: marco.cantalupi.agelav
P Save a tree - please do not print this email unless you really need to
AVVISO DI RISERVATEZZA
"Informazioni riservate possono essere contenute nel messaggio o nei suoi allegati. Se non siete i destinatari indicati nel messaggio, o responsabili per la sua consegna alla persona, o se avete ricevuto il messaggio per errore, non dovete trascriverlo, copiarlo o inviarlo a nessuno. In tal caso, dovete cancellare/distruggere il messaggio ed i suoi allegati. Grazie."
CONFIDENTIALITY NOTICE
"Confidential information may be contained in this message or in its attachments. If you are not the addressee indicated in this message, or responsible for message delivering to that person, or if you have received this message in error, you may not transcribe, copy or deliver this message to anyone. In that case, you should delete/destroy this message and its attachments. Thank you".
