Cem July 2020

cem@lists.gking.harvard.edu

3 participants
3 discussions

by Blom, drs. M. - WODC/CRS

Hi All, I can't get why the analyzes below, performed in Stata / R, do not produce the same result. Please provide your suggestions for the explanation for these differences. Thanks, Martine PS (how) can I display the actual cutpoints in Stata (imb)? -- STATA . import delimited "G:\Recidive\Projecten\Verkeer\Algemeen\Artikel effectiviteit ASP tijdens\Data\Data voor R.csv", clear (8 vars, 8,627 obs) . imb d_sekse lftbegrec d_addinfo_ind_best lftinsz1 vgalgexuz vgverkexuz vgrijoiuz, treatment (treat) (using the scott break method for L1 distance) Multivariate L1 distance: .68243298 Univariate imbalance: L1 mean min 25% 50% 75% max d_sekse .01557 -.01557 0 0 0 0 0 lftbegrec .0687 -1.5626 0 -1 -2 -1 -5 d_addinfo_ind_best .00497 -.00497 0 0 0 0 0 lftinsz1 .08852 .91045 0 1 1 2 -12 vgalgexuz .14912 -1.2611 0 0 -1 -1 78 vgverkexuz .09852 -.41727 0 0 0 -1 4 vgrijoiuz .11859 -.36406 0 0 0 0 -4 --R > data <- read.csv(file = "G:\\Recidive\\Projecten\\Verkeer\\Algemeen\\Artikel effectiviteit ASP tijdens\\Data\\Data voor R.csv") > cov <- c("D_SEKSE", "LFTBEGREC", "D_ADDINFO_IND_BEST","LFTINSZ1", "VGALGEXUZ", "VGVERKEXUZ", "VGRIJOIUZ") > imb <- imbalance(group = data$TREAT, data = data[cov]) > imb Multivariate Imbalance Measure: L1=0.665 Percentage of local common support: LCS=11.2% Univariate Imbalance Measures: statistic type L1 min 25% 50% 75% max D_SEKSE -0.015571508 (diff) 0.015571508 0 0 0 0 0 LFTBEGREC -1.562561051 (diff) 0.025019213 0 -1 -2 -1 -5 D_ADDINFO_IND_BEST -0.004970458 (diff) 0.004970458 0 0 0 0 0 LFTINSZ1 0.910454843 (diff) 0.045867098 0 1 1 2 -12 VGALGEXUZ -1.261056065 (diff) 0.146272471 0 0 -1 -1 78 VGVERKEXUZ -0.417271627 (diff) 0.097957912 0 0 0 -1 4 VGRIJOIUZ -0.364063976 (diff) 0.118028583 0 0 0 0 -4 ________________________________ Dit bericht kan informatie bevatten die niet voor u is bestemd. Indien u niet de geadresseerde bent of dit bericht abusievelijk aan u is toegezonden, wordt u verzocht dat aan de afzender te melden en het bericht te verwijderen. De Staat aanvaardt geen aansprakelijkheid voor schade, van welke aard ook, die verband houdt met risico's verbonden aan het elektronisch verzenden van berichten. Ministerie van Justitie en Veiligheid This message may contain information that is not intended for you. If you are not the addressee or if this message was sent to you by mistake, you are requested to inform the sender and delete the message. The State accepts no liability for damage of any kind resulting from the risks inherent in the electronic transmission of messages. Ministry of Justice and Security

3 years, 8 months

Re: [cem] Cem Digest, Vol 128, Issue 1

by Oladele Campbell

Dear Victoria I have used CEM in my PhD research to realise 1:1 matched samples as shown in the attached document. So, it is possible. However, I did not use the R package but rather the SPSS add-on. I downloaded the software and the manual for installing the add-on online. I searched for the SPSS software extension bundle on my system, and i am sending a copy. But unfortunately, I could not find the manual. But if you are interested in using it, you can find it online. Best wishes Oladele Campbell On Thu, Jul 9, 2020 at 5:00 PM <cem-request(a)lists.gking.harvard.edu> wrote: > Send Cem mailing list submissions to > cem(a)lists.gking.harvard.edu > > To subscribe or unsubscribe via the World Wide Web, visit > https://lists.gking.harvard.edu/mailman/listinfo/cem > or, via email, send a message with subject or body 'help' to > cem-request(a)lists.gking.harvard.edu > > You can reach the person managing the list at > cem-owner(a)lists.gking.harvard.edu > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of Cem digest..." > > > Today's Topics: > > 1. How to get matched pairs from CEM? (Victoria Graham) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Thu, 9 Jul 2020 03:10:09 +0000 > From: Victoria Graham <victoria.graham(a)mq.edu.au> > To: "cem(a)lists.gking.harvard.edu" <cem(a)lists.gking.harvard.edu> > Subject: [cem] How to get matched pairs from CEM? > Message-ID: > < > SY2PR01MB241067CD2E5D13825698283AA5640(a)SY2PR01MB2410.ausprd01.prod.outlook.com > > > > Content-Type: text/plain; charset="utf-8" > > Hi, > > I am using Coursened Exact Matching to improve balance in my sample. The > summary matching results after applying CEM are that 11,128 pixels have > been matched: > > Sample sizes: > > Control Treated > > All 19741 8283 > > Matched 6728 4400 > > Unmatched 13013 3883 > > Discarded 0 0 > > Firstly, is it possible to specify 1:1 matching with CEM? I want to only > retain matched pairs, which requires the same amount of pixels in my > treatment & control groups. > > I am also having trouble getting a unique identifier for each matched pair > so that I can calculate the difference between the treatment and control. I > am using the code: > > pair: Produces a paired sample out of a CEM match solution > > When I execute the below command, it says there are 0 matches. > > > # we want a set of paired units > > > pairsample <- pair(match.out_cem, data) > > > > Total number of units paired in CEM strata: 0 > > Total number of units matched: 0 > > Can anyone help or point me in the right direction? > > Thanks, Tori > > Victoria Graham | MPhil., PhD candidate > Macquarie University NSW 2109 Australia > T: +61 2 9850 6298 | E: victoria.graham(a)mq.edu.au<mailto: > victoria.graham(a)mq.edu.au> > Twitter:@ToriFGraham<https://twitter.com/ToriFGraham> > > [cid:image001.png@01D5F150.30032D90] > >

3 years, 9 months

How to get matched pairs from CEM?

by Victoria Graham

Hi, I am using Coursened Exact Matching to improve balance in my sample. The summary matching results after applying CEM are that 11,128 pixels have been matched: Sample sizes: Control Treated All 19741 8283 Matched 6728 4400 Unmatched 13013 3883 Discarded 0 0 Firstly, is it possible to specify 1:1 matching with CEM? I want to only retain matched pairs, which requires the same amount of pixels in my treatment & control groups. I am also having trouble getting a unique identifier for each matched pair so that I can calculate the difference between the treatment and control. I am using the code: pair: Produces a paired sample out of a CEM match solution When I execute the below command, it says there are 0 matches. > # we want a set of paired units > pairsample <- pair(match.out_cem, data) Total number of units paired in CEM strata: 0 Total number of units matched: 0 Can anyone help or point me in the right direction? Thanks, Tori Victoria Graham | MPhil., PhD candidate Macquarie University NSW 2109 Australia T: +61 2 9850 6298 | E: victoria.graham(a)mq.edu.au<mailto:victoria.graham@mq.edu.au> Twitter:@ToriFGraham<https://twitter.com/ToriFGraham> [cid:image001.png@01D5F150.30032D90]

3 years, 9 months

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

Cem July 2020