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

 
 
 
日一二三四五六
       
       
       
       
       
       

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

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

本吧签到人数:0

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

  • 图片

  • 吧主推荐

  • 游戏

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

生信求助perl

  • 只看楼主
  • 收藏

  • 回复
  • 一粒小麦g
  • 天降软妹
    1
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
使用perl脚本逐一提取文件中的信息,却报错文件不存在

报错文件不存在

附上代码:
use strict;
use warnings;
#use File::Basename;
use XML::Simple;
#use Data::Dumper;
my @dirs=glob("*");
open(WF,">clinical.xls") or die $!;
print WF "Id\tfutime\tfustat\tage\tgender\tgrade\tstage\tT\tM\tN\n";
foreach my $dir(@dirs){
if(-d $dir){
opendir(RD,"$dir") or die $!;
while(my $xmlfile=readdir(RD)){
if($xmlfile=~/\.xml$/){
#print "$dir\\$xmlfile\n";
my $userxs = XML::Simple->new(KeyAttr => "name");
my $userxml = $userxs->XMLin("$dir\\$xmlfile");
# print output
#open(WF,">dumper.txt") or die $!;
#print WF Dumper($userxml);
#close(WF);
my $disease_code=$userxml->{'admin:admin'}{'admin:disease_code'}{'content'}; #get disease code
my $disease_code_lc=lc($disease_code);
my $patient_key=$disease_code_lc . ':patient'; #ucec:patient
my $follow_key=$disease_code_lc . ':follow_ups';
my $patient_barcode=$userxml->{$patient_key}{'shared:bcr_patient_barcode'}{'content'}; #TCGA-AX-A1CJ
if($patient_barcode eq "TCGA-AA-3521"){
print "$xmlfile\n";
}
my $gender=$userxml->{$patient_key}{'shared:gender'}{'content'}; #male/female
my $age=$userxml->{$patient_key}{'clin_shared:age_at_initial_pathologic_diagnosis'}{'content'};
my $race=$userxml->{$patient_key}{'clin_shared:race_list'}{'clin_shared:race'}{'content'}; #white/black
my $grade=$userxml->{$patient_key}{'shared:neoplasm_histologic_grade'}{'content'}; #G1/G2/G3
my $clinical_stage=$userxml->{$patient_key}{'shared_stage:stage_event'}{'shared_stage:clinical_stage'}{'content'}; #stage I
my $clinical_T=$userxml->{$patient_key}{'shared_stage:stage_event'}{'shared_stage:tnm_categories'}{'shared_stage:clinical_categories'}{'shared_stage:clinical_T'}{'content'};
my $clinical_M=$userxml->{$patient_key}{'shared_stage:stage_event'}{'shared_stage:tnm_categories'}{'shared_stage:clinical_categories'}{'shared_stage:clinical_M'}{'content'};
my $clinical_N=$userxml->{$patient_key}{'shared_stage:stage_event'}{'shared_stage:tnm_categories'}{'shared_stage:clinical_categories'}{'shared_stage:clinical_N'}{'content'};
my $pathologic_stage=$userxml->{$patient_key}{'shared_stage:stage_event'}{'shared_stage:pathologic_stage'}{'content'}; #stage I
my $pathologic_T=$userxml->{$patient_key}{'shared_stage:stage_event'}{'shared_stage:tnm_categories'}{'shared_stage:pathologic_categories'}{'shared_stage:pathologic_T'}{'content'};
my $pathologic_M=$userxml->{$patient_key}{'shared_stage:stage_event'}{'shared_stage:tnm_categories'}{'shared_stage:pathologic_categories'}{'shared_stage:pathologic_M'}{'content'};
my $pathologic_N=$userxml->{$patient_key}{'shared_stage:stage_event'}{'shared_stage:tnm_categories'}{'shared_stage:pathologic_categories'}{'shared_stage:pathologic_N'}{'content'};
$gender=(defined $gender)?$gender:"unknow";
$age=(defined $age)?$age:"unknow";
$race=(defined $race)?$race:"unknow";
$grade=(defined $grade)?$grade:"unknow";
$clinical_stage=(defined $clinical_stage)?$clinical_stage:"unknow";
$clinical_T=(defined $clinical_T)?$clinical_T:"unknow";
$clinical_M=(defined $clinical_M)?$clinical_M:"unknow";
$clinical_N=(defined $clinical_N)?$clinical_N:"unknow";
$pathologic_stage=(defined $pathologic_stage)?$pathologic_stage:"unknow";
$pathologic_T=(defined $pathologic_T)?$pathologic_T:"unknow";
$pathologic_M=(defined $pathologic_M)?$pathologic_M:"unknow";
$pathologic_N=(defined $pathologic_N)?$pathologic_N:"unknow";
my $survivalTime="";
my $vital_status=$userxml->{$patient_key}{'clin_shared:vital_status'}{'content'};
my $followup=$userxml->{$patient_key}{'clin_shared:days_to_last_followup'}{'content'};
my $death=$userxml->{$patient_key}{'clin_shared:days_to_death'}{'content'};
if($vital_status eq 'Alive'){
$survivalTime="$followup\t0";
}
else{
$survivalTime="$death\t1";
}
for my $i(keys %{$userxml->{$patient_key}{$follow_key}}){
my @survivalArr=split(/\t/,$survivalTime);
eval{
$followup=$userxml->{$patient_key}{$follow_key}{$i}{'clin_shared:days_to_last_followup'}{'content'};
$vital_status=$userxml->{$patient_key}{$follow_key}{$i}{'clin_shared:vital_status'}{'content'};
$death=$userxml->{$patient_key}{$follow_key}{$i}{'clin_shared:days_to_death'}{'content'};
};
if($@){
$followup=$userxml->{$patient_key}{$follow_key}{$i}[0]{'clin_shared:days_to_last_followup'}{'content'};
$vital_status=$userxml->{$patient_key}{$follow_key}{$i}[0]{'clin_shared:vital_status'}{'content'};
$death=$userxml->{$patient_key}{$follow_key}{$i}[0]{'clin_shared:days_to_death'}{'content'};
}
if($vital_status eq 'Alive'){
if($followup>$survivalArr[0]){
$survivalTime="$followup\t0";
}
}
else{
if($death>$survivalArr[0]){
$survivalTime="$death\t1";
}
}
}
print WF "$patient_barcode\t$survivalTime\t$age\t$gender\t$grade\t$pathologic_stage\t$pathologic_T\t$pathologic_M\t$pathologic_N\n";
}
}
close(RD);
}
}
close(WF);


  • 一粒小麦g
  • 天降软妹
    1
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
卡了好几天了,求大佬救救孩子


