50 Box Plots intensity

Load and Combine Data

folder <- "C:/Users/jeffe/OneDrive/Documents/VTPP 234/Total_MRNA"

files <- list.files(folder, pattern = "*.txt", full.names = TRUE)

all_data <- files %>%
  map_dfr(~ read.delim(.x, header = TRUE, sep = "\t", skip = 8) %>%
            filter(Raw_intensity != -9999, Normalized_intensity != -9999) %>%
            mutate(SourceFile = basename(.x)))

Reshape Data

plot_data <- all_data %>%
  select(SourceFile, Raw_intensity, Normalized_intensity) %>%
  pivot_longer(cols = c(Raw_intensity, Normalized_intensity),
               names_to = "Type",
               values_to = "Intensity")

Boxplots: Separate Panels

## Warning in transformation$transform(x): NaNs produced

## Warning in scale_y_log10(): log-10 transformation introduced infinite values.

## Warning: Removed 17868 rows containing non-finite outside the scale range
## (`stat_boxplot()`).