# NOT RUN {
x <- 1:7
GA <- c(1,3,1.5,8); DE <- c(1,3,4,2.5)
par(mfrow=c(2,2))
for (j in 1:4) {
y <- NULL
for (i in x) y <- c(y,pfakebad(x[i],h=5,Q=7,gam=GA[j],del=DE[j],p=.4))
plot(x,y,type="h",panel.first=points(x,y,pch=19),
main=paste("gamma=",GA[j]," delta=",DE[j],sep=""),ylim=c(0,.7),
ylab="Replacement probability")
}
# }
Run the code above in your browser using DataLab