construct a matrix by some given values

[epidemiologist]

Dear perl Forum,

I have no clue how to start this. I have a file with values like this:
A 2 - 4
B 3 - 5
C 1 - 6

up to the end

I need to formulate a matrix like this

ID 1 2 3 4 5 6 7 8 9 10
A 0 1 1 1 0 0 0 0 0 0
B 0 0 1 0 0 0 0 0 0 0
C 1 1 1 1 1 1 0 0 0 0

up to the end

print one in the range of the value the individual have and print 0 their before and after it. I have no clue how to start doing this.

If you guide me this will be a great help.

Thanks,
Dalila

[miller]

Greetings Dalila,

One way would be the following:

use strict;
 
# Goal:
# ID 1 2 3 4 5 6 7 8 9 10
# A 0 1 1 1 0 0 0 0 0 0
# B 0 0 1 0 0 0 0 0 0 0
# C 1 1 1 1 1 1 0 0 0 0
 
print "ID " . join(' ', (1..10)) . "\n";
 
while (<DATA>) {
    if (! /(\w+) (\d+) - (\d+)/) {
        print "Invalid line format: $_";
 
    } else {
        my $id = $1;
        my $min = $2;
        my $max = $3;
 
        print "$id ";
        print join ' ', map {$_ >= $min && $_ <= $max ? 1 : 0} (1..10);
        print "\n";
    }
}
 
__DATA__
A 2 - 4
B 3 - 5
C 1 - 6
 

Output is:

>perl scratch.pl
ID 1 2 3 4 5 6 7 8 9 10
A 0 1 1 1 0 0 0 0 0 0
B 0 0 1 1 1 0 0 0 0 0
C 1 1 1 1 1 1 0 0 0 0
 

- Miller

[epidemiologist]

Many thanks Miiler,
I will try it. If the values are irregulary incremented, I think I can put them into array and use the function map as you guided me?
I really very thankful for your guide.
Dalila

[miller]

Nothing says that you have to use map. That is simply the method that I chose. You could have just as easily used a for loop.

- Miller

[KevinADC]

There are also a bunch of matrix modules on cpan that might be worth looking into.

[epidemiologist]

Hi Kevin and Miller,
I have done the following and it was working perfectly

#!/usr/bin/perl -w
 
use strict;
 
open(FILE, "IDs.txt");
open(OUT, ">matrix.txt");
 
my @array = (1 , 2 , 3 , 4, 8 , 9);
print OUT "ID", "\t" . join("\t", (@array)) . "\n";
 
while(defined(my $line=<FILE>)) {
    chomp($line);
    my @data = split("\t", $line);
    # to make sure the IDs and intervals are seprated by tabs
    if (!($line=~ /(\w+)(\t)(\d+)(\t)(\d+)/)) {
        print "Invalid line format: $_", "\n";
    } else {
        my $id = $data[0];
        my $min = $data[1];
        my $max = $data[2];
 
        print OUT ($id , "\t");
 
        print OUT join "\t" , map {$_ >= $min && $_ <= $max ? 1 : 0} (@array);
        print "\n";
    }
}
 

I read the IDs and interval from files and reprinted them as you guided me. However, I found some IDs repeated more than once
e.g

A 1 2
A 8 9

I want to print all the intervals for each ID in one file. So, will have
ID 1 2 3 4 8 9
A 1 1 0 0 1 1

Can you help me in a clever way of doing it.
Thanks,
Dalila

[KevinADC]

is the input data different than your origianl question? Before it was ranges (n-n) , now it looks like individual values (n n) . Do you really want each id (A B C etc) as a seperate output file?

[epidemiologist]

is the input data different than your origianl question? Before it was ranges (n-n) , now it looks like individual values (n n) . Do you really want each id (A B C etc) as a seperate output file?

Hi Miller,
THe input data will be like this (I relaized I have some headers not in interval)

A 1 2
A 7 10
A 15 20
B 3 5
B 11 15
C 5 10
D 10 20

I found some intervalsc not to be equal, it is a fixed set of numbers
e.g @array = (1 , 2, 3 , 5 , 7, 10 , 11, 15 , 20)

I want all the IDs (A, B, C, D) to be in one output file, but each in only one line:
Aim

ID 1 2 3 5 7 10 11 15 20
A 1 1 0 0 1 1 0 1 1
B 0 0 1 1 0 0 1 1 0
C 0 0 0 1 1 1 0 0 0
D 0 0 0 0 0 1 1 1 1

So, all the vlaues will be "map" to the arrays but do not know how to look further to see if the Id is coming again (Id could be repeated with different values for up to 100 times! or more no fixed rule)

I tried to work it out but seems little but hard for me.
Thanks,
Dalila

[KevinADC]

hmm... I thought you had abondoned this thread. My interest has wanned during the interlude.

[epidemiologist]

Hi Kevin,
I really apologize. If any one can help this will be great. Pleaaaaaaaaaaaaaase help!
Dalila

[epidemiologist]

Thanks Guy,
I have figured it out!!

[KevinADC]

You figured it out? I think our buddy Travis on Tek-Tips figured it out for you. ;)