1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
|
makeIntersectList<-function(x,xnames=NULL) {
# if no names passed, see if x is a data frame with names
if(is.null(xnames)) xnames<-names(x)
# coerce x to a matrix
if(!is.matrix(x)) x<-as.matrix(x)
# if it was already a matrix and xnames is still NULL, try colnames
if(is.null(xnames)) xnames<-colnames(x)
dimx<-dim(x)
# if no names available, use capital letters
if(is.null(xnames)) xnames<-LETTERS[1:dimx[2]]
intersectList<-vector("list",dimx[2]+1)
for(intersect in 1:dimx[2])
intersectList[[1]][intersect]<-sum(rowSums(x) == 1 & x[,intersect])
names(intersectList[[1]])<-xnames
for(comb in 2:dimx[2]) {
nn<-choose(dimx[2],comb)
intersectList[[comb]]<-rep(0,nn)
currentnames<-names(intersectList[[comb]])<-pasteCols(combn(xnames,comb))
currentcombs<-combn(1:dimx[2],comb,simplify=TRUE)
for(intersect in 1:nn) {
combvec<-rep(0,dimx[2])
combvec[currentcombs[,intersect]]<-1
intersectList[[comb]][intersect]<-sum(colSums(apply(x,1,"==",combvec))==dimx[2])
}
}
intersectList[[dimx[2]+1]]<-dimx[1]
names(intersectList[[dimx[2]+1]])<-"Total"
class(intersectList)<-"intersectList"
return(intersectList)
}
|
