网页
资讯
视频
图片
知道
文库
贴吧
地图
采购
进入贴吧
全吧搜索
吧内搜索
搜贴
搜人
进吧
搜标签
日
一
二
三
四
五
六
签到排名:今日本吧第
个签到,
本吧因你更精彩,明天继续来努力!
本吧签到人数: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
看高清直播、视频!
贴吧页面意见反馈
违规贴吧举报反馈通道
贴吧违规信息处理公示