data(HorseBeginners)
# There is one outlier
with(HorseBeginners,plot(paired(Actual,Imaginary),type="profile"))
# This outlier has a great influence
# on the non robust Pitman-Morgan test of variances
with(HorseBeginners,var.test(paired(Actual,Imaginary)))
with(HorseBeginners[-1,],var.test(paired(Actual,Imaginary)))
with(HorseBeginners,grambsch.var.test(paired(Actual,Imaginary)))
with(HorseBeginners,bonettseier.var.test(paired(Actual,Imaginary)))
Run the code above in your browser using DataLab