Cacher is the code snippet organizer for pro developers

ls /home/laramie/CLEANED_DATA/SAFR_PTSD_wave1 | grep .fam | awk '{print "/home/laramie/CLEANED_DATA/SAFR_PTSD_wave1/"$1}'   >> famlist.txt

Rscript familycounter.Rscript famlist.txt

args <- commandArgs(trailingOnly = TRUE)

family_files <- args[1]

dat <- read.table(family_files)


for (i in dat$V1)
{
  options(stringsAsFactors=F)
  dat <- read.table(i,header=F)
  names(dat) <- c("FIDL","IID","M","F","G","P")
  print(i)
  print(table(dat$P))
}