MyoclonusGraphs

This document contains the code to generate time-to-event graphs for each Myoclonus experiment where myoclonus was observed. Five trials of this experiment were conducted to optimize the Myoclonus protocol to elicit Myoclonic seizures in CSTB mice to study ULD. Data is showcased for Myclonus 2, 3, & 5. Myoclonus 1 is omitted because no seizures occured, myoclonus 4 is missing because I cannot locate the recording or data in onenote or find the designated SD card.

This Graph is for Myoclonus #2. The experiment objective was to determine whether adding ONLY nesting material from the mouse’s respective cage impacted the number of observed Myoclonic Seizures during a 2-hour long recording.

mice3 <- read_xlsx('Myoclonus2090426.xlsx')
datatable(data = mice3)

ggplot(mice3, aes(x = Time, y = Mouse_ID, color = Genotype)) +
  geom_point(size = 3) +
  labs(
    title = "Myoclonus Events Over Time by Mouse",
    x = "Time (minutes)",
    y = "Mouse ID",
    color = "Genotype"
  ) +
  theme_minimal()

This Graph is for Myoclonus 3. The experiment objective was to determine whether the age of the mice impacted the number of Myoclonic Seizures observed during a 60-minute long recording.

mice2 <- read_xlsx('Myoclonus3090426.xlsx')
datatable(data = mice2)

ggplot(mice2, aes(x = Time, y = Mouse_ID, color = Age_wks, shape = Genotype)) +
  geom_point(size = 3) +
  scale_color_gradient(low = "lightblue", high = "darkblue") +
  labs(
    title = "Myoclonus Events Over Time by Mouse",
    x = "Time (hours)",
    y = "Mouse ID",
    shape = "Genotype", color = "Age (weeks)"
  ) +
  theme_minimal()

This Graph is for Myoclonus #6. The experiment objective was to determine whether glass vs. the plastic beaker impacted the number of myoclonic seizures recorded during a 60 minute long recording period- after a 30-minute acclimation period.

mice <- read_excel("Myoclonus5090426.xlsx")
datatable(data = mice)

mice <- mice %>%
  mutate( time_to_event = as_hms(Myoclonus),       time_seconds = as.numeric(time_to_event), time_minutes = time_seconds / 60
  )
mice <- mice %>%
  mutate(had_event = !is.na(time_minutes))  
ggplot(mice, aes(x =time_minutes, y = Beaker_Type, color = Mouse_ID, shape = Genotype)) +
  geom_point(size = 3) +
  labs(
    title = "Myoclonus Events Glass vs. Plastic Beaker",
    x = "Time to Event (Minutes)",
    y = "Beaker Type",
    color = "Mouse_ID"
  ) +
  theme_minimal()

## Warning: Removed 3 rows containing missing values or values outside the scale range
## (`geom_point()`).