問題タブ [lm]

からlmオブジェクトを取得し、標準エラーを抽出したい

関数の概要、名前、係数を知っています。

ただし、要約は、標準エラーに手動でアクセスする唯一の方法のようです。

どうすればseを出力できるかわかりますか？

OpenEarsを使用しようとしていますが、迷惑なエラーが発生します。プロジェクトにlanguagemodelファイルを追加した後、アプリを実行しようとすると、すぐにそのエラーが発生します。ここにエラーがあります：

これが言語モデルファイルです：

私はあなたが私を助けてくれることを本当に望んでいます。よろしく、アディ。

Rの線形モデルの要約テーブルがある場合、行番号を数えることなく、交互作用の推定値またはグループの切片などに関連付けられたp値を取得するにはどうすればよいですか？

たとえば、連続およびカテゴリlm(y ~ x + group)のようなモデルの場合、オブジェクトのサマリーテーブルには次の推定値があります。xgrouplm

1. 切片
2. x、すべてのグループにわたる勾配
3. 全体的なインターセプトからのグループ内の5の違い
4. 全体的な勾配からのグループ内の5の違い。

グループの数やモデルの式が変わっても、これらをそれぞれp値のグループとして取得する方法を見つけたいと思います。要約テーブルが行をグループ化するために何らかの方法で使用する情報があるのではないでしょうか。

以下は、2つの異なるモデルを使用したデータセットの例です。最初のモデルには、個別に取得したい4つの異なるp値のセットがありますが、2番目のモデルには2つのp値のセットしかありません。

subset私はRの関数の引数がどのように機能するかを理解しようとしてきましたlm()。特に、次のコードは私には疑わしいようです。

いずれの場合も、回帰には32個の観測値があります

それでも、係数は（R²とともに）変化します。ヘルプは、この問題についてあまり多くの情報を提供していません。

サブセットフィッティングプロセスで使用される観測値のサブセットを指定するオプションのベクトル

問題：高次のパラメーター（つまり、交互作用）がモデルに残っている限り、モデルの低次のパラメーター（たとえば、主効果パラメーター）を削除できません。その場合でも、モデルはリファクタリングされ、新しいモデルは上位モデルにネストされません。
次の例を参照してください（私が使用しているANOVAから来ているためcontr.sum）：

ご覧のとおり、1つのパラメーター（A）を削除しましたが、新しいモデル（m2）はリファクタリングされ、より大きなモデル（ ）にネストされていませんm1。手あたりの因子を数値コントラスト変数に変換すると、目的の結果を得ることができますが、Rの因子機能を使用してどのように得るのですか？

質問： Rの低次因子を削除して、このパラメーターを実際に見逃し、リファクタリングされていないモデルを取得するにはどうすればよいですか（つまり、小さいモデルのパラメーターの数を少なくする必要があります）。

しかし、なぜ？パッケージの関数lmerを使用して、モデルのp値のような「タイプ3」を取得したいと思います。したがって、この例は実際には単なる例です。KRmodcomppbkrtest

なぜCrossValidatedではないのですか？これは実際には統計の質問よりもRの方が多いと感じています（つまり、モデルを相互作用に適合させるべきではないが、主要な効果の1つを持たないようにする必要があることはわかっていますが、それでもやりたいと思います）。

私はRを初めて使用し、1つの問題に悩まされています。例を挙げて説明します。

私はPHPでRを使用しています。線形回帰を計算するRスクリプトが1つあります。

そして、私は次のような回帰結果reg_resultを使用していますprediction：

それは非常にうまく機能します。

しかし、私の問題は、1時間ごとなど、いくつかの時間間隔で使用predictする必要があることです。reg_result

それで、reg_result後でそれを使用するために回帰の結果を永続化できる方法はありますか？

私はRMySQLを使用しましたが、回帰の結果をmysqlに保存する方法がわかりません。

どんな助けでも大歓迎です。

また、必要に応じて詳細を提供できます。

編集：MySQL blob RMySQLでRオブジェクトを保存することは可能ですか？私がグーグルで見つけたようにCurrently there are no facilities to import/export BLOBs. in RMySQL

Rで単純な線形モデル（Y〜X）を実行しています。ここで、予測子はカテゴリ変数（0〜10）です。ただし、この変数は正規分布ではなく、使用可能な変換手法 (log、sq など) はどれも役に立ちません。データが負/正に歪んでいるのではなく、いたるところに歪んでいるからです。lm の場合、結果変数 (Y) は正規分布する必要があることは承知していますが、これは予測子にも必要ですか? はいの場合、これを行う方法の提案は大歓迎です。

また、私が見ているデータには患者とコントロールの 2 つのグループがあるため (ご想像のとおり、グループの違いに興味があります)、データが 2 つのグループ内で正規に分布しているか、全体的に分布しているかを確認する必要がありますか? 2つのグループ？
ありがとう。

lmRIでは、多くの場合、タイプオブジェクト（、、glmなどlme）のリストを循環clmし、新しい式または他の引数でそれらを更新する必要があります。lapply()名前をそのままにして同じリスト形式で返すため、これを使用してこれを行います。これにより、プロセスを続行できます。

私が渡す引数は動的であり、モデルごとに異なる場合があります。例えば、

ここjjで、はリスト内の各モデルに追加される用語であり、モデルごとに他のリストで定義されます。

特定のモデルの意図された結果式がであると想定されたとしましょうy~a+b+c。その代わりに、結果のオブジェクト内の呼び出しオブジェクトには。lmが含まれますy~a+b+jj。そのため、状況によっては、結果のlmオブジェクトは正常に動作しますが、数式を評価する必要があるときはいつでも、jj長い間使用されていないためにエラーになります。update()出力にオブジェクトを書き込むときに引数内のすべての変数を強制的に展開するか、オブジェクトcallを使用する関数にグローバルコンテキストではなくオブジェクトの内部コンテキストでcall呼び出しを評価するように強制する推奨される方法はありますか？少なくとも、スタイルオブジェクトに通常含まれているまたはオブジェクトlmに新しい変数を詰め込むことができますか？modeldatalm

また、動的変数で更新される部分がデータ引数である場合はどうなりますか？例えば：

...kk現在のモデルに固有の場所（おそらく、sharedByAllModelsデータフレーム内の特定の列のある種の集約）。

データフレームの操作

私はこのエラーを取得しています：

私は何が間違っているのですか？

これがdputの出力です：

対数関係にあると思われるデータポイント (x と y) がいくつかあります。

3ここで、グラフに適合し、他のデータポイント (または)を推測できる基になる関数を見つけたいと思います82。私はについて読みましたがlm、nls実際にはどこにも行きません。

最初に、プロットに最も似ていると思われる関数を作成しました。

その後、以下を使用してフィッティング モデルを生成しようとしましたnls。

ここから何をすべきかについて、誰かが私を正しい方向に向けることができますか?

ファローアップ

あなたのコメントを読み、もう少しグーグルで調べた後、 の開始パラメータを調整するとa、突然モデルが収束しました。bc