## -----------------------------------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>",
  eval = FALSE
)

## -----------------------------------------------------------------------------
# library(ggseg)
# library(ggplot2)
# library(dplyr)

## -----------------------------------------------------------------------------
# subjects_dir <- Sys.getenv("SUBJECTS_DIR")
# stats_file <- file.path(subjects_dir, "bert/stats/lh.aparc.stats")
# 
# data <- read_freesurfer_stats(stats_file)
# data

## -----------------------------------------------------------------------------
# ggplot(data |> mutate(label = paste0("lh_", label))) +
#   geom_brain(atlas = dk(), mapping = aes(fill = ThickAvg))

## -----------------------------------------------------------------------------
# dat <- read_atlas_files(subjects_dir, "aparc.stats$")
# dat

## -----------------------------------------------------------------------------
# ggplot(dat) +
#   geom_brain(atlas = dk(), mapping = aes(fill = ThickStd))

## -----------------------------------------------------------------------------
# library(tidyr)
# 
# dat |>
#   pivot_longer(-c(subject, label), names_to = "stat", values_to = "val") |>
#   ggplot() +
#   geom_brain(atlas = dk(), mapping = aes(fill = val)) +
#   facet_wrap(~stat)

## -----------------------------------------------------------------------------
# table_path <- "path/to/aparc.volume.table"
# 
# read_freesurfer_table(table_path)

## -----------------------------------------------------------------------------
# dat <- read_freesurfer_table(table_path, measure = "volume")
# dat

## -----------------------------------------------------------------------------
# ggplot(dat |> filter(grepl("lh|rh", label))) +
#   geom_brain(atlas = dk(), mapping = aes(fill = volume))