Rのkde2d関数の多次元版…スケーリングとか、bandwidth.nrd関数からの出力の補正とかの点でちょっと怪しい… kdeNd<-function (x,h, n = 25) { max<-apply(x,2,max) min<-apply(x,2,min) nx <- ncol(x)#No. axis ns <- nrow(x)#No. samples g<-mapply(seq.int,…

data<-read.table("data.txt",sep="\t") #多次元ベクトルデータの読み込み(たとえば400人分、数百マーカー分) Cdata<-cov(data,use="pairwise") # covariance matrix Cdatacor<-cov2cor(Cdata) # correlation matrixへの変換(ただしエラーが･･･)?? answer<-p…

多次元ベクトルデータを持つサンプルがある サンプル間の関係を相関係数で表した、相関係数行列(サンプル数xサンプル数の正方対象行列、対角成分は1)ができる この相関係数行列は、サンプル数次元の単位超球面上の点を表現しているとみなせる(A Spherical Re…