Nothing
# ApproxHypergeo.R
# This demo compares a Wallenius' and a Fisher's noncentral hypergeometric
# distribution with the same mean rather than the same odds in order to
# make them approximate each other better.
require(BiasedUrn)
require(stats)
ApproxHypPlot <- function(m1, m2, n, w.odds) {
xmin <- minHypergeo(m1, m2, n)
xmax <- maxHypergeo(m1, m2, n)
x <- xmin : xmax
w.mean <- meanWNCHypergeo(m1, m2, n, w.odds)
f.odds <- oddsFNCHypergeo(w.mean, m1, m2, n)
wnc <- dWNCHypergeo(x, m1, m2, n, w.odds)
fnc <- dFNCHypergeo(x, m1, m2, n, f.odds)
fnc0 <- dFNCHypergeo(x, m1, m2, n, w.odds)
plot (x, fnc, type="l", col="red",
main = "Hypergeometric distributions",
sub = "Blue = Wallenius, Red = Fisher w. same mean,\n Green = Fisher w. same odds",
xlab = "", ylab = "Probability")
points (x, wnc, type="l", col="blue")
points (x, fnc0, type="l", col="green", lty="dashed")
}
ApproxHypPlot(80, 60, 100, 0.5)
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.