Skip to content
Permalink
Branch: master
Find file Copy path
Find file Copy path
Fetching contributors…
Cannot retrieve contributors at this time
59 lines (46 sloc) 1.59 KB
#
# common-dis-distrib.R, 6 Aug 16
# Data from:
#
# Example from:
# Evidence-based Software Engineering: based on the publicly available data
# Derek M. Jones
#
# TAG example
source("ESEUR_config.r")
plot_layout(4, 1)
par(mar=c(1.2, 2.0, 1, 0.8))
brew_col=rainbow(3)
plot(c(0, 1, 2, 3, 4, 5), c(0, 0.25, 0.25, 0.25, 0.25, 0), type="h",
bty="n", yaxt="n",
col=brew_col[1],
xlab="", ylab="")
x_geom=0:10
plot(x_geom, dgeom(x_geom, prob=0.2),
bty="n", yaxt="n",
col=brew_col[1],
xlim=range(x_geom), ylim=c(0, 0.5),
xlab="", ylab="")
points(x_geom, dgeom(x_geom, prob=0.4), col=brew_col[2])
points(x_geom, dgeom(x_geom, prob=0.6), col=brew_col[3])
legend(x="topright", legend=c("prob=0.2", "prob=0.4", "prob=0.6"), bty="n", fill=brew_col, cex=1.1)
x_binom=0:10
plot(x_binom, dbinom(x_binom, size=10, prob=0.2),
bty="n", yaxt="n",
col=brew_col[1],
xlim=range(x_binom), ylim=c(0, 0.3),
xlab="", ylab="")
points(x_binom, dbinom(x_binom, size=10, prob=0.4), col=brew_col[2])
points(x_binom, dbinom(x_binom, size=25, prob=0.2), col=brew_col[3])
legend(x="topright", legend=c("size=10\nprob=0.2\n", "size=10\nprob=0.4\n", "size=25\nprob=0.2"), bty="n", fill=brew_col, cex=1.1)
lambda_str=function(num) as.expression(substitute(lambda == num))
x_pois=0:10
plot(x_pois, dpois(x_pois, lambda=2),
bty="n", yaxt="n",
col=brew_col[1],
xlim=range(x_pois), ylim=c(0, 0.3),
xlab="", ylab="")
points(x_pois, dpois(x_pois, lambda=4), col=brew_col[2])
points(x_pois, dpois(x_pois, lambda=7), col=brew_col[3])
legend(x="topright", legend=c(lambda_str(2), lambda_str(4), lambda_str(7)),
bty="n", fill=brew_col, cex=1.1)
You can’t perform that action at this time.