Hi all,

I´m re-starting (as my name indicates) my little knowlegde of CRAN R.

I made this function time before but I don´t know where is the error
because nothing appears as an error but the histogram plot doesn´t appear.
Should I install some special library to run sapply?


pru<-function(){

randz<-matrix(rnorm(200000),100,2000)

H<-matrix(0,100,2000)


for (j in 2:2000){
for (i in 2:100){
H[1,]<-randz[1,]
H[i,j]<-H[i-1,j]+randz[i,j]


}}

hy<-nrow(H)-1

estima<-H[2:nrow(H),]

estima2<-H[1:hy,]


a<-estima
b<-estima2


mycoef <- function(x, y) coefficients( lm(y ~ x-1) )

rest <- sapply(2:2000, function(i){
               y <- a[,i]
               x <- b[,i]
              mycoef(x,y)
               }
            )



print(summary(rest))

hist(rest,col="blue",breaks=seq(0.6,1.05,0.01),prob=TRUE)
lines(density(rest,bw=0.03))
rug(rest)



}

Many Thanks¡¡¡

	[[alternative HTML version deleted]]

