Calculating entropy as suggested by Wand and King (2007) is straightforward. The entropy() function simply calls anchors(), mrop(), and predict() for every permutation of vignettes and calculates the minimum known entropy and estimated entropy for each permutation.
data(poleff)
ent <- entropy(xsayself ~ xsay1 + xsay2 + xsay3 + xsay4 + xsay5,
covar = cbind(Cs, Ce) ~ as.factor(china) + age + male + educyrs,
data = poleff)
summary(ent)
plot(ent)
For more details, please see help(entropy) in R. The plot() method for class entropy plots minimum known entropy
against estimated entropy, so you must specify a formula in covar for this to work.