Skip to content

Commit bce60af

Browse files
clarepaciniclarepacini
authored
Initial commit Common essentials calculation
1 parent 5b98688 commit bce60af

1 file changed

Lines changed: 166 additions & 0 deletions

File tree

‎CalcCommonEssentials.R‎

Lines changed: 166 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,166 @@
1+
2+
library(here)
3+
source("EssNonEssfunctions.R")
4+
5+
#load a reference set of essential genes
6+
library(CRISPRcleanR)
7+
data(BAGEL_essential)
8+
data(BAGEL_nonEssential)
9+
10+
source("Combat_HKfunctions.R")
11+
library(ADaM2)
12+
library(ggplot2)
13+
library(gridExtra)
14+
library(magrittr)
15+
library(tidyverse)
16+
library(dplyr)
17+
library(magrittr)
18+
19+
dir.MergeFile<-"./Results"
20+
dir.Results<-"./ResultsFilter"
21+
dir.Input<-"/path/to/downloaded/figshare/"
22+
23+
cmp<-read.csv(paste0(dir.Input,"/model_list_latest.csv"),header=T,stringsAsFactors = F)
24+
cmp2<-cmp
25+
cmp2$model_id<-cmp2$BROAD_ID
26+
cmp<-rbind(cmp,cmp2)
27+
###use Broad 90th depletion method ###
28+
CCR_correctedPC2<-readRDS(file=paste0(dir.MergeFile,"/CCR_SQ_Combat_PC2_All_merge_F.Rds"))
29+
CERES_correctedPC2<-readRDS(file=paste0(dir.MergeFile,"/CERES_SQ_Combat_PC2_All_merge_F.Rds"))
30+
CCRJ_correctedPC2<-readRDS(file=paste0(dir.MergeFile,"/CCR_SQ_Combat_PC2_All_JACKS_merge_F.Rds"))
31+
32+
CCR_correctedPC1<-readRDS(file=paste0(dir.MergeFile,"/CCR_SQ_Combat_PC1_All_merge_F.Rds"))
33+
CERES_correctedPC1<-readRDS(file=paste0(dir.MergeFile,"/CERES_SQ_Combat_PC1_All_merge_F.Rds"))
34+
CCRJ_correctedPC1<-readRDS(file=paste0(dir.MergeFile,"/CCR_SQ_Combat_PC1_All_JACKS_merge_F.Rds"))
35+
36+
CCR_correctedQN<-readRDS(file=paste0(dir.MergeFile,"/CCR_SQ_Combat_All_merge_F.Rds"))
37+
CERES_correctedQN<-readRDS(file=paste0(dir.MergeFile,"/CERES_SQ_Combat_All_merge_F.Rds"))
38+
CCRJ_correctedQN<-readRDS(file=paste0(dir.MergeFile,"/CCR_SQ_Combat_All_JACKS_merge_F.Rds"))
39+
40+
CCR_corrected<-readRDS(file=paste0(dir.MergeFile,"/CCR_SQ_Combat_All_NoNorm_merge_F.Rds"))
41+
CERES_corrected<-readRDS(file=paste0(dir.MergeFile,"/CERES_SQ_Combat_All_NoNorm_merge_F.Rds"))
42+
CCRJ_corrected<-readRDS(file=paste0(dir.MergeFile,"/CCR_SQ_Combat_All_NoNorm_JACKS_merge_F.Rds"))
43+
44+
B90CCR<-Broad90(CCR_corrected,BAGEL_essential,BAGEL_nonEssential)
45+
B90CCRJ<-Broad90(CCRJ_corrected,BAGEL_essential,BAGEL_nonEssential)
46+
B90CERES<-Broad90(CERES_corrected,BAGEL_essential,BAGEL_nonEssential)
47+
48+
B90CCRQN<-Broad90(CCR_correctedQN,BAGEL_essential,BAGEL_nonEssential)
49+
B90CCRJQN<-Broad90(CCRJ_correctedQN,BAGEL_essential,BAGEL_nonEssential)
50+
B90CERESQN<-Broad90(CERES_correctedQN,BAGEL_essential,BAGEL_nonEssential)
51+
52+
B90CCRPC1<-Broad90(CCR_correctedPC1,BAGEL_essential,BAGEL_nonEssential)
53+
B90CCRJPC1<-Broad90(CCRJ_correctedPC1,BAGEL_essential,BAGEL_nonEssential)
54+
B90CERESPC1<-Broad90(CERES_correctedPC1,BAGEL_essential,BAGEL_nonEssential)
55+
56+
B90CCRPC2<-Broad90(CCR_correctedPC2,BAGEL_essential,BAGEL_nonEssential)
57+
B90CCRJPC2<-Broad90(CCRJ_correctedPC2,BAGEL_essential,BAGEL_nonEssential)
58+
B90CERESPC2<-Broad90(CERES_correctedPC2,BAGEL_essential,BAGEL_nonEssential)
59+
60+
save(B90CCR,file=paste0(dir.Results,"/CE_IntCCR90.Rdata"))
61+
save(B90CCRJ,file=paste0(dir.Results,"/CE_IntCCRJ90.Rdata"))
62+
save(B90CERES,file=paste0(dir.Results,"/CE_IntCERES90.Rdata"))
63+
64+
save(B90CCRQN,file=paste0(dir.Results,"/CE_IntCCRQN90.Rdata"))
65+
save(B90CCRJQN,file=paste0(dir.Results,"/CE_IntCCRJQN90.Rdata"))
66+
save(B90CERESQN,file=paste0(dir.Results,"/CE_IntCERESQN90.Rdata"))
67+
68+
save(B90CCRPC1,file=paste0(dir.Results,"/CE_IntCCRPC190.Rdata"))
69+
save(B90CCRJPC1,file=paste0(dir.Results,"/CE_IntCCRJPC190.Rdata"))
70+
save(B90CERESPC1,file=paste0(dir.Results,"/CE_IntCERESPC190.Rdata"))
71+
72+
save(B90CCRPC2,file=paste0(dir.Results,"/CE_IntCCRPC290.Rdata"))
73+
save(B90CCRJPC2,file=paste0(dir.Results,"/CE_IntCCRJPC290.Rdata"))
74+
save(B90CERESPC2,file=paste0(dir.Results,"/CE_IntCERESPC290.Rdata"))
75+
76+
load(paste0(dir.Results,"/BinaryCCRJ4m.RData"))
77+
load(paste0(dir.Results,"/BinaryCCR4m.RData"))
78+
load(paste0(dir.Results,"/BinaryCERES4m.RData"))
79+
80+
PCcoreCERESPC2<-AdamTissue(BinaryCERES4,cmp)
81+
PCcoreCCRPC2<-AdamTissue(BinaryCCR4,cmp)
82+
PCcoreCCRJPC2<-AdamTissue(BinaryCCRJ4,cmp)
83+
84+
save(PCcoreCCRJPC2,file=paste0(dir.Results,"/AdamCoreCCRJPC2.RData"))
85+
save(PCcoreCCRPC2,file=paste0(dir.Results,"/AdamCoreCCRPC2.RData"))
86+
save(PCcoreCERESPC2,file=paste0(dir.Results,"/AdamCoreCERESPC2.RData"))
87+
88+
89+
load(paste0(dir.Results,"/BinaryCCRJ3m.RData"))
90+
load(paste0(dir.Results,"/BinaryCCR3m.RData"))
91+
load(paste0(dir.Results,"/BinaryCERES3m.RData"))
92+
93+
PCcoreCERESPC1<-AdamTissue(BinaryCERES3,cmp)
94+
PCcoreCCRPC1<-AdamTissue(BinaryCCR3,cmp)
95+
PCcoreCCRJPC1<-AdamTissue(BinaryCCRJ3,cmp)
96+
97+
save(PCcoreCCRJPC1,file=paste0(dir.Results,"/AdamCoreCCRJPC1.RData"))
98+
save(PCcoreCCRPC1,file=paste0(dir.Results,"/AdamCoreCCRPC1.RData"))
99+
save(PCcoreCERESPC1,file=paste0(dir.Results,"/AdamCoreCERESPC1.RData"))
100+
101+
load(paste0(dir.Results,"/BinaryCCRJ2m.RData"))
102+
load(paste0(dir.Results,"/BinaryCCR2m.RData"))
103+
load(paste0(dir.Results,"/BinaryCERES2m.RData"))
104+
105+
PCcoreCERESQN<-AdamTissue(BinaryCERES2,cmp)
106+
PCcoreCCRQN<-AdamTissue(BinaryCCR2,cmp)
107+
PCcoreCCRJQN<-AdamTissue(BinaryCCRJ2,cmp)
108+
109+
save(PCcoreCCRJQN,file=paste0(dir.Results,"/AdamCoreCCRJQN.RData"))
110+
save(PCcoreCCRQN,file=paste0(dir.Results,"/AdamCoreCCRQN.RData"))
111+
save(PCcoreCERESQN,file=paste0(dir.Results,"/AdamCoreCERESQN.RData"))
112+
113+
load(paste0(dir.Results,"/BinaryCCRJ1m.RData"))
114+
load(paste0(dir.Results,"/BinaryCCR1m.RData"))
115+
load(paste0(dir.Results,"/BinaryCERES1m.RData"))
116+
117+
PCcoreCERES<-AdamTissue(BinaryCERES1,cmp)
118+
PCcoreCCR<-AdamTissue(BinaryCCR1,cmp)
119+
PCcoreCCRJ<-AdamTissue(BinaryCCRJ1,cmp)
120+
121+
save(PCcoreCCRJ,file=paste0(dir.Results,"/AdamCoreCCRJ.RData"))
122+
save(PCcoreCCR,file=paste0(dir.Results,"/AdamCoreCCR.RData"))
123+
save(PCcoreCERES,file=paste0(dir.Results,"/AdamCoreCERES.RData"))
124+
125+
#load in original Broad and Sanger datasets:
126+
load(paste0(dir.MergeFile,"/BroadDataCERES.Rdata"))
127+
load(paste0(dir.MergeFile,"/SangerDataCERES.Rdata"))
128+
load(paste0(dir.MergeFile,"/BroadData.Rdata"))
129+
load(paste0(dir.MergeFile,"/SangerData.Rdata"))
130+
load(paste0(dir.MergeFile,"/JACKSsanger.Rdata"))
131+
load(paste0(dir.MergeFile,"/JACKSbroad.Rdata"))
132+
133+
dn<-dimnames(BroadDataCERES)
134+
BroadDataCERES<-normalize.quantiles(BroadDataCERES)
135+
dimnames(BroadDataCERES)<-dn
136+
137+
dn<-dimnames(SangerDataCERES)
138+
SangerDataCERES<-normalize.quantiles(SangerDataCERES)
139+
dimnames(SangerDataCERES)<-dn
140+
141+
142+
B90CCR_sanger<-Broad90(SangerData,BAGEL_essential,BAGEL_nonEssential)
143+
B90CCR_broad<-Broad90(BroadData,BAGEL_essential,BAGEL_nonEssential)
144+
B90CERES_sanger<-Broad90(SangerDataCERES,BAGEL_essential,BAGEL_nonEssential)
145+
B90CERES_broad<-Broad90(BroadDataCERES,BAGEL_essential,BAGEL_nonEssential)
146+
147+
save(B90CCR_sanger,file=paste0(dir.Results,"/B90CCR_sanger.RData"))
148+
save(B90CCR_broad,file=paste0(dir.Results,"/B90CCR_broad.RData"))
149+
save(B90CERES_sanger,file=paste0(dir.Results,"/B90CERES_sanger.RData"))
150+
save(B90CERES_broad,file=paste0(dir.Results,"/B90CERES_broad.RData"))
151+
152+
load(paste0(dir.Results,"/BinaryCCRSanger.RData"))
153+
load(paste0(dir.Results,"/BinaryCCRBroad.RData"))
154+
load(paste0(dir.Results,"/BinaryCERESSanger.RData"))
155+
load(paste0(dir.Results,"/BinaryCERESBroad.RData"))
156+
157+
158+
PCcoreCERES_Sanger<-AdamTissue(BinaryCERESSanger,cmp)
159+
PCcoreCERES_Broad<-AdamTissue(BinaryCERESBroad,cmp)
160+
PCcoreCCR_Sanger<-AdamTissue(BinaryCCRSanger,cmp)
161+
PCcoreCCR_Broad<-AdamTissue(BinaryCCRBroad,cmp)
162+
163+
save(PCcoreCERES_Sanger,file=paste0(dir.Results,"/AdamCoreCERES_Sanger.RData"))
164+
save(PCcoreCERES_Broad,file=paste0(dir.Results,"/AdamCoreCERES_Broad.RData"))
165+
save(PCcoreCCR_Sanger,file=paste0(dir.Results,"/AdamCoreCCR_Sanger.RData"))
166+
save(PCcoreCCR_Broad,file=paste0(dir.Results,"/AdamCoreCCR_Broad.RData"))

0 commit comments

Comments
 (0)