从品牌网站建设到网络营销策划,从策略到执行的一站式服务
本篇文章给大家分享的是有关如何使用R语言实现对SSR数据做主成分分析,小编觉得挺实用的,因此分享给大家学习,希望大家阅读完这篇文章后可以有所收获,话不多说,跟着小编一起来看看吧。
创新互联建站专注为客户提供全方位的互联网综合服务,包含不限于成都做网站、成都网站建设、成都外贸网站建设、桥东网络推广、小程序开发、桥东网络营销、桥东企业策划、桥东品牌公关、搜索引擎seo、人物专访、企业宣传片、企业代运营等,从售前售中售后,我们都将竭诚为您服务,您的肯定,是我们最大的嘉奖;创新互联建站为所有大学生创业者提供桥东建站搭建服务,24小时服务热线:18980820575,官方网址:www.cdcxhl.com
有朋友留言问如何对0,1矩阵做主成分分析(PCA),查了一下参考资料找到了一个办法
示例数据来自于R语言包poppr
,csv文件存储,数据格式如下
使用到的是R语言的poppr
包中的read.genalex()
函数poppr
第一次使用需要先安装
install.packages('poppr')
读入数据
library(poppr)
mydf<-read.genalex("rootrot.csv")
mydf
mydf1<-genclone2genind(mydf)
读入数据直接是 genclone object,使用函数genclone2genind()
将其转换成genind object,接下来使用ade4
包中的dudi.pca()
函数做主成分分析
mydf1.pca<-dudi.pca(mydf1,scannf = F,nf=2)
mydf1.pca.scores<-mydf1.pca$li
主成分的结果存储在li中
还是认为的分个组,然后做散点图
mydf1.pca.scores$Pop<-ifelse(mydf1.pca.scores$Axis1>0,"pop1","pop2")
library(ggplot2)
ggplot(mydf1.pca.scores,aes(x=Axis1,y=Axis2))+
geom_point(aes(color=Pop))+
stat_ellipse(aes(x=Axis1,y=Axis2,fill=Pop),
geom="polygon",
alpha=0.2,
lty="dashed",
color="black")+
theme_bw()
以上是第一种方法,刚刚看到昨天提到的参考资料,原来和昨天的方法一样也可以读入数据
以上就是如何使用R语言实现对SSR数据做主成分分析,小编相信有部分知识点可能是我们日常工作会见到或用到的。希望你能通过这篇文章学到更多知识。更多详情敬请关注创新互联行业资讯频道。
成都网站建设公司地址:成都市青羊区太升南路288号锦天国际A座10层 建设咨询028-86922220
成都快上网科技有限公司-四川网站建设设计公司 | 蜀ICP备19037934号 Copyright 2020,ALL Rights Reserved cdkjz.cn | 成都网站建设 | © Copyright 2020版权所有.
专家团队为您提供成都网站建设,成都网站设计,成都品牌网站设计,成都营销型网站制作等服务,成都建网站就找快上网! | 成都网站建设哪家好? | 网站建设地图