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]]
try write.table(H, file = "datad,txt") On Sun, May 6, 2012 at 5:02 PM, Trying To learn again <tryingtolearnagain at 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 at 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. >-- Jim Holtman Data Munger Guru What is the problem that you are trying to solve? Tell me what you want to do, not how you want to do it.
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
<tryingtolearnagain at 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 at 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/