# Les dades tenen una distribució parabòlica
# Codi original: Vincent Zoonekynd (http://zoonek2.free.fr/UNIX/48_R/09.html#3)
N = 100
X = runif(N,-1,1)
Y = rnorm(N, X * X, .1)
DADES_RLS_PAR = data.frame(Y, X)                       # Es crea una taula de dades
save(DADES_RLS_PAR, file="dades/dades_rls_par.RData")  # Se desa la taula de dades anàlisi posteriors
load("dades/dades_rls_par.RData")                      # Es carregen les dades (no cal)
svg("gra/reg_ls_par.svg") 
plot(Y ~ X, data = DADES_RLS_PAR, main = "Relació parabòlica", xlab="X", ylab="Y")
abline(lm(Y ~ X, data = DADES_RLS_PAR), col="red")
mtext("La recta de regressió (en vermell) és de la Y sobre la X") # Subtítol
dev.off()