問題タブ [r-car]

の x 軸と y 軸を変更する方法はありmmpsますか?

正規分布に変換する必要があるデータセットがあります。

まず、再現可能なデータセットを生成します。

次に、関数を定義します。この関数は、P > 0.05 になるまで df を変換し続けます。変換された df が生成され、y という名前が付けられます。

第三に、自由度をテストする

P < 0.05 であるため、df を変換します。

次に、次のエラーメッセージが表示されます。

私は何を間違えましたか？

バージョン 0.99.467 – © 2009-2015 for R Studio を使用していますが、Mac OS X 10.8.5 に CAR パッケージをインストールできません。

コンソールで install.packages("car") を試し、CRAN Web サイトから手動で zip ファイルをインポートしました: car_2.0-26.tgz / car_2.0-25.tgz

また、車のアーカイブから古いソースをインストールしようとしましたが、どれも機能しません。

統計の割り当てにこのパッケージが必要ですが、うまくロードできません。

キャレットパッケージをインストールできません!!

「車」パッケージを個別にインストールすると、次のように表示されます。

私は使っている

互換性の問題が発生していたため、R 3.2 から R 3.1.2 にダウングレードしました。
「car」パッケージはR 3.2以降のものであるため、R 3.2にアップグレードせずにキャレットを使用する方法はありますか?

次の形式で boxTidwell 関数を実行しようとしています:

次のエラーが表示されます。

これを解釈できる唯一の方法は、x 値の一部が <0 であることを示しているということです。しかし、これは真実ではありません - 私は確認しました。

何かアドバイス？

ありがとう！

主な質問

被験者間因子を持たないデータセットに対して、R で一方向反復測定 MANOVA を設定するためのヘルプを探しています。

バックグラウンド

被験者間因子を使用して RM MANOVA を設定するための優れたガイドはたくさんありますが、完全に被験者内のデザインを使用している場合、私はまだ見つけることができませんでした. この問題はかなり単純明快なように思えますが、私は MANOVA を初めて使用するので、問題に正しく取り組んでいるかどうかはわかりません。私は主carに R でパッケージを使用してきましたが、これを別の方法で行う方法についての提案は受け付けています。

この問題を説明するために、OBrienKaiser データ セットのサブセットを使用し、Hours被験者内因子の各レベルが異なる従属変数の測定値を表すと仮定します。次に、preandpost条件を、単一の被験者内独立変数の 2 つのレベルとします。簡潔にするために、 の最初の 3 つのレベルだけを見ていきますHours。

したがって、私のデータ セットには、3 つの異なる従属変数 ( 、、および) で 2 つの異なる条件 (preおよび) で測定された 16 人の被験者がいます。post123

私の目標は、3 つの異なる従属変数の組み合わせの違いpreを探すことです。post私はこのガイドに大きく依存してきました...

http://socserv.mcmaster.ca/jfox/Books/Companion/appendix/Appendix-Multivariate-Linear-Models.pdf

...しかし、そこで使用されているケースはまったく同じではなく、私が持っていない被験者間の条件が含まれています。これまでのところ、次のアプローチで潜在的に正しい結果を生成することができましたが、私の問題は、問題を正しく設定したかどうか完全に確信が持てないことです。一般に、私のアプローチは上記の手順に従うことでしたが、被験者間の条件は単純に省略しました。

私のアプローチ

まず、条件と従属変数を 2 つの異なる被験者内因子として扱うことによりidata、関数呼び出しの行列を設計します。Anova()

次に、被験者間の要因を無視して、データセットで多変量線形モデルを構築します。

次に、以前に定義した を使用しAnova()て線形モデル オブジェクトを呼び出し、被験者内計画を に設定します。mod.mlmidata~Condition

これにより、次の出力が得られます...

私には、このプロセスとこの結果は合理的に思えます。私を一時停止させるものは2つあります。

1. 線形モデルを設定するとき、予測変数がないのは奇妙に思えますが、Anova()withの呼び出しで定義されることになると思いますidesign。それがAnovaを設定する方法である場合car、素晴らしいです。明示的に関心のある予測変数があるときに、予測変数なしで線形モデルを構築するのは奇妙に思えます。

2. 出来上がっsummary(an.out)たモデルをよく見てみると、デザインのコントラストが見えてきます。上記のアプローチ コードpreで生成されるコントラストは 1 およびpost-1 です。これが私がやろうとしていることの適切な対比であると信じていますが、完全には確信が持てません。imatrixまたはcontrastsの呼び出しでカスタム コントラストを渡すことができるので、テストしようとしているもの (つまり、3 つの従属変数間の差) が実際にテストAnova()しているものであることを確認したいと思います。prepost

このコンテキストで一般的に反復測定MANOVAを理解する方法に関するヘルプおよび/またはアドバイス、およびRでこれを実装する方法に関する具体的なアドバイスをいただければ幸いです。

ボーナス

私はMatlabでも同じことをしたいので、誰かがそれについて具体的なアドバイスを持っているなら、それをいただければ幸いです(ただし、これには独自の質問が必要になる場合があります)。

aov() 用に書かれたコードを car-package の Anova() 関数に変換したい。

このコードにより、次の出力が得られます。

ただし、 car-package の Anova() 関数を使用しようとすると、次のようになります。

これにより、この出力が得られます。

aov() に使用するコードを Anova() に合わせて作り直す方法を知っている人はいますか。https://gribblelab.wordpress.com/2009/03/09/repeated-measures-anova-using-r/のチュートリアルに従ってみました 。

Anova() を正しく取得しようとします。しかし、似たような出力は得られません。また、ウェブページから、Mauchlys と Greenhouse が提供されていると思われることがわかりますが、これはわかりません。また、Anovaの結果でイータ二乗を取得する方法を知っている人はいますか? または、別の関数を使用して eta (etaSquared()) を計算する必要があります。

以下のデータはテストに使用されました。時間、性別、グループ間の「値」と、因子間の交互作用効果に有意差があるかどうかをテストしようとしています。