Column_Manipulation

library(tidyverse)

## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr     1.1.2     ✔ readr     2.1.4
## ✔ forcats   1.0.0     ✔ stringr   1.5.0
## ✔ ggplot2   3.4.2     ✔ tibble    3.2.1
## ✔ lubridate 1.9.2     ✔ tidyr     1.3.0
## ✔ purrr     1.0.1     
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors

library(dplyr)
library(infer)

msdlabs <- read.csv("msd_labs.csv")

###msdlabs$gender <- msdlabs$FEMALE
msdlabs$male <- ifelse (msdlabs$FEMALE == 0,msdlabs$HOUSEHOLD, NA)
msdlabs$women <- ifelse(msdlabs$FEMALE == 1, msdlabs$HOUSEHOLD, NA)

var.test (msdlabs$male,  msdlabs$women,
          paired = FALSE,
          alternative = "two.sided",
          sigma.squared = 100)

## 
##  F test to compare two variances
## 
## data:  msdlabs$male and msdlabs$women
## F = 0.65992, num df = 709, denom df = 845, p-value = 1.116e-08
## alternative hypothesis: true ratio of variances is not equal to 1
## 95 percent confidence interval:
##  0.5732264 0.7603869
## sample estimates:
## ratio of variances 
##          0.6599153

var.test (msdlabs$HOUSEHOLD ~ msdlabs$FEMALE,
          paired = FALSE,
          alternative = "two.sided",
          sigma.squared = 100)

## 
##  F test to compare two variances
## 
## data:  msdlabs$HOUSEHOLD by msdlabs$FEMALE
## F = 0.65992, num df = 709, denom df = 845, p-value = 1.116e-08
## alternative hypothesis: true ratio of variances is not equal to 1
## 95 percent confidence interval:
##  0.5732264 0.7603869
## sample estimates:
## ratio of variances 
##          0.6599153