R Notebook

Inter-system correlations - TOTAL

Conditions in this set of graphs: - 1. ‘Total’ cells only - 2. Density measures - 3. Correlations use “complete.obs” - 4. Ordering by “hclust” - 5. Removed bad freezers from young group

# source("LoadLibs.R")
library(ggpubr)
library(grid)
library(gridExtra)
cluster.method = "hclust"
use.method = "complete.obs"
corr.type="pearson"

Define functions for getting correlation tests, making heatmaps

# DEFINE FUNCTION FOR GETTING CORRELATION TEST
get_corrtest_sys = function(df) {
  data = df %>%
    # mutate(grain2 = paste(CHANNEL, System)) %>%
    ddply(.(CHANNEL, ID), transform, z = scale(weighted)) %>%
    dcast(ID ~ System, value.var="z", fun.aggregate = mean) 
  
  # IMPUTE MISSING VALUES 
  data = complete(mice(data, m=5, seed=500))
  
  # rcorr = cor(data[,2:length(data)], data[,2:length(data)], use=use.method, method=corr.type)
    rcorr = rcorr(as.matrix(data[,2:length(data)]), type=corr.type)
    # return(list(data=data, rcorr=rcorr))
    return(rcorr)
}
# DEFINE FUNCTION FOR MAKING HEATMAP
make_heatmap = function(test.results, title) {
  # inputs the results from 'get_corrtest'
  # matrix = test.results$r
  heatmap.2(test.results, symm=T, 
            col=rev(brewer.pal(11,"RdBu")), main=title, 
            trace="none", revC=T, margins=c(10,10),
            key=T, scale="none"
            )
}
# DEFINE FUNCTION FOR MAKING HEATMAP FOR NUMBER ANIMALS
make_heatmap2 = function(test.results, title) {
  # inputs the results from 'get_corrtest'
  # matrix = test.results$r
  # heatmap(test.results, symm=T, col=heat.colors, main=title)
  
  corrplot(test.results$r, title=title, p.mat=test.results$P, insig="label_sig",pch="*", sig.level=0.05,
           order=cluster.method, method="color", mar=c(0,0,1,0), pch.col="white",pch.cex = 2,
           diag=T, tl.pos = "l")
}
compare_matrix_r(arcodd, arceven, "odd", "even")

[1] "n1=" "24" 
[1] "n2=" "27" 
Tests of correlation matrices 
Call:cortest.normal(R1 = x$r, R2 = y$r, n1 = min(x$n), n2 = min(y$n))
 Chi Square value 111.84  with df =  105   with probability < 0.31 

Define Corr. Comparisons

flattenCorrMatrix <- function(cormat, pmat) {
  ut <- upper.tri(cormat)
  data.frame(
    row = rownames(cormat)[row(cormat)[ut]],
    column = rownames(cormat)[col(cormat)[ut]],
    cor  =(cormat)[ut],
    p = pmat[ut]
  )
}
compare_sigs = function(mat1, mat2) {
}
# DEFINE MATRIX FOR MAKING GENERAL COMPARISONS
compare_matrix_r = function(x,y, namex, namey) {
  
  # MAKE COR MATS
  make_heatmap2(x, namex)
  make_heatmap2(y, namey)
  # print n
  print(paste(c("n1=", min(x$n))))
  print(paste(c("n2=", min(y$n))))
  
  # STEIGER TEST STATS
  print(
    cortest.normal(R1 = x$r,
                 R2 = y$r,
                 n1 = min(x$n),
                 n2 = min(y$n))
  )
  
  # print(
  #   cortest.mat(R1 = x$r,
  #                R2 = y$r,
  #                n1 = min(x$n),
  #                n2 = min(y$n))
  # )
    
  # PRINT DIFFERENCE MATRIX
  yx = y$r - x$r
  title = paste("Difference (", namey, " - ", namex, ")")
  corrplot(yx, title=title, order=cluster.method, method="color", 
           mar=c(0,0,1,0), diag=T, tl.pos = "l") 
  
  
  # PRINT T-TEST PAIRED
  x.vector = as.vector(x$r[lower.tri(x$r, diag=F)])
  y.vector = as.vector(y$r[lower.tri(y$r, diag=F)])
  # print(t.test(fisherz(x.vector), 
  #            fisherz(y.vector), 
  #            paired=T))
  
    # PLOT DISTRIBUTION OF R VALUES
  df = data.frame(var = c(rep(namex, length(x.vector)), rep(namey, length(y.vector))), 
               values = c(x.vector, y.vector)) 
  print(
    df %>%
      ggplot() + aes(values, group=var, fill=var, alpha=0.5) + geom_density() + 
      theme_minimal() + xlim(-1,1) + xlab("Pearson's R") +
      ylab("Frequency") + ggtitle("Distribution of R values")
  )
    print(
    df %>%
      mutate(z = fisherz(as.vector(values))) %>%
      ggplot() + aes(var, z) +
      stat_summary(fun.y = "mean", geom="bar") +
      stat_summary(fun.data = "mean_se", geom="errorbar") +
      theme_minimal())
  
  flat1 = flattenCorrMatrix(x$r, x$P) %>% mutate(comp = paste(row, column))
  flat2 = flattenCorrMatrix(y$r, y$P) %>% mutate(comp = paste(row, column))
  flat.both = merge(flat1, flat2 %>% select(p, comp), by="comp")
  flat.both = flat.both %>% mutate(p.delta = case_when(
    p.x < 0.05 & p.y < 0.05 ~ "s-s",
    p.x < 0.05 & p.y > 0.05 ~ "s-ns",
    p.x > 0.05 & p.y < 0.05 ~ "ns-s",
    p.x > 0.05 & p.y > 0.05 ~ "ns-ns")
  )
  results = summary(as.factor(flat.both$p.delta))
  results
  
  barplot(results, main=paste("Change: ", namex, "to", namey))
  
}

