Is there any heteroskedasticity test in the package? Something that would flag a sample like x <- c(rnorm(1000), rnorm(1000, 0, 1.2)) Alberto Monteiro