[R] Cbind() on the right-side of a formula in xYplot()

Dylan Beaudette debeaudette at ucdavis.edu
Mon Oct 26 18:55:38 CET 2009


Using the latest rms package I am able to make nice plots of model predictions 
+/- desired confidence intervals like this:

# need this

# setup data
d <- data.frame(x=rnorm(100), y=rnorm(100))
dd <- datadist(d)

# fit model
l <- ols(y ~ rcs(x), data=d)

# predict along original limits of data
l.pred <- Predict(l)

# plot of fit and conf. int.
xYplot(Cbind(yhat, lower, upper) ~ x, data=l.pred, method='filled', 
col.fill=grey(0.9), col=1, type='l')

Is there any way in which I can turn this figure on its side, by plotting x ~ 
y... something like this:

# doesn't work
xYplot(x ~ Cbind(yhat, lower, upper), data=l.pred, method='filled', 
lable.curves=FALSE, col.fill=grey(0.9), col=1, type='l')

# standard lattice, but without the fancy filled area
xyplot(x ~ yhat + lower + upper, data=l.pred, type='l', lty=c(1,2,2), col=1)

Any suggestions? If it is not possible, then I will try and manually make the 
figure with basic lattice functions.


Dylan Beaudette
Soil Resource Laboratory
University of California at Davis

More information about the R-help mailing list