#!/usr/bin/perl use strict; # https://perlmonks.org/?node_id=11146766 use warnings; use List::AllUtils qw( zip_by ); my @matrix = ( [qw(t1 t1 t2 t2 t2)], # <- transcripts [qw(a1 a2 a1 a1 a2)], # <- alleles [qw(intron intron UTR_CG UTR UTR)], # <- locations ); use Data::Dump 'dd'; dd 'original matrix', \@matrix; my @t = zip_by { [ @_ ] } @matrix; use Data::Dump 'dd'; dd '@t', \@t; my %cg = map { $_->[0] . "\n" . $_->[1] => 1 } grep { $_->[2] eq 'UTR_CG' } @t; use Data::Dump 'dd'; dd '%cg', \%cg; @t = grep { $_->[2] ne 'UTR' or not $cg{ $_->[0] . "\n" . $_->[1] } } @t; use Data::Dump 'dd'; dd 'modified @t', \@t; my @finalmatrix = zip_by { [ @_ ] } @t; use Data::Dump 'dd'; dd 'final matrix', \@finalmatrix;