# multi-figure output
op <- par(no.readonly = TRUE)
par(mfrow = c(2, 2))
# defaults
plotDAF('C')
plotDAF('G')
plotDAF('B')
plotDAF('D')
# adjust intercept term for DAF 'C'
plotDAF('C', intercept = 0.1)
title(sub = 'intercept = 0.1', font.sub = 3)
plotDAF('C', intercept = 0.3)
title(sub = 'intercept = 0.3', font.sub = 3)
# adjust threshold term for DAF 'G'
plotDAF('G', threshold = 0.5)
title(sub = 'threshold = 0.5', font.sub = 3)
plotDAF('G', threshold = 0.9)
title(sub = 'threshold = 0.9', font.sub = 3)
# reset output device options
par(op)
Run the code above in your browser using DataLab