Współpracownik analizuje niektóre dane biologiczne pod kątem swojej pracy doktorskiej z pewną nieprzyjemną heteroscedastycznością (rysunek poniżej). Analizuje to za pomocą modelu mieszanego, ale nadal ma problemy z resztkami. Logarytmiczne przekształcanie zmiennych odpowiedzi czyści wszystko i na...
Wiemy, że sparowany test t jest tylko specjalnym przypadkiem jednokierunkowej ANOVA z powtarzanymi pomiarami (lub wewnątrz podmiotu), a także liniowym modelem mieszanego efektu, który można zademonstrować za pomocą funkcji lme () pakiet nlme w języku R jak pokazano niżej. #response data from 10...
Użyłem wielokrotnej imputacji, aby uzyskać liczbę kompletnych zestawów danych. Użyłem metod bayesowskich na każdym z kompletnych zestawów danych, aby uzyskać rozkłady tylne dla parametru (efekt losowy). Jak mogę połączyć / połączyć wyniki dla tego parametru? Więcej kontekstu: Mój model jest...
Jeśli f1,…,fkf1,…,fkf_1,\ldots,f_k są znanymi gęstościami, z których mogę symulować, tj. Dla których dostępny jest algorytm. a jeśli iloczyn ∏i=1kfi(x)αiα1,…,αk>0∏i=1kfi(x)αiα1,…,αk>0\prod_{i=1}^k f_i(x)^{\alpha_i}\qquad \alpha_1,\ldots,\alpha_k>0 jest do zabudowy, istnieje ogólne podejście...
Czy przekształcenie wielowymiarowego modelu regresji liniowej jako wielokrotnej regresji liniowej jest całkowicie równoważne? Ja nie odnosząc się po prostu działa ttt oddzielnych regresji. Przeczytałem o tym w kilku miejscach (Bayesian Data Analysis - Gelman i wsp. Oraz Multivariate Old School -...
Używam R i analizowałem moje dane za pomocą GLM z łączem dwumianowym. Chcę wiedzieć, jakie jest znaczenie przecięcia w tabeli wyjściowej. Punkt przecięcia dla jednego z moich modeli jest znacząco inny, jednak zmienna tak nie jest. Co to znaczy? Co to jest przechwycenie. Nie wiem, czy po prostu...
Mam pewne dane w [0,1], które chciałbym przeanalizować za pomocą regresji beta. Oczywiście należy coś zrobić, aby uwzględnić wartości 0,1. Nie lubię modyfikować danych, aby pasowały do modelu. również nie uważam, aby inflacja zero i 1 była dobrym pomysłem, ponieważ uważam, że w tym przypadku...
To wydaje się takie elementarne, ale zawsze utknąłem w tym momencie… Większość danych, z którymi mam do czynienia, jest nienormalna, a większość analiz opartych na strukturze GLM. Do mojej obecnej analizy mam zmienną odpowiedzi, która jest „prędkością marszu” (metry / minutę). Łatwo jest mi...
Rozważ następujące dane z dwustronnego projektowania przedmiotów: df <- "http://personality-project.org/r/datasets/R.appendix4.data" df <- read.table(df,header=T) head(df) Observation Subject Task Valence Recall 1 1 Jim Free Neg 8 2 2 Jim Free Neu 9 3 3 Jim Free Pos 5 4 4 Jim Cued Neg 7 5 5...
Próbuję użyć lmez nlmepakietu do replikacji wyników aovdla ANOVA z powtarzanymi pomiarami. Zrobiłem to dla eksperymentu z powtarzanymi pomiarami z jednym czynnikiem i dla eksperymentu z dwoma czynnikami z jednym czynnikiem między poddanymi i jednym czynnikiem wewnątrz poddanych, ale mam problem z...
Powiedzmy, że mamy GLMM mod1 <- glmer(y ~ x + A + (1|g), data = dat) mod2 <- glmer(y ~ x + B + (1|g), data = dat) Te modele nie są zagnieżdżone w zwykłym sensie: a <- glmer(y ~ x + A + (1|g), data = dat) b <- glmer(y ~ x + A + B + (1|g), data = dat) więc nie możemy zrobić tego,...
Mam arkusz wyników niektórych drużyn. Pierwsza drużyna z 10 punktami wygrywa. W każdej drużynie jest 2 graczy. Gracze cały czas grają z różnymi członkami drużyny, chociaż nie są wybierani idealnie losowo. Żadne indywidualne wyniki nie są przechowywane. Więc w zasadzie mamy Billa i Boba pokonali...
Potrzebuję porady dotyczącej dwóch głównych dylematów w moich badaniach, które są studium przypadku 3 dużych farmaceutyków i innowacji. Liczba patentów rocznie jest zmienną zależną. Moje pytania są Jakie są najważniejsze kryteria dobrego modelu? Co jest ważniejsze / mniej ważne? Czy to, że...
Załóżmy, że pracujemy nad modelem efektów losowych niektórych danych zliczanych w czasie i chcemy kontrolować niektóre trendy. Zwykle zrobiłbyś coś takiego: lmer(counts ~ dependent_variable + (1+t+I(t^2)|ID), family="poisson") zawierać kwadratowy kształt dla t. Czy można zastosować bardziej...
Jest takie rozróżnienie, które wprawia mnie w zakłopotanie w przypadku modeli mieszanych i zastanawiam się, czy mógłbym uzyskać trochę jasności. Załóżmy, że masz mieszany model danych zliczania. Istnieje zmienna, o której wiesz, że chcesz jako efekt stały (A), oraz inna zmienna czasu (T),...
Zasadniczo zastanawiam się, w jaki sposób wymuszane są różne struktury kowariancji i jak obliczane są wartości w tych macierzach. Funkcje takie jak lme () pozwalają nam wybrać, którą strukturę chcielibyśmy, ale chciałbym wiedzieć, jak są szacowane. Rozważ liniowy model efektów mieszanych...
Próbuję stworzyć model, dla którego mam zmienną odpowiedzi, która jest proporcją między 0 a 1, obejmuje to całkiem sporo zer i 1, ale także wiele wartości pomiędzy nimi. Myślę o próbie regresji beta. Pakiet, który znalazłem dla R (betareg), dopuszcza tylko wartości z zakresu od 0 do 1, ale bez...
Jeśli naprawię wartości zaobserwowanych węzłów MRF, czy stanie się to
Interesuję się uczeniem się, jak oceniać i klasyfikować osoby w grupie, które wchodzą w interakcje / rywalizują tylko w parach (tj. Systemy takie jak system oceny ELO dla szachów). Czy są dostępne jakieś metody przejścia lub bardziej dokładne i zaawansowane metody? Czy są jakieś pakiety R, które...
Mam wykres wartości resztkowych modelu liniowego w funkcji dopasowanych wartości, w których heteroscedastyczność jest bardzo wyraźna. Jednak nie jestem pewien, jak powinienem postępować teraz, ponieważ o ile rozumiem ta heteroscedastyczność powoduje, że mój model liniowy jest nieważny. (Czy to...