溫馨提示×

溫馨提示×

您好，登錄后才能下訂單哦！

密碼登錄×
登錄注冊×
其他方式登錄
點擊 登錄注冊 即表示同意《億速云用戶服務(wù)條款》

如何使用R語言實現(xiàn)對SSR數(shù)據(jù)做主成分分析

發(fā)布時間：2021-11-22 15:50:47 來源：億速云 閱讀：761 作者：柒染 欄目：大數(shù)據(jù)

本篇文章給大家分享的是有關(guān)如何使用R語言實現(xiàn)對SSR數(shù)據(jù)做主成分分析，小編覺得挺實用的，因此分享給大家學(xué)習(xí)，希望大家閱讀完這篇文章后可以有所收獲，話不多說，跟著小編一起來看看吧。

有朋友留言問如何對0,1矩陣做主成分分析（PCA），查了一下參考資料找到了一個辦法

 首先是準(zhǔn)備數(shù)據(jù)

示例數(shù)據(jù)來自于R語言包poppr,csv文件存儲，數(shù)據(jù)格式如下

如何使用R語言實現(xiàn)對SSR數(shù)據(jù)做主成分分析  
image.png
  • 1 位點總數(shù)
  • 2 樣本總數(shù)
  • 3 群體總數(shù)
  • 4、5、6 每個群體中 的樣本數(shù)以及群體名
  • 7 是樣本編號
  • 8 是群體名稱
  • 9 是位點編號
  • 10 是0,1型的SSR數(shù)據(jù)
 R語言讀入數(shù)據(jù)

使用到的是R語言的poppr包中的read.genalex()函數(shù)poppr第一次使用需要先安裝

install.packages('poppr')
 

讀入數(shù)據(jù)

library(poppr)
mydf<-read.genalex("rootrot.csv")
mydf
mydf1<-genclone2genind(mydf)
 

讀入數(shù)據(jù)直接是 genclone object，使用函數(shù)genclone2genind()將其轉(zhuǎn)換成genind object，接下來使用ade4包中的dudi.pca()函數(shù)做主成分分析

mydf1.pca<-dudi.pca(mydf1,scannf = F,nf=2)
mydf1.pca.scores<-mydf1.pca$li
 

主成分的結(jié)果存儲在li中

還是認(rèn)為的分個組，然后做散點圖

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語言實現(xiàn)對SSR數(shù)據(jù)做主成分分析  
image.png

以上是第一種方法，剛剛看到昨天提到的參考資料，原來和昨天的方法一樣也可以讀入數(shù)據(jù)

如何使用R語言實現(xiàn)對SSR數(shù)據(jù)做主成分分析  
image.png
如何使用R語言實現(xiàn)對SSR數(shù)據(jù)做主成分分析  

以上就是如何使用R語言實現(xiàn)對SSR數(shù)據(jù)做主成分分析，小編相信有部分知識點可能是我們?nèi)粘９ぷ鲿姷交蛴玫降?。希望你能通過這篇文章學(xué)到更多知識。更多詳情敬請關(guān)注億速云行業(yè)資訊頻道。

向AI問一下細(xì)節(jié)

免責(zé)聲明：本站發(fā)布的內(nèi)容（圖片、視頻和文字）以原創(chuàng)、轉(zhuǎn)載和分享為主，文章觀點不代表本網(wǎng)站立場，如果涉及侵權(quán)請聯(lián)系站長郵箱：is@yisu.com進(jìn)行舉報，并提供相關(guān)證據(jù)，一經(jīng)查實，將立刻刪除涉嫌侵權(quán)內(nèi)容。

AI
助
手