ryamada22 https://blog.hatena.ne.jp/ryamada22/ ryamadaの遺伝学・遺伝統計学メモ https://ryamada22.hatenablog.jp/ 特異値分解 固有値分解 PCA R 参考こちら が特異値分解。 変形して 。 これを解いて、S,V->Uが得られる。 今、Xを中心化すると はXの分散共分散行列に比例した値になるので、中心かした特異値分解と固有値分解は、同じようなもの。 #構造化集団をシミュレート Nm<-1000 #マーカー数 Npop<-4 #亜集団数 Ns<-c(100,150,200,250) #集団別人数 M<-NULL #全ジェノタイプデータを納める行列 #亜集団別にアレル頻度を振ってシミュレーション for(j in 1:Npop){ tmpM<-matrix(rep(0,Nm*Ns[j]),nrow=Nm) for(i in 1:Nm){ af<… 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Fryamada22.hatenablog.jp%2Fentry%2F20100223%2F1266883282" title="中心化後特異値分解と固有値分解 - ryamadaの遺伝学・遺伝統計学メモ" class="embed-card embed-blogcard" scrolling="no" frameborder="0" style="display: block; width: 100%; height: 190px; max-width: 500px; margin: 10px 0px;"></iframe> Hatena Blog https://hatena.blog 2010-02-23 09:01:22 rich https://ryamada22.hatenablog.jp/entry/20100223/1266883282 1.0 100%