setwd("D:\\Documents\\uoa\\EngSci\\Teaching\\ENGSCI355\\2017\\Labs\\Pages\\SHCScheduledAppointments") #Lab2data = read.table("Lab2_Separate.dat", sep="\t", Lab2data = read.table("Lab2_Final.dat", sep="\t", col.names=c("RunNumber", "WalkUp.WaitingTime", "WalkUp.TimeInSystem", "WalkUp.Queue", "Doctor.Utilisation", "Appointment.WaitingTime", "Appointment.TimeInSystem", "Appointment.Queue", "Doctor2.Utilisation" ), skip=2) aggregate(.~RunNumber,data=Lab2data, mean) aggregate(.~RunNumber,data=Lab2data, FUN = function(x) { tTest = try(t.test(x), silent=TRUE) if (is(tTest, "try-error")) { interval = c(mean(x), mean(x)) return(interval) } else { if (is.na(tTest$conf.int[1])) { return(c(mean(x), mean(x))) } else { return(tTest$conf.int) } } }) Lab2scen = Lab2data[Lab2data$RunNumber == 4,] t.test(Lab2scen$WalkUp.WaitingTime, Lab2scen$Appointment.WaitingTime) t.test(Lab2scen$WalkUp.TimeInSystem, Lab2scen$Appointment.TimeInSystem) t.test(Lab2scen$WalkUp.Queue, Lab2scen$Appointment.Queue) t.test(Lab2scen$Doctor.Utilisation, Lab2scen$Doctor2.Utilisation)