The Brown-Forsyth test is included in the HH library in function hov(). hov is "homogeneity of variance". Use hov(y ~ x, data=mydata) and plot.hov(y ~ x, data=mydata) Rich