ggbiplot设置分组_prcomp和ggbiplot:无效的’rot’值
我正在尝试使用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’值相关推荐
- ggbiplot设置分组_比PCA更好用的监督排序—LDA分析、作图及添加置信-ggord
线性判别分析LDA 线性判别分析,英文Linear Discriminant Analysis, 以下简称LDA.LDA在模式识别领域(比如人脸识别,舰艇识别等图形图像识别领域)中有非常广泛的应用,在 ...
- ggbiplot设置分组_R语言安装ggbiplot
ggbiplot画主成分图的案例,让人印象很深,但是用起来好像没那么容易,需要在github上下载安装.但是,ggbiplot在安装的时候经常遇到问题.按照网上的安装流程: install.packa ...
- pandas使用groupby函数对dataframe进行分组统计、使用as_index参数设置分组聚合的结果中分组变量不是dataframe的索引(index)
pandas使用groupby函数对dataframe进行分组统计.使用as_index参数设置分组聚合的结果中分组变量不是dataframe的索引(index) 目录
- R语言ggplot2可视化分组散点图、使用scale_shape_manual函数、scale_color_manual函数、scale_size_manual函自定义设置分组散点的形状、大小、颜色
R语言ggplot2可视化分组散点图.使用scale_shape_manual函数.scale_color_manual函数.scale_size_manual函数自定义设置分组散点的形状.大小.颜色 ...
- SSRS 动态设置分组依据及行组个数
动态设置行组及行组级数 业务场景: A公司想在同一张BI报告中通过切换参数,分别从不同角度分析统计产品销售情况, 同时又想同时从两个角度或三个角度统计分析. 本场景主要涉及两个技能点: 1) 动态设 ...
- 调试远程服务器上的代码时报错:调试设置中的Python路径无效
错误 当在本地调试远程服务器上的代码时,报错: 调试设置中的Python路径无效 踩坑 根据提示去launch.json文件里 以为是"python": "${comma ...
- 设置CFLAGS后,编译无效
问题: 1)设置CFLAGS后,编译无效 2)修改Makefile.am,编译无效 解决方法:./configure重新生成Makefile <script>window._bd_shar ...
- DialogFragment中通过dataBinding绑定View,设置点击事件无效,通过getWindow设置dialog位置和大小无效。
1.问题描述:DialogFragment中通过dataBinding绑定View,设置点击事件无效. private val binding: CoreDialogEditContentBindin ...
- 设置 核心内存转储后无效_从内存转储的第1部分(共3部分)调查无效的程序异常
设置 核心内存转储后无效 Datadog automated instrumentation for .NET works by rewriting the IL of interesting met ...
最新文章
- 4000字超干货!《统计学习方法》啃书指南(1)
- SpringMVC 原理和流程
- Java数据校验(Bean Validation / JSR303)
- 关于报表在移动端展现需你需要知道哪些?
- Drupal Working with nodes, content types and fields
- php java session共享_PHP实现session共享
- 我也想对广大程序员说一些话
- 安卓查看php文件是否存在,Android_Android编程判断SD卡是否存在及使用容量查询实现方法,本文实例讲述了Android编程判断 - phpStudy...
- 笨办法学 Python · 续 第七部分:大作业
- 需求跟踪矩阵模板_需求可追溯性矩阵(RTM)生成程序使用说明
- 第14题 计算机网络故障排除的第二步是,计算机网络故障诊断与排除第1章网络故障和网络诊断测试工具(习题)(ok)...
- Pandas系列(十五)stack和pivot实现数据透视
- Unity 2D角色复活点与复活等待时间设置
- adapthisteq和histeq
- 前端基础知识复习之html
- python词频统计并按词频排序
- SeaweedFS使用小结
- 如何导出微信好友?通讯录好友保存
- JavaScript绘制矢量图
- 【芯片前端】Filelist -f/-F的文件解析方式探究
热门文章
- Bailian4097 报站【暴力】
- Bailian3659 判断是否为C语言的合法标识符【文本处理】
- Bailian4013 中位数【中位数】
- 【数学】十万个为什么(一) —— 为什么乘法会分为左乘和右乘,除法会分为左除和右除?
- 强悍的 Linux —— tar
- C/C++基础::sizeof() 与 sizeof(string)
- np.random.RandomState、np.random.rand、np.random.random、np.random_sample
- linux buffer cache 过高_你真的理解Linux的内存监控吗?
- vsan Linux 网络磁盘,理解vSAN中的磁盘空间占用
- python编程100个小程序-100个Python练手小程序