臨床統計もおもしろいですよ、その3

18 ：卵の名無しさん：2020/03/09(月) 19:01:23.07 ID:Ost2Ton5.net

Sys.setenv(lang='en')
rm(list=ls())

library(gtools)
n=4
a=permutations(n,n) # ｎの順列
r=nrow(a)
f<-function(x){ # x=c(1,2,3) -> rbind(a[1],a[2],a[3])
n=length(x)
ans=NULL
for(i in 1:n){
ans=rbind(ans,a[i,])
}
return(ans)
}
(b=combn(r,n,f)) # すべての行列の配列 b[,,i]
(c=dim(b)[3])　 # その個数

diag2 <- function(x){ # 右上から左下への対角線の配列を返す
n=nrow(x)
ans=numeric(n)
for(i in 1:n) ans[i]=x[i,n+1-i]
return(ans)
}