# load example data:
data(coal)
data(coalobs)
param<-coal[1:50000,2]
simstats<-coal[1:50000,4:6]
# use matrix below just in case to preserve dimensions.
obsstats<-matrix(coalobs[1,4:6],nrow=1)
obsparam<-matrix(coalobs[1,1])
# example of entropy minimization algorithm:
tmp <-mincrit(obsstats, param, simstats, tol=.01, method="rejection",
do.crit=TRUE)
tmp$critvals
Run the code above in your browser using DataLab