問題タブ [anova]

使用しようとしましTukeyHSD(my_anova$aov)たが、エラーが発生します：

Google は、「aovlist」で事後処理する方法はないと言っています。しかし、おそらく、ezANOVA 出力を使用した事後処理について何らかの考えがあるでしょう。

例：

multcomp を使用してみてください:

使用してみてくださいlme：

年ごとの人口を見る時系列の単純化された例が与えられた場合

年間の変化/どの年が互いに有意に異なるかという点で有意性をテストするための最良の方法は何ですか?

私はRを学んでおり、lm()因子変数を処理する方法とANOVAテーブルを理解する方法を理解しようとしています. 統計初心者なので、お手柔らかにお願いします。

Rotten Tomatoes のムービー データを次に示します。G、PG、PG-13、および R と評価された 4 つのグループのすべての映画の平均スコアに基づいて、各映画のスコアをモデル化しようとしています。

および ANOVA 出力:

このテーブルのすべての数値を取得する方法を理解しています。例外Sum SqとMean Sqfor as.factor(movies$rating). Sum Sq誰かが私のデータからどのように計算されるか説明してもらえますか? で割ったMean Sqだけだと分かります。Sum SqDf

目標: 反復測定 ANOVA から eta 値の 95% 信頼区間を計算します。

計画は 2 因子計画 (因子 1 には 3 つの水準、因子 2 には 7 つの水準) であり、30 人の被験者と完全に交差しています。ただし、イータの信頼区間を計算するためにブートストラップを正しく行っているかどうかはわかりません。boot() 関数を使用して被験者ごとにサンプルをブートストラップしますが、因子と因子内のレベルはそのままにします。したがって、これを正しく行っているかどうかはわかりません-ファクター/レベルでリサンプリングするより複雑なブートストラップ/リサンプリングを行う必要がありますか、それとも被験者のレベルで行っても大丈夫ですか. 私のコードは妥当な結果をもたらすようです...

オブジェクトからデータを抽出するための、それほど複雑でない方法を見つけようとしていaovます。以下に示すようなデータセットがありa、そのデータに基づいてANOVAを実行した結果、という名前のオブジェクトが作成されたとしますa.model。を使用してデータを検索しようとしましたstr(a.model)が、見つかりませんでした。lmオブジェクトからデータを抽出する方法を知っているので、私が行ったことはlm(a.model)$model$score、を使用することでした。しかし、最初にオブジェクトをオブジェクトにa.model変換せずに、から直接データを抽出することは可能ですか？-変換の「余分な」ステップはそれほど多くの作業ではないため、これは何よりも好奇心から外れていると思います。aovlm

lme4 ライブラリを使用して anova モデルを実行しようとしています。実際のデータがないので、MASS から Cars93 データをヤンクしました。lmer(...) を機能させることはできましたが、anova テーブルを印刷できないようです。

何が起こったかは次のとおりです。

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

まず、「ANOVA」というタグですみません、MANOVAについてです（まだタグにはなりません...）

私が見つけたチュートリアルから、すべての例は小さな行列を使用しています。多くの研究の場合と同様に、大きな行列の場合、それらに従うことは現実的ではありません。

14 のサンプリング ポイントに対して 2 つのマトリックスを取得しました。1 つは生物 ID (4493 ID) 用、もう 1 つは化学プロファイル (190 変数) 用です。

2 つの行列は、spearman によって相関され、相関に基づいて、4 つのクラスター (平方ユークリッド クラスタリング値に関する k 平均)、行の ID、および化学プロファイルの行に分割されます。

それらの間の違いはある程度明らかですが、より堅牢な方法でクラスター間およびクラスター内の違いを示すためにMANOVAを実行したいと思います-もちろん、それは結論の重要な要素です.

問題は、8時間試行した後、分析に受け入れられる形式でデータを入力することさえできなかったことです.

私が見つけたチュートリアルは、非常に少数の変数に合わせて設計されており、それを克服したと思っても、プログラムは行列を長さの違いで比較できないと言っています。

各クラスターには、すべて同じ変数セットを共有する独自の ID セットがあります。

私は何をすべきか？

前もって感謝します。

小川ディオゴ

MATLAB の一連のデータに対して双方向ランダム効果 ANOVA を実行しようとしています。これは私が使用しているコードです:

ただし、これを行うとエラーが発生します。

このエラーの意味はわかりますが、なぜエラーが発生するのかよくわかりません。コードを次のいずれかに変更すると、エラーは発生しませんが、必要なデータが得られません。

「相互作用」モデルを使用しようとすると、うまくいきません。なんで？本当にデータが多すぎますか? 各入力変数 (modelfit、RATEES、RATERS) のサイズはすべて 14879 x 1 です。つまり、14879 の値を持つ 1 列の配列です。

これが本当に多すぎる場合、必要な相互作用効果を得るにはどうすればよいですか?

編集:関連する場合、私は64ビットのMac OS X 10.6.8(Snow Leopard)とMATLAB R2008aを使用しています。

EDIT2：anovanコマンドの「ランダム」要素の使用方法もわかりません。インデックスを入力する必要がありますが、単一の変数をランダム効果として選択することしかできませんか? RATEES と RATERS の両方をランダムにしたいのですが、選択できるのは 1 つだけです。MATLAB がこの種の解析に柔軟性がないのはなぜですか?

aov()私はループの中を走っていました。各反復の出力は、リスト項目として保存されました。以前model.frame()に各aovオブジェクトのデータ フレームを確認したところ、それらはすべて同一でした。実際、それらはすべて、基本的に最後の反復のデータ フレームでした。さらに興味深いことに、これaovは誤差項を含むモデルでのみ発生するようです。 これを説明するために、R の ChickWeight データを使用します (以下を参照)。

対照的に、これはエラー項がない場合には起こりません。

何が起こっているのか誰かが知っているのだろうか。

私は不在-存在データに関するモデルを当てはめています.ランダム因子が有意であるかどうかを確認したいと思います. これを行うには、glmm と glm を比較し、どちらが最も重要であるかを LR テストで確認する必要があります。

しかし、ANOVA(glm,glmm) を実行すると、Deviance Table の分析が得られ、モデルを比較する出力は得られません。

両方のモデルを比較して、希望する出力を得るにはどうすればよいですか?

前もってありがとう、コーエン