Define PCA function

get_data_matrix = function(df) {
  data = df %>%
    # mutate(grain2 = paste(CHANNEL, System)) %>%
    ddply(.(System), transform, z = scale(weighted)) %>%
    dcast(ID ~ System, value.var="z", fun.aggregate = mean, .drop=T) 
  data
}
get_pca = function(df, dfname) {
  
  adm = df %>% get_data_matrix()
  adm[is.na(adm)] = 0
  rownames(adm) = adm$ID # make ID the rownames
  adm = adm %>% select(2:16) # remove ID
  pr = prcomp(na.omit(adm), scale=T)
  print(summary(pr))
  print(biplot(pr))
  
  
  plot(pr$x[,1], pr$x[,2]) # 
  pr.var = pr$sdev^2 # get variance by squaring eigenvalues
  pr.var.pct = round(pr.var/(sum(pr.var)) * 100, 1) # get variance %
  
  print(
  data.frame(name = paste(1:length(pr.var.pct)), percent = pr.var.pct) %>%
    ggplot() + geom_bar(stat="identity", show.legend = F) + aes(sort(as.numeric(name), decreasing=F), percent, fill=name) + 
    ggtitle(paste(dfname, "- Scree plot")) + 
    xlab("Principal Component") + ylab("% variance explained") + 
    theme_minimal()
)
  
  
  pr.data = data.frame(Sample=rownames(pr$x), X = pr$x[,1], Y = pr$x[,2])
  pr.data$ID = pr.data$Sample
  
  pr.data=merge(pr.data, animal_details[2:length(animal_details)])
  pr.data$GRAIN = paste(pr.data$GROUP, pr.data$VALENCE)
  
  
  # PRINT SOME GRAPHS
  
  print(
    ggplot(pr.data) + aes(X, Y, label=GRAIN,color=GRAIN) + geom_point() + theme_minimal()  + stat_ellipse() + facet_wrap(~EXPT) +
    xlab("PC1") + ylab("PC2") + ggtitle(paste(dfname,"PCA by experimental group"))
  )
  
  print(
    ggplot(pr.data) + aes(X, Y, label=FREEZE,color=VALENCE) + geom_text() + theme_minimal() + stat_ellipse() +
      xlab("PC1") + ylab("PC2") + ggtitle(paste(dfname, "Component space, by shock/context, labeled freezing score"))
  )
    
  print(
    ggplot(pr.data) + aes(X, Y, label=FREEZE,color=VALENCE) + geom_text() + theme_minimal() + stat_ellipse() + 
      xlab("PC1") + ylab("PC2") + ggtitle(paste(dfname, "Component space, by young/old, labeled freezing score")) + 
      facet_wrap(~EXPT)
  )
  
  pr.load = pr$rotation[,1] # PC1
  sys.scores = abs(pr.load)
  sys.scores.ranked = sort(sys.scores, decreasing=T)
  as.matrix(pr.load[names(sys.scores.ranked)])
}

ARCIHC

# SUBSET FOR ARC ONLY
arc = sys.weighted %>% filter(CHANNEL=='ArcIHC') 
get_pca(arc, "ArcIHC")

Importance of components:
                          PC1     PC2     PC3     PC4     PC5     PC6     PC7     PC8     PC9    PC10    PC11    PC12
