找到约 15 条结果

问题 请教:基因结构分析时,输出的gene_exon_info.gff文件为空。该如何解决?

... ../script/get_gene_exon_from_gff.pl -in1 Carb_anhydrase_removed_redundant_and_confirmed_IDlist.txt -in2 ../Phaeodactylum_tricornutum.ASM15095v2.47.gff3 -out gene_exon_info.gff 图1 图2

文章 人类癌症数据库突变注释

cosmic70: the latest COSMIC database with somatic mutations from cancer and the frequency of occurence in each subtype of cancer. For more updated cosmic, see instructions below on how to make them.icgc21: International Cancer Genome Consortium version 21 mutations.nci60: NCI-60 human tumor cell lin...

问题 GSVA Error: useNames = NA is defunct. Instead, specify either useNames = TRUE or useNames = FALSE.

...rixStats 1.2 that released a couple of days ago. Can you downgrade to 1.1 and try?remotes::install_version("matrixStats", version="1.1.0") # restart your session and run previous scripts👍3

问题 单细胞测序数据FindAllMaker查找不同Cluster的Marker基因

...st one zero, cannot compute log geometric means data frame with 0 columns and 0 rows Error in `group_by_prepare()`: ! Must group by variables found in `.data`. • Column `cluster` is not found. Backtrace:     ▆  1. ├─obj.markers %>% group_by(cluster) %>% ...  2. ├─dplyr:...

文章 如何处理矫正GEO数据批次效应(batch effect)

...上面的数据利用李 dat = edata, batch = pheno$batchlibrary("limma")design <- model.matrix(~pheno$batch+pheno$hasCancer)limma.edata <- removeBatchEffect(edata,                                batch = pheno$batch,                                design = desi...

文章 运行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...

文章 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可能导致权限中断,可以尝试使用虚...

文章 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==...

文章 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 ...