201130_Math217_FP

Intro Stuff

library(tidyverse)

## Warning: package 'tidyverse' was built under R version 4.0.2

## -- Attaching packages --------------------------------------- tidyverse 1.3.0 --

## v ggplot2 3.3.2     v purrr   0.3.4
## v tibble  3.0.4     v dplyr   1.0.2
## v tidyr   1.1.2     v stringr 1.4.0
## v readr   1.4.0     v forcats 0.5.0

## Warning: package 'tibble' was built under R version 4.0.3

## Warning: package 'tidyr' was built under R version 4.0.2

## Warning: package 'dplyr' was built under R version 4.0.2

## -- Conflicts ------------------------------------------ tidyverse_conflicts() --
## x dplyr::filter() masks stats::filter()
## x dplyr::lag()    masks stats::lag()

library(infer)

## Warning: package 'infer' was built under R version 4.0.2

library(pastecs)

## Warning: package 'pastecs' was built under R version 4.0.3

## 
## Attaching package: 'pastecs'

## The following objects are masked from 'package:dplyr':
## 
##     first, last

## The following object is masked from 'package:tidyr':
## 
##     extract

library(mosaic)

## Warning: package 'mosaic' was built under R version 4.0.2

## Registered S3 method overwritten by 'mosaic':
##   method                           from   
##   fortify.SpatialPolygonsDataFrame ggplot2

## 
## The 'mosaic' package masks several functions from core packages in order to add 
## additional features.  The original behavior of these functions should not be affected by this.

## 
## Attaching package: 'mosaic'

## The following object is masked from 'package:Matrix':
## 
##     mean

## The following objects are masked from 'package:infer':
## 
##     prop_test, t_test

## The following objects are masked from 'package:dplyr':
## 
##     count, do, tally

## The following object is masked from 'package:purrr':
## 
##     cross

## The following object is masked from 'package:ggplot2':
## 
##     stat

## The following objects are masked from 'package:stats':
## 
##     binom.test, cor, cor.test, cov, fivenum, IQR, median, prop.test,
##     quantile, sd, t.test, var

## The following objects are masked from 'package:base':
## 
##     max, mean, min, prod, range, sample, sum

library(Hmisc)

## Warning: package 'Hmisc' was built under R version 4.0.3

## Loading required package: survival

## Loading required package: Formula

## Warning: package 'Formula' was built under R version 4.0.3

## 
## Attaching package: 'Hmisc'

## The following objects are masked from 'package:dplyr':
## 
##     src, summarize

## The following objects are masked from 'package:base':
## 
##     format.pval, units

library(gmodels)

## Warning: package 'gmodels' was built under R version 4.0.3

getwd()

## [1] "C:/Users/Jerome/Documents/From_Toshiba_HD_Work_Files/0000_Montgomery_College/Math_217/Final_Project/Working_Folder"

family_z34 <- read.csv ("family_z34.csv")

Examine Variables

table(family_z34$FM_SIZE)

## 
##  1  2  3  4  5  6  7  8  9 10 
## 70 88 61 56 30 17  5  1  1  1

table(family_z34$FM_TYPE)

## 
##   1   2   3   4 
##  70 112  36 112

table(family_z34$FM_STRP)

## 
## 11 12 21 22 23 31 32 41 42 43 44 45 
## 63  7 41 18 53 32  4 30 10 16 48  8

table(family_z34$FM_KIDS)

## 
##   0   1   2   3   4   5   6 
## 182  49  57  26  11   4   1

table(family_z34$FM_ELDR)

## 
##   0   1   2 
## 305  20   5

table(family_z34$FM_EDUC1)

## 
##  1  2  3  4  5  6  7  8  9 
##  1 15  5 66 98 48 16 62 19

table(family_z34$F10DVYN)

## 
##   1   2 
## 138 192

table(family_z34$FDMEDYN)

## 
##   1   2 
## 190 140

table(family_z34$FNMEDYN)

## 
##   1   2 
## 172 158

table(family_z34$FLIADLCT)

## 
##   0   1   2 
## 279  49   2

table(family_z34$FWKLIMCT)

## 
##   0   1   2   3 
## 189 117  19   5

table(family_z34$FHSTATEX)

## 
##   0   1   2   3   4   5   6 
## 221  47  30  11  12   5   4

family_z35 <- mutate(family_z34, health_ins = ifelse(FHIPRVCT > 0 | FHISINCT > 0  |FHICARCT > 0 | FHICADCT >0 |FHICHPCT >0 | FHIMILCT >0  | FHIPUBCT > 0 | FHIOGVCT > 0 | FHIIHSCT> 0 | FHIEXCT > 0, 1,0))

table(family_z35$health_ins)

## 
##   1 
## 330

table(family_z34$FHIPRVCT)

## 
##   0   1   2   3   4   5   6   7 
##   2 165  81  39  19  15   6   3

table(family_z34$FHISINCT)

## 
##   0   1   2   3   4   5   6   7 
## 139  89  45  23  13  14   5   2

table(family_z34$FHICARCT)

## 
##   0   1   2 
## 277  45   8