Standard deviation     3.4199 0.95121 0.94046 0.69454 0.57504 0.41544 0.37900 0.35300 0.27889 0.23505 0.20277 0.18284
Proportion of Variance 0.7797 0.06032 0.05896 0.03216 0.02205 0.01151 0.00958 0.00831 0.00519 0.00368 0.00274 0.00223
Cumulative Proportion  0.7797 0.84002 0.89898 0.93114 0.95319 0.96469 0.97427 0.98258 0.98776 0.99144 0.99418 0.99641
                          PC13    PC14    PC15
Standard deviation     0.16274 0.13273 0.09848
Proportion of Variance 0.00177 0.00117 0.00065
Cumulative Proportion  0.99818 0.99935 1.00000
NULL

                     [,1]
HYPOTHAL       -0.2799714
THALAMIC       -0.2789033
TEMP ASSOC CTX -0.2769579
OLFACTORY      -0.2758282
VISC CTX       -0.2704783
SOMAT CTX      -0.2703665
CORTICO-AMY    -0.2698390
NUCLEAR AMY    -0.2635804
POST ASSOC CTX -0.2622890
MOTOR CTX      -0.2619307
AUD CTX        -0.2550255
SUBIC-HPC      -0.2320983
MIDBRAIN       -0.2230450
PFC            -0.2217221
STRIATUM       -0.2174334

Is there a difference between odd/even animals? That should be a control. ### Control

arcodd = arc %>%
  filter((ID %%2) == T) %>%
  get_corrtest_sys()

 iter imp variable
  1   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  1   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  1   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  1   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  1   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  2   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  2   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  2   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  2   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  2   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  3   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  3   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  3   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  3   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  3   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  4   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  4   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  4   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  4   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  4   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  5   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  5   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  5   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  5   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  5   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC

arceven = arc %>%
  filter((ID%%2) ==F) %>%
  get_corrtest_sys()

 iter imp variable
  1   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  1   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  1   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  1   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  1   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  2   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  2   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  2   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  2   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  2   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  3   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  3   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  3   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  3   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  3   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  4   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  4   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  4   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  4   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  4   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  5   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  5   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  5   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  5   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  5   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX

compare_matrix_r(arcodd, arceven, "odd", "even")

[1] "n1=" "24" 
[1] "n2=" "27" 
Tests of correlation matrices 
Call:cortest.normal(R1 = x$r, R2 = y$r, n1 = min(x$n), n2 = min(y$n))
 Chi Square value 111.84  with df =  105   with probability < 0.31 

Arc::VALENCE stats

av1 = arc %>%
  filter(VALENCE == 'CONTEXT') %>%
  get_corrtest_sys()

 iter imp variable
  1   1  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  1   2  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  1   3  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  1   4  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  1   5  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  2   1  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  2   2  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  2   3  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  2   4  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  2   5  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  3   1  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  3   2  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  3   3  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  3   4  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  3   5  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  4   1  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  4   2  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  4   3  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  4   4  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  4   5  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  5   1  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  5   2  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  5   3  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  5   4  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  5   5  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC

av2 = arc %>%
  filter(VALENCE == 'SHOCK') %>%
  get_corrtest_sys()

 iter imp variable
  1   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX

compare_matrix_r(av1, av2, "Context", "Shock")

[1] "n1=" "26" 
[1] "n2=" "25" 
Tests of correlation matrices 
Call:cortest.normal(R1 = x$r, R2 = y$r, n1 = min(x$n), n2 = min(y$n))
 Chi Square value 182.18  with df =  105   with probability < 4.5e-06 

Arc::EXPT

# Get Z dist of CONTEXT
arc.age1 = arc %>%
  filter(EXPT == 'YOUNG') %>%
  get_corrtest_sys() 

 iter imp variable
  1   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX

arc.age2 = arc %>%
  filter(EXPT == 'OLD') %>%
  get_corrtest_sys() 

 iter imp variable
  1   1  AUD CTX  MIDBRAIN
  1   2  AUD CTX  MIDBRAIN
  1   3  AUD CTX  MIDBRAIN
  1   4  AUD CTX  MIDBRAIN
  1   5  AUD CTX  MIDBRAIN
  2   1  AUD CTX  MIDBRAIN
  2   2  AUD CTX  MIDBRAIN
  2   3  AUD CTX  MIDBRAIN
  2   4  AUD CTX  MIDBRAIN
  2   5  AUD CTX  MIDBRAIN
  3   1  AUD CTX  MIDBRAIN
  3   2  AUD CTX  MIDBRAIN
  3   3  AUD CTX  MIDBRAIN
  3   4  AUD CTX  MIDBRAIN
  3   5  AUD CTX  MIDBRAIN
  4   1  AUD CTX  MIDBRAIN
  4   2  AUD CTX  MIDBRAIN
  4   3  AUD CTX  MIDBRAIN
  4   4  AUD CTX  MIDBRAIN
  4   5  AUD CTX  MIDBRAIN
  5   1  AUD CTX  MIDBRAIN
  5   2  AUD CTX  MIDBRAIN
  5   3  AUD CTX  MIDBRAIN
  5   4  AUD CTX  MIDBRAIN
  5   5  AUD CTX  MIDBRAIN

