Sinh số giả ngẫu nhiên

Lấy mẫu trong R

Richie Cotton

Data Evangelist at DataCamp

“Ngẫu nhiên” nghĩa là gì?

{adjective} được tạo ra, thực hiện, xảy ra hoặc chọn mà không theo phương pháp hay quyết định có ý thức.

  • Oxford Languages
Lấy mẫu trong R

Số ngẫu nhiên thật

  • Sinh từ quá trình vật lý, như tung đồng xu.
  • Hotbits dùng phân rã phóng xạ.
  • RANDOM.ORG dùng nhiễu khí quyển.
    • Có sẵn trong R qua gói random.
  • Ngẫu nhiên thật sự tốn kém.
1 https://www.fourmilab.ch/hotbits 2 https://www.random.org
Lấy mẫu trong R

Sinh số giả ngẫu nhiên

  • Số “ngẫu nhiên” kế tiếp được tính từ số “ngẫu nhiên” trước đó.
  • Số “ngẫu nhiên” đầu tiên được tính từ một seed.
  • Nếu bắt đầu cùng seed, mọi số ngẫu nhiên về sau sẽ giống hệt.
seed <- 1
calc_next_random(seed)

3

calc_next_random(3)

2

calc_next_random(2)

6
Lấy mẫu trong R

Hàm sinh số ngẫu nhiên

function distribution function distribution function distribution
rbeta Beta rgeom Geometric rsignrank Wilcoxon signed rank
rbinom Binomial rhyper Hypergeometric rt t
rcauchy Cauchy rlnorm Lognormal runif Uniform
rchisq Chi-squared rlogis Logistic rweibull Weibull
rexp Exponential rnbinom Negative binomial rwilcox Wilcoxon rank sum
rf F rnorm Normal
rgamma Gamma rpois Poisson
Lấy mẫu trong R

Trực quan hóa số ngẫu nhiên

rbeta(5000, shape1 = 2, shape2 = 2)

[1] 0.2788 0.7495 0.6485 0.6665 0.6546 0.1575
...

[4996] 0.84719 0.35177 0.92796 0.67603 0.53960

randoms <- data.frame(
  beta = rbeta(5000, shape1 = 2, shape2 = 2)
)

ggplot(randoms, aes(beta)) +
  geom_histogram(binwidth = 0.1)

hist-beta.png

Lấy mẫu trong R

Hạt giống (seed) cho số ngẫu nhiên

set.seed(20000229)

rnorm(5)

-1.6538 -0.4028 -0.1654 -0.0734  0.5171

rnorm(5)

1.908  0.379 -1.499  1.625  0.693

set.seed(20000229)

rnorm(5)

-1.6538 -0.4028 -0.1654 -0.0734  0.5171

rnorm(5)

1.908  0.379 -1.499  1.625  0.693
Lấy mẫu trong R

Dùng seed khác

set.seed(20000229)

rnorm(5)

-1.6538 -0.4028 -0.1654 -0.0734  0.5171

rnorm(5)

1.908  0.379 -1.499  1.625  0.693

set.seed(20041004)

rnorm(5)

-0.6547 -0.7854 -0.0152  0.1514  0.5285

rnorm(5)

0.748  0.974  0.174 -0.781 -0.930
Lấy mẫu trong R

¡Vamos a practicar!

Lấy mẫu trong R

Preparing Video For Download...