mh2018

mh18 <- load("/Users/caelynsobie/Downloads/nhmss_puf_2018_r.Rdata")

mh18 <- subset(nm18, select=c(115:171))
#change all values of non-responses to NA
mh18[mh18 == -6 |  mh18 == -1 | mh18 == -3 | mh18 == -2 | mh18 == 2 | mh18 == 3] <- NA

# Sum up each column, ignoring NA values
column_sums <- colSums(mh18, na.rm = TRUE)
column_sums

           UTREV           SATSUR    SMOKINGPOLICY    USEDSECLUSION   ADOPTSECLUSION            INTKE 
           10449            11082             5934             2482             9204             4561 
        SCHEDULE           ASSESS           TXPLAN         PROGRESS           DSCHRG              REF 
            6910             6392             6563             7139             6446             2952 
             LAB             DISP           MEDINT         STOREREC         SENDINFO          RECINFO 
            3051             3313             5358             4269             1219              891 
            BILL        SATSURVEY         FEESCALE          PAYASST          REVCHK1          REVCHK2 
            6214              930             6475             5854             9766             9427 
         REVCHK8          REVCHK5         REVCHK10         FUNDSMHA FUNDSTATEWELFARE     FUNDSTATEJUV 
            8073            10329             6867             6683             4867             3503 
   FUNDSTATEEDUC     FUNDOTHSTATE     FUNDLOCALGOV         FUNDCSBG         FUNDCMHG         REVCHK15 
            1996             4152             5533             2574             3658             5803 
          FUNDVA         REVCHK17         REVCHK2A          LICENMH         LICENSED          LICENPH 
            2692              944               87             8399             3944             5733 
     LICENSEDFCS         LICENHOS            JCAHO             CARF              COA              CMS 
            2457             1803             4108             2877             1177             5865 
        OTHSTATE           OTHFAC           FACNUM           IPSERV          IPTOTAL        IPSEXTOTM 
             421            10036              343             1863             5380             3286 
       IPSEXPERM        IPSEXTOTF        IPSEXPERF 
            2740             2753             2222 

# Example: Summing the first 10 columns (adjust as needed for your variables)
v1 <- rowSums(mh18[, 1:10], na.rm = TRUE)
hist(v1)

v2 <- rowSums(mh18[, 11:20], na.rm = TRUE)
hist(v2)

v3 <- rowSums(mh18[, 21:30], na.rm = TRUE)
hist(v3)

v4 <- rowSums(mh18[, 31:40], na.rm = TRUE)
hist(v4)

v5 <- rowSums(mh18[, 41:50], na.rm = TRUE)
hist(v5)

v5_1 <- rowSums(mh18[, c(41,42  :50)], na.rm = TRUE)
hist(v5_1)

v6 <- rowSums(mh18[, 51:57], na.rm = TRUE)
hist(v6)