2010-04-20から1日間の記事一覧

等分するベクトル

この頃の話題(記事はこちら)のベクトルをRで作る カテゴリ数nc(自由度=次元=df) 実行はこちらでも CategoryVector<-function(nc=3){ df<-nc-1 d<-df+1 diagval<-1:d diagval<-sqrt((df+1)/df)*sqrt((df-diagval+1)/(df-diagval+2)) others<--diagval/(df-(0:…