lme4 i nlme są pakietami R stosowanymi do dopasowania liniowych, uogólnionych liniowych i nieliniowych modeli efektów mieszanych. W przypadku ogólnych pytań dotyczących modeli mieszanych użyj znacznika [mieszany model].
Moje pytanie opiera się na tej odpowiedzi, która wykazała, który lme4::lmermodel odpowiada dwukierunkowej ANOVA z powtarzanymi pomiarami: require(lme4) set.seed(1234) d <- data.frame( y = rnorm(96), subject = factor(rep(1:12, 4)), a = factor(rep(1:2, each=24)), b = factor(rep(rep(1:2, each=12))), c = factor(rep(rep(1:2, each=48)))) # standard two-way repeated measures ANOVA: summary(aov(y~a*b+Error(subject/(a*b)), d[d$c == …
Jak uzyskać znormalizowane wagi regresji (efekt stały) z regresji wielopoziomowej? I jako „dodatek”: Jaki jest najłatwiejszy sposób na uzyskanie tych standardowych wag z obiektu mer-object (z lmerfunkcji lme4pakietu w R)?
mgcvOpakowanie Rposiada dwie funkcje montowania interakcji produktów napinacz: te()i ti(). Rozumiem podstawowy podział pracy między nimi (dopasowanie interakcji nieliniowej vs. rozkładanie tej interakcji na główne efekty i interakcję). To, czego nie rozumiem, to dlaczego te(x1, x2)i ti(x1) + ti(x2) + ti(x1, x2)może powodować (nieznacznie) różne wyniki. MWE (dostosowany z ?ti): …
Pracuję nad zestawem danych w celu oceny wpływu suszenia na aktywność mikrobiologiczną osadów. Celem jest ustalenie, czy wpływ suszenia różni się w zależności od rodzaju osadu i / lub głębokości osadu. Projekt eksperymentalny jest następujący: Pierwszy czynnik Osad odpowiada trzem rodzajom osadu (kodowany Sed1, Sed2, Sed3). Dla każdego rodzaju osadu …
Mam nadzieję, że ktoś może pomóc w tym, co uważam za stosunkowo proste pytanie, i myślę, że znam odpowiedź, ale bez potwierdzenia stała się ona czymś, czego po prostu nie mogę być pewien. Mam dane zliczania jako zmienną odpowiedzi i chcę zmierzyć, jak ta zmienna zmienia się wraz z proporcjonalną …
Próbuję odtworzyć kilka testów interakcji między obiema lmi lmerpowtarzanymi pomiarami (2x2x2). Powodem, dla którego chcę porównać obie metody, jest to, że GLM SPSS dla powtarzanych pomiarów daje dokładnie takie same wyniki, jak lmprzedstawione tutaj podejście, więc na koniec chcę porównać SPSS z R-lmer. Do tej pory udało mi się tylko …
Staram się ręcznie wyliczyć przewidywania efektu losowego z liniowego modelu mieszanego i używając notacji dostarczonej przez Wooda w uogólnionych modelach addytywnych: wprowadzenie do R (str. 294 / str. 307 z pdf), mylę się co do każdego parametru reprezentuje. Poniżej znajduje się podsumowanie z Wood. Zdefiniuj liniowy model mieszany Y= Xβ+ …
Jako odpowiednik tego postu pracowałem nad symulacją danych ze zmiennymi ciągłymi, nadając się do skorelowanych przechwyceń i nachyleń. Chociaż na stronie i poza nią są świetne posty na ten temat , miałem trudności ze znalezieniem od początku do końca przykładu z symulowanymi danymi, który byłby podobny do prostego scenariusza z …
W lmerfunkcji w lme4w Ristnieje wezwanie do skonstruowania matrycy modelu efektów przypadkowych, ZZZ , jak opisano tutaj , na stronach 7 - 9. Obliczanie obejmuje produkty KhatriRao i / lub Kronecker dwóch matryc, i . J i X iZZZjotjaJiJ_iXjaXiX_i Macierz to kęs: „Macierz wskaźników wskaźników współczynnika grupowania”, ale wydaje się, …
Przeprowadzam testy post-hoc na liniowym modelu mieszanych efektów w R( lme4pakiecie). Używam multcomppakietu ( glht()funkcji) do przeprowadzania testów post-hoc. Mój plan eksperymentalny to powtarzane pomiary z przypadkowym efektem blokowania. Modele są określone jako: mymod <- lmer(variable ~ treatment * time + (1|block), data = mydata, REML = TRUE) Zamiast załączać …
Podane trzy zmienne, yi x, które są pozytywne i ciągły z, co jest kategoryczny, mam dwa modele kandydujące dana przez: fit.me <- lmer( y ~ 1 + x + ( 1 + x | factor(z) ) ) i fit.fe <- lm( y ~ 1 + x ) Mam nadzieję porównać …
Potrzeby jeden wzór, aby określić za szkolenie wielopoziomowego modelu (korzystając lmerz lme4 Rbiblioteki) zawsze mnie trafia. Czytałem niezliczone podręczniki i samouczki, ale nigdy nie zrozumiałem tego poprawnie. Oto przykład z tego samouczka , który chciałbym zobaczyć sformułowany w równaniu. Staramy się modelować częstotliwość głosu jako funkcję płci (kobiety mają wyższy …
Analizuję wyniki eksperymentu czasu reakcji w R. Przeprowadziłem ANOVA z powtarzanymi pomiarami (1 czynnik wewnątrz podmiotu z 2 poziomami i 1 czynnik między podmiotem z 2 poziomami). Uruchomiłem podobny liniowy model mieszany i chciałem podsumować wyniki lmera w formie tabeli ANOVA lmerTest::anova. Nie zrozumcie mnie źle: nie spodziewałem się identycznych …
Dla prostego przykładu załóżmy, że istnieją dwa modele regresji liniowej 1 Model posiada trzy czynniki prognostyczne, x1a, x2b, ix2c Model 2 ma trzy predyktory z modelu 1 i dwa dodatkowe predyktory x2aorazx2b Istnieje równanie regresji populacji, w którym wyjaśniona wariancja populacji wynosi ρ2(1)ρ(1)2\rho^2_{(1)} dla Modelu 1 i ρ2(2)ρ(2)2\rho^2_{(2)} dla Modelu …
Używamy plików cookie i innych technologii śledzenia w celu poprawy komfortu przeglądania naszej witryny, aby wyświetlać spersonalizowane treści i ukierunkowane reklamy, analizować ruch w naszej witrynie, i zrozumieć, skąd pochodzą nasi goście.
Kontynuując, wyrażasz zgodę na korzystanie z plików cookie i innych technologii śledzenia oraz potwierdzasz, że masz co najmniej 16 lat lub zgodę rodzica lub opiekuna.