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

pythonをRstudioで使うためにはRnotebookを使わないとだめかと思いましたが(こちら)、Rmarkdownで行けるようです こちらを真似してやってみる --- title: "pythontest" author: "ryamada" date: "2016年10月15日" output: html_document: default pdf_docume…