ASSES/ASS3/qu14-16.R

In agracy2246/MATH4753grac0009: Adam Gracy MATH4753

# Q 14

library(readxl)
dird = ""

myread = function(xls){
  require(readxl)
read_xls(paste0(dird, xls))
}

leadcop = read.csv("LEADCOPP.csv")
#solrad = myread("SOLARAD.xls")
#diazinon = myread("DIAZINON.xls")

t.test(leadcop$LEAD,conf.level = 0.99)

# Or use a fomula
x = leadcop$LEAD
alpha  =0.01
n = length(x)

mean(x) + mp*qt(1-alpha/2,n-1 )*sd(x)/sqrt(n)

#names(solrad)
#names(diazinon)

#with(solrad, t.test(STJOS, IOWA, paired=TRUE, conf.level = 0.80))

#with( diazinon, t.test(DAY, NIGHT, paired =TRUE))