ex 22

RNA-sequencing intro

Examine count data

d <- read_csv(here("data","unfilt_counts.csv.gz")) |> as.matrix()


df <- tibble(variance=??,
                mean=??)

ggplot(??) +
        geom_point(aes(x=??, y=??)) + 
        scale_y_log10(limits = c(1,1e9)) +
        scale_x_log10(limits = c(1,1e9)) +
        geom_abline(intercept = 0, slope = 1, color="red") +
  theme_cowplot()

estimateSizeFactors

d <- read_csv(here("data","unfilt_counts.csv.gz")) |> as.matrix()

# estimate size factors