~ubuntu-branches/debian/sid/rkward/sid

local({

## Prepare

## Compute

## Print result

# parameters:

m <- 12; n <- 12; k <- 15;

# mean and variances of the distribution of sample averages:

avg.exp <- k*m/(m+n);

avg.var <- (k*m*n*(m+n-k)/((m+n)^2*(m+n-1)))/10;

# generate the entire data:

data <- matrix(rhyper(nn=10000, m=m, n=n, k=k), nrow=10);

# get the sample averages:

avg <- colMeans(data);

# generate random normal samples:

normX <- seq(from=min(avg), to=max(avg), length=1000);

normY <- dnorm (normX, mean = avg.exp, sd = sqrt(avg.var));

dist.hist <- hist(avg, plot=FALSE, breaks="Sturges");

# calculate the ylims appropriately:

ylim <- c(0,max(c(dist.hist$density, normY)));

rk.graph.on ()

try ({

        plot(dist.hist, ylim=ylim, freq=FALSE, lty="solid", density=-1, xlab="Sample Averages", main="Hypergeometric")

        lines (x=normX, y=normY, type="l", col="red")

})

rk.graph.off ()

})

.rk.rerun.plugin.link(plugin="rkward::plot_hypergeometric_clt", settings="drawnorm.state=1\nfunction.string=hist\nhistogram_opt.addtoplot.state=\nhistogram_opt.barlabels.state=\nhistogram_opt.density.real=-1.00\nhistogram_opt.doborder.state=1\nhistogram_opt.freq.state=0\nhistogram_opt.histbordercol.color.string=\nhistogram_opt.histbreaksFunction.string=Sturges\nhistogram_opt.histlinetype.string=solid\nhistogram_opt.rightclosed.state=1\nhistogram_opt.usefillcol.state=\nk.real=15.00\nm.real=12.00\nn.real=12.00\nnAvg.real=10.00\nnDist.real=1000.00\nnormlinecol.color.string=red\nnormpointtype.string=l\nplotoptions.add_grid.state=0\nplotoptions.asp.real=0.00\nplotoptions.main.text=\nplotoptions.pointcolor.color.string=\nplotoptions.pointtype.string=\nplotoptions.sub.text=\nplotoptions.xaxt.state=\nplotoptions.xlab.text=\nplotoptions.xlog.state=\nplotoptions.xmaxvalue.text=\nplotoptions.xminvalue.text=\nplotoptions.yaxt.state=\nplotoptions.ylab.text=\nplotoptions.ylog.state=\nplotoptions.ymaxvalue.text=\nplotoptions.yminvalue.text=\nscalenorm.state=0", label="Run again")

.rk.make.hr()