Hi Trying,
Jim already suggested youuse write.table(), which I think is really
what you want. I also wanted to point out that your outer loop is
unnecessary. The following yields identical results and is *much*
faster.
randz <- matrix(rnorm(1000000), 500, 2000)
H <- matrix(0, 500, 2000)
H[1, ] <- randz[1, ]
for (i in 2:500){
if(i < 251) {
H[i, ] <- 0.6 * H[i-1, ] + randz[i, ]
} else {
H[i, ] <- H[i-1, ] + randz[i, ]
}
}
write.table(H, file = "datad.txt")
There may be ways to optimize (or remove) the remaining loop, but at
least this first pass should move things along considerably.
Cheers,
Josh
On Sun, May 6, 2012 at 2:02 PM, Trying To learn again
<tryingtolearnag...@gmail.com> wrote:
> Hi all,
>
> I´m trying to use write function to save the output of a program (my
> constructed "H" matrix)
>
>
> randz<-matrix(rnorm(1000000),500,2000)
>
> H<-matrix(0,500,2000)
>
> H[1,]<-randz[1,]
>
> for (j in 1:2000){
> for (i in 2:500){
> if(i<251)
> H[i,j]<-0.6*H[i-1,j]+randz[i,j]
>
> else H[i,j]<-H[i-1,j]+randz[i,j]
>
> }}
>
> write(H, file = "datad.txt",2000)
>
> If I ommit the 2000 on write function it only puts 5 columns.
>
> The problem is that if I use this it seems it is not saving the same data I
> have simulated....or this seems to me.
>
> You see if I type H[,1] it is not the same that includen on the firs column
> on datad.txt?
>
> I feel very slow witted....
>
> Many thanks in advance for all
>
> [[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.
>
--
Joshua Wiley
Ph.D. Student, Health Psychology
Programmer Analyst II, Statistical Consulting Group
University of California, Los Angeles
https://joshuawiley.com/
______________________________________________
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.
