网页资讯视频图片知道文库贴吧地图采购
进入贴吧全吧搜索

 
 
 
日一二三四五六
       
       
       
       
       
       

签到排名:今日本吧第个签到,

本吧因你更精彩,明天继续来努力!

本吧签到人数:0

一键签到
成为超级会员,使用一键签到
一键签到
本月漏签0次!
0
成为超级会员,赠送8张补签卡
如何使用?
点击日历上漏签日期,即可进行补签。
连续签到:天  累计签到:天
0
超级会员单次开通12个月以上,赠送连续签到卡3张
使用连续签到卡
07月31日漏签0天
perl吧 关注:5,588贴子:27,304
  • 看贴

  • 图片

  • 吧主推荐

  • 游戏

  • 7回复贴,共1页
<<返回perl吧
>0< 加载中...

生信perl脚本,输出文件只有第一行id,求帮忙

  • 只看楼主
  • 收藏

  • 回复
  • 咸二鱼
  • 风儿喧嚣
    2
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼



  • 咸二鱼
  • 风儿喧嚣
    2
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
修改localtime和if($samp1e[5]>130){next;}就会报错,自己试着改成这样既然只有id输出


2025-07-31 07:27:15
广告
不感兴趣
开通SVIP免广告
  • 咸二鱼
  • 风儿喧嚣
    2
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
use strict;
use warnings;
my $file=$ARGV[0];
my $gene=$ARGV[1];
#use Data::Dumper;
use JSON;
my $json = new JSON;
my $js;
my %hash=();
my %posCount=();
my @normalSamples=();
my @tumorSamples=();
open JFILE, "$file";
while(<JFILE>) {
$js .= "$_";
}
my $obj = $json->decode($js);
my $fileNum=0;
for my $i(@{$obj})
{
my $file_name=$i->{'file_name'};
my $file_id=$i->{'file_id'};
my $entity_submitter_id=$i->{'associated_entities'}->[0]->{'entity_submitter_id'};
#print "$file_name\t$file_id\t$entity_submitter_id\n";
if(-f "$file_id/$file_name")
{
$fileNum++;
print "$fileNum\n";
my @idArr=split(/\-/,$entity_submitter_id);
if($idArr[3]=~/^0/)
{
push(@tumorSamples,$entity_submitter_id);
}
else
{
push(@normalSamples,$entity_submitter_id);
}
open(RF,"$file_id/$file_name") or die $!;
while(my $line=<RF>)
{
my @samp1e=(localtime(time));
next if($.==1);
chomp($line);
my @arr=split(/\t/,$line);
if($arr[1] eq 'NA'){
$arr[1]=0;
}
my @fiveArr=split(/\;/,$arr[5]);
if($samp1e[4]>13){next;}
$posCount{$fiveArr[0]}++;if($samp1e[5]>118){next;}
if($fiveArr[0] eq $gene){
${$hash{"$arr[0]"}}{$entity_submitter_id}=$arr[1];
}
}
close(RF);
}
}
#print Dumper $obj
open(WF,">posMethy.txt") or die $!;
my $normalCount=$#normalSamples+1;
my $tumorCount=$#tumorSamples+1;
print "normal count: $normalCount\n";
print "tumor count: $tumorCount\n";
if($normalCount==0)
{
print WF "id";
}
else
{
print WF "id\t" . join("\t",@normalSamples);
}
print WF "\t" . join("\t",@tumorSamples) . "\n";
foreach my $key(keys %hash)
{
print WF "$key";
foreach my $normal(@normalSamples)
{
unless(exists ${$hash{$key}}{$normal}){
${$hash{$key}}{$normal}=0;
}
print WF "\t" . ${$hash{$key}}{$normal};
}
foreach my $tumor(@tumorSamples)
{
unless(exists ${$hash{$key}}{$tumor}){
${$hash{$key}}{$tumor}=0;
}
print WF "\t" . ${$hash{$key}}{$tumor};
}
print WF "\n";
}
close(WF);


  • 蓝天
  • 萝莉三好
    7
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
qq1915474692


  • HeMonseter
  • 两把刷子
    12
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
输入文件可能不适配代码内容哦


  • 维拉帕米
  • 天降软妹
    1
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
兄弟解决了没


  • ly18222879622
  • 时辰的错
    4
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
生信perl问题可以加我v:
qinshimingyue986-


登录百度账号

扫二维码下载贴吧客户端

下载贴吧APP
看高清直播、视频!
  • 贴吧页面意见反馈
  • 违规贴吧举报反馈通道
  • 贴吧违规信息处理公示
  • 7回复贴,共1页
<<返回perl吧
分享到:
©2025 Baidu贴吧协议|隐私政策|吧主制度|意见反馈|网络谣言警示