lme4

0熱

1答えて

glmerの機能がlme4から別のマシンで異なる結果を生ずる理由はありますか？マシン内のハードウェアは実質的に異なりますが、すべて同じOS、R、およびパッケージのバージョンを実行しています（これは実際には当てはまりません）。この式には、グループ化された2項応答変数と22個の連続した固定効果があり、それらはすべて同じスケールであり、いくつかのランダムな効果があり、文字列であり、logitリンク機能

1熱

1答えて

ランダム効果を有する一方向ANOVAの事後検定

私は連続応答変数yldとカテゴリープレディラーcheck（3つのレベル）を持っています。私は、一方向ANOVAと事後検定を行って、どのレベルがお互いに異なるかを調べました。 mdl<-aov(sqrt(var$yld) ~ var$check); summary(mdl);TukeyHSD(mdl) Df Sum Sq Mean Sq F value Pr(>F) var$c

1熱

1答えて

私のランダムエフェクトが別のランダムエフェクトの下にネストされている場合はRコード

Rで私の研究のために混合線形モデルを実行しようとしています。私のコードが正しいかどうかを知りたいと思います。私のデザイン - 5つのサイト、各サイト内に2つのサブサイト、各サイト内に2つの永続的な四角形があります。私は5つのサイト、10のサブサイト、20のクアドラットを持っています。私はすべての四角でコロニーの大きさ（サンゴ）を測定しました。私の質問は、サイズの構造はサイトによって異なります

0熱

1答えて

lmer時間変数の標準誤差はR

です。したがって、私は時間の経過とともに何か（MS_TOT）の変化を予測しようとしています。私は自分の時間変数（ExamStage：1,3,6、および12ヶ月）と私の薬物使用変数（acstatus、ファクタw 3レベル）を持っています。私は（他の変数を無視する）薬物使用は、私は私のlmerモデルに入れた時間ワットの相互作用、としてこの変更に影響を与えることを考える：あなたのよう > summary

4熱

1答えて

dplyr :: failwithはlme4 :: lmerでは動作しませんが、lmで動作します

do、failwith、lmerで動作していたワークフローがありました。今それは壊れているようです。私が行っていることはlmと動作するので、バグか何かを変更する必要があるのだろうかと思います。 library(dplyr) library(lme4) lmer( Yield ~ 1|Batch, Dyestuff, REML = FALSE) lm( Yield ~ Batch, Dye

2熱

1答えて

レベル3のランダム切片をR

に指定する 120人の被験者を6回測定した縦断研究を分析するために、Rにlmer()の機能（lme4パッケージ）を使用しています。最初のインスタンスでは、私はこのようなモデルを指定： library(lme4) model1 = lmer(DV ~ 1 + X1*X2 + (1+X1|SubjectID), REML="false") X1は、時間的に変化する変数（レベル1）であり、X2は、

4熱

1答えて

ランダムな効果変数は、Rのlmer（またはlme）の要素として自動的に取り込まれますか？

ランダム効果混合効果モデルで連続変数または数値変数を持つことは意味がありません（たとえば、hereを参照）。しかし、私は思ったんだけどRでlme4::lmerまたはnlme::lmeが意図的にそうすることからあなたを防ぐ場合ですが... 具体的には、私が求めていることである：私はlmer（またはlme）を供給する場合には任意の非ランダムな効果としての因子（非カテゴリー）変数は自動的にそれを因子と

0熱

1答えて

混合エフェクトモデルによる縦方向データの分析R

混合エフェクトモデル（lme4パッケージ）を使用して、Rでシミュレートされたいくつかの縦方向データを分析しようとしています。模擬データ：25人の被験者は5つの連続する時点で2つのタスクを実行しなければならない。 #Simulate longitudinal data N <- 25 t <- 5 x <- rep(1:t,N) #task1 beta1 <- 4 e1 <- rno