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.
Właśnie natknąłem się na ten artykuł , który opisuje, jak obliczyć powtarzalność (aka niezawodność, aka korelacja wewnątrzklasowa) pomiaru za pomocą modelowania efektów mieszanych. Kod R byłby następujący: #fit the model fit = lmer(dv~(1|unit),data=my_data) #obtain the variance estimates vc = VarCorr(fit) residual_var = attr(vc,'sc')^2 intercept_var = attr(vc$id,'stddev')[1]^2 #compute the unadjusted repeatability …
W modelu wielopoziomowym, jakie są praktyczne i związane z interpretacją implikacje oszacowania w porównaniu z niedoszacowaniem parametrów korelacji efektu losowego? Praktycznym powodem pytania jest to, że w ramce Lmer w R nie ma zaimplementowanej metody szacowania wartości p za pomocą technik MCMC, gdy dokonuje się szacunków w modelu korelacji między …
Ponieważ ogólną zgodą wydaje się być stosowanie modeli mieszanych lmer()w R zamiast klasycznej ANOVA (z często cytowanych powodów, takich jak niezrównoważone projekty, skrzyżowane efekty losowe itp.), Chciałbym spróbować z moimi danymi. Martwię się jednak, że będę w stanie „sprzedać” to podejście mojemu przełożonemu (który w końcu oczekuje klasycznej analizy z …
Problem: Przeczytałem w innych postach, które predictnie są dostępne dla lmermodeli z efektami mieszanymi {lme4} w [R]. Próbowałem zgłębić ten temat za pomocą zestawu danych o zabawkach ... Tło: Zestaw danych jest dostosowany z tego źródła i dostępny jako ... require(gsheet) data <- read.csv(text = gsheet2text('https://docs.google.com/spreadsheets/d/1QgtDcGJebyfW7TJsB8n6rAmsyAnlz1xkT3RuPFICTdk/edit?usp=sharing', format ='csv')) Oto pierwsze …
Korzystam z GEE na 3-poziomowych niezrównoważonych danych, używając łącza logit. Czym to się różni (pod względem wniosków, które mogę wyciągnąć i znaczenia współczynników) od GLM z efektami mieszanymi (GLMM) i linkiem logit? Więcej szczegółów: Obserwacje są pojedynczymi próbami bernoulli. Są one pogrupowane w klasy i szkoły. Za pomocą R. Przypadkowe …
Mam następujące dane wyjściowe: Generalized linear mixed model fit by the Laplace approximation Formula: aph.remain ~ sMFS2 +sAG2 +sSHDI2 +sbare +season +crop +(1|landscape) AIC BIC logLik deviance 4062 4093 -2022 4044 Random effects: Groups Name Variance Std.Dev. landscape (Intercept) 0.82453 0.90804 Number of obs: 239, groups: landscape, 45 Fixed effects: …
Używam modelu mieszanego w R( lme4) do analizy niektórych danych z powtarzanymi pomiarami. Mam zmienną odpowiedzi (zawartość włókna w kale) i 3 stałe efekty (masa ciała itp.). Moje badanie ma tylko 6 uczestników, z 16 powtarzanymi pomiarami dla każdego z nich (chociaż dwóch ma tylko 12 powtórzeń). Podmiotami są jaszczurki, …
Przeszukałem wiele stron pomocy i wciąż nie jestem pewien, jak określić bardziej skomplikowane terminy zagnieżdżone w modelu mieszanym. Jestem również mylić jak wykorzystanie :i /oraz |w określaniu i gniazdowania interakcje z czynnikami losowymi wykorzystujących lmer()w lme4opakowaniu w R. Na potrzeby tego pytania załóżmy, że dokładnie przedstawiłem moje dane za pomocą …
W Internecie wiele znalazłem na temat interpretacji efektów losowych i stałych. Jednak nie udało mi się uzyskać źródła określającego: Jaka jest matematyczna różnica między efektami losowymi i stałymi? Rozumiem przez to matematyczne sformułowanie modelu i sposób szacowania parametrów.
Cytat blokowy poniżej, od liderów w dziedzinie modelowania efektów mieszanych, twierdzi, że koordynacja przesunięć w modelach z zerową korelacją między efektami losowymi (modele „ZCP”) zmienia prognozy modeli. Ale czy ktoś może rozwinąć lub uzasadnić swoje roszczenia? Sprawozdania te zostały wydane z Bates wsp za 2015 na papierze lme4, Montaż liniowe …
lmerTestPakiet zawiera anova()funkcję liniowe modele mieszane z ewentualnie zbliżenia Satterthwaite'a (domyślnie) Kenward-Roger ze stopniami swobody (df). Jaka jest różnica między tymi dwoma podejściami? Kiedy wybrać, który?
Obecnie korzystam z modeli liniowych z mieszanym efektem. Korzystam z pakietu „lme4” w języku R. Moje modele mają postać: model <- lmer(response ~ predictor1 + predictor2 + (1 | random effect)) Przed uruchomieniem moich modeli sprawdziłem możliwą wielokoliniowość między predyktorami. Zrobiłem to przez: Utwórz ramkę danych predyktorów dummy_df <- data.frame(predictor1, …
Przeprowadziłem powtarzający się projekt, w którym przetestowałem 30 mężczyzn i 30 kobiet w trzech różnych zadaniach. Chcę zrozumieć, jak różni się zachowanie mężczyzn i kobiet i jak to zależy od zadania. Użyłem zarówno pakietu lmer, jak i lme4, aby to zbadać, jednak utknąłem przy próbie sprawdzenia założeń dla każdej z …
Jestem trochę zdezorientowany co do zalet mieszanych modeli w zakresie modelowania predykcyjnego. Ponieważ modele predykcyjne mają zwykle przewidywać wartości wcześniej nieznanych obserwacji, wydaje mi się oczywiste, że jedynym sposobem, w jaki model mieszany może być użyteczny, jest jego zdolność do przewidywania na poziomie populacji (to znaczy bez dodawania żadnych efektów …
Mam problem z uznaniem korzyści oznaczania czynnika modelowego za losowy z kilku powodów. Wydaje mi się, że prawie we wszystkich przypadkach optymalnym rozwiązaniem jest traktowanie wszystkich czynników jako ustalonych. Po pierwsze, rozróżnienie między ustalonym a losowym jest dość arbitralne. Standardowe wyjaśnienie jest takie, że jeśli ktoś interesuje się konkretnymi jednostkami …
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.