モデルの結果の種類

Rで学ぶGeneralized Additive Models (GAMs) による非線形モデリング

Noam Ross

Senior Research Scientist, EcoHealth Alliance

結果の種類

連続的な結果

バイクの速度（mph）
車の燃費（mpg）
土壌中の汚染物質の濃度（g/kg）

二値的な結果

ある場所における生物の有無
購入の有無
アンケートのはい/いいえの回答

確率とログオッズ：ロジスティック関数

確率とログオッズ：ロジット関数

Rにおけるロジスティック関数とロジット関数

plogis() # Logistic
qlogis() # Logit

qlogis(plogis(0.5))

0.5

qlogis(0.25) == log(1/3)

TRUE

mgcvによるロジスティックGAM

gam(y ~ x1 + s(x2),
   data = dat,
   family = binomial,
   method = "REML")

Family: binomial 
Link function: logit 
Formula:
y ~ s(x1) + s(x2)
Parametric coefficients:
            Estimate Std. Error z value Pr(>|z|)    
(Intercept)   0.7330     0.1208    6.07 1.28e-09 ***
 ---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Approximate significance of smooth terms:
        edf Ref.df Chi.sq  p-value    
s(x1) 1.367  1.646  25.83 1.23e-05 ***
s(x2) 5.754  6.890  51.37 8.12e-09 ***

plogis(0.733)

0.6754633

head(csale)

  purchase n_acts bal_crdt_ratio avg_prem_balance retail_crdt_ratio
1        0     11        0.00000         2494.414           0.00000
2        0      0       36.09506         2494.414          11.49123
3        0      6       17.60000         2494.414           0.00000
4        0      8       12.50000         2494.414           0.80000
5        0      8       59.10000         2494.414          20.80000
6        0      1       90.10000         2494.414          11.49123
  avg_fin_balance mortgage_age cred_limit
1        1767.197     182.0000      12500
2        1767.197     138.9601          0
3           0.000     138.9601          0
4        1021.000     138.9601          0
5         797.000      93.0000          0
6        4953.000     138.9601          0

では、練習しましょう！

Rで学ぶGeneralized Additive Models (GAMs) による非線形モデリング