Login or Register to Ask a Question and Join Our Community

Sponsored Content
Full Discussion: diagonal matrix to square matrix
Top Forums Shell Programming and Scripting diagonal matrix to square matrix Post 302352541 by durden_tyler on Friday 11th of September 2009 06:22:45 PM
Old 09-11-2009
Here's a Perl solution for the type of data you posted:

Code:
$
$ cat diagmtx.txt
"244901_AT" 1.000 0.234 0.435 0.123 0.012 0.102 0.325 0.412 0.087 0.098
"243903_AT" 1.000 0.111 0.412 0.115 0.058 0.091 0.190 0.045 0.058
"244501_AT" 1.000 0.205 0.542 0.335 0.054 0.117 0.203 0.125
"254902_AT" 1.000 0.587 0.159 0.357 0.258 0.654 0.341
"247906_AT" 1.000 0.269 0.369 0.687 0.145 0.125
"242901_AT" 1.000 0.222 0.451 0.134 0.333
"243906_AT" 1.000 0.112 0.217 0.095
"244908_AT" 1.000 0.508 0.701
"294902_AT" 1.000 0.663
"245902_AT" 1.000
$
$ # show the Perl program
$ cat convert_2.pl
#!/usr/bin/perl -w
@mtx = ();
$file = $ARGV[0];
open (F,$file) or die "Can't open $file: $!";
while(<F>){
  @x = split;
  $id = shift @x;
  $therest = "@x\n";
  shift @x;
  push @mtx, [@x];
  if ($.==1) {
   print;
  } else {
    print $id," ";
    for($i=0;$i<=$.-2;$i++) {
      print $mtx[$i][0]," ";
      shift @{$mtx[$i]};
    }
    print $therest;
  }
}
close (F) or die "Can't close $file: $!";
$
$ perl convert_2.pl diagmtx.txt
"244901_AT" 1.000 0.234 0.435 0.123 0.012 0.102 0.325 0.412 0.087 0.098
"243903_AT" 0.234 1.000 0.111 0.412 0.115 0.058 0.091 0.190 0.045 0.058
"244501_AT" 0.435 0.111 1.000 0.205 0.542 0.335 0.054 0.117 0.203 0.125
"254902_AT" 0.123 0.412 0.205 1.000 0.587 0.159 0.357 0.258 0.654 0.341
"247906_AT" 0.012 0.115 0.542 0.587 1.000 0.269 0.369 0.687 0.145 0.125
"242901_AT" 0.102 0.058 0.335 0.159 0.269 1.000 0.222 0.451 0.134 0.333
"243906_AT" 0.325 0.091 0.054 0.357 0.369 0.222 1.000 0.112 0.217 0.095
"244908_AT" 0.412 0.190 0.117 0.258 0.687 0.451 0.112 1.000 0.508 0.701
"294902_AT" 0.087 0.045 0.203 0.654 0.145 0.134 0.217 0.508 1.000 0.663
"245902_AT" 0.098 0.058 0.125 0.341 0.125 0.333 0.095 0.701 0.663 1.000
$
$

I think summer_cherry's program is a much more optimized version. It -
(i) does not store the entire N X N square matrix in any data structure.
(ii) stores only the information present in the file, since that is sufficient to generate the other half of the square matrix.
(iii) uses hashes for fast access.

My second version uses a multi-dimensional array to store only the necessary information i.e. everything after the "id" and "1.000" per element. It also keeps on chopping the array element right after it is printed, by using the "shift" operator. So while the run time would be higher for this, the memory consumption should be lesser.

HTH,
tyler_durden
This User Gave Thanks to durden_tyler For This Post:
yifangt
 

9 More Discussions You Might Find Interesting

1. Programming

matrix pointer

Can anyone tell me what the following statements do? float (*tab); tab=(float (*)) calloc(MAXCLASS, (MAXCLASS+1)*sizeof(float)); (3 Replies)
Discussion started by: littleboyblu
3 Replies

2. Shell Programming and Scripting

Matrix

Hi All I would like to merge multiple files with the same row and column size into a matrix format In a folder I have multiple files in the following format vi 12.txt a 1 b 5 c 7 d 0 vi 45.txt a 3 b 6 c 9 d 2 vi 9.txt a 4 (7 Replies)
Discussion started by: Lucky Ali
7 Replies

3. Ubuntu

How to convert full data matrix to linearised left data matrix?

Hi all, Is there a way to convert full data matrix to linearised left data matrix? e.g full data matrix Bh1 Bh2 Bh3 Bh4 Bh5 Bh6 Bh7 Bh1 0 0.241058 0.236129 0.244397 0.237479 0.240767 0.245245 Bh2 0.241058 0 0.240594 0.241931 0.241975 ... (8 Replies)
Discussion started by: evoll
8 Replies

4. Shell Programming and Scripting

awk? adjacency matrix to adjacency list / correlation matrix to list

