找到约 15 条结果

问题 xpclr 软件不能运行

...p --maxsnps 200 --minsnps 5donemkdir: cannot create directory ‘xpclr’: File exists-bash: xpclr: command not found-bash: xpclr: command not found-bash: xpclr: command not found-bash: xpclr: command not found-bash: xpclr: command not found-bash: xpclr: command not found-bash: xpclr: command not fo...

文章 rgb函数选择颜色(并控制深浅)

...alpha, names = NULL, maxColorValue = 1) red, blue, green, alphanumeric vectors with values in [0, M] where M is maxColorValue. When this is255, the red, blue, green, and alpha values are coerced to integers in0:255 and the result is computed most efficiently.namescharacter vector. The na...

文章 外源插入位置查找好方法

... Esher SK, Granek JA, Alspaugh JA. Rapid mapping of insertional mutations to probe cell wall regulation in Cryptococcus neoformans. Fungal genetics and biology : FG & B. 2015;82:9-21. doi:10.1016/j.fgb.2015.06.003.

文章 awk正则表达式使用外部变量

...号,再加单引号即可,如: for i in read write print; do cat file.txt | awk '{ if ($7~/'"$i"'/) {print $0} }'| wc -l ; done 在写shell脚本时,经常会使用到awk程序。但是有些复杂的逻辑,可能需要在awk中使用在shell中定义的变量,而且awk程序处...

文章 corrplot包绘制相关性图,简单又好看!

...deg") #读入数据 data<-read.table("fpkm.xls", header=TRUE,sep = "\t",check.names=FALSE,row.names=1) head(data) 基因表达量数据格式如下: ID A1-1 A1-2 A1-3 B1-1 B1-2 B1-3 A2-1 A2-2 A2-3 B2-1 B2-2 B2-3 FtPinG0000001800.01 15.1995160901025 ...

文章 NCBI 中搜索某类基因蛋白序列技巧

...filter] 筛选到可靠的结果之后,我们就可以,点击send to 把想要的基因的蛋白序列批量下载下来,同样的方法我们可以再搜索一下水稻的, 更多生物信息课程: 1. 文章越来越难发?是你没发现新思路,基因家族分析发2-4...

文章 Uniprot蛋白数据库介绍及使用详解!

...mapping,如下图,查询蛋白列表可直接粘贴在下图1. Provide your identifiers文本框中,也可以将蛋白ID单列粘贴于TXT文本中提交到网站。另外该页面2. Select options 还可提供ID转换功能,支持多种数据库间的ID转换。 提交好蛋白列表后...

问题 我在call snp,使用GATK导入db,得到的结果中没有callset.json文件,我一共15条染色体,db文件夹里只有两个染色体的文件夹内容,麻烦老师解答

...$tmpdir  -R $REF --batch-size 5 \ --reader-threads 5 --max-num-intervals-to-import-in-parallel 5 \ --genomicsdb-workspace-path db --sample-name-map cohort.sample_map

问题 用TBtools绘制两个物种的进化树+motif+基因结构三连图时报错问题

...========= IDs that you used as input do       NOT present in gff/gtf file.      Examples of valid IDs is showed as bellowedGuess FileType is:GTFtranscript:AT1G03987.1transcript:AT1G01020.6transcript:AT1G01010.1transcript:AT1G01020.3transcript:AT1G01020.2transcript:AT1G01020.5transcript:AT1...

文章 影响GWAS关联分析好坏的因素

...on in Biotechnology, 2006, 17(2):155–160. 5.Mcmullen M D, Stephen K, Hector Sanchez V, et al. Genetic Properties of the Maize Nested Association Mapping Population[J]. Science, 2009, 325(5941):737-. 6.Buckler E S, Holland J B, Bradbury P J, et al. The genetic architecture of maize flowering time...

文章 rainbow函数取色

...行设置: With rainbow, the parameters start and end can be used to specify particular subranges of hues. The following values can be used when generating such a subrange: red = 0, yellow = 1/6, green = 2/6, cyan = 3/6, blue = 4/6 and magenta = 5/6. 譬如取色由红色到蓝色: ...

文章 SNP指纹图谱软件SNPT

.../cgi-bin/snpt 这个软件可以做的分析: Indentify unique SNPs profiles and assign sequence types(MST);Calculate allele frequencies;Find informative SNPs;Calculate D (Index of the Discrininatory Ability);Find the set(s) of minimum SNPs;Output MST data file in mega and interleaved mega forma...

文章 R语言统计数据范围并用ggplot2绘制饼图

...ggplot2绘图饼图代码如下: library(reshape2) local({r <- getOption("repos") ;r["CRAN"] <- "http://mirrors.tuna.tsinghua.edu.cn/CRAN/" ;options(repos=r)}) library(ggplot2) library(cowplot) library(RColorBrewer) pairjoinC <- function(x,y){   ran1=paste(paste(head(x,-1),"%"...

文章 R 语言colors()657种颜色

...个颜色名称具体对应的颜色,可以通过绘图呈现出来 pdf(file = "RColors.pdf",height = 20,width = 30)par(mar=c(0,0,0,0))plot(0, 0, type = "n", xlim = c(0, 1), ylim = c(0, 1), axes = FALSE, xlab = "", ylab = "")line=30col=22rect(  rep((0:(col - 1)/col)...

问题 RNAseq有参转录组数据自主分析课程 ,运行建立基因组索引命令时,报错,无法获取gene_length.txt文件,请问是什么原因?如何处理呢?

...a_sinensis.gene.gtf -p gene_length Traceback (most recent call last):   File "/work/TS/scripts/get_gene_length_from_gtf.py", line 53, in <module>     if kvs['gene_id'] in geneL and kvs['transcript_id'] in geneL[kvs['gene_id']]: KeyError: 'gene_id' 基因组gff文件如下: gff...