inst/CaseStudies/tmprv.R

In matloff/revisit: "Statistical Audit": Verify and Revisit Statistical Analysis

# RV history start
# original code
# RV history end
# original code
### data(pima)
datapath <- system.file("data", package="revisit")                              
pima <- read.table(paste0(datapath, "/pima.txt"), header = TRUE)                
# divide into diabetic, non-diabetics
d <- which(pima$Diab == 1)
diab <- pima[d,]
nondiab <- pima[-d,]
# form a confidence interval for each variable, difference between
# diabetics and non-diabetics
for (i in 1:8)  {
   tmp <- revisit::t.test.rv(diab[,i],nondiab[,i],bonf=8)$conf.int
   cat(names(pima)[i],'  ',tmp[1],tmp[2],'\n')
}