問題タブ [anova]

を使用して ANCOVA を実行したいのcar::Anovaですが、主効果としてのみ共変量を追加する方法があるかどうかを確認できません (つまり、何とも相互作用しないでください)。

私が ANCOVA を理解している限り、共変量はモデルに追加されたもう 1 つの主効果 (つまり、もう 1 つの効果) であり、それによってこの共変量の全体的な相加的影響を制御します。したがって、共変量は他の要因と相互作用しません。ただし、Anova被験者内要因と相互作用しない変数を追加することはできません (つまり、私の最終モデルは ANCOVA に適合していないようです)。

の例で私の問題を説明しましょう?Anova。OBrienKaiserデータセットには、( と ) の間に 2 個、(およびtreatment)gender内に 2個phaseありますhour。ここageで、参加者の も記録し、それを共変量として任意の分析に追加したいとします。

結果が示すように、結果には、共変量age、つまり被験者内 (または反復測定) 因子との相互作用phase、hourおよびそれらの相互作用が含まれていphase:hourます。

私の質問は:で ANCOVA を実行できますか?car::Anovaもしそうなら、この ANCOVA を相互作用なしで指定する方法はありますageか?

更新(2012 年 7 月 22 日): R-helpでこの質問をしましたが、今のところ回答がありません。ニュースがあれば、ここに投稿します。

Rで分散が等しくないANOVAを実行する方法があるかどうか疑問に思っていましたか?

次の例を想像してください。

この例からわかるように、R で ANOVA を実行するさまざまな方法は、2 つのグループのみの場合、等分散の t.test によって取得されたものと常に同じ p 値になります。繰り返しますが、分散が等しくない ANOVA を実行する方法はありますか?

タイプ III エラーと種と生息地内のペアワイズ比較を使用しながら、ホームレンジのサイズに対する種と生息地の影響を比較したいデータセットがあります。
データのサブセットを次に示します。

現在、私は 3 つの種のデータを分割し、それぞれについて別々の ANOVA を実行していますが、1 つの ANOVA で同時に種と生息地について尋ねる方が理にかなっていると思います。ある種に対して実行した ANOVA の例を次に示します。

aov() はタイプ I エラーを使用しているようです。. . 適切ではないと思います。さらに、テューキーの検定はペアワイズ比較のアプローチとして保守的すぎる可能性があると思います。ホームレンジに対する種と生息地の両方の影響を考慮した1つのANOVAを実行できるアプローチを誰かが助けてくれますか?

4つのグループx、y、z、bを比較して、どのグループが大幅に異なるかを取得したいと思います。

ありがとう！

一元配置分散分析テストを使用して、適合度、つまり、実験データと予測した他のデータとの差を測定します。Minitab ソフトウェアでテストを実施し、F 値と P 値の両方で許容値を取得しましたが、R2 値は以下に示すように 0.30 % と小さかったです。

一元配置分散分析: ANN Sr、Exp.Sr

出典: DF=1 , SS= 0.824 , MS= 0.824 , F =1.89 , P=0.170

S = 0.6605、R-Sq = 0.30%、R-Sq(adj) = 0.14%

fとpの値（それぞれ1.89と0.17）によると、予測データは正確であり、実験データと予測データの間に違いはありません。

私の質問は：

R-sq=0.30% というこの小さな値が私の決定に影響を与えることはありますか???私が予測したデータを却下し、f&p 値が許容可能であっても実験データとの有意差があると言うことはできますか???.

ありがとうございました

ezANOVAパッケージからezループ内の複数の従属変数に対してを実行し、結果をいくつかの変数に保存したいと思います。各従属変数は、同じデータフレームの個別の列内にあります。

これはデータフレームの例です。

このアプローチの問題は、Rがインデックスdependent.variableをデータフレーム内の列の指定子として解釈するaov.dataため、次のエラーが返されることです。

「dependent.variable」は、提供されたデータフレーム内の変数ではありません。

インデックスをeval()またはでラップしようとしましprint()たが、役に立ちませんでした。

私は多くの ANOVA テーブルを実行していますが、最終的には xtable パッケージの 'xtable' 関数を使用してそれらを実行し、LaTeX にエクスポートしたいと考えています。ただし、テーブル オブジェクトを使用して xtable を提供する必要があります。ANOVA オブジェクトは受け入れられません。私は基本的にANOVAテーブルをテーブルオブジェクトとして作りたいと思っています。再現可能なコードを次に示します。

通常の容疑者（as.table、print、およびxtable（summary（npk.aoE）））をすべて試しましたが、成功しませんでした。どんな助けでも大歓迎です

純粋なエラーとを使用してANOVA表を作成しようとしていますが、関数を使用しようとするとlack-of-fit errors、Rから次のメッセージが表示されますpure.error.anova。

この関数が見つからないのはなぜですか？

私は問題があります。scipy.stats f_oneway（）ANOVAを、次のようなnumpy配列のグループを含むデータアーカイブをロードするスクリプトで実行したいと思います。

今私の問題は、グループの数が異なるデータアーカイブに対して固定されていないことです。つまり、アーカイブにグループがいくつあるか（また、必ずしもそれらの名前が何であるかはわかりません）、事前にはわかりません。

一元配置分散分析のscipy実装は、次のように入力としてコンマ区切りの配列のみを受け入れます。

リスト、タプル、多次元配列をすべて提供しようとしましたが、成功しませんでした。したがって、現在、このANOVA実装を使用できる唯一の方法は、毎回グループ変数を手動で入力することです。これにより、スクリプトでこれを実行することが事実上不可能になります。これを解決する方法や、f_oneway（）のこれらの非常に具体的なデータ形式の要件を回避する方法を知っている人がいるかどうか疑問に思います。

lmer()個々の効果を考慮に入れて、生物学的データに対する治療効果を経時的にモデル化したいと思います。通常使用される手順は、固定効果と相互作用の項を段階的に削除していくつかのモデルを構築し、anova(model1,model2)有意でない差が見つかった場合に最も偏狭なモデルを使用して維持することで最適なモデルを見つけることです。summary(model1)ただし、との出力の間に異なるモデル適合度値（AIC、BIC ..）が見つかりましたanova (model1,model2)。コードは次のとおりです。

出力を比較すると、2つのモデルの間に14のAICユニットの違いがありsummary()ますが、の中には9つしかありませんanova()。この違いはどこから来るのですか？前もって感謝します。