Hi friends
I plotted this box-wisker for TCGA HTSeq data in R.
I want to have harf of them as red and half as blue (control vs treatment groups). or is there any better way for boxplot?
How can I do that?
I just used simple code:

    data <- read.table("mydata.csv", header=TRUE, row.names = 1)
    library(DESeq2)
    mat <- as.matrix(data)
    log <- rlog(mat)
    boxplot(log)

my plot:
https://ibb.co/xXsHtC0
maybe this one? A dodged box plot?

I used  "ToothGrowth" public data for this demonstration:

    data("ToothGrowth") 
    library(ggplot2)
    
    ggplot(ToothGrowth, aes(x = as.factor(dose), y = len)) + 
    geom_boxplot(aes(fill = supp), position = position_dodge(0.9)) +
    scale_fill_manual(values = c("#09E359", "#E31009")) + 
    theme_bw()

![enter image description here][1]

Also, a dogged violin plot would be a very informative way to present expression data. Personally, I would rather use a split violin plot for expression data when I have many cases to be plotted.  The codes for creating a split violin plot in python can be found on my [GitHub,][2] which I used to plot top DE genes from RNA-seq.



  [1]: https://raw.githubusercontent.com/hamidghaedi/file_sharing/main/box_demo.JPG
  [2]: https://github.com/hamidghaedi/RNA-seq-differential-expression/blob/master/violin_plot_DE.ipynb
  [3]: https://i.stack.imgur.com/tYQWX.png