Za chwilę wypróbuję środowisko w stylu BUGS do szacowania modeli bayesowskich. Czy są jakieś ważne zalety, które należy wziąć pod uwagę, wybierając pomiędzy OpenBugs lub JAGS? Czy jedna może zastąpić drugą w dającej się przewidzieć przyszłości? Będę używać wybranego Gibbsa Samplera z R. Nie mam jeszcze konkretnej aplikacji, ale raczej …
Znalazłem kilka rozkładów, dla których BŁĘDY i R mają różne parametryzacje: Normalny, log-Normalny i Weibull. Dla każdego z nich zbieram, że drugi parametr użyty przez R musi zostać odwrócony (1 / parametr) przed użyciem w BŁĘDACH (lub JAGS w moim przypadku). Czy ktoś wie o wyczerpującej liście tych transformacji, która …
Wyjaśnię mój problem na przykładzie. Załóżmy, że chcesz przewidzieć dochód danej osoby na podstawie niektórych atrybutów: {Wiek, płeć, kraj, region, miasto}. Masz taki zestaw danych szkoleniowych train <- data.frame(CountryID=c(1,1,1,1, 2,2,2,2, 3,3,3,3), RegionID=c(1,1,1,2, 3,3,4,4, 5,5,5,5), CityID=c(1,1,2,3, 4,5,6,6, 7,7,7,8), Age=c(23,48,62,63, 25,41,45,19, 37,41,31,50), Gender=factor(c("M","F","M","F", "M","F","M","F", "F","F","F","M")), Income=c(31,42,71,65, 50,51,101,38, 47,50,55,23)) train CountryID RegionID CityID …
Witaj guru statystyczni i kreatorzy programowania R, Interesuje mnie modelowanie chwytów zwierząt jako funkcji warunków środowiskowych i dnia w roku. W ramach innego badania mam liczbę przechwyceń przez ~ 160 dni w ciągu trzech lat. Na każdy z tych dni mam temperaturę, opady, prędkość wiatru, wilgotność względną itp. Ponieważ dane …
Zastanawiałem się, jaki pakiet statystyczny oprogramowania polecacie do przeprowadzania wnioskowania bayesowskiego. Na przykład wiem, że możesz uruchamiać openBUGS lub winBUGS jako standalones lub też możesz wywoływać je z R. Ale R ma również kilka własnych pakietów (MCMCPack, BACCO), które mogą wykonywać analizę bayesowską. Czy ktoś ma jakieś sugestie co do …
Zamknięte. To pytanie jest nie na temat . Obecnie nie przyjmuje odpowiedzi. Chcesz poprawić to pytanie? Zaktualizuj pytanie, aby było tematem dotyczącym weryfikacji krzyżowej. Zamknięte w zeszłym roku . Podążam kursem statystyki bayesowskiej używając BŁĘDÓW i R. Teraz już znam BŁĘDY, to świetnie, ale tak naprawdę nie lubię używać osobnego …
Mam zestaw danych szeregów czasowych, do którego próbuję dopasować ukryty model Markowa (HMM) w celu oszacowania liczby stanów ukrytych w danych. Mój pseudo-kod do tego jest następujący: for( i in 2 : max_number_of_states ){ ... calculate HMM with i states ... optimal_number_of_states = "model with smallest BIC" ... } Teraz, …
Wydaje się, że pełne warunki warunkowe są często dość trudne do uzyskania, ale programy takie jak JAGS i BŁĘDY wyprowadzają je automatycznie. Czy ktoś może wyjaśnić, w jaki sposób algorytmicznie generuje pełne warunki warunkowe dla dowolnej specyfikacji modelu?
W R„uprzedniej wadze” możemy glmregresję za pomocą parametru wag . Na przykład: glm.D93 <- glm(counts ~ outcome + treatment, family = poisson(), weights=w) Jak można tego dokonać w modelu JAGSlub BUGS? Znalazłem jakiś artykuł na ten temat, ale żaden z nich nie stanowi przykładu. Interesują mnie głównie przykłady Poissona i …
Próbuję zbudować model, w którym odpowiedź jest proporcjonalna (w rzeczywistości jest to liczba głosów, jaką partia zdobywa w okręgach wyborczych). Jego rozkład nie jest normalny, więc postanowiłem modelować go z rozkładem beta. Mam też kilka predyktorów. Nie wiem jednak, jak napisać to w BŁĘDACH / JAGACH / STANU (JAGS byłby …
Próbuję dopasować model czasu dyskretnego do R, ale nie jestem pewien, jak to zrobić. Czytałem, że możesz zorganizować zmienną zależną w różnych wierszach, po jednym dla każdej obserwacji czasu, i użyć glmfunkcji z łączem logit lub cloglog. W tym sensie, mam trzy kolumny: ID, Event(1 lub 0, w każdym okresie …
Co dzieje się, gdy nie masz pojęcia o rozkładzie parametrów? Jakie podejście powinniśmy zastosować? Przez większość czasu dążymy do podkreślenia, czy dana zmienna ma jakikolwiek wpływ na obecność / nieobecność określonego gatunku, a zmienna jest akceptowana lub nie, zgodnie ze znaczeniem zmiennej. Oznacza to, że przez większość czasu nie myślimy …
Jestem nowym użytkownikiem WinBUGS i mam jedno pytanie do twojej pomocy. Po uruchomieniu następującego kodu uzyskałem parametry beta0through beta4(statystyki, gęstość), ale nie wiem, jak uzyskać prognozę ostatniej wartości h, którą ustawiłem NAdo modelowania w kodzie. Czy ktoś może mi podpowiedzieć? Wszelkie porady będą mile widziane. model { for(i in 1: …
Gelman i Hill (2006) mówią: W Bugs brakujące wyniki w regresji można łatwo rozwiązać, po prostu włączając wektor danych, NA i wszystkie. Błędy jawnie modelują zmienną wynikową, dlatego użycie tego modelu jest banalne, aby w efekcie przypisywać brakujące wartości przy każdej iteracji. Brzmi to jak prosty sposób na wykorzystanie JAGS …
Zamknięte. To pytanie jest nie na temat . Obecnie nie przyjmuje odpowiedzi. Chcesz poprawić to pytanie? Zaktualizuj pytanie, aby było tematem dotyczącym weryfikacji krzyżowej. Zamknięte 2 lata temu . R wydaje się być w stanie generować ładne wykresy podsumowań z obiektów bugsi jagsgenerowanych przez funkcje R2WinBUGS :: bugs i R2jags: …
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.