利用perl或者R实现不规则文本的行列转置

最新推荐文章于 2023-04-26 04:19:15 发布

原创最新推荐文章于 2023-04-26 04:19:15 发布 · 2.8k 阅读

8 ·

CC 4.0 BY-SA版权

文章标签：

#perl #r语言

perl 同时被 2 个专栏收录

14 篇文章

订阅专栏

1 篇文章

订阅专栏

本文介绍如何利用R语言内置函数实现文本数据的行列转置操作，包括使用Perl代码与R语言代码对比演示，并展示转置后的结果。

如下是一个不规则的文本（行列中有空缺位），将其行列转置，求新的文本。

Array.lst：

1 aa bb cc dd

2 re uy oo jj

3 uh njk kjh jll

4 bgjkh klllj hhgj kk

5 hjj kk fjj

Perl代码：

#! /usr/bin/perl -w
use strict;
die "perl $0 \n" unless @ARGV==1;
my $lst=shift;
open IN,$lst;
my (@a,@b);
my $len;
my $max=0;
while(<IN>){
        chomp;
        @b=split/\t/,$_;
        $len=@b;
        $max=$max > $len ? $max:$len;
        push @a,[@b];
}
close IN;
for my $i(0..$max-1){
        for(@a){
                @$_[$i]||="";
                print "@$_[$i]\t";
        }
        print "\n";
}
__END__

命令行： perl transposition.array.pl array.lst > b

显然，单纯用perl实现起来比较繁琐，当然可以借助perl的相关模块实现快速转化，此处略去不讲。

接下来在R语言中，这个例子就很容易实现了，请看代码：

> a<-read.table(file="array.lst",fill=T,sep="\t")
> write.table(t(a),file="./b",row.names=F,col.names=F,quote=F,sep="\t")

结果为：b.

1 2 3 4 5

aa re uh bgjkh hjj

bb uy njk klllj kk

cc oo kjh hhgj

dd jj jll kk fjj