Commit 0200b0f89b2d3d0e78e84b671aff59b1a2dad4c4
1 parent
671da0d26a
Exists in
master
This code is to calculate the entropy of any 3 variable.
Showing
1 changed file
with
25 additions
and
0 deletions
Show diff stats
entropySorted.R
| ... | ... | @@ -0,0 +1,25 @@ |
| 1 | +#Please read README file. | |
| 2 | + | |
| 3 | +EntropySorted = | |
| 4 | + function (filename = "/home/zgong001/Documents/SprinklerDataset/sprinkler1001.txt", alphalevel = 0.05, methodu="x2"){ | |
| 5 | + library(entropy) | |
| 6 | + library(bnlearn) | |
| 7 | + result = data.frame() | |
| 8 | + sprinkler <- read.table(filename, header = TRUE) | |
| 9 | + columname = colnames(sprinkler) | |
| 10 | + m=length(columname) | |
| 11 | + x=seq(0,m-1,by=1) | |
| 12 | + colnames(sprinkler)=x | |
| 13 | + columname = colnames(sprinkler) | |
| 14 | + ThreeVComb = combn(columname,3) | |
| 15 | + for(jt in 1:ncol(ThreeVComb)){ | |
| 16 | + ThV = c(ThreeVComb[1,jt], ThreeVComb[2,jt], ThreeVComb[3,jt]) | |
| 17 | + tbT = table(sprinkler[,ThV[1]], sprinkler[,ThV[2]],sprinkler[,ThV[3]]) | |
| 18 | + en = entropy(tbT) | |
| 19 | + lable = paste(ThV[1], ThV[2], ThV[3], sep=", ") | |
| 20 | + CIlist= data.frame("Index" = lable, "entropy" = en) | |
| 21 | + result = rbind(result, CIlist) | |
| 22 | + } | |
| 23 | + result2 = result[order(-result[,2]), ] | |
| 24 | + return(result2) | |
| 25 | + } | |
| 0 | 26 | \ No newline at end of file |