-
Notifications
You must be signed in to change notification settings - Fork 0
/
get_atac_bed.R
19 lines (16 loc) · 1.03 KB
/
get_atac_bed.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
data <- read.csv("chromatin_access/Human_atac_dif.csv",header=T,row.names=1)
bed <- matrix(unlist(strsplit(rownames(data),"[:-]")),nrow(data),3,byrow=T)
bed <- cbind(bed,rownames(data))
write.table(bed,"Human_atac.bed",sep="\t",quote=F,row.names=F,col.names=F)
data <- read.csv("chromatin_access/Macaque_atac_dif.csv",header=T,row.names=1)
bed <- matrix(unlist(strsplit(rownames(data),"[:-]")),nrow(data),3,byrow=T)
bed <- cbind(bed,rownames(data))
write.table(bed,"Macaque_atac.bed",sep="\t",quote=F,row.names=F,col.names=F)
data <- read.csv("chromatin_access/Marmoset_atac_dif.csv",header=T,row.names=1)
bed <- matrix(unlist(strsplit(rownames(data),"[:-]")),nrow(data),3,byrow=T)
bed <- cbind(bed,rownames(data))
write.table(bed,"Marmoset_atac.bed",sep="\t",quote=F,row.names=F,col.names=F)
data <- read.csv("chromatin_access/Mouse_atac_dif.csv",header=T,row.names=1)
bed <- matrix(unlist(strsplit(rownames(data),"[:-]")),nrow(data),3,byrow=T)
bed <- cbind(bed,rownames(data))
write.table(bed,"Mouse_atac.bed",sep="\t",quote=F,row.names=F,col.names=F)