久保先生の本「データ解析のための統計データモデリング」のRコードをpythonで。