# The three Models
intercept <- c(1.3,1.1,2)
slope <- c(2,2.2,0.9)
age <- c(18:90)

# simulate errors for each model
t=1
z=t*NROW(age)
error1 <- rnorm(100*z,0,1)
error2 <- rnorm(50*z,0,1.4)
error3 <- rnorm(10*z,0,2)

m1 <- intercept[1]+slope[1]*age
m2 <- intercept[2]+slope[2]*age
m3 <- intercept[3]+slope[3]*age

m1 <- data.frame(time=rep(m1,t)+error1,age)
m2 <- data.frame(time=rep(m2,t)+error2,age)
m3 <- data.frame(time=rep(m3,t)+error3,age)

Agereg <- rbind(m1,m2,m3)

new.age <- lm(time~age,data=Agereg)
summary(new.age)

# Plot
plot(m1$age,m1$time,type="p")
points(m2$age,m2$time,col="blue")
points(m3$age,m3$time,col="red")
abline(new.age,col="green")