Próbkowanie z wymianą ma dwie zalety w porównaniu z próbkowaniem bez wymiany, tak jak to widzę: 1) Nie musisz się martwić o skończoną korektę populacji. 2) Istnieje szansa, że elementy z populacji zostaną narysowane wiele razy - wtedy możesz przetworzyć pomiary i zaoszczędzić czas. Oczywiście...
Jeśli wszystko, co robisz, to ponowne próbkowanie z rozkładu empirycznego, dlaczego nie po prostu przestudiować rozkład empiryczny? Na przykład zamiast badać zmienność poprzez powtarzanie prób, dlaczego nie po prostu skwantyfikować zmienność z rozkładu
Mam prosty problem z próbkowaniem, w którym moja wewnętrzna pętla wygląda następująco: v = sample_gamma(k, a) gdzie sample_gammapróbki z rozkładu gamma tworzą próbkę Dirichleta. Działa dobrze, ale w przypadku niektórych wartości k / a niektóre z niższych obliczeń są niedopełnione. Dostosowałem...
Chcę próbkować z gęstości jednowymiarowej ale znam tylko związek:faXfXf_X faX( x ) = ∫faX| Y( x | y) fY( y) dy.fX(x)=∫fX|Y(x|y)fY(y)dy.f_X(x) = \int f_{X\vert Y}(x\vert y)f_Y(y) dy. Chcę uniknąć używania MCMC (bezpośrednio na reprezentacji całkowej), a ponieważ i są łatwe do próbkowania, myślałem...
Mam ważoną próbkę, dla której chcę obliczyć kwantyle. 1 Najlepiej, przy czym masy są takie same (zarówno = 1 lub inaczej), wyniki mogą być zgodne z tymi, scipy.stats.scoreatpercentile()i R: quantile(...,type=7). Jednym prostym podejściem byłoby „pomnożenie” próbki przy użyciu podanych wag. To...
Czy w R, jeśli ustawię set.seed (), a następnie użyję przykładowej funkcji do losowej listy, czy mogę zagwarantować, że nie wygeneruję tej samej permutacji? to znaczy... set.seed(25) limit <- 3 myindex <- seq(0,limit) for (x in seq(1,factorial(limit))) { permutations <- sample(myindex)...
Uwaga: to pytanie jest repost, ponieważ moje poprzednie pytanie musiało zostać usunięte ze względów prawnych. Porównując PROC MIXED z SAS z funkcją lmez nlmepakietu w R, natknąłem się na pewne dość mylące różnice. Mówiąc dokładniej, stopnie swobody w różnych testach różnią się między PROC MIXEDi...
Obserwuję czasy przetwarzania procesu przed i po zmianie, aby dowiedzieć się, czy proces poprawił się przez zmianę. Proces poprawił się, jeśli czas przetwarzania został skrócony. Rozkład czasu przetwarzania jest gruby, dlatego porównanie na podstawie średniej nie jest rozsądne. Zamiast tego...
Zawsze uczono mnie, że CLT działa, gdy powtarzasz próbkowanie, a każda próbka jest wystarczająco duża. Wyobraź sobie na przykład, że mam kraj 1 000 000 obywateli. Rozumiem, że CLT jest taki, że nawet jeśli rozkład ich wysokości nie był normalny, gdybym pobrał 1000 próbek 50 osób (tj. Przeprowadził...
Moje wykształcenie to informatyka. Jestem całkiem nowy w metodach próbkowania Monte Carlo i chociaż rozumiem matematykę, trudno mi znaleźć intuicyjne przykłady ważnego próbkowania. Dokładniej, czy ktoś mógłby podać przykłady: oryginalny rozkład, z którego nie można próbkować, ale można go...
W metodach MCMC wciąż czytam o burn-inczasie lub liczbie próbek "burn". Co to dokładnie jest i dlaczego jest potrzebne? Aktualizacja: Czy po stabilizacji MCMC pozostaje stabilny? W jaki sposób pojęcie burn-inczasu jest powiązane z czasem
Testy takie jak Z, t i kilka innych zakładają, że dane są oparte na losowym próbkowaniu. Dlaczego? Załóżmy, że prowadzę badania eksperymentalne, w których bardziej zależy mi na wewnętrznej ważności niż na zewnętrznej. Tak więc, jeśli moja próbka może być nieco stronnicza, to dobrze, ponieważ...
Powiedzmy, że mamy zmienną losową z zakresu wartości ograniczonego a , gdzie jest wartością minimalną i wartości maksymalnej.aaabbbaaabbb Powiedziano mi, że jako , gdzie jest naszą wielkością próby, rozkład próbkowania naszych średnich próbek jest rozkładem normalnym. Oznacza to, że wraz ze...
Zastanawiałem się, jak CI bootstrap (i BCa w układzie dwubiegunowym) działają na normalnie dystrybuowanych danych. Wydaje się, że dużo pracy analizuje ich wydajność w różnych typach dystrybucji, ale nie można znaleźć niczego w normalnie dystrybuowanych danych. Ponieważ najpierw wydaje się rzeczą...
Mam bardzo duży zestaw danych i brakuje około 5% wartości losowych. Te zmienne są ze sobą skorelowane. Poniższy przykładowy zestaw danych R jest tylko zabawkowym przykładem z fałszywymi skorelowanymi danymi. set.seed(123) # matrix of X variable xmat <- matrix(sample(-1:1, 2000000, replace =...
Moje dane są następujące. Mam dwie grupy pacjentów. Pacjenci w każdej grupie mieli inny rodzaj operacji oka. 5 zmiennych zmierzono u pacjentów w każdej grupie. Chcę porównać te zmienne między dwiema grupami za pomocą testu permutacji lub MANOVA. Oko, na którym przeprowadzono operację, nie ma tak...
Próbuję skonstruować dowód na problem, nad którym pracuję, a jednym z założeń, które robię, jest to, że zbiór punktów, z których próbuję, jest gęsty na całej przestrzeni. Praktycznie używam łacińskiego próbkowania hipersześcianu, aby uzyskać punkty w całej przestrzeni próbki. Chciałbym wiedzieć,...
Implementacja randomForest nie pozwala na pobieranie próbek poza liczbę obserwacji, nawet w przypadku pobierania próbek z wymianą. Dlaczego to? Działa w porządku: rf <- randomForest(Species ~ ., iris, sampsize=c(1, 1, 1), replace=TRUE) rf <- randomForest(Species ~ ., iris, sampsize=3,...
Ciągle czytam o potrzebie sprawdzenia autokorelacji w MCMC. Dlaczego ważne jest, aby autokorelacja była niska? Co mierzy w kontekście
Jestem kompletnym nowicjuszem :) Robię badanie na próbie 10 000 z populacji około 745,000. Każda próbka reprezentuje „procentowe podobieństwo”. Zdecydowana większość próbek wynosi około 97% -98%, ale kilka ma od 60% do 90%, co oznacza, że rozkład jest mocno negatywnie wypaczony. Około 0,6%...