library(ggplot2)
library(qplot)
lsfit.a<-array(0,dim=c(100,2))
for(i in 1:100)
{
a1<-20
a2<-0.6
n<-25
x<-c(1:n)
u<-runif(n,min=0,max=1)
y<-a1+a2*x+u
lsfit.a[,i]<-c(lsfit(x,y)$coefficients)
}
lsfit.amean<-apply(lsfit.a,2,mean)
lsfit.amean
library(qplot)
lsfit.a<-array(0,dim=c(100,2))
for(i in 1:100)
{
a1<-20
a2<-0.6
n<-25
x<-c(1:n)
u<-runif(n,min=0,max=1)
y<-a1+a2*x+u
lsfit.a[,i]<-c(lsfit(x,y)$coefficients)
}
lsfit.amean<-apply(lsfit.a,2,mean)
lsfit.amean
