Dlaczego generatory liczb losowych, takie jak runif()w R, nie generują za każdym razem tego samego wyniku? Na przykład: X <- runif(100) X generuje różne wyniki za każdym razem. Jaki jest powód generowania różnych wyników za każdym razem? Jakie funkcje działa w tle, aby to zrobić?
Mega Millions to dziś ponad 500 milionów dolarów. Pamiętam, jak czytałem artykuł JSTOR o niektórych liczbach, których wybranie jest mało prawdopodobne. Na przykład wiele osób wybiera 7, ponieważ to ich szczęśliwa liczba, a ja chcę odwrotnie. Jednak skończyło się moje członkostwo w JSTOR. Które liczby będą najmniej wybierane przez ludzi …
Chciałbym wygenerować losową macierz korelacji, tak aby rozkład jej elementów nie przekątnych wyglądał w przybliżeniu normalnie. Jak mogę to zrobić? Motywacja jest taka. Dla zestawu danych szeregów czasowych rozkład korelacji często wygląda dość normalnie. Chciałbym wygenerować wiele „normalnych” macierzy korelacji w celu przedstawienia ogólnej sytuacji i użyć ich do obliczenia …
Wysłałem poprzednie pytanie , jest to powiązane, ale myślę, że lepiej jest rozpocząć inny wątek. Tym razem zastanawiam się, jak wygenerować równomiernie rozmieszczone punkty w sferze jednostki 3-d oraz jak sprawdzić rozkład wizualnie i statystycznie? Nie widzę strategii tam zamieszczonych, które można by bezpośrednio przenieść na tę sytuację.
Zastanawiam się, jak mogę zasymulować próbkę n okresów życia rozkładu Weibulla, które obejmują obserwacje z cenzurą po prawej stronie typu I. Na przykład niech n = 3, kształt = 3, skala = 1, a szybkość cenzury = 0,15, a czas cenzury = 0,88. Wiem, jak wygenerować próbkę Weibulla, ale nie …
W tym dokumencie , który dotyczy polecenia „ustaw ziarno”, ludzie Stata omawiają kwestie związane z ustawieniem nasion podczas generowania liczb pseudolosowych. Znamienne „nie” to „nie używaj szeregowo liczb naturalnych jako nasion, ponieważ ma to wzór i zagraża pseudolosowości”. Jedyną jedną czwartą żartobliwie godną uwagi „do” jest ustawienie tylko jednego ziarna …
Muszę narysować liczby losowe z rozkładu log-cauchy'ego, który ma gęstość: Czy ktoś może mi pomóc lub wskazać mi książkę / artykuł, który mógłby mi pokazać, w jaki sposób?fa( x ; μ , σ) = 1x πσ[ 1 + ( l n ( x ) - μσ)2)].fa(x;μ,σ)=1xπσ[1+(ln(x)-μσ)2)].f(x;\mu,\sigma)=\frac{1}{x\pi\sigma\left[1+\left(\frac{ln(x)-\mu}{\sigma}\right)^2\right]}.
Jak powinienem efektywnie próbkować z następującej dystrybucji? x∼B(α,β), x>kx∼B(α,β), x>k x \sim B(\alpha, \beta),\space x > k Jeśli kkk nie jest zbyt duże, próbowanie odrzucenia może być najlepszym podejściem, ale nie jestem pewien, jak postępować, gdy kkk jest duże. Być może istnieje jakieś asymptotyczne przybliżenie, które można zastosować?
Jeśli X∈Rn, X∼N(0–,σ2I)X∈Rn, X∼N(0_,σ2I)X\in\mathbb{R}^n,~X\sim \mathcal{N}(\underline{0},\sigma^2\mathbf{I}) tj. fX(x)=1(2πσ2)n/2exp(−||x||22σ2)fX(x)=1(2πσ2)n/2exp(−||x||22σ2) f_X(x) = \frac{1}{{(2\pi\sigma^2)}^{n/2}} \exp\left(-\frac{||x||^2}{2\sigma^2}\right) Chcę analogicznej wersji skróconego rozkładu normalnego w przypadku wielu odmian. Dokładniej, chcę wygenerować ograniczony normą (do wartości ≥a≥a\geq a ) wielowymiarowy Gaussowski YYY st fY(y)={c.fX(y), if ||y||≥a0, otherwise .fY(y)={c.fX(y), if ||y||≥a0, otherwise . f_Y(y) = \begin{cases} c.f_X(y), \text{ if …
Biorąc pod uwagę monetę o nieznanym nastawieniu , jak mogę generować zmienne - tak wydajnie, jak to możliwe - które są dystrybuowane według Bernoulliego z prawdopodobieństwem 0,5? Oznacza to, że użycie minimalnej liczby rzutów na wygenerowaną zmienną jest zmienne.ppp
Mam GLMM w postaci: lmer(present? ~ factor1 + factor2 + continuous + factor1*continuous + (1 | factor3), family=binomial) Kiedy używam drop1(model, test="Chi"), otrzymuję inne wyniki niż w przypadku korzystania Anova(model, type="III")z pakietu samochodowego lub summary(model). Te dwa ostatnie dają te same odpowiedzi. Korzystając z wielu sfabrykowanych danych, odkryłem, że te …
Muszę utworzyć losowe wektory liczb rzeczywistych, spełniające następujące ograniczenia: abs(a_i) < c_i; sum(a_i)< A; # sum of elements smaller than A sum(b_i * a_i) < B; # weighted sum is smaller than B aT*A*a < D # quadratic multiplication with A smaller than D where c_i, b_i, A, B, D …
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 . Mój cel: Chciałbym mieć funkcję, która pobiera adres e-mail i wyświetla quasi-losową liczbę 1, 2, 3 lub 4. Mały szczegół: Przez …
Korzystam z symulacji na R i klastrze komputerów i mam następujący problem. Na każdym z X komputerów uruchamiam: fxT2 <- function(i) runif(10) nessay <- 100 c(mclapply(1:nessay, fxT2), recursive=TRUE) Jest 32 komputerów, każdy z 16 rdzeniami. Jednak około 2% liczb losowych jest identycznych. Jakie strategie zastosowałbyś, aby tego uniknąć? Udało mi …
W moim programie muszę uruchomić N osobnych wątków, każdy z własnym RNG, który służy do próbkowania dużego zestawu danych. Muszę być w stanie zaszczepić cały ten proces jedną wartością, aby móc odtwarzać wyniki. Czy wystarczy po prostu sekwencyjnie zwiększać ziarno dla każdego indeksu? Obecnie używam numpytych, RandomStatektóre korzystają z generatora …
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.