|
|
|
||||
|
|
|||||
The hsb2.sdd data frame to be used in SPLUS
The hsb2.csv data frame to be used in R.
Suppose that we would like to overlay two different plots in one single trellis graph. For example, we would like to overlay two plots consisting of a scatter plot with a regression line but in each graph we would like to have a different dependent and independent variable. In one trellis graph we would like to have a scatter plot and regression line of read regressed on write and a scatter plot and regression line of science regression on math both conditional on the variable ses.
Reading in the hsb2 data frame.
#In SPLUS only
hsb2 <- openData('c:/hsb2.sdd')
#In R only
hsb2 <- read.table('c:/hsb2.csv', header=T, sep=",")
#In R we need to download the lattice package in order to generate trellis graphs
library(lattice)
Let's first look at how we would generate each trellis graph individually.
#creating a factor variable ses.f from the variable ses
hsb2$ses.f <- factor(hsb2$ses, label=c("low", "medium", "high"))
xyplot(write~read | ses.f, hsb2,
panel=function(x, y){
panel.xyplot(x, y, pch=16)
panel.lmline(x, y, lty=4)
}, as.table=T)
xyplot(math~science | ses.f, hsb2,
panel=function(x, y){
panel.xyplot(x, y, pch=3)
panel.lmline(x, y)
}, as.table=T)
In order to combine the two graphs we will make use of the subscript parameter which is a parameter in every panel function. The subscript parameter is an indicator of which indices corresponds to which panel in the trellis graph. Thus, we include the subscript parameter in the function we create for the panel argument in the trellis graph. Since we want to have two scatter plots overlaid we must have two panel.xyplot functions, likewise we want to have two regression lines in each plot and therefore we must include two panel.lmline functions. One of the panel.xyplot function will use the variables specified in the formula argument (write~read | ses.f) of the xyplot function which indicates which observation will be used in each of the panels. In the other panel.xyplot function we need to use the subscript parameter to indicate which observations will be used in each of the panels because we are no longer using the variables in the formula given in the formula argument in the xyplot function.
xyplot(write~read | ses.f, hsb2,
panel=function(x, y, subscripts){
panel.xyplot(x, y, pch=16)
panel.lmline(x, y, lty=4)
panel.xyplot(hsb2$science[subscripts], hsb2$math[subscripts], pch=3)
panel.lmline(hsb2$science[subscripts], hsb2$math[subscripts])
}, as.table=T, subscripts=T)
UCLA Researchers are invited to our Statistical Consulting Services
We recommend others to our list of Other Resources for Statistical Computing Help
These pages are Copyrighted (c) by UCLA Academic Technology Services
