Isotonic regressionを扱っている。 制約条件行列の作成で、行列が大きくなると時間がかかるので、for loopをなんとか回避して関数の高速化を行った。 library(quadprog) #for loopの関数 f1 <- function(Mat, multiple=3){ N <- nrow(Mat) # 観測用量点の数…