角座標・多次元球・カイ自乗分布・非心カイ自乗分布・αエラー・パワーなど library(sphere) ThetaAngle<-function(cos,sin){ t<-acos(cos) t[which(sin<0)]<-2*pi-t[which(sin<0)] t } Simplex<-function(Ts,a,r){ A<-Ts%*%t(a) R<-Ts%*%r L<-sqrt(A^2+R^2) …

geomSeries<-function(stval,endval,n){ exp(seq(from=log(stval),to=log(endval),length.out=n)) }> stval<-runif(1) > endval<-runif(1) > n<-10 > gsout<-geomSeries(stval,endval,n) > for( i in 1:(length(gsout)-1)){ print(gsout[i]/gsout[i+1])} [1]…

こちらがＰＤＦへのリンク先 こちらが紹介記事

２つのベクトルがあるとする。 その線形結合としてのようにパラメタ表示することがある この場合は、tの値によって、を結んだ直線上に「等間隔」で点を取ることができる 今回は、角座標表示で考えて、がなす角から、「角度として等間隔」で点をとりたい 掲載…