@@ -154,27 +154,29 @@ write.table(merge.fusion.df, paste(outdir,pheno,"_fusion_gene_table.txt", sep =
154
154
155
155
156
156
# ## extract expression of oncogene,tsg and protein kinase genes
157
- png(paste(outdir , pheno ," _fusion_gene_plot.png" , sep = " " ), width = 800 , height = 700 )
157
+ png(paste(outdir , pheno ," _fusion_gene_plot.png" , sep = " " ), width = 1600 , height = 1400 )
158
+ # pdf(paste(outdir, pheno,"_fusion_gene_plot.pdf", sep = ""), width = 8, height = 7)
158
159
ggplot(merge.fusion.df , aes(x = Phenotype , y = log10(Expression + 1 ))) +
159
160
geom_violin(trim = TRUE ) +
160
- geom_jitter(aes(color = Type ),shape = 16 # , position=position_jitter(0.15)
161
+ geom_jitter(aes(color = Type ),shape = 16 , size = 4 # , position=position_jitter(0.15)
161
162
)+
162
163
theme_bw()+
163
164
scale_colour_manual(name = " Type" , values = c(Oncogene = " #a6bddb" ,Tumor_Suppressor = " #99d8c9" , Both = " #4DAF4A" ,Cancer_Driver = " grey" ))+
164
165
new_scale_color() +
165
- geom_label_repel(data = subset(merge.fusion.df , Type != " Cancer_Driver " & Type != " Others" ),
166
+ geom_label_repel(data = subset(merge.fusion.df , Type != " Others" ),
166
167
aes(label = Gene , fill = Type ), # alpha = 0.5,#fill = "white",
167
168
fontface = ' bold' ,label.size = 0.15 , family = " serif"
168
169
)+
169
170
scale_fill_manual(values = c(Oncogene = " #a6bddb" ,Tumor_Suppressor = " #99d8c9" , Both = " #4DAF4A" ,Cancer_Driver = " grey" ),guide = FALSE )+
170
171
scale_colour_manual(name = " Target" , values = c(left_target = " #636363" ,right_target = " #e6550d" ))+
171
172
scale_x_discrete(name = " " , labels = c(" APre" = " Pre" )) +
172
- theme(axis.text.x = element_text(angle = 0 ,size = 12 ,face = " bold" ,hjust = 0.5 ),
173
- axis.text.y = element_text(size = 12 ,face = " bold" ,hjust = 1 ),
173
+ theme(axis.text.x = element_text(angle = 0 ,size = 20 ,face = " bold" ,hjust = 0.5 ),
174
+ axis.text.y = element_text(size = 20 ,face = " bold" ,hjust = 1 ),
174
175
axis.title.x = element_blank(),
175
- axis.title.y = element_text(size = 12 ,face = " bold" ),
176
- legend.text = element_text(size = 12 ),
177
- strip.text = element_text(size = 12 , face = " bold" ),
176
+ axis.title.y = element_text(size = 20 ,face = " bold" ),
177
+ legend.text = element_text(size = 20 ),
178
+ legend.title = element_text(size = 20 ),
179
+ strip.text = element_text(size = 20 , face = " bold" ),
178
180
legend.position = " right" )
179
181
dev.off()
180
182
0 commit comments