【R言語】統計解析フリーソフトＲ 第４章【GNU R】
レス数が950を超えています。1000を超えると書き込みができなくなります。
0001１３２人目の素数さん
2011/01/09(日) 17:31:49統計計算で重宝するデータ型や、複数要素を処理する演算や関数、
解析結果を表示するグラフィックなど、多彩な機能を提供します。
●関連サイト
The R Project
http://www.r-project.org/
RjpWiki
http://www.okada.jp.org/RWiki/
リンク集
http://www.okada.jp.org/RWiki/?%A5%EA%A5%F3%A5%AF%BD%B8
0968１３２人目の素数さん
2013/08/23(金) NY:AN:NY.AN下記データ( r t ）があります。
r <- c( 0.727 , 0.634 , 0.500 , 0.443 , 0.416 , 0.328 , 0.272 , 0.225 , 0.133 )
t <- 1:9
このとき、以下のようにロジスティック回帰を行うと、
glm( r ~ x , family=binomial )
>In eval(expr, envir, enclos) : 二項 glm で整数でない成功数がありました!
とエラーになります。
ここで質問です。
１．小数データを目的変数に、ロジスティック回帰はできないのですか？
できるのであれば、どうすればよいのでしょうか？
２．ふと疑問に思ったのですが、目的変数が０、１って、そのままでは、ロジット変換できないと思うのですが、
内部的にどのように計算されているのでしょうか？（微小数を加減算している？0+0.0000001と1-0.0000001 ）
0969１３２人目の素数さん
2013/08/23(金) NY:AN:NY.AN> エラーになります
エラーになりません。「警告」と「エラー」は違います。
> 小数データを目的変数に、ロジスティック回帰はできないのですか？
できますが、多分あなたが思っているものと違います。
glm()でロジスティック回帰をやるときには、
目的変数の設定仕方が３通りあります。
>>968を見ると、目的変数は、二値でも、説明変数カテゴリごとの計数(整数)でもないから、
説明変数カテゴリごとの割合(小数)だよね。
その場合は、分母に相当するweightsを指定しないと駄目だよ。
何を言っているのか、分からなければ、Rの解説書を読んだ方がよいよ。
Rで学ぶデータサイエンスシリーズ「一般線形モデル」なら、pp.53-56のあたり。
2. の質問は、目的変数をロジット変換するという思いこみによる誤解。
実測値(0,1)とモデル値exp(beta_1x_1+beta_0)/(1-exp(beta_1x_1+beta_0))を比較して、
当てはまりを見る話じゃないの？
0970１３２人目の素数さん
2013/08/23(金) NY:AN:NY.ANレスありがとうございます。
確かに「警告」であって「エラー｣ではないですね。
私が提示した例は、説明変数カテゴリごとの割合(小数)です。
その上で、明示的に分母に相当するweightsを指定しなくても
同じウエイトが設定されるのでは？と思いました。如何でしょう？
そうだとしたら、「警告」を出力しなくても良いのではと思いました。
2.
> 実測値(0,1)とモデル値exp(beta_1x_1+beta_0)/(1-exp(beta_1x_1+beta_0))を比較
なるほど。すっきりしました！
実測値(0,1)とモデル値が、完全に一致しない（限りなく近くはなるが・・・）ところに
私は何か引っかかっていたようです。
レス数が950を超えています。1000を超えると書き込みができなくなります。