compare_matrix_r(arc.age1, arc.age2, "Young", "Old")

[1] "n1=" "24" 
[1] "n2=" "27" 
Tests of correlation matrices 
Call:cortest.normal(R1 = x$r, R2 = y$r, n1 = min(x$n), n2 = min(y$n))
 Chi Square value 163.65  with df =  105   with probability < 0.00022 

Arc::CONTEXT

# Get Z dist of CONTEXT
arc.c1 = arc %>%
  filter(CONTEXT == 'AB') %>%
  get_corrtest_sys() 

 iter imp variable
  1   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  1   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  1   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  1   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  1   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  2   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  2   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  2   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  2   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  2   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  3   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  3   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  3   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  3   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  3   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  4   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  4   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  4   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  4   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  4   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  5   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  5   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  5   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  5   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  5   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC

arc.c2 = arc %>%
  filter(CONTEXT == 'AA') %>%
  get_corrtest_sys() 

 iter imp variable
  1   1  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  SOMAT CTX  STRIATUM  VISC CTX
  1   2  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  SOMAT CTX  STRIATUM  VISC CTX
  1   3  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  SOMAT CTX  STRIATUM  VISC CTX
  1   4  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  SOMAT CTX  STRIATUM  VISC CTX
  1   5  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  SOMAT CTX  STRIATUM  VISC CTX
  2   1  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  SOMAT CTX  STRIATUM  VISC CTX
  2   2  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  SOMAT CTX  STRIATUM  VISC CTX
  2   3  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  SOMAT CTX  STRIATUM  VISC CTX
  2   4  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  SOMAT CTX  STRIATUM  VISC CTX
  2   5  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  SOMAT CTX  STRIATUM  VISC CTX
  3   1  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  SOMAT CTX  STRIATUM  VISC CTX
  3   2  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  SOMAT CTX  STRIATUM  VISC CTX
  3   3  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  SOMAT CTX  STRIATUM  VISC CTX
  3   4  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  SOMAT CTX  STRIATUM  VISC CTX
  3   5  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  SOMAT CTX  STRIATUM  VISC CTX
  4   1  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  SOMAT CTX  STRIATUM  VISC CTX
  4   2  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  SOMAT CTX  STRIATUM  VISC CTX
  4   3  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  SOMAT CTX  STRIATUM  VISC CTX
  4   4  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  SOMAT CTX  STRIATUM  VISC CTX
  4   5  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  SOMAT CTX  STRIATUM  VISC CTX
  5   1  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  SOMAT CTX  STRIATUM  VISC CTX
  5   2  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  SOMAT CTX  STRIATUM  VISC CTX
  5   3  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  SOMAT CTX  STRIATUM  VISC CTX
  5   4  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  SOMAT CTX  STRIATUM  VISC CTX
  5   5  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  SOMAT CTX  STRIATUM  VISC CTX

compare_matrix_r(arc.c1, arc.c2, "AB", "AA")

[1] "n1=" "25" 
[1] "n2=" "26" 
Tests of correlation matrices 
Call:cortest.normal(R1 = x$r, R2 = y$r, n1 = min(x$n), n2 = min(y$n))
 Chi Square value 136.32  with df =  105   with probability < 0.022 

Arc - EXPT:VALENCE

arc.ys = arc %>%
  filter(EXPT=='YOUNG') %>%
  filter(VALENCE=='SHOCK') %>%
  get_corrtest_sys()

 iter imp variable
  1   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX

arc.yc = arc %>%
  filter(EXPT=='YOUNG') %>%
  filter(VALENCE=='CONTEXT') %>%
  get_corrtest_sys()

 iter imp variable
  1   1  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  1   2  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  1   3  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  1   4  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  1   5  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  2   1  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  2   2  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  2   3  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  2   4  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  2   5  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  3   1  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  3   2  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  3   3  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  3   4  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  3   5  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  4   1  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  4   2  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  4   3  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  4   4  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  4   5  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  5   1  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  5   2  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  5   3  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  5   4  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  5   5  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC

