Biostar Beta. Not for public use.
Perl Script To Make Statistics Of Mirna Abundances For Many Samples
0
Entering edit mode
2.2 years ago
biolab ♦ 1.1k

HI everyone,

I am making statistics of mirna abundances for many samples. Below is an example for input file.

SAMPLE    MIR    ABUNDANCE
sample1   mir1   30
sample1   mir3   100
sample1   mir4   120
sample2   mir1   40
sample2   mir2   200
sample3   mir1   190

......

I need to change the directions of tthe above matrix, and an ideal output is like below.

          sample1    sample2    sample3
mir1      30           40         190
mir2      0            200         0
mir3      190          0           0
mir4      120          0           0
......

i tried to write perl hash of hash to sort out the problem (see below). However, I am new with this relatively complex hash. Could anyone provide suggestions. I believe it's a good stuff to learn for other perl beginners too. Thank you very much!

open FH, '<', $ARGV[0] or die "open failed:$!";
my %h;
while (<>){
        my ($sample, $mir, $abun) = /(.+?)\t(.+)\t(.+)/;
        $h{$sample}{$mir} = $abun; 
}
foreach my $sample (keys %h){
        foreach my $mir (keys %{h{$sample}})
                print "   "      # i am stuck here. Need your help!
}
perl mirna • 2.4k views
ADD COMMENTlink
5
Entering edit mode

May I ask, is perl an absolute requirement? What you are looking for can be accomplished by a few commands using Pandas in Python, and I am happy to provide code for you if you would like to.

ADD REPLYlink
2
Entering edit mode

I agree with eric. Though it can be done with perl, it can more easily be accomplished by other means. In my experience, knowing how to work with R is very useful for numerical analyses. Have a look into the aggregate() function.

ADD REPLYlink
0
Entering edit mode

Thanks Eric and Irsan, I agree. Learning R is really necessary.

ADD REPLYlink
7
Entering edit mode
20 months ago
Neilfws 48k
Sydney, Australia

I realise that you asked for a Perl solution but as others noted in the comments, sometimes it's good to know that the right tool for the job exists already.

With that in mind, here is _acast()_ from the R/reshape2 package.

library(reshape2)
mi <- read.table("mi.txt", header = T, stringsAsFactors = F)
acast(mi, MIR ~ SAMPLE)

#      sample1 sample2 sample3
# mir1      30      40     190
# mir2      NA     200      NA
# mir3     100      NA      NA
# mir4     120      NA      NA
ADD COMMENTlink
0
Entering edit mode

Thanks a lot for all solutions. They are really helpful!

ADD REPLYlink
4
Entering edit mode
2.5 years ago
Glasgow, UK

Okay quickly wrote this perl one-liner for you assuming that your input file is tab-delimited (just redirect the output to another file):

Best Wishes, Umer

ADD COMMENTlink
3
Entering edit mode

This rather pushes the definition of "one liner" :)

ADD REPLYlink
0
Entering edit mode

saves u from one "wget" though ;)

ADD REPLYlink
0
Entering edit mode

This answer doesn't make sense.

ADD REPLYlink
0
Entering edit mode

Looks like something got lost in the transition to the new biostars.

ADD REPLYlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.1