找到约 15 条结果

文章 统计不同样品-生物学重复基因表达量数据合并和绘图统计

...的PCA分布图,查看样品之间的相互关系,要求不同的处理相图的形状区分,不同的stage不同的颜色区分,相同的stage颜色相同;   4 绘制Venn图,看看不同的stage 基因的表达变化等; 生物学重复按平均值,合并后的数据: ...

文章 hclust_analysis.r 转录组数据层次聚类分析

...类数量选择: 探究模型确定聚成几类合适 ,NbClust 包的 NbClust 函数确定聚类K值 ;  sum of squared error (SSE)组内平方误差和来确定最佳聚类数目,随着聚类数目增多,每一个类别数量越来越少,距离越来越近,因此WSS值肯...

问题 老师你好。请问共享目录挂载不成功,重装也一样,而且vbox屏幕不能全屏,no such device

图二为导入贵公司的linux后,打开设置的共享文件夹时,预置的两个共享文件夹。但在电脑没有查询到,我就把它删了。重新创建了共享文件夹。不删而启动Linux将出现图三

问题 perl脚本提取基因上游的1500bp(指定位置)序列时提示出现Incorrect parameters to subseq的错误

... 基因位置信息的文件:附件有上传(gen-weizhi.txt),大概格式如下: SETIT_034190mgIX4350641943510019+ SETIT_033939mgIX4555147145555562- SETIT_034271mgIX4764394847648195+ SETIT_034139mgIX5103987451043619+ SETIT_038909mgIX531353465...

文章 R脚本绘制多层环形饼图

“ 饼图(扇形图)时常于对统计数据的可视化,在 R,很多R包于绘制饼图,例如ggplot2、ggsunburst就能绘制出好看的饼图,操作如下” 1. 加载所需要的包 library(RColorBrewer)  #调色模板library(ggplot2)  #绘图library(dplyr)  #数据...

问题 xpclr报错AttributeError: 'NoneType' object has no attribute 'keys'

有哪位大佬过xpclr,我的命令:xpclr --out /hpc/home/t2024096/FST/out --format vcf --input /hpc/home/t2024096/data/cxx.vcf.gz --samplesA /txtlist/nong.txt --samplesB /txtlist/19.txt --size 200000 --step 100000 --chr chr1A出现以下报错 2025-02-18 17:34:41 : INFO : r...

文章 R-常见的字符串处理方法

...e",B,invert=T,value=T)[1] "abcdgegh" ##grepl字符串起始匹配也可以"^"锚定 > grepl("ab",A)[1] TRUE> grepl("ab",B)[1] TRUE TRUE 匹配替换 sub(pattern, replacement, x, ignore.case = FALSE, perl = FALSE, fixed = FALSE, useBytes = FALSE)gsub(pattern, replacement, x, ignore.cas...

文章 linux 命令当出现特殊字符无法传递参数怎么办

linux命令特殊符号包括:*、\、$、[、]、+、-、&、%、#、!、~、`、空格等,有特殊的含义,命令行要注意这些符号的使,不然会报错的。 举个例子,我这里有个目录:目录的名字间有个空格(aa bb  是一个目录,不...

文章 vcftools随机抽取指定个样品

...的vcf文件,到的选项为  --max-indv  ,指定要从vcf文件随机抽取指定个样品。 具体法如下: 下图为原始vcf文件。 随机抽取5个样品,执行以下代码: vcftools  --vcf  snp.vcf  --max-indv  5  --remove-indels  --recode  --out outfil...

文章 lm()多元回归分析结果解读

...F统计量和对应的p值。 数据存在一个叫myData的data.frame,数据前几列长下面这样:      ID        X        Y   Ele     Temp1 50353 126.6333 51.73333 173.9 0.0550242 50468 127.4500 50.25000 166.4 1.6190743 50564 127.3500 49.43333 234.5 1.4689174 50674 ...

文章 linux shell 将多行文件转换为一行

 1. 采awk       awk BEGIN{RS=EOF}'{gsub(/\n/," ");print}' file       说明:awk默认将记录分隔符(record separator即RS)设置为\n,此行代码将RS设置为EOF(文件结束),也就是把文件视为一个记录,然后通过gsub函数将\n替换成空格,最...

问题 各位老师好,我要分析的基因家族没有fa. pep. gff3. 和cds文件,请问还有其他办法进行分析吗?

NCBI只有fna.和gbff,其他的网址没有这个物种,真的不知道怎么做了,球球各位老师帮帮忙。

文章 2019年国自然结果分析,生命科学领域难度又又又增加了吗?

...达1194项,项目总金额达到6.4亿元,毫无疑问位列第一! 山大学立项数高达987项,位居第二位,项目总金额达到4.6亿。 浙江大学立项数达到838项,数目也很惊人。 立项数前十的高校还有华科技大学、复旦大学、北京大学、...

文章 IF=8.86|肝癌免疫相关lncRNA预后标志物文献解读

... 3.专注于肿瘤免疫浸润的特征在癌症的诊断、评估和治疗显示出有希望的预测和预后价值 数据分析思路 从TCGA数据库获得LIHC基因表达谱数据,ensemble获得gtf文件,根据注释文件得到lncRNA的表达数据;从ImmPort获得免疫相关基...

文章 Python模块导入的 路径索引

...Python到底是从哪里、以什么样的顺序调的呢? 在Python有内建函数(built-in)、第三方库(site-packages)以及自义库三种可以 调的模块。在 import时,Python的搜索顺序是先built-in模块,然后搜索 sys.path这个路径列表的模块,其...