plotCumProp <- function(x, line.col="blue", xlims=c(-1,1), add=FALSE, xlab="score", res=1000, invert=FALSE, ...) {
if((min(x) < xlims[1]) | (max(x) > xlims[2])) warning("scores outside plotting range!")
breaks <- seq(xlims[1], xlims[2], length.out=res)
mean.cut <- cut(x, breaks, right=F)
mean.freq <- table(mean.cut)
if(invert) {cumfreq0 = c(0, cumsum(mean.freq)) / length(x)
} else {
cumfreq0 = 1 - c(0, cumsum(mean.freq)) / length(x)
}
if(add) {
lines(breaks, cumfreq0 , xlab=xlab, ylab="cumulative proportion", type="l", lwd=2, col=line.col)
} else {
#plot.new()
plot(breaks, cumfreq0 , xlab=xlab, ylab="cumulative proportion", type="l", lwd=2, col=line.col, xlim=xlims, ylim=c(0,1), ...)
}
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.