library(poliscidata)  # run every time

# indicator variable created for each abortion variable
# this is going to run without console output
gss$abdefect.n    = as.numeric(gss$abdefect=="YES")
gss$abhlth.n 	= as.numeric(gss$abhlth  =="YES")
gss$abnomore.n	= as.numeric(gss$abnomore=="YES")
gss$abpoor.n 	= as.numeric(gss$abpoor  =="YES")
gss$abrape.n 	= as.numeric(gss$abrape  =="YES")
gss$absingle.n	= as.numeric(gss$absingle=="YES")
gss$abany.n 	= as.numeric(gss$abany   =="YES")

# individual indicators added together to create 0-7 scale
# this does not create console output
gss$abortion.scale = gss$abdefect.n + gss$abhlth.n + gss$abnomore.n +  
  gss$abpoor.n   + gss$abrape.n + gss$absingle.n + 
  gss$abany.n

# describe the additive index to see results
cat("Frequency Distribution Table:")
freqC(gss$abortion.scale, gss$wtss)