問題タブ [glm]

私はドイツの信用データの分析に従っていますが、ROCRパッケージを既にインストールしているため、エラーを修正できませんでした。以下は、ROCRを使用するコードです。

を実行するpredと、次のエラーが表示されます。

これは、次の行を実行しようとしたときにも発生します。

パッケージ「betareg」（バージョン3.0-3）を使用して、Rでベータ回帰を行っています。モンテカルロ シミュレーションの場合、分散をモデル化したいと思います - パラメータ "phi" をロジット モデルとして。残念ながら、パッケージは phi のリンク関数「logit」をサポートしていません。自己記述の「glm - リンク」オブジェクトを提供できることは確認しましたが、これを行う方法がわかりません。どうすれば私の問題を解決し、「ファイ」をロジスティックモデルとしてモデル化できるか、誰にもわかりません。

link.phi:
精密モデル (phi) におけるリンク関数の文字指定。現在、「identity」、「log」、「sqrt」がサポートされています。数式が y ~ x 型でない限り、既定値は "log" です。既定値は "identity" です (下位互換性のため)。または、クラス「link-glm」のオブジェクトを指定できます。

glm.fit() 関数を使用して R で glm モデルを構築しました。

その後、以下を使用していくつかの予測を試みました (s を正しく選択したかどうかはわかりません)。

エラーが発生しました：

ここでhttps://stat.ethz.ch/pipermail/r-help/2004-September/058242.html問題の解決策を見つけました:

しかし、 glm.fit を使用して後で予測できないかどうかはわかりません。なぜそれが可能か不可能なのですか？また、s を正しく選択するにはどうすればよいでしょうか。

注: glm() 関数を使用すると、この問題を回避できます。しかし、 glm() 関数は数式を要求するため、場合によってはあまり便利ではありません。まだ誰かが後で glm.fit と予測を使用したい場合は、ここにいくつかの解決策があります: https://stat.ethz.ch/pipermail/r-help/2004-September/058242.html

外気温度と TOD (96 カテゴリ変数、時刻) を使用して、15 分間隔で負荷消費の前日予測を取得する次のコードがあります。以下のコードを実行すると、次のエラーが発生します。

glm() を使用して最後の行の後に次のエラーが表示されます。

そして、predict() を使用して最後の行の後に次のエラーが表示されます。

R では、回帰関数はリストglmであるクラスのオブジェクトを生成しますglm。

glmリストなので、クラスのオブジェクトをフォーマットせずにリストとして表示できるはずです。ただし、as.listこれを行うようには見えません。

たとえば、fitが関数 によって適合されたモデルである場合glm:

一方、names関連付けられたリストの 30 個の名前を生成するなど、リストに適用される他の関数は正しく機能します。

また、他のリストと同じ方法で個々の要素を表示できます。

fitフォーマットせずにリスト形式で表示できる既存の関数はありますか?

上で述べたように、リストの名前を使用して独自の関数を作成することもできますが、そのような単純なタスクには不要のようです。

predictラスター ライブラリの関数を使用して、glm一連のラスター スタック予測子でモデルを予測しています。私のモデルではいくつかの因子レベルが欠落していたので、predict関数を実行するためにそれらの層を NA に変更しました。現在、出力マップに欠損値があり、この問題を解決する方法を探しています。他の連続変数を使用してこれらの地域を予測する方法はありますか?

これらは私のサンプルコードです:

glm モデルの適合度を表す 1 つの p 値を取得する方法を探しています。lmマンページからわずかに変更された例を次に示します。

1summary(lm.D9)つを取得します。

idがglmで同じことをする場合

私は得る

lmモデル全体の要約として F 統計量がありますが、glmありません。もう一度質問します。フィットを説明する glm モデルから 1 つの p 値を取得するにはどうすればよいですか?

ありがとう

を使用して参照カテゴリを変更したい "Higher" 、 "Middle" 、および "Lower" の 3 つのレベルで構成される独立因子変数を使用して R で二項回帰を実行すると、次のrelevelエラーが発生します。

cbsnivcat3 が要因かどうかを確認しました

しかし、factor glm(tertiary ~ factor(cbsnivcat3) , family = binomial, data = data)ワークスオケと同じです。

MQ1: エラーの原因は何ですか?

ヒント、コマンド、および解決策を事前にありがとう

str(データ)

グループの回帰モデルを当てはめる場合、グループ間の不均等な分散が問題を引き起こす可能性があることはよく知られている事実です。Rでこれを説明する方法はありますか?

Stata のソリューションは、このプレゼンテーションhttp://www.stata.com/meeting/fnasug08/gutierrez.pdf (スライド 9-17) で提案されました。スライド 14 に示す実際のコードは次のとおりです。

Rに同様のソリューションはありますか?

3 方向の交互作用 (factorA、factorB、factorC) の 3 つのカテゴリ因子と、単純にモデルに追加される 4 番目の連続因子 (factorD) を持つ一般化線形モデルで分析した一連のデータがあります。モデルから一連の Tukey 文字グループ (つまり、コンパクトな文字表示) を取得しようとしていますが、相互作用をうまく含める方法が見つかりません。factorD を含めることには興味がありません。相互作用の 3 つだけです。

これでTukey調整されたペアごとの比較が得られました：

しかし、そこからコンパクトな文字表示を作成する方法がわかりませんでした。パッケージで実行multcompできますが、相互作用ではなく、そのパッケージの主効果でのみ実行する方法を見つけることができました。

それで、agricolaeこの投稿のようにパッケージを試してみました ( https://stats.stackexchange.com/questions/31547/how-to-obtain-the-results-of-a-tukey-hsd-post-hoc-test-in -a-table-showing-groupe )は、それが機能するはずであることについて説明しています。ただし、その回答の指示に従うと、HSD.test から機能しない応答が返されました。具体的には、主効果テストを正常に機能させることができました。たとえばHSD.test(my.glm,"factorA")、交互作用を機能させることができませんでした。私はこれを試しました：

しかし、HSD.test 関数が "intxns" を有効なオブジェクトとして認識しなかったことを示すエラーが表示されます。私のGLMの）：

HSD.test 関数呼び出しの factor フィールドにナンセンスを入力すると、同じエラーが発生します。inxtns オブジェクトをチェックしたところ、良さそうに見え、行数が residu の数と一致していました。agricolaeノートでは、実際には HSD.test での相互作用の使用については説明していませんが、機能すると思います。

HSD.test をインタラクションで動作させる方法を知っている人はいますか? または、相互作用のある glm 用のコンパクトな文字表示を作成するために機能するようになった他の機能はありますか?

私はこれに何日も取り組んできましたが、解決策を見つけることができませんでした。

ありがとう！