こちら simplex <- function(n) { qr.Q(qr(matrix(1,nrow=n)),complete=T)[,-1] } QR分解は行列を直交行列と上三角行列の積に分解する 今、行数n、列数1の行列をQR分解すると、nxn直交行列と、nx1上三角行列に分解される nx1上三角行列とは、１列の行列で、…