周期データのパラメタ化

• こちらの続き
• k次元の複素数ベクトルのノルムを１にするために、なるを作りたい
• データをうまく説明するように変数を回帰推定するためには、うまく取り扱いたい
• 多次元極座標を用いて行うことにすれば、以下のように。。。

k<-3 # 次元
v<-runif(k)
sphereCoords<-function(v){
	C<-cos(v)
	S<-sin(v)
	ret<-c(1,S)

	for(i in 1:length(v)){
		for(j in i:length(v)){
			ret[i]<-ret[i]*C[j]
		}
	}
	return(ret)
}
ret<-sphereCoords(v)
# 検算
sum(ret^2)