{"categories":["R","R\u3092\u4f7f\u3044\u3053\u306a\u3059"],"description":"Isotonic regression\u3092\u6271\u3063\u3066\u3044\u308b\u3002 \u5236\u7d04\u6761\u4ef6\u884c\u5217\u306e\u4f5c\u6210\u3067\u3001\u884c\u5217\u304c\u5927\u304d\u304f\u306a\u308b\u3068\u6642\u9593\u304c\u304b\u304b\u308b\u306e\u3067\u3001for loop\u3092\u306a\u3093\u3068\u304b\u56de\u907f\u3057\u3066\u95a2\u6570\u306e\u9ad8\u901f\u5316\u3092\u884c\u3063\u305f\u3002 library(quadprog) #for loop\u306e\u95a2\u6570 f1 <- function(Mat, multiple=3){ N <- nrow(Mat) # \u89b3\u6e2c\u7528\u91cf\u70b9\u306e\u6570 M <- ncol(Mat) # \u7528\u91cf\u53cd\u5fdc\u66f2\u7dda\u306e\u6570 # \u30c7\u30d5\u30a9\u30eb\u30c8\u306ftriplicate\u3068\u3059\u308b ws <- rep(multiple, N*M) Dmat <- diag(ws) dvec <- ws*c(Mat) Q <- NULL for(j in\u2026","provider_name":"Hatena Blog","height":"190","author_url":"https://blog.hatena.ne.jp/MikuHatsune/","provider_url":"https://hatena.blog","version":"1.0","title":"Isotonic regression\u306e\u9ad8\u901f\u5316","image_url":"https://cdn-ak.d.st-hatena.com/diary/MikuHatsune/2012-03-25.png","author_name":"MikuHatsune","type":"rich","published":"2012-03-25 00:51:39","url":"https://mikuhatsune.hatenadiary.com/entry/20120325/1332604299","blog_url":"https://mikuhatsune.hatenadiary.com/","html":"<iframe src=\"https://hatenablog-parts.com/embed?url=https%3A%2F%2Fmikuhatsune.hatenadiary.com%2Fentry%2F20120325%2F1332604299\" title=\"Isotonic regression\u306e\u9ad8\u901f\u5316 - \u9a5a\u7570\u306e\u30a2\u30cb\u30f2\u30bf\u793e\u4f1a\u5fa9\u5e30\u306e\u4e88\u5099\" class=\"embed-card embed-blogcard\" scrolling=\"no\" frameborder=\"0\" style=\"display: block; width: 100%; height: 190px; max-width: 500px; margin: 10px 0px;\"></iframe>","width":"100%","blog_title":"\u9a5a\u7570\u306e\u30a2\u30cb\u30f2\u30bf\u793e\u4f1a\u5fa9\u5e30\u306e\u4e88\u5099"}