genes ketone ADT treatment

cargo los paquetes y la data

selecciono los genes y los datasets

hago los graficos

  
# HAGO LOS GRAFICO
genes_int = genes.quiero


for (i in 1:length(genes_int)){
  # preparo la data
  gen = genes_int[i]
  exp.gen = exp[,gen]
  exp.gen.l = log2(exp.gen + 1)
  grupo = clin$grupo
  experimento = clin$Experimento
  datitos = data.frame(grupo, experimento, exp.gen, exp.gen.l )

  # grafico
  ploteo = ggplot(data = datitos %>% filter(experimento %in% ensayos.quiero ), aes(x=grupo, y=exp.gen.l )) +
    geom_jitter(aes(color=experimento), size=3, width = 0.1)+
    xlab(NULL) +
    ylab(paste(gen, "expression \n log2 (norm counts +1)")) +
    theme(legend.position = "bottom") +
    theme_bw() +
    theme(axis.text = element_text(size = 10),
          axis.title = element_text(size = 10),
          plot.title =element_text(size = 15),
          axis.text.x = element_text(angle = 45, hjust=1)) +#,
          #legend.position = 'none') +   # para sacar la leyenda
    stat_summary(fun=mean,
                 geom="point",
                 shape= '_',
                 size=10)
  # imprimo
    print(ploteo)
    
    # lo agrego a pw
  
  # guardo el plot
  assign(paste(gen,'ploteo_DGE', sep = '_'), ploteo  )
}

componer una figura unificada, esta parte la estoy desarrollando ahora