Capstone Project

lead <- read.csv("C:/Users/YUCHIAWU/Desktop/lead.csv", header = TRUE)
boxdata <- read.csv("C:/Users/YUCHIAWU/Desktop/lead data for side by side boxplot.csv", header = TRUE)

northmean1<- mean(lead$North.1, na.rm = TRUE)
northsd1 <- sd(lead$North.1, na.rm = TRUE)
northmean1

## [1] 6.089333

northsd1

## [1] 7.319642

northmean2 <- mean(lead$North.2, na.rm = TRUE)
northsd2 <- sd(lead$North.2, na.rm = TRUE)
northmean2

## [1] 3.188396

northsd2

## [1] 3.549662

northmean3 <- mean(lead$North.3, na.rm = TRUE)
northsd3 <- sd(lead$North.3, na.rm = TRUE)
northmean3

## [1] 2.106542

northsd3

## [1] 2.957509

NEmean1 <- mean(lead$Northeast.1, na.rm = TRUE)
NEsd1 <- sd(lead$Northeast.1, na.rm = TRUE)
NEmean1

## [1] 12.17284

NEsd1

## [1] 27.82881

NEmean2 <- mean(lead$Northeast.2, na.rm = TRUE)
NEsd2 <- sd(lead$Northeast.2, na.rm = TRUE)
NEmean2

## [1] 3.024227

NEsd2

## [1] 6.999262

NEmean3 <- mean(lead$Northeast.3, na.rm = TRUE)
NEsd3 <- sd(lead$Northeast.3, na.rm = TRUE)
NEmean3

## [1] 2.296932

NEsd3

## [1] 6.304727

Southmean1 <- mean(lead$South.1, na.rm = TRUE)
southsd1 <- sd(lead$South.1, na.rm = TRUE)
Southmean1

## [1] 11.10155

southsd1

## [1] 23.28593

Southmean2 <- mean(lead$South.2, na.rm = TRUE)
southsd2 <- sd(lead$South.2, na.rm = TRUE)
Southmean2

## [1] 6.204941

southsd2

## [1] 13.35963

Southmean3 <- mean(lead$South.3, na.rm = TRUE)
southsd3 <- sd(lead$South.3, na.rm = TRUE)
Southmean3

## [1] 4.784412

southsd3

## [1] 14.38969

centermean1 <- mean(lead$Center.1, na.rm = TRUE)
centersd1 <- sd(lead$Center.1, na.rm = TRUE)
centermean1

## [1] 10.95643

centersd1

## [1] 17.65755

centermean2 <- mean(lead$Center.2, na.rm = TRUE)
centersd2 <- sd(lead$Center.2, na.rm = TRUE)
centermean2

## [1] 5.65913

centersd2

## [1] 9.016728

centermean3 <- mean(lead$Center.3, na.rm = TRUE)
centersd3 <- sd(lead$Center.3, na.rm = TRUE)
centermean3

## [1] 3.766203

centersd3

## [1] 7.795677

NWmean1 <- mean(lead$Northwest.1, na.rm = TRUE)
NWsd1 <- sd(lead$Northwest.1, na.rm = TRUE)
NWmean1

## [1] 13.17613

NWsd1

## [1] 27.09034

NWmean2 <- mean(lead$Northwest.2, na.rm = TRUE)
NWsd2 <- sd(lead$Northwest.2, na.rm = TRUE)
NWmean2

## [1] 32.60027

NWsd2

## [1] 147.7863

NWmean3 <- mean(lead$Northwest.3, na.rm = TRUE)
NWsd3 <- sd(lead$Northwest.3, na.rm = TRUE)
NWmean3

## [1] 5.131473

NWsd3

## [1] 15.49357

boxplot(boxdata$first.draw ~ boxdata$Cardianl.direction)

boxplot(boxdata$second.draw ~ boxdata$Cardianl.direction)

boxplot(boxdata$third.draw ~ boxdata$Cardianl.direction)

hist(lead$Northwest.3)