rでは、二項()を使用し、0から1の間のpredict()値を制限する

StackOverflow https://stackoverflow.com//questions/9610318

  •  09-12-2019
  •  | 
  •  

質問

i GLM、Family= Binomial(link= logit)モデルは、以下に見られるPredict()関数内に適用されるモデルです。予測値はゼロを超えて1を超えていますが、それらを確率として保持したいと思います。そのため、適用機能で使用できるBinomial()$逆コマンドを使用します。

これは私がそれを初めて走ったときにうまく機能しましたが、Rを閉じて再び閉じた後、私はこのエラーを得る:

     Error in get(as.character(FUN), mode = "function", envir = envir) : 
     object 'ilogit' of mode 'function' was not found"
.

このコードは通常働いていたので、私は何時間も奮闘していました。誰かが私が間違っていることについてアイデアを持っていますか?これをするより良い方法はありますか?

マイコードは以下の通りです。私は他の変動も試しましたが、仕事に入れることができません。

    ## predicted probabilities 
    pp <- predict(logit_model,
            newdata=data,
            type="link",
            se.fit=T)

   ilogit <- binomial()$inverse
   yhat_prob <- lapply(pp,ilogit) #converts to probabilities
.

役に立ちましたか?

解決

確率を望む場合は、ドキュメントのtype="response"

あなたが得るエラーメッセージの場合は、おそらく?pregict.glmが必要です。

> str( binomial() )
List of 12
 $ family    : chr "binomial"
 $ link      : chr "logit"
 $ linkfun   :function (mu)  
 $ linkinv   :function (eta)  
 $ variance  :function (mu)  
 ...
.

エラー不足はおそらくあなたがロードしたいくつかのパッケージによるもので、binomial()$linkinv関数を定義しました。

ライセンス: CC-BY-SA帰属
所属していません StackOverflow
scroll top