arc.os = arc %>%
  filter(EXPT=='OLD') %>%
  filter(VALENCE=='SHOCK') %>%
  get_corrtest_sys()

 iter imp variable
  1   1  AUD CTX  MIDBRAIN
  1   2  AUD CTX  MIDBRAIN
  1   3  AUD CTX  MIDBRAIN
  1   4  AUD CTX  MIDBRAIN
  1   5  AUD CTX  MIDBRAIN
  2   1  AUD CTX  MIDBRAIN
  2   2  AUD CTX  MIDBRAIN
  2   3  AUD CTX  MIDBRAIN
  2   4  AUD CTX  MIDBRAIN
  2   5  AUD CTX  MIDBRAIN
  3   1  AUD CTX  MIDBRAIN
  3   2  AUD CTX  MIDBRAIN
  3   3  AUD CTX  MIDBRAIN
  3   4  AUD CTX  MIDBRAIN
  3   5  AUD CTX  MIDBRAIN
  4   1  AUD CTX  MIDBRAIN
  4   2  AUD CTX  MIDBRAIN
  4   3  AUD CTX  MIDBRAIN
  4   4  AUD CTX  MIDBRAIN
  4   5  AUD CTX  MIDBRAIN
  5   1  AUD CTX  MIDBRAIN
  5   2  AUD CTX  MIDBRAIN
  5   3  AUD CTX  MIDBRAIN
  5   4  AUD CTX  MIDBRAIN
  5   5  AUD CTX  MIDBRAIN

arc.oc = arc %>%
  filter(EXPT=='OLD') %>%
  filter(VALENCE=='CONTEXT') %>%
  get_corrtest_sys()

 iter imp variable
  1   1  MIDBRAIN
  1   2  MIDBRAIN
  1   3  MIDBRAIN
  1   4  MIDBRAIN
  1   5  MIDBRAIN
  2   1  MIDBRAIN
  2   2  MIDBRAIN
  2   3  MIDBRAIN
  2   4  MIDBRAIN
  2   5  MIDBRAIN
  3   1  MIDBRAIN
  3   2  MIDBRAIN
  3   3  MIDBRAIN
  3   4  MIDBRAIN
  3   5  MIDBRAIN
  4   1  MIDBRAIN
  4   2  MIDBRAIN
  4   3  MIDBRAIN
  4   4  MIDBRAIN
  4   5  MIDBRAIN
  5   1  MIDBRAIN
  5   2  MIDBRAIN
  5   3  MIDBRAIN
  5   4  MIDBRAIN
  5   5  MIDBRAIN

Arc Young CTX vs SHK

compare_matrix_r(arc.yc, arc.ys, "youngCTX", "youngSHK")

[1] "n1=" "13" 
[1] "n2=" "11" 
Tests of correlation matrices 
Call:cortest.normal(R1 = x$r, R2 = y$r, n1 = min(x$n), n2 = min(y$n))
 Chi Square value 83.86  with df =  105   with probability < 0.94 

