我正在尝试使用R对我的数据进行PCA分析,我使用prcomp和ggbiplot找到了

this nice guide.我的数据是两种样品类型,每种样品具有三个生物重复(即6行)和约20000个基因(即变量).首先,使用指南中描述的代码获取PCA模型不起作用:

>pca=prcomp(data,center=T,scale.=T)

Error in prcomp.default(data, center = T, scale. = T) :

cannot rescale a constant/zero column to unit variance

但是,如果我删除了比例. = T部分,它工作正常,我得到一个模型.这是为什么,这是下面错误的原因?

> summary(pca)

Importance of components:

PC1 PC2 PC3 PC4 PC5

Standard deviation 4662.8657 3570.7164 2717.8351 1419.3137 819.15844

Proportion of Variance 0.4879 0.2861 0.1658 0.0452 0.01506

Cumulative Proportion 0.4879 0.7740 0.9397 0.9849 1.00000

其次,绘制PCA.即使只是使用基本代码,我得到一个错误和一个空的情节图像:

> ggbiplot(pca)

Error: invalid 'rot' value

这意味着什么,我该如何解决?是否与制作PCA的(非)规模有关,还是有所不同?我认为它必须与我的数据有关,因为如果我使用标准示例代码(下面),我会得到一个非常好的PCA图.

> data(wine)

> wine.pca=prcomp(wine,scale.=T)

> print(ggbiplot(wine.pca, obs.scale = 1, var.scale = 1, groups = wine.class,

ellipse = TRUE, circle = TRUE))

[编辑1]我尝试以两种方式对我的数据进行子集化:1)删除所有行都是0的所有列,以及2)删除所有列,任何行都是0.第一个子集仍然给出了比例误差,但不是已删除任何0的列的那些.为什么是这样?这对我的PCA有何影响?

此外,我尝试对原始数据(非缩放)和上面的子集化数据使用普通的biplot命令,并且它在两种情况下都有效.那么这与ggbiplot有关吗?

[编辑2]我已经上传了我的数据子集,当我不删除所有零并且在我执行时工作时,会给我错误.我之前没有用过吉斯特,但我认为this就是这样.或者this ……

ggbiplot设置分组_prcomp和ggbiplot:无效的’rot’值相关推荐

  1. ggbiplot设置分组_比PCA更好用的监督排序—LDA分析、作图及添加置信-ggord

    线性判别分析LDA 线性判别分析,英文Linear Discriminant Analysis, 以下简称LDA.LDA在模式识别领域(比如人脸识别,舰艇识别等图形图像识别领域)中有非常广泛的应用,在 ...

  2. ggbiplot设置分组_R语言安装ggbiplot

    ggbiplot画主成分图的案例,让人印象很深,但是用起来好像没那么容易,需要在github上下载安装.但是,ggbiplot在安装的时候经常遇到问题.按照网上的安装流程: install.packa ...

  3. pandas使用groupby函数对dataframe进行分组统计、使用as_index参数设置分组聚合的结果中分组变量不是dataframe的索引(index)

    pandas使用groupby函数对dataframe进行分组统计.使用as_index参数设置分组聚合的结果中分组变量不是dataframe的索引(index) 目录

  4. R语言ggplot2可视化分组散点图、使用scale_shape_manual函数、scale_color_manual函数、scale_size_manual函自定义设置分组散点的形状、大小、颜色

    R语言ggplot2可视化分组散点图.使用scale_shape_manual函数.scale_color_manual函数.scale_size_manual函数自定义设置分组散点的形状.大小.颜色 ...

  5. SSRS 动态设置分组依据及行组个数

    动态设置行组及行组级数 业务场景: A公司想在同一张BI报告中通过切换参数,分别从不同角度分析统计产品销售情况, 同时又想同时从两个角度或三个角度统计分析. 本场景主要涉及两个技能点: 1)  动态设 ...

  6. 调试远程服务器上的代码时报错:调试设置中的Python路径无效

    错误 当在本地调试远程服务器上的代码时,报错: 调试设置中的Python路径无效 踩坑 根据提示去launch.json文件里 以为是"python": "${comma ...

  7. 设置CFLAGS后,编译无效

    问题: 1)设置CFLAGS后,编译无效 2)修改Makefile.am,编译无效 解决方法:./configure重新生成Makefile <script>window._bd_shar ...

  8. DialogFragment中通过dataBinding绑定View,设置点击事件无效,通过getWindow设置dialog位置和大小无效。

    1.问题描述:DialogFragment中通过dataBinding绑定View,设置点击事件无效. private val binding: CoreDialogEditContentBindin ...

  9. 设置 核心内存转储后无效_从内存转储的第1部分(共3部分)调查无效的程序异常

    设置 核心内存转储后无效 Datadog automated instrumentation for .NET works by rewriting the IL of interesting met ...

最新文章

  1. 4000字超干货!《统计学习方法》啃书指南(1)
  2. SpringMVC 原理和流程
  3. Java数据校验(Bean Validation / JSR303)
  4. 关于报表在移动端展现需你需要知道哪些?
  5. Drupal Working with nodes, content types and fields
  6. php java session共享_PHP实现session共享
  7. 我也想对广大程序员说一些话
  8. 安卓查看php文件是否存在,Android_Android编程判断SD卡是否存在及使用容量查询实现方法,本文实例讲述了Android编程判断 - phpStudy...
  9. 笨办法学 Python · 续 第七部分:大作业
  10. 需求跟踪矩阵模板_需求可追溯性矩阵(RTM)生成程序使用说明
  11. 第14题 计算机网络故障排除的第二步是,计算机网络故障诊断与排除第1章网络故障和网络诊断测试工具(习题)(ok)...
  12. Pandas系列(十五)stack和pivot实现数据透视
  13. Unity 2D角色复活点与复活等待时间设置
  14. adapthisteq和histeq
  15. 前端基础知识复习之html
  16. python词频统计并按词频排序
  17. SeaweedFS使用小结
  18. 如何导出微信好友?通讯录好友保存
  19. JavaScript绘制矢量图
  20. 【芯片前端】Filelist -f/-F的文件解析方式探究

热门文章

  1. Bailian4097 报站【暴力】
  2. Bailian3659 判断是否为C语言的合法标识符【文本处理】
  3. Bailian4013 中位数【中位数】
  4. 【数学】十万个为什么(一) —— 为什么乘法会分为左乘和右乘,除法会分为左除和右除?
  5. 强悍的 Linux —— tar
  6. C/C++基础::sizeof() 与 sizeof(string)
  7. np.random.RandomState、np.random.rand、np.random.random、np.random_sample
  8. linux buffer cache 过高_你真的理解Linux的内存监控吗?
  9. vsan Linux 网络磁盘,理解vSAN中的磁盘空间占用
  10. python编程100个小程序-100个Python练手小程序