David Winsemius
On Mar 5, 2009, at 9:18 AM, Michael Friendly wrote:
> In the following, I'm fitting a logistic regression model, and using
> car:::influencePlot. When I run the latter with
> output to the screen, it calls identify() that lets me label
> observations with large CookD.
>
> However, if I use postscript() to get .eps output, identify() seems
> not to be called at all. If instead, I
> use dev.copy2eps() after getting output to the screen, the point
> labels do not appear in the resulting .eps
> graph. Why? Is there a workaround?
>
> library(vcd)
> data(Arthritis)
> # define Better
> Arthritis$Better <- Arthritis$Improved > 'None'
> arth.mod1 <- glm(Better ~ Age + Sex + Treatment , data=Arthritis,
> family='binomial')
> library(car)
> # influence plot
> #postscript(file="arthritis-diag2.eps", paper="special", height=6,
> width=6, onefile=TRUE, horizontal=FALSE)
> influencePlot(arth.mod1, main="Arthritis data: influencePlot")
> #dev.off()
>
> -Michael
>
> --
> Michael Friendly Email: friendly AT yorku DOT ca Professor,
> Psychology Dept.
> York University Voice: 416 736-5115 x66249 Fax: 416 736-5814
> 4700 Keele Street http://www.math.yorku.ca/SCS/friendly.html
> Toronto, ONT M3J 1P3 CANADA
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.