1360 tissue specific MHL in tissue cluster

From ZhangLabWiki
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)