scripts/check_surf_clusters.R

In bbuchsbaum/FROIAtlas: FROIAtlas

#! /usr/bin/env Rscript

options(warn=-1)

suppressMessages(library(neuroim))

args <- commandArgs(trailingOnly=TRUE)
surf <- args[1]
hemi <- args[2]

source("~/.voltosurf")

if (hemi == "lh") {
  cmd <- paste("SurfClust -spec",
               LH_STD_SPEC,
               "-surf_A smoothwm -input", surf,  
               "4 -rmm -2 -no_cent -thresh_col 4 -thresh .1 -amm2 10")
  system(cmd)
} else {
  cmd <- paste("SurfClust -spec",
               RH_STD_SPEC,
               "-surf_A smoothwm -input", surf,  
               "4 -rmm -2 -no_cent -thresh_col 4 -thresh .1 -amm2 10")
  system(cmd)
}