找到约 15 条结果

文章 运行python STAMP程序报错

...don't work.You need to also specify an X-11 server on the host machine to handle the launch of GUId programs. If the client is a Windows machine install Xming. If the client is a Linux machine install XQuartz.Now suppose this is Windows connecting to Linux. In order to be able to launch X11 programs...

文章 pennCNV软件linux系统安装

... -c -o khmm_wrap.o khmm_wrap.c khmm_wrap.c: In function ‘SWIG_AsCharPtrAndSize’: khmm_wrap.c:1576: error: ‘bool’ undeclared (first use in this function) khmm_wrap.c:1576: error: (Each undeclared identifier is reported only once khmm_wrap.c:1576: error: for each function it appears in.) khm...

问题 使用vcftools过滤vcf文件缺失位点时报错

#代码 vcftools --gzvcf xxx.vcf.gz --max-missing 1.0 --recode --stdout | gzip -c > xxx.recode.vcf.gz #报错 Error: Require Genotypes in variant file to filter by frequency and/or call rate

文章 借助perl脚本进行文件整理 | 将按第一列信息整理的文件改成按照第二列信息整理

...有perl脚本如下:print "perl $0  <id_list> <out_file>\n" and die unless(@ARGV==2);open IN,"$ARGV[0]" or die "$!; can't open file $ARGV[0]\n";open OUT,">$ARGV[1]" or die "$!; can't open file $ARGV[1]\n";my%id;my%hash;while(<IN>){        chomp;        my@a=split(/\s/...

文章 Root microbiota assembly and adaptive differentiation among European Arabidopsis populations

...动欧洲各地的细菌菌群组成趋同。3. Root endosphere bacteria and A. thaliana exhibit stable associations(根内细菌和拟南芥种群存在稳定联系)首先作者将在80%以上地方都存在的微生物OTUs定义为地理广泛存在的OTUs,20%~80%定义为普遍存在的OTUs...

文章 GWAS分析哪个模型最优

...并绘制QQ图。CW trait (H = 80% 下图A) 关联结果看, ANOVA, GLM, and SUPER假阳性高;MLM, CMLM, and ECMLM这几个模型只是尾部少量SNP翘起说明假阳性减少但可能存在假阴性,因为尾部翘起的点较少。而FarmCPU 尾部高高翘起而且多说明这个模型假...

文章 筛选编码蛋白的CDS序列

...args.out_dir)    dout=os.path.abspath(args.out_dir)count=0total=0output_handle = open(dout+'/'+args.prefix+'.fa', "w")for rec in SeqIO.parse(args.fasta, "fasta"):    seq=str(rec.seq)    #seq=seq.upper()    total+=1    if(len(seq)>=args.len and re.search("^ATG",seq,re.I) and len(seq)%3==...

文章 pip报错

...WARNING: Running pip as the 'root' user can result in broken permissions and conflicting behaviour with the system package manager. It is recommended to use a virtual environment instead: https://pip.pypa.io/warnings/venv 以“root”用户运行pip可能导致权限中断,可以尝试使用虚...

文章 exp_immu_cor_plot.r 基因表达量与免疫侵润相关性绘图分析

...e/1488 optional arguments: -h, --help show this help message and exit -e expset, --expset expset set event column name [required] -i immu, --immu immu input data file path[required] -m method, --method method a character string ind...

文章 multi_cox.r 基因表达量做多因素cox分析

...onoptional arguments:  -h, --help            show this help message and exit  -i data, --data data  input data file path[required]  -t time, --time time  set suvival time column name [required]  -e event, --event event                        set event column name must 0 or 1 ...

文章 基因家族的定义和形成方式(加倍与扩张)

...f several similar genes, formed by duplication of a single original gene, and generally with similar biochemical functions.  基因家族的形成方式(加倍与扩张): 1.染色体組多倍化( Whole genome duplication or  polyploidization):大部分植物都经历过古老的全...

问题 安装github上的包报错:

...) Using github PAT from envvar GITHUB_PAT. Use `gitcreds::gitcreds_set()` and unset GITHUB_PAT in .Renviron (or elsewhere) if you want to use the more secure git credential store instead. Error: Failed to install 'unknown package' from GitHub:   HTTP error 401.   Bad credentials   Rate lim...

问题 TCGAbiolinks包

你好,我遇到了跟你一样的问题Downloading genome information (try:0) Using: Human genes (GRCh38.p13) 错误: $ operator is invalid for atomic vectors 请问怎么更新R包解决呢?

问题 TCGA的RNA-seq count annotation过滤,有些不知道是否需要删除

...ion:General:This subject was duplicated by TCGA-24-0981; 0981 was redacted and this one retained 2. patient_annotation:Notification:Item in special subset:This is part of the ovarian triplet set. 3. aliquot_annotation:CenterNotification:Item flagged DNU:QC Warning- Expression quality metrics ...

问题 Error running system command:

...ant这几个大小为0,运行提示如下NOTICE: Running with system command <convert2annovar.pl  -includeinfo -allsample -withfreq -format vcf4 /work/my_reseq/4.snp_indel/all.clean.snp.vcf.gz > /work/my_reseq/5.var_ann/snp.avinput> NOTICE: Finished reading 5099177 lines from VCF file ...