Error in corrplot(yx, title = title, order = cluster.method, method = "color",  : 
  The matrix is not in [-1, 1]!

Arc Old CTX vs SHK = ns

Arc OLD CTX vs YOUNG CTX

Arc OLD SHK vs YOUNG SHK

H2BGFP

# SUBSET FOR ARC ONLY
gfp = sys.weighted %>% filter(CHANNEL=='H2BGFP') 
get_pca(gfp, "H2BGFP")

Importance of components:
                          PC1    PC2    PC3     PC4     PC5     PC6     PC7     PC8     PC9    PC10    PC11    PC12
Standard deviation     2.4466 1.5400 1.2608 1.01123 0.90372 0.80768 0.80353 0.69973 0.63959 0.56645 0.46090 0.45099
Proportion of Variance 0.3991 0.1581 0.1060 0.06817 0.05445 0.04349 0.04304 0.03264 0.02727 0.02139 0.01416 0.01356
Cumulative Proportion  0.3991 0.5572 0.6631 0.73131 0.78576 0.82925 0.87229 0.90494 0.93221 0.95360 0.96776 0.98132
                          PC13    PC14    PC15
Standard deviation     0.35095 0.30867 0.24851
Proportion of Variance 0.00821 0.00635 0.00412
Cumulative Proportion  0.98953 0.99588 1.00000
NULL

                      [,1]
AUD CTX         0.35647194
TEMP ASSOC CTX  0.33477023
PFC             0.33278290
VISC CTX        0.32270837
MOTOR CTX       0.31308514
SUBIC-HPC       0.29909351
POST ASSOC CTX  0.29015752
SOMAT CTX       0.28370508
OLFACTORY       0.26676408
NUCLEAR AMY     0.24880701
STRIATUM        0.19131378
MIDBRAIN       -0.12398956
CORTICO-AMY     0.07308214
THALAMIC       -0.05343342
HYPOTHAL       -0.02395417

Control

godd = gfp %>%
  filter((ID %%2) == T) %>%
  get_corrtest_sys()

 iter imp variable
  1   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  1   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  1   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  1   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  1   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  2   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  2   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  2   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  2   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  2   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  3   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  3   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  3   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  3   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  3   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  4   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  4   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  4   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  4   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  4   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  5   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  5   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  5   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  5   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  5   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC

geven = gfp %>%
  filter((ID%%2) ==F) %>%
  get_corrtest_sys()

 iter imp variable
  1   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  1   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  1   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  1   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  1   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  2   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  2   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  2   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  2   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  2   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  3   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  3   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  3   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  3   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  3   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  4   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  4   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  4   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  4   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  4   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  5   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  5   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  5   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  5   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  5   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX

compare_matrix_r(godd, geven, "odd", "even")

[1] "n1=" "24" 
[1] "n2=" "27" 
Tests of correlation matrices 
Call:cortest.normal(R1 = x$r, R2 = y$r, n1 = min(x$n), n2 = min(y$n))
 Chi Square value 106.75  with df =  105   with probability < 0.43 

GFP::VALENCE stats

gv1 = gfp %>%
  filter(VALENCE == 'CONTEXT') %>%
  get_corrtest_sys()

 iter imp variable
  1   1  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  1   2  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  1   3  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  1   4  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  1   5  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  2   1  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  2   2  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  2   3  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  2   4  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  2   5  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  3   1  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  3   2  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  3   3  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  3   4  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  3   5  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  4   1  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  4   2  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  4   3  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  4   4  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  4   5  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  5   1  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  5   2  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  5   3  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  5   4  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  5   5  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC

gv2 = gfp %>%
  filter(VALENCE == 'SHOCK') %>%
  get_corrtest_sys()

 iter imp variable
  1   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX

compare_matrix_r(gv1, gv2, "context", "shock")

[1] "n1=" "26" 
[1] "n2=" "25" 
Tests of correlation matrices 
Call:cortest.normal(R1 = x$r, R2 = y$r, n1 = min(x$n), n2 = min(y$n))
 Chi Square value 130.58  with df =  105   with probability < 0.046 

GFP::EXPT

gfp.age1 = gfp %>%
  filter(EXPT == 'YOUNG') %>%
  get_corrtest_sys() 

 iter imp variable
  1   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX

gfp.age2 = gfp %>%
  filter(EXPT == 'OLD') %>%
  get_corrtest_sys()

 iter imp variable
  1   1  AUD CTX  MIDBRAIN
  1   2  AUD CTX  MIDBRAIN
  1   3  AUD CTX  MIDBRAIN
  1   4  AUD CTX  MIDBRAIN
  1   5  AUD CTX  MIDBRAIN
  2   1  AUD CTX  MIDBRAIN
  2   2  AUD CTX  MIDBRAIN
  2   3  AUD CTX  MIDBRAIN
  2   4  AUD CTX  MIDBRAIN
  2   5  AUD CTX  MIDBRAIN
  3   1  AUD CTX  MIDBRAIN
  3   2  AUD CTX  MIDBRAIN
  3   3  AUD CTX  MIDBRAIN
  3   4  AUD CTX  MIDBRAIN
  3   5  AUD CTX  MIDBRAIN
  4   1  AUD CTX  MIDBRAIN
  4   2  AUD CTX  MIDBRAIN
  4   3  AUD CTX  MIDBRAIN
  4   4  AUD CTX  MIDBRAIN
  4   5  AUD CTX  MIDBRAIN
  5   1  AUD CTX  MIDBRAIN
  5   2  AUD CTX  MIDBRAIN
  5   3  AUD CTX  MIDBRAIN
  5   4  AUD CTX  MIDBRAIN
  5   5  AUD CTX  MIDBRAIN

compare_matrix_r(gfp.age1, gfp.age2, "young", "old")

[1] "n1=" "24" 
[1] "n2=" "27" 
Tests of correlation matrices 
Call:cortest.normal(R1 = x$r, R2 = y$r, n1 = min(x$n), n2 = min(y$n))
 Chi Square value 91.09  with df =  105   with probability < 0.83 

GFP - EXPT:VALENCE

gfp.ys = gfp %>%
  filter(EXPT=='YOUNG') %>%
  filter(VALENCE=='SHOCK') %>%
  get_corrtest_sys()

 iter imp variable
  1   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX

gfp.yc = gfp %>%
  filter(EXPT=='YOUNG') %>%
  filter(VALENCE=='CONTEXT') %>%
  get_corrtest_sys()

 iter imp variable
  1   1  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  1   2  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  1   3  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  1   4  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  1   5  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  2   1  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  2   2  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  2   3  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  2   4  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  2   5  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  3   1  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  3   2  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  3   3  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  3   4  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  3   5  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  4   1  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  4   2  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  4   3  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  4   4  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  4   5  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  5   1  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  5   2  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  5   3  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  5   4  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  5   5  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC

gfp.os = gfp %>%
  filter(EXPT=='OLD') %>%
  filter(VALENCE=='SHOCK') %>%
  get_corrtest_sys()

 iter imp variable
  1   1  AUD CTX  MIDBRAIN
  1   2  AUD CTX  MIDBRAIN
  1   3  AUD CTX  MIDBRAIN
  1   4  AUD CTX  MIDBRAIN
  1   5  AUD CTX  MIDBRAIN
  2   1  AUD CTX  MIDBRAIN
  2   2  AUD CTX  MIDBRAIN
  2   3  AUD CTX  MIDBRAIN
  2   4  AUD CTX  MIDBRAIN
  2   5  AUD CTX  MIDBRAIN
  3   1  AUD CTX  MIDBRAIN
  3   2  AUD CTX  MIDBRAIN
  3   3  AUD CTX  MIDBRAIN
  3   4  AUD CTX  MIDBRAIN
  3   5  AUD CTX  MIDBRAIN
  4   1  AUD CTX  MIDBRAIN
  4   2  AUD CTX  MIDBRAIN
  4   3  AUD CTX  MIDBRAIN
  4   4  AUD CTX  MIDBRAIN
  4   5  AUD CTX  MIDBRAIN
  5   1  AUD CTX  MIDBRAIN
  5   2  AUD CTX  MIDBRAIN
  5   3  AUD CTX  MIDBRAIN
  5   4  AUD CTX  MIDBRAIN
  5   5  AUD CTX  MIDBRAIN

gfp.oc = gfp %>%
  filter(EXPT=='OLD') %>%
  filter(VALENCE=='CONTEXT') %>%
  get_corrtest_sys()

 iter imp variable
  1   1  MIDBRAIN
  1   2  MIDBRAIN
  1   3  MIDBRAIN
  1   4  MIDBRAIN
  1   5  MIDBRAIN
  2   1  MIDBRAIN
  2   2  MIDBRAIN
  2   3  MIDBRAIN
  2   4  MIDBRAIN
  2   5  MIDBRAIN
  3   1  MIDBRAIN
  3   2  MIDBRAIN
  3   3  MIDBRAIN
  3   4  MIDBRAIN
  3   5  MIDBRAIN
  4   1  MIDBRAIN
  4   2  MIDBRAIN
  4   3  MIDBRAIN
  4   4  MIDBRAIN
  4   5  MIDBRAIN
  5   1  MIDBRAIN
  5   2  MIDBRAIN
  5   3  MIDBRAIN
  5   4  MIDBRAIN
  5   5  MIDBRAIN

gfp Young CTX vs SHK

gfp Old CTX vs SHK

NormOL

# SUBSET FOR ARC ONLY
nol = sys.weighted %>% filter(CHANNEL=='NormOL')
get_pca(nol, "NormOL")

Importance of components:
                          PC1   PC2    PC3    PC4     PC5     PC6    PC7     PC8     PC9    PC10    PC11    PC12   PC13
Standard deviation     2.4589 1.534 1.3703 1.0356 0.89848 0.82603 0.7046 0.61775 0.53374 0.51935 0.48467 0.42834 0.3854
Proportion of Variance 0.4031 0.157 0.1252 0.0715 0.05382 0.04549 0.0331 0.02544 0.01899 0.01798 0.01566 0.01223 0.0099
Cumulative Proportion  0.4031 0.560 0.6852 0.7567 0.81054 0.85603 0.8891 0.91457 0.93356 0.95154 0.96720 0.97943 0.9893
                          PC14    PC15
Standard deviation     0.33139 0.22399
Proportion of Variance 0.00732 0.00334
Cumulative Proportion  0.99666 1.00000
NULL

                      [,1]
TEMP ASSOC CTX -0.35223944
PFC            -0.34622668
MOTOR CTX      -0.34204612
AUD CTX        -0.29739292
OLFACTORY      -0.29499301
POST ASSOC CTX -0.29149819
STRIATUM       -0.28938061
SOMAT CTX      -0.28931878
VISC CTX       -0.28403330
SUBIC-HPC      -0.23968086
CORTICO-AMY    -0.22219523
THALAMIC       -0.10285461
MIDBRAIN        0.08365725
NUCLEAR AMY    -0.07368919
HYPOTHAL        0.02602978

Control

nodd = nol %>%
  filter((ID %%2) == T) %>%
  get_corrtest_sys()

 iter imp variable
  1   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  1   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  1   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  1   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  1   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  2   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  2   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  2   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  2   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  2   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  3   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  3   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  3   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  3   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  3   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  4   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  4   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  4   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  4   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  4   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  5   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  5   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  5   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  5   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC
  5   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SUBIC-HPC  THALAMIC

neven = nol %>%
  filter((ID%%2) ==F) %>%
  get_corrtest_sys()

 iter imp variable
  1   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  1   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  1   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  1   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  1   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  2   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  2   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  2   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  2   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  2   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  3   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  3   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  3   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  3   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  3   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  4   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  4   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  4   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  4   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  4   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  5   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  5   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  5   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  5   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX
  5   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  VISC CTX

compare_matrix_r(nodd, neven, "odd", "even")

[1] "n1=" "24" 
[1] "n2=" "27" 
Tests of correlation matrices 
Call:cortest.normal(R1 = x$r, R2 = y$r, n1 = min(x$n), n2 = min(y$n))
 Chi Square value 101.99  with df =  105   with probability < 0.57 

NormOL::VALENCE stats

nv1 = nol %>%
  filter(VALENCE == 'CONTEXT') %>%
  get_corrtest_sys() 

 iter imp variable
  1   1  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  1   2  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  1   3  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  1   4  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  1   5  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  2   1  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  2   2  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  2   3  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  2   4  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  2   5  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  3   1  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  3   2  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  3   3  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  3   4  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  3   5  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  4   1  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  4   2  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  4   3  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  4   4  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  4   5  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  5   1  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  5   2  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  5   3  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  5   4  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC
  5   5  AUD CTX  MIDBRAIN  POST ASSOC CTX  SUBIC-HPC

nv2 = nol %>%
  filter(VALENCE == 'SHOCK') %>%
  get_corrtest_sys() 

 iter imp variable
  1   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX

compare_matrix_r(nv1, nv2, "context", "shock")

[1] "n1=" "26" 
[1] "n2=" "25" 
Tests of correlation matrices 
Call:cortest.normal(R1 = x$r, R2 = y$r, n1 = min(x$n), n2 = min(y$n))
 Chi Square value 165.65  with df =  105   with probability < 0.00015 

NormOL::EXPT

# Get Z dist of CONTEXT
nol.age1 = nol %>%
  filter(EXPT == 'YOUNG') %>%
  get_corrtest_sys()

 iter imp variable
  1   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  1   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  2   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  3   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  4   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   1  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   2  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   3  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   4  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX
  5   5  AUD CTX  CORTICO-AMY  HYPOTHAL  MIDBRAIN  NUCLEAR AMY  POST ASSOC CTX  SOMAT CTX  STRIATUM  SUBIC-HPC  THALAMIC  VISC CTX

nol.age2 = nol %>%
  filter(EXPT == 'OLD') %>%
  get_corrtest_sys()

 iter imp variable
  1   1  AUD CTX  MIDBRAIN
  1   2  AUD CTX  MIDBRAIN
  1   3  AUD CTX  MIDBRAIN
  1   4  AUD CTX  MIDBRAIN
  1   5  AUD CTX  MIDBRAIN
  2   1  AUD CTX  MIDBRAIN
  2   2  AUD CTX  MIDBRAIN
  2   3  AUD CTX  MIDBRAIN
  2   4  AUD CTX  MIDBRAIN
  2   5  AUD CTX  MIDBRAIN
  3   1  AUD CTX  MIDBRAIN
  3   2  AUD CTX  MIDBRAIN
  3   3  AUD CTX  MIDBRAIN
  3   4  AUD CTX  MIDBRAIN
  3   5  AUD CTX  MIDBRAIN
  4   1  AUD CTX  MIDBRAIN
  4   2  AUD CTX  MIDBRAIN
  4   3  AUD CTX  MIDBRAIN
  4   4  AUD CTX  MIDBRAIN
  4   5  AUD CTX  MIDBRAIN
  5   1  AUD CTX  MIDBRAIN
  5   2  AUD CTX  MIDBRAIN
  5   3  AUD CTX  MIDBRAIN
  5   4  AUD CTX  MIDBRAIN
  5   5  AUD CTX  MIDBRAIN

compare_matrix_r(nol.age1, nol.age2, "young", "old")

[1] "n1=" "24" 
[1] "n2=" "27" 
Tests of correlation matrices 
Call:cortest.normal(R1 = x$r, R2 = y$r, n1 = min(x$n), n2 = min(y$n))
 Chi Square value 160.8  with df =  105   with probability < 0.00038 

NormOL - EXPT:VALENCE

NormOL Young CTX vs SHK

NormOL Old CTX vs SHK

GFAP

CONTROL

GFAP - VALENCE

GFAP - EXPT