2025-05-17 16:10:40
广告
  • HeMonseter
  • 两把刷子
    12
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
贴吧里面有你的同学。你们俩的脚步一模一样。前几天她才问我了


  • HeMonseter
  • 两把刷子
    12
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
哦,我知道问题了,这个脚本他写的是用在win环境下的,win环境下的目录结构是用\\
你这个mac所以目录结构是用要用 /


  • 保护小青铜
  • 天降软妹
    1
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
你是怎么安装Perl的?我安装的只有PerlCritic


  • wangyang7611
  • 天降软妹
    1
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
请问楼主用这个脚本提取的临床信息有没有unknow的,我的TMN分期都是unknow


  • 丨依然轮回
  • 天降软妹
    1
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
同问,解决了吗


  • 吃糕的大笨钟
  • 风儿喧嚣
    2
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
同生信 运行merge脚本出现文件夹为空 不知道怎么回事


2025-05-17 16:04:40
广告
  • 灬沐小涵
  • 天降软妹
    1
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
楼楼可以把安装包给我发一份吗,我下载了好久了都没下好马上就是作业ddl了呜呜


  • 那年夜夜蒲
  • 天降软妹
    1
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
兄弟加个微信啊,339672430


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


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


  • 18892112835cy
  • 天降软妹
    1
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
友友们,生信刚入门,你们都是哪里 获取的perl的脚本啊


登录百度账号

扫二维码下载贴吧客户端

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