# heteroscedasticitat
# Codi original de Vincent Zoonekynd (http://zoonek2.free.fr/UNIX/48_R/11.html#9)
X             = runif(100)
Y             = 1 - 2 * X + .3 * X * rnorm(100)
DADES_RLS_CED = data.frame(Y, X)                       # Es crea taula de dades
save(DADES_RLS_CED, file="dades/dades_rls_ced.RData")     # Se desa la taula de dades anàlisi posteriors
load("dades/dades_rls_ced.RData")                         # Carregar dades (no cal)
svg("gra/reg_ls_ced.svg")
plot(Y ~ X, data = DADES_RLS_CED, main = "Variància no constant (heteroscedasticitat)", xlab="X", ylab="Y")
abline(lm(Y ~ X, data = DADES_RLS_CED), col='red')                # Recta de regressió
mtext("La recta de regressió (en vermell) és de la Y sobre la X") # Subtítol
title(sub="Codi original de Vincent Zoonekynd (http://zoonek2.free.fr/UNIX/48_R/11.html#9)")
dev.off()