24 Apr
2006
24 Apr
'06
7:20 p.m.
Hi all,
I'm seeking advice about code I've put together to test a counterfactual.
Specifically, I'm trying to ask: What if all observations for our treatment
variable, the infant mortality rate (infantlag), were set to the 25th
quantile? How factual is this counterfactual? I'm using the code from page 4
of the Section 9 handout as my guide. My code reads:
library(WhatIf)
fact <- na.omit(Civil.Data)
counterfact <- na.omit(Civil.Data)
## I also tried running these prior two lines without using na.omit
counterfact$infantlag <- quantile(Civil.Data$infantlag, .25, na.rm=T)
whif <- whatif(data=fact,cfact=counterfact)
The last line of code returns the error message:
"Error in whatif(data = fact, cfact = counterfact) :
observed covariate data not coercable to numeric matrix due to
character column(s)"
Interestingly, I get the same error message when I try to replicate Dan and
Ian's code from that handout using the coalition data. I'd appreciate any
thoughts on what's causing this error message.
Thanks,
Janet
24 Apr
24 Apr
7:24 p.m.
New subject: [gov2001-l] ternarypoints()
I'm having some trouble adding ternarypoints to an existing ternary
plot. The plot is created, very pretty actually, but the ternarypoints
doesn't see it for some reason. Suggestions? --Doru
ternaryplot(s.out$qi$ev, pch = ".", col =
"blue",
dimnames=c("Dem.","Ind.","Rep."),
main = "Individual effects of Prop.
209") ### <- This works and it creates the plot
ternarypoints(ev.hi, pch = ".", col =
"red") ### <-ev.hi was
correctly defined earlier
Error in plot.xy(xy.coords(x, y), type = type, pch = pch, col = col, bg
= bg, :
plot.new has not been called yet
9:51 p.m.
New subject: [gov2001-l] ternarypoints()
Hi Doru,
You aren't shutting the plot window in between the two commands, right?
Usually you get that error message when you try to plot points wihout
using the plot command first. But if you are just entering the two lines
of code right after each other, it seems like it should work.
Best,
Ian
On Mon, 24 Apr 2006, Doru Gavril wrote:
I'm having some trouble adding ternarypoints to
an existing ternary
plot. The plot is created, very pretty actually, but the ternarypoints
doesn't see it for some reason. Suggestions? --Doru
ternaryplot(s.out$qi$ev, pch = ".", col
= "blue",
dimnames=c("Dem.","Ind.","Rep."),
main = "Individual effects of Prop.
209") ### <- This works and it creates the plot
ternarypoints(ev.hi, pch = ".", col =
"red") ### <-ev.hi was
correctly defined earlier
Error in plot.xy(xy.coords(x, y), type = type, pch = pch, col = col, bg
= bg, :
plot.new has not been called yet
_______________________________________________
gov2001-l mailing list
gov2001-l(a)lists.fas.harvard.edu
http://lists.fas.harvard.edu/mailman/listinfo/gov2001-l
25 Apr
25 Apr
6:46 a.m.
New subject: [gov2001-l] ternarypoints()
So here's the interesting thing. Ternaryplot exists in both vcd and
Zelig. If Zelig is open and you load vcd ternarypoints (only defined in
Zelig) does not recognize the new vcd function. If you leave out vcd it
works.
On 4/25/2006 1:51 AM, Ian Brett Yohai wrote:
Hi Doru,
You aren't shutting the plot window in between the two commands, right?
Usually you get that error message when you try to plot points wihout
using the plot command first. But if you are just entering the two lines
of code right after each other, it seems like it should work.
Best,
Ian
On Mon, 24 Apr 2006, Doru Gavril wrote:
I'm having some trouble adding ternarypoints
to an existing ternary
plot. The plot is created, very pretty actually, but the ternarypoints
doesn't see it for some reason. Suggestions? --Doru
_______________________________________________
gov2001-l mailing list
gov2001-l(a)lists.fas.harvard.edu
http://lists.fas.harvard.edu/mailman/listinfo/gov2001-l
ternaryplot(s.out$qi$ev, pch = ".", col
= "blue",
dimnames=c("Dem.","Ind.","Rep."),
main = "Individual effects of Prop.
209") ### <- This works and it creates the plot
ternarypoints(ev.hi, pch = ".", col =
"red") ### <-ev.hi was
correctly defined earlier
Error in plot.xy(xy.coords(x, y), type = type, pch = pch, col = col, bg
= bg, :
plot.new has not been called yet
_______________________________________________
gov2001-l mailing list
gov2001-l(a)lists.fas.harvard.edu
http://lists.fas.harvard.edu/mailman/listinfo/gov2001-l
24 Apr
24 Apr
7:48 p.m.
Janet, not familiar with your particular problem, but when I have seen error
messages like the one you have, it usually is because I'm trying to apply a
matrix operation to a data frame. My guess is that you have some kind of
data frame containing strings which has to be "coerced" into a matrix, which
may only take on numeric (?possibly NA??) values. Try creating a subsetted
data frame with only your numeric variables of interest.
On 4/24/06, Janet Ingram Lewis <jilewis(a)fas.harvard.edu> wrote:
Hi all,
I'm seeking advice about code I've put together to test a counterfactual.
Specifically, I'm trying to ask: What if all observations for our treatment
variable, the infant mortality rate (infantlag), were set to the 25th
quantile? How factual is this counterfactual? I'm using the code from page 4
of the Section 9 handout as my guide. My code reads:
library(WhatIf)
fact <- na.omit(Civil.Data)
counterfact <- na.omit(Civil.Data)
## I also tried running these prior two lines without using na.omit
counterfact$infantlag <- quantile(Civil.Data$infantlag , .25, na.rm=T)
whif <- whatif(data=fact,cfact=counterfact)
The last line of code returns the error message:
"Error in whatif(data = fact, cfact = counterfact) :
observed covariate data not coercable to numeric matrix due to
character column(s)"
Interestingly, I get the same error message when I try to replicate Dan
and Ian's code from that handout using the coalition data. I'd appreciate
any thoughts on what's causing this error message.
Thanks,
Janet
_______________________________________________
gov2001-l mailing list
gov2001-l(a)lists.fas.harvard.edu
http://lists.fas.harvard.edu/mailman/listinfo/gov2001-l
6596
days inactive
6597
days old
gov2001@lists.gking.harvard.edu
4 comments
4 participants
participants (4)
-
dgavril＠fas.harvard.edu -
jilewis＠fas.harvard.edu -
nall＠fas.harvard.edu -
yohai＠fas.harvard.edu