example-script.R

# This needs to already have been done
# install.packages("psych")
library(psych)

SCRATCH <- Sys.getenv("SCRATCH")
setwd(SCRATCH)

# This command assumes the file exists at /scratch/$USER/emotionzs.csv
alldata = read.csv("emotionzs.csv", header = FALSE)

ev <- eigen(cor(alldata)) # get eigenvalues
ap <- parallel(subject=nrow(alldata),var=ncol(alldata),

rep=100,cent=.05)
nS <- nScree(x=ev$values, aparallel=ap$eigen$qevpea)

# This will save to /scratch/$USER
save(nS)