11 Jan
2006
11 Jan
'06
8:52 a.m.
thanks for your note. yes, relogit in Zelig will do everything you're
interested in and so you don't need to calculate C or the var-cov matrix
yourself; all the corrections are made for you.
the key is to recognize that there are 2 separate forms of bias
correction, one for case-control studies, for which you'd use the
case.correct option, and the other for rare events bias (which is
appropriate whether or not you use case control designs), for which you'd
use the bias.correct option.
Best of luck with your research,
Gary
---
Gary King,
Institute for Quantitative Social Science
Harvard University, 1737 Cambridge St, Cambridge, MA 02138
http://GKing.Harvard.Edu, King(a)Harvard.Edu
Direct 617-495-2027, Assistant 495-9271, eFax 812-8581
On Tue, 10 Jan 2006, Dave Armstrong wrote:
Dear Prof King,
I read your paper "Logistic Regression in Rare Events Data" and have just a
couple of quick questions about it. For consistency of notation, I've put
them in the attached PDF. If you would have some time to look at them or
point them and me in the right direction, it would be very much appreciated.
Thanks for your help,
Dave.
--
Dave Armstrong
University of Maryland
Dept of Government and Politics
3140 Tydings Hall
College Park, MD 20742
Office: 2103L Cole Field House
Phone: 301-405-9735
e-mail: darmstrong(a)gvpt.umd.edu
web: www.davearmstrong-ps.com
Facts are meaningless. You can use facts to prove anything that's even
remotely true. - Homer Simpson
To this day, philosophers suffer from Plato's disease: the assumption that
reality fundamentally consists of
abstract essences best described by words or geometry. (In truth, reality is
largely a probabilistic affair best
described by statistics) - Steve Sailer "The Unexpected Uselessness of
Philosophy"
-
Zelig Mailing List, served by Harvard-MIT Data Center
Send messages: zelig(a)latte.harvard.edu
[un]subscribe Options: http://lists.hmdc.harvard.edu/?info=zelig
Zelig program information: http://gking.harvard.edu/zelig/
30 Mar
30 Mar
4:04 p.m.
New subject: Cbind errors
I'm using R-2.2.1 (win) + MatchIt 2.2-5
Now, after matching and verifying that object m.out2 exists, invoking
R$ matches <- match.data(m.out2)
yields the following error message "Error in cbind(data, object$distance) :
cannot coerce type closure to list vector".
The match.data() function worked fine until last week, and I know I updated
some packages this week, but MatchIt certainly wasn't one of them.
thanks!
Ani
Anirudh V. S. Ruhil, Ph.D.
Sr. Research Associate
Voinovich Center for Leadership and Public Affairs
Ohio University
Building 21, The Ridges
Athens, OH 45701-2979
Tel: 740.597.1949 | Fax: 740.597.3057
-
Zelig Mailing List, served by Harvard-MIT Data Center
Send messages: zelig(a)latte.harvard.edu
[un]subscribe Options: http://lists.hmdc.harvard.edu/?info=zelig
Zelig program information: http://gking.harvard.edu/zelig/
8:29 p.m.
New subject: Cbind errors
Dear Ani,
Did the exact same command worked last week? If so, you might want to
reinstall the latest version of MatchIt and try again. If you still get an
error, could you also send the code that generated "m.out2" to the list?
Thanks
Kosuke
-----------------------------------------------------
Kosuke Imai Office: Corwin Hall 041
Assistant Professor Phone: 609-258-6601
Department of Politics eFax: 973-556-1929
Princeton University Email: kimai(a)Princeton.Edu
Princeton, NJ 08544-1012 http://imai.princeton.edu
-----------------------------------------------------
On Thu, 30 Mar 2006, Anirudh V. S. Ruhil wrote:
I'm using R-2.2.1 (win) + MatchIt 2.2-5
Now, after matching and verifying that object m.out2 exists, invoking
R$ matches <- match.data(m.out2)
yields the following error message "Error in cbind(data, object$distance) :
cannot coerce type closure to list vector".
The match.data() function worked fine until last week, and I know I updated
some packages this week, but MatchIt certainly wasn't one of them.
thanks!
Ani
Anirudh V. S. Ruhil, Ph.D.
Sr. Research Associate
Voinovich Center for Leadership and Public Affairs
Ohio University
Building 21, The Ridges
Athens, OH 45701-2979
Tel: 740.597.1949 | Fax: 740.597.3057
-
Zelig Mailing List, served by Harvard-MIT Data Center
Send messages: zelig(a)latte.harvard.edu
[un]subscribe Options: http://lists.hmdc.harvard.edu/?info=zelig
Zelig program information: http://gking.harvard.edu/zelig/
-
Zelig Mailing List, served by Harvard-MIT Data Center
Send messages: zelig(a)latte.harvard.edu
[un]subscribe Options: http://lists.hmdc.harvard.edu/?info=zelig
Zelig program information: http://gking.harvard.edu/zelig/
31 Mar
31 Mar
2:05 p.m.
New subject: Cbind errors
Dear Kosuke;
I reinstalled matchIt 2.2-5 and the error persists. The relevant code and
associated output (truncated to save bandwidth but show m.out2 exists)
follow ...
------------------------------------
library(MatchIt)
m.out2 <- matchit(scale ~ population + urbanization
+ + pctownocc + mdyearstruc
+ + mdhomevalue + pcincome
+ + pctjuvenile + pctjuvespoor,
+ data = scale, method = "nearest")
summary(m.out2)
Call:
matchit(formula = scale ~ population + urbanization + pctownocc +
mdyearstruc + mdhomevalue + pcincome + pctjuvenile + pctjuvespoor,
data = scale, method = "nearest")
{....}
{....}
Percent Balance Improvement:
Mean Diff. eQQ Med eQQ Mean eQQ Max
distance 98.45 96.71 94.96 89.44
population 48.63 24.30 62.40 86.64
urbanization 95.06 57.10 94.97 96.52
pctownocc 84.82 81.91 82.74 65.41
mdyearstruc 86.78 83.33 76.47 73.08
mdhomevalue 92.66 87.10 87.44 93.22
pcincome 96.39 87.71 91.88 92.71
pctjuvenile 85.76 80.44 84.25 86.66
pctjuvespoor 86.45 80.75 81.74 67.84
Sample sizes:
Control Treated
All 160 13
Matched 13 13
Unmatched 147 0
Discarded 0 0
m.out2$match.matrix
1
1 "5"
31 "59"
35 "44"
53 "34"
70 "64"
81 "62"
87 "146"
100 "102"
137 "127"
148 "154"
158 "30"
171 "84"
173 "114"
matches <- match.data(m.out2)
Error in
cbind(data, object$distance) : cannot coerce type closure to list
vector
----------------------------------
best
Ani
Anirudh V. S. Ruhil, Ph.D.
Sr. Research Associate
Voinovich Center for Leadership and Public Affairs
Ohio University
Building 21, The Ridges
Athens, OH 45701-2979
Tel: 740.597.1949 | Fax: 740.597.3057
-
Zelig Mailing List, served by Harvard-MIT Data Center
Send messages: zelig(a)latte.harvard.edu
[un]subscribe Options: http://lists.hmdc.harvard.edu/?info=zelig
Zelig program information: http://gking.harvard.edu/zelig/
4 Apr
4 Apr
9:18 p.m.
New subject: Cbind errors
Dear Ani,
I just released MatchIt version 2.2-6 which includes this bug fix. You
can install this latest version via
install.packages("MatchIt", repos="http://gking.harvard.edu")
Or you can wait for a few days until it appears on CRAN.
Cheers,
Kosuke
-----------------------------------------------------
Kosuke Imai Office: Corwin Hall 041
Assistant Professor Phone: 609-258-6601
Department of Politics eFax: 973-556-1929
Princeton University Email: kimai(a)Princeton.Edu
Princeton, NJ 08544-1012 http://imai.princeton.edu
-----------------------------------------------------
On Fri, 31 Mar 2006, Anirudh V. S. Ruhil wrote:
Dear Kosuke;
I reinstalled matchIt 2.2-5 and the error persists. The relevant code and
associated output (truncated to save bandwidth but show m.out2 exists)
follow ...
------------------------------------
-
Zelig Mailing List, served by Harvard-MIT Data Center
Send messages: zelig(a)latte.harvard.edu
[un]subscribe Options: http://lists.hmdc.harvard.edu/?info=zelig
Zelig program information: http://gking.harvard.edu/zelig/
library(MatchIt)
m.out2 <- matchit(scale ~ population +
urbanization
+ + pctownocc + mdyearstruc
+ + mdhomevalue + pcincome
+ + pctjuvenile + pctjuvespoor,
+ data = scale, method = "nearest")
summary(m.out2)
Call:
matchit(formula = scale ~ population + urbanization + pctownocc +
mdyearstruc + mdhomevalue + pcincome + pctjuvenile + pctjuvespoor,
data = scale, method = "nearest")
{....}
{....}
Percent Balance Improvement:
Mean Diff. eQQ Med eQQ Mean eQQ Max
distance 98.45 96.71 94.96 89.44
population 48.63 24.30 62.40 86.64
urbanization 95.06 57.10 94.97 96.52
pctownocc 84.82 81.91 82.74 65.41
mdyearstruc 86.78 83.33 76.47 73.08
mdhomevalue 92.66 87.10 87.44 93.22
pcincome 96.39 87.71 91.88 92.71
pctjuvenile 85.76 80.44 84.25 86.66
pctjuvespoor 86.45 80.75 81.74 67.84
Sample sizes:
Control Treated
All 160 13
Matched 13 13
Unmatched 147 0
Discarded 0 0
m.out2$match.matrix
1
1 "5"
31 "59"
35 "44"
53 "34"
70 "64"
81 "62"
87 "146"
100 "102"
137 "127"
148 "154"
158 "30"
171 "84"
173 "114"
matches <- match.data(m.out2)
Error in
cbind(data, object$distance) : cannot coerce type closure to list
vector
----------------------------------
best
Ani
Anirudh V. S. Ruhil, Ph.D.
Sr. Research Associate
Voinovich Center for Leadership and Public Affairs
Ohio University
Building 21, The Ridges
Athens, OH 45701-2979
Tel: 740.597.1949 | Fax: 740.597.3057
6619
days inactive
6703
days old
4 comments
3 participants
participants (3)
-
Anirudh V. S. Ruhil -
Gary King -
Kosuke Imai