問題タブ [lsmeans]

私は R に非常に慣れていないので、理解するために最善を尽くしていますが、現時点では使用するのは簡単だと思うので、あなたの助けを求めます.

を使用してモデルを適合させます

を使用して ls-means (または調整された平均) を計算します。

グループ化変数の平均 (エラーを含むバー) をプロットし、fixed2 どのグループが統計的に異なるかを確認するにはどうすればよいですか?

ありがとう！

MASSの関数を使用して、序数データに比例オッズ累積ロジット モデルを当てはめましたpolr(この場合、さまざまな種類のチーズを優先するデータに) :

モデルの予測をプロットするために、次を使用して効果プロットを作成しました

パッケージによって報告された予測平均から、effectsチーズの種類ごとの平均的な好みのようなものを、これに関する 95% conf 間隔と共にプロットできるかどうか疑問に思っていましたか?

編集: もともと、テューキー ポストホック テストを取得する方法についても尋ねましたが、その間に、それらを使用して取得できることがわかりました

またはパッケージlsmeansを次のように使用する

nnetの関数を使用して多項式モデルを適合させましたmultinom(この場合、オスとメスの食事の好みと、さまざまな湖のワニのさまざまなサイズのクラスを与えるデータに基づいています) :

使用して取得できる因子の全体的な有意性

そして、たとえば「湖」を使用して得た効果プロット

全体的な Anova テストに加えて、ペアワイズ テューキー ポストホック テストも実行して、獲物が食べられる多項分布の全体的な違いをテストしたいと思います。

私は最初glhtにパッケージで関数を使用することを考えましmultcompたが、これは機能していないようです。たとえば、 factor の場合lake:

lsmeans代替手段は、これにパッケージを使用することでした。

ただし、これは特定の種類の食品の割合の違いをテストします。

全体的な多項分布をさまざまな湖で比較するテューキー事後検定を取得することも何らかの方法で可能かどうか疑問に思っていました。で試しました

しかし、それはうまくいかないようです:

何かご意見は？

または、モデルglhtで機能させる方法を知っている人はいますか?multinom

線形混合モデルでの事後テストの実行について質問があります。

lme43 つのグループ、グループごとに 5 匹のヘビ、各グループの換気率が異なり ( Vent)、異なる時点で測定値が取得され ( Time)、ヘビが変量効果として指定されている ( ID)線形混合モデルを実行しています。

以下のデータサブセット:

コード：

ベントの効果のテスト:

ベントと比較:

時間の影響をテストする:

したがって、事後テストを試してください。

これは、次のエラーが表示される場所です。

以下を使用して、修正を伴うペアワイズ テストを実行できます。

しかし、これは他の変数が相互作用する場合 (私の他のデータの場合のように) には機能しないことを知っておいてください。これは可能lsmeans()ですか？

ご意見ありがとうございます。尤度比検定自体が物議をかもしていることは承知しています。混合効果の ANOVA を検討しましたが、これを不可能にするいくつかの欠落データ ポイントがあります。データは以前に別の学生によって二元配置分散分析として分析され、繰り返し測定はありませんでしたが、各ヘビが繰り返し測定されたため、これは不適切であると感じました

モデルをプロットするためにlsmipfrom lsmeansを使用しています。

しかし、凡例の位置、サイズ、タイトルなどを処理する方法について、インターネットで提案を得ることができませんでした。以前trellis.par.get()はパラメーターを確認していましたが、私の問題に関連するものを見つけることができませんでした。グラフからわかるように、凡例は "Photoperiod*Ws" のはずですが、Ws は表示されません。

グループごとにLSmean値でデータフレームを要約する方法はありますか? この例では、3 つの異なるスタディがあり、それぞれが RCB 設計です

単純な平均値は次のとおりです。

最初のスタディの LSmeans 値は次のとおりです。

すべてを一度に行うことができるのか、それとも 1 つずつ行うことができるのか、よくわかりません... よろしくお願いします!

実行したい分析があり、対比を計画しました (事後比較ではありません!) 治療群間で行いたいと考えています。治療グループ変数には (k =) 4 つのレベルがあります。合計 3 つの異なる比較を行う予定です。そのため、私の理解が正しければ、比較は k-1 であるため、計算される p 値を調整する必要はありません。

これを行うには、Rのmultcomporパッケージを使用したいと思います。lsmeans私の質問は次のとおりです。信頼区間（およびp値）を調整せずに、この計画された比較を行うことが可能かどうかを誰かが知っていますか? 私が見たビネットと私が見た例からわかる限り、summary.glht()関数はデフォルトとして調整を行い、どのオプションがこれを元に戻すかは明確ではありません.

再現可能な例が必要な場合は、 http : //www.ats.ucla.edu/stat/r/faq/testing_contrasts.htmで見つけた次の例を使用できます。

orangeslsmeans で提供されるデータを使用しています。

信頼区間は をプロットすることで視覚化できますがplot(contrast(days.lsm, "trt.vs.ctrl", ref = c(5,6)))、days_contr.lsm

信頼区間を data.frame に抽出するにはどうすればよいですか?