Nothing
######################################################################
## File: test4.R
## Author: Ray Griner
## Date: November 11, 2011 (111111!)
## Desc: Short examples of predicted interval plots (PIPs)
######################################################################
set.seed(12345)
library(PIPS)
#################################################################################
# EXAMPLE 4: Overriding default graphic options. These are only examples of the
# customizations that can be made. Users should be able to customize almost any
# feature of the graph. (Type help(par) in an R session for a list of R graphics
# parameters. Any of these parameters can be passed to the plot function.)
#################################################################################
myY<-c(rep(1,times=20),rep(0,times=80),rep(1,times=25),rep(0,times=25))
myGroup<-c(rep('A',100),rep('B',50))
pips <- pred.int(y=myY, group=myGroup, N=c(400,400), data.type="binary", iters=100)
print(pips, pi.count=50)
plot(pips, # The object containing the pips data needs to be first
xlab="My x axis label", # Change default x label
main="My graph for #BY#", # Remember, #BY# is replaced with "B vs A". Not really useful
# here with one group but would be useful if multiple groups
xlim=c(.12,.6), # Change limits of x-axis
ylab="My y axis label",
cex.sub=.8, # Change size of footnote text
sub="This example shows how to customize the appearance of graphics elements", # Footnote
col.main="blue", # Change color of main title
axes=FALSE) # Do not create default axes. Will use custom axes with axis statements
###########################################
# Custom axes
###########################################
axis(1, at=seq(from=0,to=1,by=.04),las=1)
axis(2, at=seq(from=0,to=100,by=10), las=1)
###########################################
# Create a legend
###########################################
legend("topright",c("Actual","Predicted"), col=c("red","black"), lwd=c(2,1))
###########################################
# Create a couple vertical lines
###########################################
abline(v=c(.3,.4), lty=2, lwd=2)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.