What may cause my z-score-FDR distribution
0
0
Entering edit mode
4 months ago
Pallondyle • 0

I'm trying to recurrent the FDR-Zscore dot plot in this article: Genome-wide CRISPR Screens Reveal Host Factors Critical for SARS-CoV-2 Infection, figure 1.B enter image description here I used the MaGeck sgRNA test data as my input (DIFF):

plot_dat <- DIFF %>%mutate(z_score = (LFC- mean(LFC)) / sd(LFC))
plot_dat <- plot_dat %>% mutate(fdr=ifelse(`pos.fdr`<`neg.fdr`,`pos.fdr`,`neg.fdr`))
plot_nont<-plot_dat[which(plot_dat$Gene=="Non-Targeting Control"),]
ggplot(plot_dat,aes(x=z_score,y=-log10(FDR)))+geom_point()+
  geom_point(data=plot_nont,color="red")+
  scale_y_continuous(limits=c(0,5))

and my plot is like this. enter image description here enter image description here I want to know what causes my data distributed like this. Is there any problem with my raw data?

FDR Z-score CRISPR_screen • 248 views
ADD COMMENT
0
Entering edit mode

Nobody can tell you. Did you run the exact same code they published? Same normalization? I don't see how people could help here, we cannot magically see your screen or the code from whatever paper this is.

ADD REPLY

Login before adding your answer.

Traffic: 2234 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6