閉じる
閉じる
×
【追記】
今ニコ生で放送中なので、本記事について放送内で質問していただければ可能な限り答えます。
なんかやる-ニコニコミュニティ
回帰分析を用いて、大型艦建造における建造率の推定を行いました。サンプルとしてしたらばにある艦これ専用掲示板の【報告】大型艦/新型艦建造スレ25【検証】 - 艦これ専用掲示板にて報告された39855件のデータを使用しました。
したらばにて報告のあった1187種類のレシピを、確率が大きく変わらないと推定された、いくつかのグループにまとめ、その中でもっとも資源投入量の少ないレシピの結果として併合しています。レシピごとの報告数が少なくても、グループごとに報告をまとめることで精度を高めています。見たいレシピが表に載っていない場合は、それよりも投入量の少ないレシピを参照してください。
表中にある95%CIとは95%信頼区間の略で、実際の建造率は95%の確率でこの区間の中に納まっていると推定されたことを表しています。
今回の推定では、開発資材・空きドックともに建造率への影響は見られませんでした。しかし、今後さらに細かく推定していけば、開発資材・空きドックの影響が見られるようになる可能性は残されています。
大鳳
大鳳は3600/2000/5300/5200が最も出やすく、大鳳率11.06%と推定されました。ボーキサイトを7000まで盛ると、確率はかえって下がると推定されました。それ以外では、資源をこれ以上持っても確率は変わらないと推定されました。なお、このレシピでのあきつ丸率は7.27%と推定されました。
大鳳/ 計 確率 (95%CI 下限-上限)
後空0/ 資材 1/ 3000/1500/3500/1000 "大鳳 1/ 7359 0.00% ( 0.00%- 0.01%)"
後空0/ 資材 1/ 3000/1500/3500/3000 "大鳳 32/ 6208 0.59% ( 0.43%- 0.80%)"
後空0/ 資材 1/ 3000/1500/3500/3200 "大鳳 70/ 1625 4.25% ( 3.46%- 5.20%)"
後空0/ 資材 1/ 3000/2000/5000/7000 "大鳳 0/ 2 2.65% ( 1.95%- 3.60%)"
後空0/ 資材 1/ 3500/2000/4100/5200 "大鳳 174/ 3643 4.70% ( 4.09%- 5.41%)"
後空0/ 資材 1/ 3500/2000/5000/5000 "大鳳 22/ 317 6.11% ( 4.85%- 7.68%)"
後空0/ 資材 1/ 3500/2000/5000/5200 "大鳳 896/13296 6.75% ( 6.35%- 7.19%)"
後空0/ 資材 1/ 3500/2000/5000/7000 "大鳳 75/ 1771 4.26% ( 3.45%- 5.26%)"
後空0/ 資材 1/ 3600/2000/5300/5200 "大鳳 74/ 677 11.06% ( 9.00%- 13.52%)"
後空0/ 資材 1/ 3600/2000/5300/7000 "大鳳 9/ 114 7.10% ( 5.33%- 9.41%)"
大和
大和は4010/5000/6000/2500が最も出やすく、大和率13.25%と推定されました。また、資源をこれ以上増やしても、確率は変わらないと推定されました。しかし、燃料4010の報告件数が少ないため、燃料の境目が甘く出ている可能性があります。余裕をもって、報告が多くあがっている燃料4500程まで増やした方がいいかもしれません。なお、このレシピでのあきつ丸率は3.56%と推定されました。
大和/ 計 確率 (95%CI 下限-上限)
後空0/ 資材 1/ 3000/4000/4000/1000 "大和 0/ 90 0.22% ( 0.08%- 0.63%)"
後空0/ 資材 1/ 3000/4000/4000/1500 "大和 64/ 1484 4.39% ( 3.46%- 5.55%)"
後空0/ 資材 1/ 3500/4000/4000/2000 "大和 296/ 6198 4.70% ( 4.21%- 5.25%)"
後空0/ 資材 1/ 4000/5000/6000/2500 "大和 272/ 3811 7.22% ( 6.46%- 8.07%)"
後空0/ 資材 1/ 4010/5000/6000/1000 "大和 2/ 161 0.43% ( 0.15%- 1.26%)"
後空0/ 資材 1/ 4010/5000/6000/1500 "大和 0/ 0 8.27% ( 5.78%- 11.70%)"
後空0/ 資材 1/ 4010/5000/6000/2000 "大和 20/ 277 8.83% ( 6.71%- 11.55%)"
後空0/ 資材 1/ 4010/5000/6000/2500 "大和 30/ 209 13.25% ( 10.20%- 17.05%)"
あきつ丸
あきつ丸は4000/2000/6000/7000が最も出やすくあきつ丸率17.67%、3500/4500/4000/4500が次に出やすくあきつ丸率14.27%と推定されました。しかしどちらも報告件数が少ないため、推定が甘く出ている可能性があります。報告件数が十分多い中では3500/4500/4000/2000がもっとも出やすく9.41%、3500/2000/4000/4500が次に多く7.27%と推定されました。
あきつ丸/ 計 確率 (95%CI 下限-上限)
後空0/ 資材 1/ 3000/1500/3500/1000 "あきつ丸 11/ 740 1.49% ( 0.83%- 2.66%)"
後空0/ 資材 1/ 3000/2000/4000/2000 "あきつ丸 12/ 776 3.04% ( 2.38%- 3.87%)"
後空0/ 資材 1/ 3000/2000/4000/4500 "あきつ丸 9/ 171 4.79% ( 3.79%- 6.03%)"
後空0/ 資材 1/ 3000/4500/4000/2000 "あきつ丸 94/ 1340 6.24% ( 5.19%- 7.50%)"
後空0/ 資材 1/ 3000/4500/4000/4500 "あきつ丸 0/ 1 9.64% ( 7.66%- 12.08%)"
後空0/ 資材 1/ 3000/4500/5500/2000 "あきつ丸 1/ 19 2.31% ( 1.77%- 3.01%)"
後空0/ 資材 1/ 3000/4500/5500/4500 "あきつ丸 0/ 0 3.65% ( 2.66%- 4.99%)"
後空0/ 資材 1/ 3500/2000/4000/2000 "あきつ丸 215/ 4363 4.67% ( 4.11%- 5.30%)"
後空0/ 資材 1/ 3500/2000/4000/4500 "あきつ丸 1444/19891 7.27% ( 6.92%- 7.64%)"
後空0/ 資材 1/ 3500/4500/4000/2000 "あきつ丸 67/ 825 9.41% ( 7.84%- 11.27%)"
後空0/ 資材 1/ 3500/4500/4000/4500 "あきつ丸 6/ 39 14.27% ( 11.68%- 17.33%)"
後空0/ 資材 1/ 3500/4500/5500/2000 "あきつ丸 239/ 6719 3.56% ( 3.14%- 4.02%)"
後空0/ 資材 1/ 3500/4500/5500/4500 "あきつ丸 2/ 14 5.58% ( 4.66%- 6.67%)"
後空0/ 資材 1/ 4000/2000/6000/7000 "あきつ丸 7/ 30 17.67% ( 11.20%- 26.75%)"
後空0/ 資材 1/ 4000/4500/6000/7000 "あきつ丸 10/ 84 13.93% ( 8.79%- 21.38%)"
まるゆ
まるゆは1500/1500/2000/1000のいわゆる初期値レシピが最も出やすく、まるゆ率15.01%と推定されました。
まるゆ/ 計 確率 (95%CI 下限-上限)
後空0/ 資材 1/ 1500/1500/2000/1000 "まるゆ 546/ 3578 15.01% ( 13.92%- 16.17%)"
後空0/ 資材 1/ 3000/1500/4000/2000 "まるゆ 1784/17613 10.18% ( 9.74%- 10.63%)"
後空0/ 資材 1/ 3200/3700/4200/2000 "まるゆ 0/ 2 9.28% ( 8.08%- 10.63%)"
後空0/ 資材 1/ 3500/3000/3500/1000 "まるゆ 36/ 343 13.10% ( 11.74%- 14.58%)"
後空0/ 資材 1/ 3500/3000/4000/2000 "まるゆ 560/ 6252 8.82% ( 8.16%- 9.52%)"
後空0/ 資材 1/ 3500/3700/4200/2000 "まるゆ 306/ 3811 8.02% ( 7.20%- 8.93%)"
後空0/ 資材 1/ 4000/4500/5010/2000 "まるゆ 403/ 6802 5.92% ( 5.39%- 6.51%)"
詳細
今回は変数選択基準をBIC(ベイズ情報量基準)に戻しました。その上で、偏相関係数の2乗が0.9以上となる説明変数を除き、多重共線性を防ぎました。
最後に、ソースコード等を載せておきます。
$大鳳
Call:
glm(formula = cbind(大鳳, 大和 + あきつ丸 + まるゆ + 加賀 + 瑞鶴 +
赤城 + 蒼龍 + 飛龍 + 翔鶴 + 祥鳳 + 瑞鳳 + 隼鷹 + 飛鷹 + 龍驤 +
阿賀野 + 能代 + 矢矧 + 熊野 + 最上 + 三隈 + 筑摩 + 利根 +
鈴谷 + 伊勢 + 金剛 + 山城 + 榛名 + 長門 + 日向 + 比叡 + 扶桑 +
霧島 + 陸奥) ~ `3500/2000/4100/5200` + `2000/1500/3000/3200` +
`3000/1500/3000/3000` + `3000/1500/3500/1000` + `1500/1500/2000/3000` +
`3500/2000/5000/5000` + `2500/2000/5000/7000` + `3600/2000/5300/5200`,
family = binomial, data = data)
Deviance Residuals:
Min 1Q Median 3Q Max
-2.4063 -0.3740 -0.1086 -0.0040 3.6245
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -8.14495 1.00033 -8.142 3.88e-16 ***
`3500/2000/4100/5200` 0.10695 0.12308 0.869 0.384863
`2000/1500/3000/3200` 2.01499 0.19041 10.582 < 2e-16 ***
`3000/1500/3000/3000` 2.07254 0.58740 3.528 0.000418 ***
`3000/1500/3500/1000` -3.60095 0.36078 -9.981 < 2e-16 ***
`1500/1500/2000/3000` 4.54288 1.11420 4.077 4.56e-05 ***
`3500/2000/5000/5000` 0.38355 0.08132 4.716 2.40e-06 ***
`2500/2000/5000/7000` -0.48645 0.11646 -4.177 2.96e-05 ***
`3600/2000/5300/5200` 0.54072 0.12074 4.478 7.53e-06 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for binomial family taken to be 1)
Null deviance: 1971.0 on 1187 degrees of freedom
Residual deviance: 517.9 on 1179 degrees of freedom
AIC: 883.32
Number of Fisher Scoring iterations: 11
$大和
Call:
glm(formula = cbind(大和, 大鳳 + あきつ丸 + まるゆ + 加賀 + 瑞鶴 +
赤城 + 蒼龍 + 飛龍 + 翔鶴 + 祥鳳 + 瑞鳳 + 隼鷹 + 飛鷹 + 龍驤 +
阿賀野 + 能代 + 矢矧 + 熊野 + 最上 + 三隈 + 筑摩 + 利根 +
鈴谷 + 伊勢 + 金剛 + 山城 + 榛名 + 長門 + 日向 + 比叡 + 扶桑 +
霧島 + 陸奥) ~ `3500/4000/4000/2000` + `3000/4000/4000/1000` +
`3600/2000/4100/4200` + `3000/3000/3000/1500` + `4000/5000/6000/2500` +
`4010/5000/6000/1000` + `3500/1500/3500/5200`, family = binomial,
data = data)
Deviance Residuals:
Min 1Q Median 3Q Max
-1.9871 -0.3081 -0.0634 -0.0259 3.2424
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -7.31136 0.50582 -14.454 < 2e-16 ***
`3500/4000/4000/2000` 0.07235 0.13902 0.520 0.60275
`3000/4000/4000/1000` 1.18928 0.26023 4.570 4.87e-06 ***
`3600/2000/4100/4200` -3.10995 1.05996 -2.934 0.00335 **
`3000/3000/3000/1500` 3.04057 0.53218 5.713 1.11e-08 ***
`4000/5000/6000/2500` 0.45575 0.08315 5.481 4.22e-08 ***
`4010/5000/6000/1000` 0.67479 0.15203 4.438 9.06e-06 ***
`3500/1500/3500/5200` -1.38248 0.45924 -3.010 0.00261 **
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for binomial family taken to be 1)
Null deviance: 1859.23 on 1187 degrees of freedom
Residual deviance: 306.02 on 1180 degrees of freedom
AIC: 528.52
Number of Fisher Scoring iterations: 11
$あきつ丸
Call:
glm(formula = cbind(あきつ丸, 大鳳 + 大和 + まるゆ + 加賀 + 瑞鶴 +
赤城 + 蒼龍 + 飛龍 + 翔鶴 + 祥鳳 + 瑞鳳 + 隼鷹 + 飛鷹 + 龍驤 +
阿賀野 + 能代 + 矢矧 + 熊野 + 最上 + 三隈 + 筑摩 + 利根 +
鈴谷 + 伊勢 + 金剛 + 山城 + 榛名 + 長門 + 日向 + 比叡 + 扶桑 +
霧島 + 陸奥) ~ `3000/2000/4000/4500` + `3500/2000/4000/2000` +
`3000/2000/4000/2000` + `3000/1500/3500/1000` + `3000/4500/5500/2000` +
`3000/4500/4000/2000` + `4000/2000/6000/7000`, family = binomial,
data = data)
Deviance Residuals:
Min 1Q Median 3Q Max
-2.8472 -0.5354 -0.3461 -0.1384 3.3598
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -7.02879 0.57761 -12.169 < 2e-16 ***
`3000/2000/4000/4500` 0.47155 0.07291 6.468 9.94e-11 ***
`3500/2000/4000/2000` 0.44472 0.12250 3.630 0.000283 ***
`3000/2000/4000/2000` 0.73180 0.32946 2.221 0.026338 *
`3000/1500/3500/1000` 2.83501 0.65262 4.344 1.40e-05 ***
`3000/4500/5500/2000` -1.03544 0.12119 -8.544 < 2e-16 ***
`3000/4500/4000/2000` 0.75304 0.11445 6.580 4.72e-11 ***
`4000/2000/6000/7000` 1.00680 0.27246 3.695 0.000220 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for binomial family taken to be 1)
Null deviance: 1370.74 on 1187 degrees of freedom
Residual deviance: 764.59 on 1180 degrees of freedom
AIC: 1261.2
Number of Fisher Scoring iterations: 9
-
次の記事これより新しい記事はありません。
-
前の記事2014-02-28 07:52:03[艦これ][大型艦建造] 建造率の推定 (2月28日更新)
広告
コメントコメントを書く
他80件のコメントを表示