プロビット回帰分析
プロビット回帰分析の場合
glsyi = [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
glsxij = [
[1, 24],
[1, 18],
[0, 15],
[1, 16],
[0, 10],
[1, 26],
[1, 2],
[0, 24],
[1, 18],
[1, 22],
[1, 3],
[1, 6],
[0, 15],
[0, 12],
[1, 6],
[0, 6],
[1, 12],
[0, 12],
[1, 18],
[1, 3],
[1, 8],
[0, 9],
[0, 12],
[0, 6],
[0, 8],
[1, 12],
]
regana = Num4RegAnaLib::ProBitRegAnaLib.new
regana.non_line_reg_ana(glsyi, glsxij)
AIC(赤池の情報基準)の場合
reg = {
:intercept=> -5.0497, # 定数項
:slope=> [2.2379, 0.2973], # 回帰係数
}
xij = [
[1, 24],
[1, 18],
[0, 15],
[1, 16],
[0, 10],
[1, 26],
[1, 2],
[0, 24],
[1, 18],
[1, 22],
[1, 3],
[1, 6],
[0, 15],
[0, 12],
[1, 6],
[0, 6],
[1, 12],
[0, 12],
[1, 18],
[1, 3],
[1, 8],
[0, 9],
[0, 12],
[0, 6],
[0, 8],
[1, 12],
]
regana = Num4RegAnaLib::ProBitRegAnaLib.new
regana.get_aic(reg, xij)