找到约 15 条结果

文章 核苷酸替换模型-GTR是所有模型中考虑参数最多最复杂的模型,其他所有模型是GTR模型的特例

...r DNA site substitution. The models are listed here from the least complex to the most parameter rich. Jukes-Cantor (JC, nst=1): equal base frequencies, all substitutions equally likely (PAUP* rate classification: aaaaaa, PAML: aaaaaa) (Jukes and Cantor 1969) Felsenstein 1981 (F81, nst=1): variabl...

问题 老师您好,我在做mcscanx共线性分析基因家族加倍实操课程时,安装mcscanx报错

根据咱们已经回答的问题,把报错make: *** [mcscanx] Error 1,进行了修改,只需要给MCScanX 目录下的 msa.h, dissect_multiple_alignment.h, and detect_collinear_tandem_arrays.h 这三个文件内容 最前面添加#include <unistd.h>。 但是又出现了另一种报...

问题 老师,我在用比较基因组学方法做两个物种之间的共线性分析时,是可以筛到某一种基因家族在两个物种中对应的共线性的,但是用下面这种方法做是做不出来某一基因家族在两物种之间的共线性关系的

... /biosoft/MCScanX/MCScanX/MCScanX ST_AT 得到以下结果 Reading BLAST file and pre-processingGenerating BLAST list170 matches imported (1 discarded)44 pairwise comparisons0 alignments generatedPairwise collinear blocks written to ST_AT.collinearity [0.013 seconds elapsed]Writing multiple synten...

文章 提取基因组所有基因的启动子序列

...nfig::General;use Cwd qw(abs_path getcwd);use FindBin qw($Bin $Script);use File::Basename qw(basename dirname);use Bio::SeqIO;use Bio::Seq;my $version = "1.3"; ## prepare parameters ######################################################################### --------------------------------------------...

问题 ggplot2画图

...F VS P-F " ##è¯»å–æ•°æ® demo_go=read.table(gofile,sep="\t",check.names = F,header = T,comment.char = "") head(demo_go) 我输入需要做图的enrichment .txt 报错如图 请问哪里出了问题

文章 perl对文件夹处理

...opendir ( DIR, $dirname ) || die "Error in opening dir $dirnamen";while(($filename = readdir(DIR))){print("$filenamen");       #循环输出该目录下的文件。### do something ###}closedir(DIR);   例 把每个文件复制一次: #!/usr/bin/perl use warnings;use stric...

问题 老师,您好,找不到临床数据该如何解决?

> # 读取输入的参数 > expr_file <- "/Volumes/HUBGENE/TCGA/methylation/lab/GDC/TCGA_OV_methylation.txt" > clin_file <- "/Volumes/HUBGENE/TCGA/methylation/lab/Clinical/GDC/TCGA_OV_clinical.txt" > work_dir <- "/Volumes/HUBGENE/TCGA/methylation/lab" > p_cutoff <- 0.0...

文章 stacks RAD 分析方法

... processed by Stacks version 2.6[1].We used the process_radtags function to clean the tags, discarding reads of low quality (-q), removing reads with uncalled bases (-c) and rescuing barcodes and radtags (-r). SNP calling for each individual was done using the program denovo_map.pl of Stacks. Th...

文章 TCGA临床信息如何正确的保存格式

...如,采用星号(*)去分割。 # 将数据保存到文件 clinical_file <- paste0(DataDirectory, "_","clinical",".txt") write.table(clinical, file = clinical_file, row.names = F, col.names=T,quote = T,sep='*',na = "NA") 之后打开excel 时,对数据进行分列,同样选用星号...

问题 单倍型报错

...cter",default=getwd(),                     help="input file directory [default %(default)s]",                     metavar="filepath")parser$add_argument( "-v", "--vcf", type="character",required=T,                     help=" vcf files[requi...

问题 在GDC数据库下载数据出错

...GDCdownload(query = query, +             directory = DataDirectory,files.per.chunk=6, method='client') Downloading data for project TCGA-CHOL 试开URL’https://gdc.cancer.gov/files/public/file/gdc-client_v1.6.0_Windows_x64-py3.7_0.zip' Content type 'application/zip' length 16207515 byte...

问题 老师,我利用脚本得到对应基因的蛋白序列: perl script/get_fa_by_id.pl WRKY_removed_redundant_IDlist.txt Arabidopsis_thaliana.TAIR10.pe p.all.fa WRKY_pep_need_to_confirm.fa 时,程序出现错误。 ,这是怎么回事啊?

这是输出的结果 这是脚本

问题 GEOquery 下载不了数据

...   AnnotGPL = F,                 getGPL = F)  save(eSet,file='GSE112093_eSet.Rdata')}load('GSE112093_eSet.Rdata')b = eSet[[1]]raw_exprSet=exprs(b) phe=pData(b)library(stringr)group_list= str_split(as.character(phe$title),' ',simplify = T)[,1]save(raw_exprSet,group_list,     fi...

问题 向管理员要了metagenomics:v1.0的本地镜像,但是输入 docker load -i metagenomics-v1.0.tar.gz后一直报错 open /var/lib/docker/tmp/docker-import-2588968250/metagenomics/json: no such file or directory

问题 TCGAbiolinks使用求助 Downloading genome information (try:0) Using: Human genes (GRCh38.p13) Error: $ operator is invalid for atomic vectors

...pare函数的时候,出现报错,请各位给予分析或者 Starting to add information to samples  => Add clinical information to samples Add FFPE information. More information at: => https://cancergenome.nih.gov/cancersselected/biospeccriteria => http://gdac.broadinstitute.org/...