# Load needed objects
chcemat <- Quant_13B_problem_dataList$chcemat
index <- Quant_13B_problem_parmList$index
noption <- matrix(5,24,1)
# compute mutual entropies for all pairs of the first 6 items
Mvec <- 1:6
Mlen <- length(Mvec)
Hmutual <- matrix(0,Mlen,Mlen)
for (i1 in 1:Mlen) {
for (i2 in 1:i1) {
Result <- entropies(index, Mvec[i1], Mvec[i2], chcemat, noption)
Hmutual[i1,i2] = Result$Hmutual
Hmutual[i2,i1] = Result$Hmutual
}
}
print("Matrix of mutual entries (off-digagonal) and self-entropies (diagonal)")
print(round(Hmutual,3))
Run the code above in your browser using DataLab