On Jun 26, 2012, at 11:29 AM, Sarah Goslee wrote:
You could use points() instead of plot() for the second command.
Ummm. Maybe not. I think think that plot.Predict uses lattice
graphics. You may need to use trellis.focus() followed by lpoints().
Or use the "+" operation with suitable objects.
--
David.
Sarah
On Tue, Jun 26, 2012 at 8:37 AM, achaumont <agnes.chaum...@live.be>
wrote:
Hello,
I have a question about the “plot.predict” function in Frank
Harrell's rms
package.
Do you know how to superpose in the same graph the prediction curve
of ols
and raw data points?
Put most simply, I would like to combine these two graphs:
fit_linear <- ols (y4 ~ rcs(x2,c(5,10,15,20,60,80,90)), x=TRUE,
y=TRUE)
p <- Predict(fit_linear,x2,conf.int=FALSE)
plot (p, ylim =c(-2,0.5), xlim=c(0,100)) # graph n°1
z <- plot
(x2,y4,ylim=c(-2,0.5),xlim=c(0,100),type="p",lwd=6,col="blue")
# graph n°2
Thanks all,
Agnès
--
Sarah Goslee
http://www.functionaldiversity.org
______________________________________________
R-help@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.
David Winsemius, MD
West Hartford, CT
______________________________________________
R-help@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.