ℚuantStudio Real-Time PCR software 是我们经常使用的 RT-PCR 软件, 它上面的可视化只能简单看看, 不满足发论文的需求. 如果需要得到发表级的图片, 还是需要 用 ggplot 大法加持. 为了能够使用这些数据, 首先需要导出文件. 为了方便操作, 文件导出时, 选择 *.txt 格式, 每个面板导出成一个单独文件. 将文件放在 data 文件夹中. … 阅读全文 →

上次提及 PCA 分析的方法有很多种。那不同方法之间的得到的结果会有差异吗？ 最近采用 PCA 分析 RNA-seq 样本之间的差异，得到了下面的结果。 […] 生成一个含有 1000 个基因, 27 个样品的数据集. 这 27 个样品来自于 3 个基因型(WT, Mutant1, Mutant2), 3 种处理(CK, Trt1, Trt2), 共分为 \(3 * 3 = 9\) … 阅读全文 →

有一个数据分析的小项目，使用Git记录分析过程。 在添加了一些文件、删除了一些文件之后，在RStudio的Git标签页下面点击复选框Stage，界面失去响应。反复重启，仍然如此。 上网查了一下，确实很多人吐槽类似的问题。 不过，在当前Git目录下打开Git GUI使用，一切正常。 RStudio的Git GUI性能确实比较差。 阅读全文 →

给自己的网站添加SSL证书其实并不复杂，但是关键一点就是首先要拥有一个SSL证书。由于SSL证书价格不菲，很多个人站长会选择放弃使用Https。但是，自从开源、免费的Let’s Encrypt证书出现后，我觉得SSL也是我们草根站长可以玩的了。 https://blog.csdn.net/ithomer/article/details/78075006 […] … 阅读全文 →

相信大家在做微生物多样性研究时经常听到PCA分析、PCoA分析，NMDS分析，CCA分析，RDA分析。它们对物种（或基因、功能）的分析具有重要作用，因而频频出现在16S测序及宏基因组测序中。以上分析本质上都属于排序分析（Ordination analysis）。 排序分析（ordination analysis），最早是生态学（ecology）中研究群落（communities）的一大类多元分析手 … 阅读全文 →

𝕝ibrary(vegan) ## Loading required package: permute ## Loading required package: lattice ## This is vegan 2.6-4 data("varespec") pca <- rda(varespec) 首先看一下结果： summary(pca) ## ## Call: ## … 阅读全文 →

本文介绍 Bray Curtis Dissimilarity 的概念和计算方法。 Bray Curtis Dissimilarity（Bray-Curtis 相异度）是生态学中用来衡量不同样地物种组成差异的参数。 其定义和计算公式为： $$BC_{ij}=1-2C_{ij}/(S_{i}+S_{j})$$ 其中： […] 有两个水族箱： […] … 阅读全文 →

本分析方法来在2018年在 ISME J 发表的一篇论文1，原来的分析是用 Mathematica 做的，我尝试在 R 里面实现它。 […] Mallon CA et al. (2018). The impact of failure: unsuccessful bacterial invasions steer the soil microbial community away … 阅读全文 →

设置 blogdown 主要参考了 这里，设置 git 认证参考了 这里. 软件的安装不提，简单的步骤如下： 阅读全文 →

我的工作环境是Windows，但要时常用到Linux，于是在VirtualBox中部署了Ubuntu Linux系统。 以前都是打开虚拟机上使用，但是虚拟机毕竟有桌面，占用资源比较多。窗口切来切去也不方便。 设置虚拟机SSH登录可以让虚拟机作为一个进程运行在后台，很Linux。 […] 要安装SSH server。 sudo apt-get install … 阅读全文 →