Take Your Time

仕事や研究、コンピューターとの付き合い方

2013-11-30から1日間の記事一覧

Goodness-of-fitの取り出し方

R

lmの場合、summary(fit)$r.sqで決定係数、summary(fit)$adj.rで修正済み決定係数。観察数はnobs(fit)。 glmの場合、summary(fit)["null.deviance"]でNull Deviance、summary(fit)["deviance"]でResidual Deviance これとcreateTexregを合わせると相当自由に…

交差検定

R

英:Cross-validation パッケージcvToolsが使える。 線形モデルの場合、 library("cvTools") fit cv y = fit$model["y"], cost = mspe, K = nrow(fit$model)) これでone-leave-out cross-validationができる。 Kに他の数字を入れるとk-fold cross-validation…