Modele mieszane (inaczej wielopoziomowe lub hierarchiczne) to modele liniowe, które obejmują zarówno efekty stałe, jak i efekty losowe. Służą do modelowania danych podłużnych lub zagnieżdżonych.
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 …
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ć …
Zanim zadałem to pytanie, przeszukałem naszą stronę i znalazłem wiele podobnych pytań (jak tutaj , tutaj i tutaj ). Ale wydaje mi się, że na te powiązane pytania nie udzielono odpowiedzi lub nie omówiono ich, dlatego chciałbym ponownie zadać to pytanie. Uważam, że powinna istnieć duża liczba odbiorców, którzy chcieliby, …
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 …
W modelu mieszanym zakładamy, że efekty losowe (parametry) są zmiennymi losowymi, które mają rozkład normalny. Wygląda bardzo podobnie do metody bayesowskiej, w której zakłada się, że wszystkie parametry są losowe. Czy zatem model efektu losowego jest rodzajem specjalnego przypadku metody bayesowskiej?
Zawsze uczono mnie, że efekty losowe wpływają tylko na wariancję (błąd), a efekty stałe wpływają tylko na średnią. Ale znalazłem przykład, w którym efekty losowe wpływają również na średnią - oszacowanie współczynnika: require(nlme) set.seed(128) n <- 100 k <- 5 cat <- as.factor(rep(1:k, each = n)) cat_i <- 1:k # …
Istnieje kilka wątków na temat R-sig-ME na temat uzyskiwania przedziałów ufności dla prognoz przy użyciu lme4iw nlmeR. Na przykład tutaj i tutaj w 2010 r., W tym komentarz Dougalsa Batesa, jednego z autorów obu pakietów. Waham się, by zacytować go dosłownie, z obawy, że zostaną wyjęci z kontekstu, ale tak …
Pracuję nad zestawem danych. Po zastosowaniu niektórych technik identyfikacji modelu, wyszłam z modelem ARIMA (0,2,1). Użyłem detectIOfunkcji w pakiecie TSAw R do wykrycia innowacyjnej wartości odstającej (IO) przy 48. obserwacji mojego oryginalnego zestawu danych. Jak włączyć tę wartość odstającą do mojego modelu, aby móc jej używać do celów prognozowania? Nie …
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.