-
Notifications
You must be signed in to change notification settings - Fork 1
/
1_modulizer_output_to_rds.R
21 lines (17 loc) · 1.46 KB
/
1_modulizer_output_to_rds.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
library(dplyr)
library(tidyr)
binaryEvents <- list(cassette=read.table("data/modulize/cassette.tsv", header=T, sep="\t"),
alternative_intron=read.table("data/modulize/alternative_intron.tsv", header=T, sep="\t"),
alt5prime=read.table("data/modulize/alt5prime.tsv", header=T, sep="\t"),
alt3prime=read.table("data/modulize/alt3prime.tsv", header=T, sep="\t"),
alt3and5prime=read.table("data/modulize/alt3and5prime.tsv", header=T, sep="\t"),
alternate_first_exon=read.table("data/modulize/alternate_first_exon.tsv", header=T, sep="\t"),
alternate_last_exon=read.table("data/modulize/alternate_last_exon.tsv", header=T, sep="\t"),
p_alt5prime=read.table("data/modulize/p_alt5prime.tsv", header=T, sep="\t"),
p_alt3prime=read.table("data/modulize/p_alt3prime.tsv", header=T, sep="\t"),
p_alternate_first_exon=read.table("data/modulize/p_alternate_first_exon.tsv", header=T, sep="\t"),
p_alternate_last_exon=read.table("data/modulize/p_alternate_last_exon.tsv", header=T, sep="\t"),
mutually_exclusive=read.table("data/modulize/mutually_exclusive.tsv", header=T, sep="\t"),
tandem_cassette=read.table("data/modulize/tandem_cassette.tsv", header=T, sep="\t"),
multi_exon_spanning=read.table("data/modulize/multi_exon_spanning.tsv", header=T, sep="\t"))
saveRDS(binaryEvents, "rds_files/binaryEventsByModulizer.rds")