Hi everyone I am very new at awk but think that that might be the best strategy for this. I have a matrix very similar to a correlation matrix and in practical terms I need to convert it into a list containing the values from the matrix (one value per line) with the first field of the line (row... (5 Replies)
Discussion started by: stonemonkey
5 Replies

5. Shell Programming and Scripting

Table to Matrix

Hi, I have a table in the format: 1 0 -1 1 0 2 0 1 -1 0 0 0 3 0 1 1 0 0 0 0 0 0 etc. I am trying to input this to a program, however it is complaining about the fact that it is not in matrix format. How do I add 0's to end of the rows to make them even? Thanks in advance! (2 Replies)
Discussion started by: Rhavin
2 Replies

6. Shell Programming and Scripting

Square matrix to columns

Hello all, I am quite new in this but I need some help to keep going with my analysis. I am struggling with a short script to read a square matrix and convert it in two collumns. A B C D A 0.00 0.06 0.51 0.03 B 0.06 0.00 0.72 0.48 C 0.51 0.72 0.00 ... (7 Replies)
Discussion started by: EvaAM
7 Replies

7. Shell Programming and Scripting

Maybe by AWK: printing help diagonal matrix characters into line

Hi Experts, I want to print this charts diagonal data into straight lines. This is a matrix 24X24 Horizontal and vertical. - I want to print all the diagonal cutting characters into straight line: Data: E F S S A H A L L A T M C N O T S O B O D U Q H I W I B N L O C N I L N L A N S I N... (9 Replies)
Discussion started by: rveri
9 Replies

8. Shell Programming and Scripting

MATRIX to CSV

Hello friends, A big question for the UNIX INTELLIGENCE I have a CSV file as follows: VALUE,USER1,relatedUSER1,relatedUSER2 -1,userA,userB,userC 1,userN,userD,userB 0,userF,userH,userG 0,userT,userH,userB 1,userN,userB,userA -1,userA,userF,userC 0,userF,userH,userG... (15 Replies)
Discussion started by: kraterions
15 Replies

9. UNIX for Beginners Questions & Answers

Matrix multiplication

I have two files. Row id in File1 matches the column id in file2 (starting from column7 )except the last 2 characters. File1 has 50 rows and File 2 has 56 columns. If the id matches I want to multiply the value in column3 of File1 to the entire column in File2. and in the final output print only... (11 Replies)
Discussion started by: Akang
11 Replies

LEARN ABOUT DEBIAN

oakdecode

oakdecode(1)						      General Commands Manual						      oakdecode(1)

NAME

       oakdecode - Decode an OAKT printer stream into human readable form.

SYNOPSIS

       oakdecode [options] <OAKT-file

DESCRIPTION

       oakdecode decodes an OAKT printer stream into human readable form.

       An OAKT printer stream is the printer langauge used by the HP Color LaserJet 1500 and other printers.

COMMAND LINE OPTIONS

       These are the options that can appear on the command line.

       -d basename
	      Basename of .pbm file for saving decompressed planes.

       -r basename
	      Basename of .jbg file for saving raw planes

       -i     Suppress display of image records.

       -o     Print file offsets.

       -D level
	      Set Debug level [0].

EXAMPLES

       Decode an OAKT file created by foo2oak.

       $ oakdecode < testpage.oak
       0d(80) 1 OTHER
       0c(64) Wed Nov 05 16:30:50 2003        a07d3   100005  32001e
       0a(80) testpage.pdf
       14(16) (no args)
       28(16) Source=Tray1
       29(80) PaperType=0 UNK8=2,0,0,0, blanks(63)
       2a(32) Copies=1        UNK=0
       2b(32) papercode=25    xwid=4648       ywid=9000       UNK=0
       33(64)
	       u0      u1      w       h       resx    resy    nBits
	       x0      x0      2128    4300    600     600     x1
       15(16) (no args)
	    bih0     w	   h	l0     bih5   dlen  plen unk yOff P subP
       3c(64) 00010000 2176  256  256 58030020  1050  1056 000   64 3 0
		       DL = 0, D = 0, P = 1, - = 0, XY = 2176 x 256
		       L0 = 256, MX = 32, MY = 0
		       Order   = 3  ILEAVE SMID
		       Options = 88  LRLTWO TPDON TPBON
		       1 stripes, 0 layers, 1 planes
       3c(64) 00010000 2176  256  256 58030020  3668  3680 000  320 3 0
       3c(64) 00010000 2176  256  256 58030020  1463  1472 000  640 3 0
       3c(64) 00010000 2176  256  256 58030020  1975  1984 000  896 3 0
       3c(64) 00010000 2176  224  224 58030020  2744  2752 000 1152 3 0
       3c(64) 00010000 2176  256  256 58030020   988	992 000 1440 3 0
       3c(64) 00010000 2176  256  256 58030020  2892  2896 000 1696 3 0
       3c(64) 00010000 2176  256  256 58030020  3634  3648 000 1952 3 0
       3c(64) 00010000 2176  256  256 58030020  3236  3248 000 2208 3 0
       3c(64) 00010000 2176  256  256 58030020  2279  2288 000 2464 3 0
       3c(64) 00010000 2176  256  256 58030020  3746  3760 000 2720 3 0
       3c(64) 00010000 2176  200  200 58030020  2404  2416 000 2976 3 0
       3c(64) 00010000 2176  256  256 58030020  3114  3120 000 3240 3 0
       3c(64) 00010000 2176   96   96 58030020  1142  1152 000 3496 3 0
       3c(64) 00010000 2176  256  256 58030020  2094  2112 000 3752 3 0
       3c(64) 00010000 2176  256  256 58030020  1319  1328 000 4008 3 0
       3c(64) 00010000 2176   36   36 58030020   208	224 000 4264 3 0
       17(16) (no args)
       18(16) UNK=0
       0b(16) (no args)

FILES

       /usr/bin/oakdecode

SEE ALSO

       foo2oak-wrapper(1), foo2oak(1), jbg2pbm(1)

AUTHOR

       Rick Richardson <rick.richardson@comcast.net>
       http://foo2oak.rkkda.com/
All times are GMT -4. The time now is 10:47 AM.
Unix & Linux Forums Content Copyright 1993-2022. All Rights Reserved.
Privacy Policy