如何分析R中的挥发物?

时间:2019-06-14 08:51:21

标签: regression analysis

我混合了老植物和幼植物(o,y)和3种处理A,B,C(老植物和幼植物相同)的挥发性混合物,现在我想分析年龄和处理对品质和效果的影响挥发物的数量。我可以通过哪种统计测试解决该问题?进一步如何确定旧的/治疗C和年轻的/治疗C的有趣挥发物?

我尝试进行PCA,但随后恐怕并非所有具有影响力的挥发物都能得到实际解决。此外,我进行了MVA比较组。现在我不知道这是否有意义以及如何进行。另外,我尝试使用mixomics软件包创建一个热图,但是出现错误。

# First pair of components

plotDiablo(DIABLO,ncomp=1)

# Second pair of components

plotDiablo(DIABLO,ncomp=2)

class(DIABLO)[1] <- "block.splsda" # R-related step
circosPlot(DIABLO,cutoff=0.2,line=FALSE)

cimDiablo(DIABLO)
network(DIABLO,cutoff=0.8)
> class(DIABLO)[1] <- "block.splsda" # R-related step
> circosPlot(DIABLO,line=FALSE)
Error: 'cutoff' is missing
> circosPlot(DIABLO,cutoff=0.2,line=FALSE)
Error: Column names `Compound_1`, `Octane`, `3-Hexen-1-ol, (Z)`, `Cyclopentane, 1,2,3,4,5-pentamethyl`, `Nonanal`, ... (and 4 more) must not be duplicated.
Use .name_repair to specify repair.
Call `rlang::last_error()` to see a backtrace

0 个答案:

没有答案