Create and histogram for transcripts per gene

[phil12]

Hello I have to do a project with the following instruction :
Plot the distribution of the number of transcripts per gene.
By distribution, we expect the density or binned histogram of the univariate number of transcript per gene.

For the moment I did this, but this does not work very well ...

library(dplyr)
library(tidyverse)
library(ggplot2)
library(forcats)
gencode %>%
  filter(!is.na(transcript_id)) %>%
  group_by(gene_id) %>%
  summarise(n = n_distinct(transcript_id)) %>%
  ggplot(aes(x = gene_id)) +
  stat_bin(aes(y = "count", label= "count"), geom="text", vjust=-.5) +
  geom_bar(stat="identity")
  geom_histogram(aes(x = n), 
                 col ="red", 
                 fill = "green", 
                 alpha = .2) + 
  geom_density(col = 2) + 
  labs(title = "Histogram of transcript per gene", x = "gene", y = "transcript")

I would be happy if you could help me, thanks a lot

[nirgrahamuk]

Hi, welcome!

Please have a look to our homework policy, homework inspired questions are welcome but they should not include verbatim instructions from your course.

FAQ: Homework Policy

Can I ask questions from a course I am taking here? General questions are always welcome! Please do ask general questions about things like: How to use R How to use the RStudio IDE or RStudio Cloud How to work with tidyverse packages Where to find resources to help you learn or solve problems Specific questions can be OK, if you follow these rules: Never copy-paste instructions from an assignment (even for online courses). Explicitly mention the course you are taking and use the #homework tag. Ask your question as a reproducible example (reprex), preferably prepared for posting using the reprex package. Want to maximize your chances of getting the help you need? Keep readin…

Also, please understand that your code can not be run without suitable data, which is private to you.
look into the following guide for general advice how to overcome such an issue

Finally., you write that your code 'does not work very well', but I think you should be more specific about your issues with it.

[phil12]

I'm very sorry about this actually i modified the instructions so they are not word for word the same but I will not do this anymore, the problem is that is shows mee that stat_bin() needs an astehtic x or y axe an I don't know what to do...