| 查看: 6476 | 回复: 3 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
R语言画PCA图，添加椭圆 已有2人参与
|
|||
|
大家好; 最近需要使用R语言画PCA图，而且要用椭圆把不同分组圈起来，请问椭圆怎么加啊。下面是看得网上的程序，但是不对，ellipse是加椭圆的命令，但是我画的图没有椭圆，是不是需要加载什么包啊。求教。 group = read.table("group.txt", header=F, row.names= 1, sep="\t"  otu_table_profile.sort = read.table("otu_table_profile.sort.txt", row.names= 1, header=T, sep="\t" idx = rownames(group) %in% colnames(otu_table_profile.sort) sub_group = group[idx,] ##idx组名 count = otu_table_profile.sort[, rownames(sub_group)] otu.pca <- prcomp(t(count), scale. = TRUE) ggbiplot(otu.pca, obs.scale = 1, var.scale = 1,groups = group$genotype, ellipse = TRUE,ellipse.prob = 0.68,labels =NULL, labels.size = 3, alpha = 1,var.axes = F) |
回复此楼» 猜你喜欢
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
4楼2018-09-24 08:38:57
deephill
铁杆木虫 (职业作家)
- 应助: 59 (初中生)
- 金币: 5262.3
- 散金: 100
- 红花: 16
- 帖子: 4425
- 在线: 425小时
- 虫号: 305318
- 注册: 2006-12-09
- 性别: GG
- 专业: 微量元素地球化学
2楼2018-06-24 22:21:30
3楼2018-06-25 09:16:20