table(family_z34$FHICADCT)

## 
##   0   1   2   3   4   5   6 
## 229  36  30  18  10   5   2

table(family_z34$FHICHPCT)

## 
##   0   1   2   3   4 
## 321   3   3   2   1

table(family_z34$FHIMILCT)

## 
##   0   1   2   3 
## 319   9   1   1

table(family_z34$FHIPUBCT)

## 
##   0   1   2 
## 326   2   2

table(family_z34$FHIOGVCT)

## 
##   0 
## 330

table(family_z34$FHIIHSCT)

## 
##   0   1 
## 329   1

table(family_z34$FHIEXCT)

## 
##   0   1   2 
## 322   7   1

table(family_z35$INCGRP5)

## 
##   1   2   3   4 
## 135 160  15  20

table(family_z35$RAT_CAT5)

## 
##  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16 17 18 
## 14 12 11 26 20 33 22 49 48 26 25 10  6 11  8  6  1  2

table(family_z34$COVCONF)

## 
##   1   2   3   4 
##   6  17  44 263

table(family_z34$FHICOST)

## 
##   0   1   2   3   4   5 
##  13  56 133  47  38  43

table(family_z34$FMEDBNOP)

## 
##   1   2 
## 227 103

table(family_z34$FSAF)

## 
##   1   2   9 
##  49 278   3

table(family_z34$FHICOVCT)

## 
##  1  2  3  4  5  6  7  8 10 
## 86 90 61 48 25 13  4  2  1

table(family_z34$FHICOVYN)

## 
##   1 
## 330

table(family_z34$FHIEBCCT)

## 
##   1   2   3   4   5   6   7 
## 175  76  36  17  16   5   3

table(family_z34$FMEDBILL)

## 
##   1 
## 330

write.csv(family_z35, file = "family_z35.csv", row.names = FALSE)
family_z35 <- read.csv("family_z35.csv")

family_z31 <- read.csv("family_z31.csv")
family_z50 <- mutate(family_z31, food_sec = ifelse(FSRUNOUT == 1 & FSLAST == 1 & FSSKIP == 1 & FSBALANC == 1 & FSLESS == 1 & FSHUNGRY == 1 & FSWEIGHT == 1 & FSNOTEAT == 1, 1,0))

table(family_z50$food_sec)

## 
##   0   1 
## 321   9

family_z51 <- mutate(family_z31, health_ins = ifelse(FHIPRVCT > 0 & FHISINCT > 0  & FHICARCT > 0 & FHICADCT >0 |FHICHPCT >0 & FHIMILCT >0  & FHIPUBCT > 0 & FHIOGVCT > 0 & FHIIHSCT> 0 & FHIEXCT > 0, 1,0))

table(family_z51$health_ins)

## 
##   0   1 
## 318  12

family_z17 <- read.csv("family_z17.csv")
family_z52 <- mutate(family_z17, health_ins = ifelse(FHIPRVCT > 0 & FHISINCT > 0  & FHICARCT > 0 & FHICADCT >0 & FHICHPCT >0 & FHIMILCT >0  & FHIPUBCT > 0 & FHIOGVCT > 0 & FHIIHSCT> 0 & FHIEXCT > 0, 1,0))

table(family_z52$health_ins)

## 
##    0 
## 2357

family_z52 <- mutate(family_z17, health_ins = ifelse(FHIPRVCT > 0 | FHISINCT > 0  | FHICARCT > 0 | FHICADCT >0 | FHICHPCT >0 | FHIMILCT >0  | FHIPUBCT > 0 | FHIOGVCT > 0 | FHIIHSCT> 0 | FHIEXCT > 0, 1,0))
table(family_z52$health_ins)

## 
##    0    1 
##  230 2127

family_0 <- read.csv("family_0.csv")
private_ins <- filter(family_0, FHIPRVCT > 0)

family_z53 <- mutate(family_0, health_ins = ifelse(FHIPRVCT > 0 | FHISINCT > 0  | FHICARCT > 0 | FHICADCT >0 | FHICHPCT >0 | FHIMILCT >0  | FHIPUBCT > 0 | FHIOGVCT > 0 | FHIIHSCT> 0 | FHIEXCT > 0, 1,0))
table(family_z53$health_ins)

## 
##     0     1 
##  1591 31566

family_z54 <- mutate(family_z53, food_sec = ifelse(FSRUNOUT == 1 | FSLAST == 1 | FSSKIP == 1 | FSBALANC == 1 | FSLESS == 1 | FSHUNGRY == 1 | FSWEIGHT == 1 | FSNOTEAT == 1, 1,0))
table(family_z54$food_sec)

## 
##    1 
## 2967

table(family_z54$health_ins)

## 
##     0     1 
##  1591 31566

write.csv(family_z54, file = "family_z54.csv", row.names = FALSE)

family_z54 <- read.csv("family_z54.csv")

table(family_z54$food_sec, family_z54$health_ins)

##    
##        0    1
##   1  273 2694

201130_Math217_FP_Analysis

Jerome

11/30/2020

Intro Stuff

Examine Variables