setwd("D:\\Documents\\uoa\\EngSci\\Teaching\\ENGSCI355\\2017\\Labs\\Pages\\ExtendedHealthClinic") #{ [Simulation].RunIndex(1) } { [WaitingTriage].SampleAverage } { [WaitingTest].SampleAverage } { [WaitingTreat].SampleAverage } { [WaitingTime].SampleAverage } { [TimeInSystem].SampleAverage } { [WalkupQueue].QueueLengthAverage } { [Doctor].Utilisation } { [WaitingTime2].SampleAverage } { [TimeInSystem2].SampleAverage } { [AppointmentQueue].QueueLengthAverage } { '[Doctor2Appointment].Utilisation + [Doctor2Walkup].Utilisation' } #Lab3data = read.table("Lab3_Statistics.dat", sep="\t", #Lab3data = read.table("Lab3_AppDelay.dat", sep="\t", Lab3data = read.table("Lab3_Final.dat", sep="\t", col.names=c("RunNumber", "WalkUp.WaitingTriage", "WalkUp.WaitingTest", "WalkUp.WaitingTreat", "WalkUp.WaitingTime", "WalkUp.TimeInSystem", "Num.WalkUps.Waiting", "Doctor.Utilisation", "Appointment.WaitingTime", "Appointment.TimeInSystem", "Appointment.Queue", "Doctor2.Utilisation" ), skip=2) Lab3scen = Lab3data[Lab3data$RunNumber == 4,2:dim(Lab3data)[2]] results = apply(Lab3scen, 2, FUN = function(x) { tTest = t.test(x) return(tTest$conf.int) }) temp = results[2,] results[2,] = apply(Lab3scen, 2, mean) results = rbind(results, temp) rownames(results) = c("ConfLower", "Mean", "ConfUpper") results t.test(Lab3scen$WalkUp.WaitingTime, Lab3scen$Appointment.WaitingTime) t.test(Lab3scen$WalkUp.TimeInSystem, Lab3scen$Appointment.TimeInSystem) t.test(Lab3scen$Num.WalkUps.Waiting, Lab3scen$Appointment.Queue) t.test(Lab3scen$Doctor.Utilisation, Lab3scen$Doctor2.Utilisation)