1360 tissue specific MHL in tissue cluster

From ZhangLabWiki
Revision as of 10:50, 22 February 2016 by >Shicheng (Created page with " # 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<...")
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
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)