Code to test the generation of hierarchies

library(knitr)
library(RMINC)
library(ggplot2)
library(dplyr)
library(tidyr)

load("~/data/Pain-Beggs/Beggs-test-data.RData")

gfBasic <- gf %>%
  filter(Neonatal != "Naive") %>%
  filter(time.to.sac == "3 days")
# mapping files
abijson <- "~/Dropbox/colourbars/allen.json"
defs <- "~/Dropbox/colourbars/Dorr_2008_Steadman_2013_Ullmann_2013_mapping_of_labels.csv"
hdefs <- makeMICeDefsHierachical(defs, abijson) 
vols <- addVolumesToHierarchy(hdefs, allvols)
hLm <- hanatLm(~Sex, gfBasic, vols)
hanatView(hLm, "tvalue.Sexmale", low=2, high=5, symmetric = T, edgeColourFromABI = T)
hanatView(hLm)


Mouse-Imaging-Centre/RMINC documentation built on Nov. 12, 2022, 1:50 p.m.