Thank you, Dieter:

I apologize to all of those who tried to help and who couldn't get my code to work.

I now know three things: 1) stop working at midnight; 2) always provide the correct package

name (I meant arm(), not car(); and 3) when testing code make sure you do it in a new R window.

Dieter has kindly provided a clean version of my code below with a resolution to my plotting

requirements.

Thank you,

Josip

----- Original Message -----

From: "Dieter Menne" <

[hidden email]>

To:

[hidden email]
Sent: Wednesday, October 1, 2008 3:09:25 AM GMT -08:00 US/Canada Pacific

Subject: Re: [R] Change color of plot points based on values of a variable

Josip Dasovic <jjd9 <at> sfu.ca> writes:

> I have run a logistic regression, used Gelman et al.'s car package to

>> simulate the parameter estimates of

"car" is by John Fox. You meant "arm".

> I would like the points from 1-24, and 976-1000 to

be one color (the default color is

> fine), and the points 25-975 to be a different color.

Please check your code in a fresh window before you submit it. There were

7 errors in your code that I almost gave up, and I am sure other people had

the same problem.

Dieter

library(arm)

set.seed(23)

y<-rbinom(100,1,0.1)

x<-rnorm(100)

fit.1<-glm(y~x, family=binomial(link="logit"))

sim.fit.1<-sim(fit.1, 1000)

plotdt =sort(invlogit(sim.fit.1$beta[,1]+sim.fit.1$beta[,2]*mean(x)))

col = rep("black",length(plotdt))

col[1:24] = "red"

plot(plotdt,col=col)

______________________________________________

[hidden email] mailing list

https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide

http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code.

______________________________________________

[hidden email] mailing list

https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide

http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code.