1360 tissue specific MHL in tissue cluster
Jump to navigation
Jump to search
# tissue specific hml, 1360 regions gsi_1<-gsi(f2[,13:58]) gsi_2<-gsi(f2) data<-subset(gsi_1,gsi_1[,3]>0.5) table(data[,2]) z<-c() for (i in names(table(data[2]))){ x<-subset(data,data[,2]==i) y<-x[order(x[,3],decreasing=T)[1:80],] z<-rbind(z,y) } z<-na.omit(z) dim(z) write.table(z,file="supplementary.1360.tsi.txt",col.names=T,row.names=F,sep="\t",quote=F)