[R] No error message no display output

2012-05-05 Thread Trying To learn again 
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(20),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.

Re: [R] No error message no display output

2012-05-05 Thread Trying To learn again 
I´m using this version of R

R version 2.13.2 (2011-09-30)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i386-pc-mingw32/i386 (32-bit)


2012/5/5 Trying To learn again tryingtolearnag...@gmail.com

 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(20),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.

Re: [R] No error message no display output

2012-05-05 Thread Uwe Ligges 

1. You need to call the function

pru()

after defining it.


2. Next excercise will be to get rid of all those loops.

3. Next excercide will be to reutn something useful from your function 
(rather than printing).


4. Upgrade to a recent version of R.

Uwe Ligges



On 05.05.2012 17:31, Trying To learn again wrote:

I´m using this version of R

R version 2.13.2 (2011-09-30)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i386-pc-mingw32/i386 (32-bit)


2012/5/5 Trying To learn againtryingtolearnag...@gmail.com


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(20),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.


__
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.