MikuHatsune https://blog.hatena.ne.jp/MikuHatsune/ 驚異のアニヲタ社会復帰の予備 https://mikuhatsune.hatenadiary.com/ R Rを使いこなす Isotonic regressionを扱っている。 制約条件行列の作成で、行列が大きくなると時間がかかるので、for loopをなんとか回避して関数の高速化を行った。 library(quadprog) #for loopの関数 f1 <- function(Mat, multiple=3){ N <- nrow(Mat) # 観測用量点の数 M <- ncol(Mat) # 用量反応曲線の数 # デフォルトはtriplicateとする ws <- rep(multiple, N*M) Dmat <- diag(ws) dvec <- ws*c(Mat) Q <- NULL for(j in… 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Fmikuhatsune.hatenadiary.com%2Fentry%2F20120325%2F1332604299" title="Isotonic regressionの高速化 - 驚異のアニヲタ社会復帰の予備" class="embed-card embed-blogcard" scrolling="no" frameborder="0" style="display: block; width: 100%; height: 190px; max-width: 500px; margin: 10px 0px;"></iframe> https://cdn-ak.d.st-hatena.com/diary/MikuHatsune/2012-03-25.png Hatena Blog https://hatena.blog 2012-03-25 00:51:39 rich https://mikuhatsune.hatenadiary.com/entry/20120325/1332604299 1.0 100%