[R] reading in a tricky computer program output

[Berwin A Turlach]

G'day Taka,

>>>>> "TM" == Taka Matzmoto <sell_mirage_ne at hotmail.com> writes:

    TM> and then assign the character vector to the numeric vector by

    TM> names<-first.10
    TM> first.10 = numeric.vector
    TM> combined.one <- cbind(names,first.10)
    TM> container <- diag(10)
    TM> for (i in 1:(10*10))
I don't really understand this loop.  If I reverse-engineer this code
correctly thenthe matrix `combined.one' is not a 2*100 matrix, so you
should get an error while exectuting this loop.

    TM> Is there any other neat way to do this?
Neat way to create those character vectors?  Or a neat way to read in
the data from a file?

If the latter, I would use the following code:

      matzmoto <- function(file, diag=TRUE){
      
        dat <- scan(file)
        if(diag){
          nvar <- sqrt(2*length(dat)+0.25) - 0.5
          nn <- nvar
        }else{
          nvar <- sqrt(2*length(dat)+0.25) + 0.5
          nn <- nvar - 1
        }
        res <- matrix(0,nvar,nvar)
        ind <- upper.tri(res, diag=diag)
      
        rind <- 1:5
        while(nn > 0){
          if( nn < 5 ){
            rind <- rind[1:nn]
            tmp <- matrix(0,nn,nvar)
          }else{
            tmp <- matrix(0,5,nvar)
          }
      
          how.many <- sum(ind[rind,])
          tmp[ind[rind,]] <- dat[1:how.many]
          res[rind,] <- tmp
      
          dat <- dat[-(1:how.many)]
      
          rind <- rind + 5
          nn <- nn - 5
        }
        t(res)
      }
      
      res <- matzmoto("matzmoto1.dat", TRUE)
      print(res)
      
      res <- matzmoto("matzmoto2.dat", FALSE)
      print(res)


After storing the two examples that you posted into the files
matzmoto1.dat and matzmoto2.dat, respectively, and removing the part
that you said you have added, I get the following result on my machine
when sourcing the above code:

      > source("matzmoto.R")
      Read 55 items
              [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]  [,8]  [,9] [,10]
       [1,]  0.000  0.000  0.000  0.000  0.000  0.000  0.000 0.000  0.00     0
       [2,]  0.001  0.000  0.000  0.000  0.000  0.000  0.000 0.000  0.00     0
       [3,] -0.002  0.019  0.000  0.000  0.000  0.000  0.000 0.000  0.00     0
       [4,]  0.012 -0.004 -0.020  0.000  0.000  0.000  0.000 0.000  0.00     0
       [5,] -0.015  0.003  0.011  0.008  0.000  0.000  0.000 0.000  0.00     0
       [6,]  0.005 -0.008 -0.005  0.002  0.005  0.000  0.000 0.000  0.00     0
       [7,]  0.008 -0.007  0.013  0.003  0.007 -0.037  0.000 0.000  0.00     0
       [8,] -0.014 -0.011 -0.010 -0.025  0.002  0.010  0.027 0.000  0.00     0
       [9,]  0.006  0.003 -0.010  0.002 -0.020  0.032 -0.004 0.008  0.00     0
      [10,]  0.006  0.010 -0.006  0.005  0.008 -0.008 -0.011 0.015 -0.02     0
      Read 45 items
              [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8]  [,9] [,10]
       [1,]  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 0.000     0
       [2,] -0.002  0.000  0.000  0.000  0.000  0.000  0.000  0.000 0.000     0
       [3,]  0.003 -0.053  0.000  0.000  0.000  0.000  0.000  0.000 0.000     0
       [4,] -0.026  0.010  0.045  0.000  0.000  0.000  0.000  0.000 0.000     0
       [5,]  0.023 -0.008 -0.025 -0.016  0.000  0.000  0.000  0.000 0.000     0
       [6,] -0.012  0.023  0.013 -0.005 -0.011  0.000  0.000  0.000 0.000     0
       [7,] -0.031  0.031 -0.054 -0.013 -0.027  0.127  0.000  0.000 0.000     0
       [8,]  0.040  0.042  0.031  0.075 -0.007 -0.035 -0.166  0.000 0.000     0
       [9,] -0.012 -0.009  0.023 -0.005  0.037 -0.083  0.015 -0.027 0.000     0
      [10,] -0.013 -0.027  0.014 -0.013 -0.020  0.021  0.047 -0.052 0.048     0

The documentation of the function is quite simple.  Just pass the name
of the file in which the output is and whether it is a file that
includes the output of the diagonal or not.  If you don't trust the
calculation of how many variables are involved, then you might want to
change the function so that this is another input paramter.

HTH.

Cheers,

        Berwin