require("MEET")
data(BackgroundOrganism)
data(iicc)
q<-iicc$q
Factortrans<-iicc$Transcriptionfactor
HXmax<-iicc$HXmax
correction<-correction.entropy(q,p=nrow(Factortrans),long=1,iicc)
Herror<-slot(correction,"Herror")
pmX<-probability(Factortrans,Prob)
HX<-switch(iicc$classentropy,"Shannon"=entropy.Shannon(pmX),"Renyi"=entropy.Renyi(pmX,q))
test<-redundancy(HX,HXmax,Herror,finite=nrow(Factortrans),correction=TRUE)
Run the code above in your browser using DataLab