問題タブ [anova]

RでN-way ANOVAを実行して、さまざまな要因間の相互依存関係を把握するのに助けが必要です。私のデータには約 100 の異なる要因があり、次のコードを使用して ANOVA を実行しています。

私が知る限り（私が間違っているかもしれません）、これは各因子のみで一元配置分散分析を実行します。何らかの理由で、x1 から x100 までの 100 グループすべての間で N-way ANOVA を実行する必要があります。次のように各要素を指定する必要がありますか、それとも簡略表記がありますか?

分散分析と構築効果テーブルを実行した後、最小値でそれらの用語のテーブルのうち5つだけをキャプチャする必要があります。要因の長いリストのため; x1からx100までの約100の因子があるため、すべてのテーブルを視覚化することはできません。

任意の用語のテーブルのラベル名は

ラベル「2」の最小値を持つこれらの用語のテーブルのうち2つだけをキャプチャする必要があります。

編集された質問：

上記の4つの表から2つの最小値のみを選択した場合、結果は-0.3008399および0.009236422になる可能性があります。

以前の投稿との対比を理解しており、正しいことをしていると思いますが、期待どおりの結果が得られていません。

これを実行すると：

ここで、私の参照はタイプ「g」であるため、typeiタイプ「g」とタイプ「i」typetの違いであり、タイプ「g」とタイプ「t」の違いです。

typei+typeg ここでさらに2つの対比、タイプ「t」の違いとタイプ「i」とタイプ「t」の違いを見たかったのです。

だからコントラスト

参照を変更して2番目のコントラストを実行しようとすると、異なる結果が得られます。コントラストのどこが悪いのかわかりません。

anova（固定または混合モデル）で単一の df 直交コントラストを実行したい。以下に例を示します。

これらのデータは、Snedecor と Cochran (1980) で分割プロット計画の例として説明されています。実験で使用された処理構造は、3 種類のアルファルファと 1943 年の 4 つの 3 回目の伐採日を含む 3\x4 完全階乗でした。実験ユニットは 6 つのブロックに配置され、それぞれが 4 つのプロットに細分化されました。アルファルファの品種 (コサック、ラダック、レンジャー) は区画にランダムに割り当てられ、3 回目の伐採日 (なし、S1 - 9 月 1 日、S20 - 9 月 20 日、O7 - 10 月 7 日) は区画にランダムに割り当てられました。各ブロックで 4 つの日付すべてが使用されました。

私はいくつかの比較 (グループ内の直交コントラスト) を実行したいと考えています。たとえば、日付の 2 つのコントラスト:

多様性係数 2 の対比の場合:

したがって、anova の出力は次のようになります。

どうすればこれを実行できますか? ...................................

サンプルが列ごとにグループ化されたデータセットがあります。次のサンプル データセットは、私のデータの形式に似ています。

上記のデータセットを使用して Excel で単一因子 ANOVA を実行すると、次の結果が得られます。

Rの典型的な形式は次のとおりです。

そして、R で ANOVA を実行するコマンドは、を使用することaov(group~measurement, data = mydata)です。 行ではなく列で整理されたサンプルを使用して、R で単一因子 ANOVA を実行するにはどうすればよいですか? つまり、R を使用して Excel の結果を複製するにはどうすればよいですか? 助けてくれてありがとう。

シーケンスの多因子分析の結果をラテックステーブル形式でエクスポートしようとしています。私はやろうとしましxtable()たが、Rはそれを##見出し##はこの機能ではそれを行うことができないと言います...誰かがアイデアを持っているでしょうか？

不均衡な設計のNestedANOVAのノンパラメトリックな対応物はありますか？さまざまな変換を行った後でも、私のデータは正常ではありません。したがって、パラメトリックテストを使用できませんでした。フリードマン検定に出くわしましたが、完全でバランスの取れた設計が必要であることを理解しています。

私が決定したいのは、AがファクターC（例：調査サイト）にネストされている場合の、B（例：キャッチ）におけるA（例：ギアタイプ）の影響です。

私はRを初めて使用し、Statの専門家でもありません。どんな助けでも大歓迎です。どうもありがとうございます。

結果変数と、カテゴリ変数であるx3 つの説明変数があります。a, b, c私の例aでは、8 つのレベル、b4 つのレベル、c35 のレベルがありますが、3 つの変数のすべての組み合わせに観測値があるわけではありません (これはおそらく重要ではありません)。

Stata で次の加法的 ANOVA モデルを実行すると、

次に、変数およびxによって調整された変数の予測を取得します。調整コマンドは、結果ウィンドウに次の表を出力し、予測が調整された変数を生成します。aby

私の問題は、変数yには と の各組み合わせの値があるのに対し、上の表にはa, bとの各組み合わせのc値しかないことです。テーブルから結果を保存するにはどうすればよいですか? 表の値と の値の関係は何ですか?aby

前もって感謝します。

更新：私はこれを見つけましたhelp adjust：

推定コマンドで使用されているが、by() 変数リストにも調整変数リストにも含まれていない変数は、観測ごとに現在の値のままになります。ここで、adjust は、by() オプションの変数によって定義された各グループ内のこれらの未指定の変数の平均を代用して、推定された予測の平均 (または対応する確率または指数化された予測) を表示します。

これは私のデータにも当てはまります。たとえば、a=75とのb=2場合c、値は 12、13、14、15、16 になります。y対応する値c=14(平均) は、まさに表に表示されているものです。しかし、値の平均がそれが取る値でない場合はどうなるでしょうか?

anova作業ディレクトリに保存されている複数のデータセットに対して実行したいと考えています。私はこれまでに思いついた：

ご覧のとおり、私はループに非常に慣れていないため、これを機能させることができません。また、すべてのサマリー出力を 1 つのファイルに追加するコードを追加する必要があることも理解しています。ディレクトリ内の複数の .csv ファイル (同じヘッダー) に対して anovas を実行し、レコードの出力を生成できる作業ループをガイドするために提供できるヘルプをいただければ幸いです。

2つのモデルをF比検定で比較するRの関数を使用する必要があるように、スクリプトを（pythonで、pypeRのR部分を使用して）作成しています。

モデルは次のようになります。

モデル1： Response ~ Predictor A + Predictor B + Predictor C.... + Predictor n
モデル2： Response ~ Predictor 1

一緒に予測子A+B+...nが構成されてPredictor 1いるので、ここでネストしても問題はありません（私を信じてください）。

Predictor A + Predictor B + Predictor C.... + Predictor n作成した関数に渡すと、1つの変数として扱われていると思います（自由度はと同じであるためModel 2）。おそらくこれは私が使用しているためですpaste()か？とにかく、モデル1の予測子の実際の数は実行ごとに変化するため（関数として必要な理由です）、を使用する以外にこれに対応する方法がわかりませんpaste()。

ここでは、ペーストが実際には問題にならない可能性があることに注意してください。問題があるのではないかと人々に知らせたかっただけです。

真の残差逸脱度と自由度を取得する方法についての提案はありますmodel 1か？それはハックである可能性があります。たとえばlength(vector of predictors) - 1、自由度を取得するために単純に減算していました。残差逸脱に対する同様のハックがどのようなものになるかはわかりません。

関数とインスタンス化の例を次に示します。