{"id":160,"date":"2015-09-05T20:43:00","date_gmt":"2015-09-05T11:43:00","guid":{"rendered":"http:\/\/arat.xyz\/wordpress\/?p=160"},"modified":"2017-01-03T11:06:09","modified_gmt":"2017-01-03T02:06:09","slug":"%e3%82%ab%e3%83%86%e3%82%b4%e3%83%aa%e3%82%ab%e3%83%ab%e3%83%87%e3%83%bc%e3%82%bf%e3%81%8b%e3%82%89%e3%83%80%e3%83%9f%e3%83%bc%e5%a4%89%e6%95%b0%e3%82%92%e4%bd%9c%e6%88%90%e3%81%99%e3%82%8br-2","status":"publish","type":"post","link":"https:\/\/www.arat.xyz\/wordpress\/?p=160","title":{"rendered":"\u30ab\u30c6\u30b4\u30ea\u30ab\u30eb\u30c7\u30fc\u30bf\u304b\u3089\u30c0\u30df\u30fc\u5909\u6570\u3092\u4f5c\u6210\u3059\u308bR\u95a2\u6570\u3092\u81ea\u4f5c\uff08\u8907\u6570\u5217\u5bfe\u5fdc\u53ef\uff09"},"content":{"rendered":"
\n

\u8ffd\u8a18\uff082017.1.3\uff09<\/h2>\n
\n

makedummies\u30d1\u30c3\u30b1\u30fc\u30b8\u304cCRAN\u306b\u767b\u9332\u3055\u308c\u307e\u3057\u305f\u3002\u30d6\u30ed\u30b0\u306e \u3053\u306e\u30da\u30fc\u30b8<\/a> \u3092\u53c2\u7167\u3057\u3066\u4e0b\u3055\u3044\u3002 <\/p>\n<\/div>\n<\/div>\n

\n

\u8ffd\u8a18\uff082015.12.26\uff09<\/h2>\n
\n

GitHub\u3067\u95a2\u6570\u3092\u516c\u958b\u3057\u307e\u3057\u305f\u3002\u30d6\u30ed\u30b0\u306e \u3053\u306e\u30da\u30fc\u30b8<\/a> \u3092\u53c2\u7167\u3057\u3066\u4e0b\u3055\u3044\u3002 <\/p>\n<\/div>\n<\/div>\n

\n

\u7d4c\u7def<\/h2>\n
\n

R\u306eglm\u95a2\u6570\u3092\u7528\u3044\u3066\u91cd\u56de\u5e30\u5206\u6790\u3092\u884c\u3046\u969b\u306b\u306f\u3001\u30ab\u30c6\u30b4\u30ea\u30ab\u30eb\u30c7\u30fc\u30bf\u3092\u81ea\u52d5\u7684\u306b\u30c0\u30df\u30fc\u5909\u6570\u306b\u5909\u63db\u3057\u3066\u304f\u308c\u308b\u305f\u3081\u3001\u3042\u307e\u308a\u82e6\u52b4\u3057\u307e\u305b\u3093\u3002 <\/p>\n

\u73fe\u5728\u3001\u500b\u4eba\u7684\u306bRStan\u3092\u52c9\u5f37\u4e2d\u3067\u3042\u308a\u3001\u540c\u3058\u89e3\u6790\u3092\u884c\u3046\u305f\u3081\u306b\u306f\u81ea\u5206\u3067\u30c0\u30df\u30fc\u5909\u6570\u3092\u4f5c\u6210\u3059\u308b\u5fc5\u8981\u304c\u3042\u308a\u307e\u3059\u3002\u3055\u3089\u306b\u30ab\u30c6\u30b4\u30ea\u30ab\u30eb\u30c7\u30fc\u30bf\u3068\u6570\u5024\u30c7\u30fc\u30bf\u304c\u6df7\u5728\u3057\u3066\u3044\u308b\u3068\u305d\u308c\u3060\u3051\u3067\u975e\u5e38\u306b\u5927\u5909\u3067\u3059\u3002 <\/p>\n

\u305d\u3053\u3067\u7fa4\u99ac\u5927\u5b66\u306e\u9752\u6728\u5148\u751f\u304c\u516c\u958b\u3055\u308c\u3066\u3044\u308bmake.dummy\u95a2\u6570\uff08 http:\/\/aoki2.si.gunma-u.ac.jp\/R\/src\/make.dummy.R<\/a> \uff09\u3092\u3082\u3068\u306b\u3057\u3066\u95a2\u6570\u3092\u4f5c\u6210\u3057\u307e\u3057\u305f\u3002 <\/p>\n<\/div>\n<\/div>\n

\n

\u4f5c\u6210\u3057\u305f\u95a2\u6570<\/h2>\n
\n

\uff082015.9.6\u8ffd\u8a18\uff09\u4ee5\u4e0b\u306e\u30d7\u30ed\u30b0\u30e9\u30e0\u3092\u4fee\u6b63\u3057\u305f\u3082\u306e\u3092\u8a18\u4e8b\u306e\u4e0b\u306b\u63b2\u8f09\u3057\u3066\u3044\u307e\u3059\u3002\u3053\u3061\u3089\u3092\u4f7f\u7528\u305b\u305a\u306b\u305d\u3061\u3089\u3092\u53c2\u8003\u306b\u3057\u3066\u4e0b\u3055\u3044\u3002 <\/p>\n

\n
# <\/span>http:\/\/aoki2.si.gunma-u.ac.jp\/R\/src\/make.dummy.R<\/span>\nmake.dummy<\/span> <-<\/span> function<\/span>(dat, basal_level= FALSE<\/span>, sep = \"_\"<\/span>) {\n    name <-<\/span> colnames(dat)\n    level <-<\/span> levels(dat[,1])\n    if<\/span> (!is.data.frame(dat))\n        dat <-<\/span> as.data.frame(dat)\n    ncat <-<\/span> ncol(dat)\n    dat[, 1:ncat] <-<\/span> lapply(dat, function<\/span>(x) {\n        if<\/span> (is.factor(x)) {\n            return<\/span>(as.integer(x))\n        } else<\/span> {\n            return<\/span>(x)\n        }\n    })\n    mx <-<\/span> sapply(dat, max)\n    start <-<\/span> c(0, cumsum(mx)[1:(ncat-1)])\n    nobe <-<\/span> sum(mx)\n    ## <\/span>modified<\/span>\n    res <-<\/span> t(apply(dat, 1, function<\/span>(obs) 1:nobe %in%<\/span> (start+obs))) + 0\n    colnames(res) <-<\/span> paste(name, level, sep = sep)\n    if<\/span> (basal_level == FALSE<\/span>) res <-<\/span> res[,-1]\n    return<\/span>(res)\n}\n\n\nmake.dummys<\/span> <-<\/span> function<\/span>(dat, ...) {\n    n <-<\/span> ncol(dat)\n    res_list <-<\/span> lapply(seq(n), function<\/span>(i) {\n        tmp <-<\/span> as.data.frame(dat[,i])\n        colnames(tmp) <-<\/span> colnames(dat)[i]\n        if<\/span> (is.factor(dat[,i])) {       # <\/span>factor or ordered<\/span>\n            make.dummy(tmp, ...)\n        } else<\/span> {\n            tmp\n        }\n    })\n\n    res <-<\/span> NULL<\/span>\n    for<\/span> (i in<\/span> seq(n)) {\n        res <-<\/span> cbind(res, res_list[[i]])\n    }\n    return<\/span>(res)\n}\n<\/pre>\n<\/div>\n<\/div>\n
\n
\u89e3\u8aac<\/h3>\n
\n<\/div>\n
\n
make.dummy\u95a2\u6570<\/h4>\n
\n
 \u9752\u6728\u5148\u751f\u304c\u516c\u958b\u3055\u308c\u3066\u3044\u308bmake.dummy\u95a2\u6570\u3092\u5c11\u3057\u3060\u3051\u6539\u5909\u3057\u307e\u3057\u305f\u3002\u30c7\u30fc\u30bf\u306f\u30c7\u30fc\u30bf\u30d5\u30ec\u30fc\u30e0\u5f62\u5f0f\u3067\u6e21\u3057\u307e\u3059\u3002 <\/p>\n
 \u5909\u66f4\u70b9\u306f\u4ee5\u4e0b\u306e\u70b9\u3067\u3059\u3002 <\/p>\n
    \n
  1. \u57fa\u6e96\u3068\u306a\u308b\u30ab\u30c6\u30b4\u30ea\u30fc\u3092\u524a\u9664\u3067\u304d\u308b\u3088\u3046\u306b\u5909\u66f4\n