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]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
