找到约 15 条结果

文章 python特殊属性和方法的详解与运用

...下划线作为变量前缀和后缀指定特殊变量 _xxx 不能用’from module import *’导入 __xxx__ 系统定义名字 __xxx 类中的私有变量名 核心风格:避免用下划线作为变量名的开始。   因为下划线对解释器有特殊的意义,而且是内建标...

文章 Cufflinks,Stringtie 合并转录本之后,如何筛选新转录本?

...y looking at the soft-masked reference sequence and applied to transcripts where at least 50% of the bases are lower case9uUnknown, intergenic transcript10xExonic overlap with reference on the opposite strand11sAn intron of the transfrag overlaps a reference intron on the opposite strand (likely due...

文章 convert 制作gif图片

...background white -extent 800x600 -colors 128 -layers Optimize *.png output.gif 逐段解释 convert:ImageMagick 主工具 -delay 100:每帧停留 1 秒(100=1 秒,数字越小越快) -loop 0:无限循环 *.png:所有 PNG 图片(你要求放最后,我会改) -resize "800x60...

文章 R语言apply批量进行t.test

...最终将所有p.value赋值给pvalue向量 pvalue=apply(data,1,function(x) t.test(x[1:3],x[4:6],paired = T)$p.value) 返回结果中针对第一第二行数据的pvalue值如下: > head(pvalue,2)TRINITY_DN100001_c0_g1 TRINITY_DN100002_c0_g1              0.7165794   ...

文章 如何看懂从TCGA下载的甲基化数据

...ociated with the genes detailed abovePosition to TSSDistance in base pairs from the CpG site to each associated transcript's start siteCGI CoordinateThe start and end coordinates of the CpG island associated with the CpG siteFeature TypeThe position of the CpG site in reference to the island: Island...

文章 BSA分析中MutMap 分析:95% 99%阈值置信区间是如何计算的

...#################################################################caluclate of genotype ratio#########################individuals_genotype<-function(number_of_total_individuals){ ratio_of_genotype<-c()for(i in 1:number_of_total_individuals){ratio_of_genotype<-c(ratio_of_genotype,genotype()) ...

问题 关于使用nohup vcftools过滤vcf文件的问题

...;ID=PID,Number=1,Type=String,Description="Physical phasing ID information, where each unique ID within a given sample (but not across samples) connects records within a phasing group">##FORMAT=<ID=PL,Number=G,Type=Integer,Description="Normalized, Phred-scaled likelihoods for genotypes as defin...

问题 braker报错

...ay, later. # Tue Jan 16 17:02:31 2024: Trying to guess $AUGUSTUS_BIN_PATH from $AUGUSTUS_CONFIG_PATH. # Tue Jan 16 17:02:31 2024: Setting $AUGUSTUS_BIN_PATH to /share/work/biosoft/augustus/latest/config//../bin # Tue Jan 16 17:02:31 2024: Did not find environment variable $AUGUSTUS_SCRIPTS_PATH (...

问题 变异结果质控过滤,去掉低质量变异结果时报错i

... ID=PID,Number=1,Type=String,Description="Physical phasing ID information, where each unique ID within a given sample (but not across samples) connects records within a phasing group"> Warning: Expected at least 2 parts in FORMAT entry: ID=PL,Number=G,Type=Integer,Description="Normalized, Phred-...

文章 ssh登陆报错

...nistrator.Add correct host key in /home/user01/.ssh/known_hosts to get rid of this message.Offending ECDSA key in /home/user01/.ssh/known_hosts:11Password authentication is disabled to avoid man-in-the-middle attacks.Keyboard-interactive authentication is disabled to avoid man-in-the-middle attacks....

问题 GWAS中gapit分析没有结果,命令行输入与屏幕输出如下

... Attaching package: ‘BiocGenerics’ The following objects are masked from ‘package:stats’:     IQR, mad, sd, var, xtabs The following objects are masked from ‘package:base’:     anyDuplicated, aperm, append, as.data.frame, basename, cbind,     colnames, dirname, do.call...

问题 合并gvcf的问题

...gvcf文件,但是提示我错误:A USER ERROR has occurred: The list of input alleles must contain <NON_REF> as an allele but that is not the case at position 763769; please use the Haplotype Caller with gVCF output to generate appropriate records。由于缺乏bam文件我只能尝试别...

问题 gwas分析群体画进化树

...istered S3 methods overwritten by 'treeio':   method              from   MRCA.phylo          tidytree   MRCA.treedata       tidytree   Nnode.treedata      tidytree   Ntip.treedata       tidytree   ancestor.phylo      tidytree   ancestor.treedata   tidytree ...

文章 图片去除大块黑色区域

...#########parser = argparse.ArgumentParser(description='This script is used to an error_bar')parser.add_argument('-i','--in_dir',help='Please input  in_dir directory path',required=True)args = parser.parse_args()pathDir =  os.listdir(args.in_dir)for allDir in pathDir:    if (re.search(r'png$',all...

文章 python版本MCScanX分析方法

...um The genomes of C. annuum, C.chinenese and C.baccatum were retrieved from PGP (http://peppergenome.snu.ac.kr/download.php).The C. annuum genome was compared to C.chinenese and C.baccatum  genomes using the MCScan toolkit (V1.1)[1]. To call synteny blocks, we performed all-against-all LAST [2]...