2012-08-12から1日間の記事一覧

メモ

CategoryVector<-function (d = 3) { if(d == 1){ return(matrix(1,ncol=1,nrow=1)) } df <- d - 1 diagval <- 1:d diagval <- sqrt((df + 1)/df) * sqrt((df - diagval + 1)/(df - diagval + 2)) others <- -diagval/(df - (0:(d - 1))) m <- matrix(rep(ot…