plot(c(0,1),c(1,0),
type="n",xlab="x",ylab="y")
xy<-locator(type="p") # marcamos los puntos en el gráfico y sus coordenadas son guardadas en xy.
xy
abline(lm(y~x,xy),col=4)#recta de cuadr. mín.
 library(MASS)
 abline(ltsreg(xy$x,xy$y),lty=2,col=2)
 

altedad <- read.table (file.choose())
 names(altedad)<-c("edad","altura")
 plot(altedad$edad,altedad$altura,
xlab="Edad (meses)",ylab="Altura (cm) ")
 identify(altedad$edad,altedad$altura)


plot(altedad$edad,altedad$altura)
 abline(lm(altedad$altura~altedad$edad))
lm(altedad$altura~altedad$edad)

aa <-lm(altedad$altura~altedad$edad)

plot (aa$fitted,aa$residuals)
abline(h=0)
identify(aa$fitted,aa$residuals)
#lines(lowess(aa$fitted,aa$residuals),col="red")
points(lowess(aa$fitted,aa$residuals),pch=19,col="blue")

resCM <-lm(altedad$altura~altedad$edad)$residuals
sort(resCM)

qqnorm(aa$residuals)
qqline(aa$residuals)
identify(qqnorm(aa$residuals))

library (MASS)
plot(altedad$edad,altedad$altura, xlab="edad", ylab="altura") 
abline(lm(altedad$altura~altedad$edad))
abline(ltsreg(altedad$edad,altedad$altura),col="red",lty= 2)
abline( ltsreg(altedad$edad,altedad$altura, quantile=12),col="blue", lty=4 )
legend(110,165,c("lm","ltsreg","ltsreg, quantile=12"),lty = c(1,2,4), col = c(1,2,4))

#legend(-1,10,c("lm","ltsreg","ltsreg, quantile=12"),lty = c(1,2,4), col = c(1,2,4))