The Brown-Forsyth test for homogeneity of variance is included in
the HH package, downloadable from CRAN.
library(HH)
x <- c(rnorm(1000), rnorm(1000, 0, 1.2))
tmp <- data.frame(x=x, group=rep(c("s1","s1.2"), c(1000,1000)))
plot.hov(x ~ group, data=tmp)
hov(x ~ group, data=tmp)