mô phỏng hai con xúc xắc

Question

Tôi đang cố mô phỏng hai con xúc xắc. Tôi đã sử dụng:

d2 <- sample(1:6, 10^6, replace = T) + sample(1:6, 10^6, replace = T) 

và nhận kết quả mong đợi. Tôi cũng đã thử

s2d <- c() 
for (i in 1:6) { 
    for (j in 1:6){ 
    s2d <- c(s2d, (i+j)) 
    } 
} 
d2 <- sample(s2d, 10^6, replace=T) 

và điều đó cũng hoạt động, nhưng chúng cảm thấy một chút "lực lượng vũ phu". Có cách nào dễ dàng hơn, thanh lịch hơn để làm điều đó không?

Nói chung, có một hàm có 2 (hoặc nhiều) sự kiện độc lập và thực hiện các thao tác trên chúng (bổ sung, nhân) không?

Answer 1

Nếu vấn đề của bạn là bạn không thể cuộn bất kỳ số lượng tùy ý các con xúc xắc, một cái gì đó như:

rowSums(replicate(2, sample(6, 10^6, replace=T))) 

sẽ được linh hoạt hơn.

Answer 2

Tôi đồng ý với David rằng dường như không có gì đặc biệt sai với tùy chọn đầu tiên của bạn. Một cách khác để đi có thể này, nếu bạn thực sự ngay sau khi tổng của hai con xúc xắc:

sample(2:12,size = 100,replace = TRUE, prob = table(outer(1:6,1:6,"+"))/36) 

Answer 3

Có một chức năng dice trong gói TeachingDemos rằng mô phỏng việc tung xúc xắc (và thậm chí còn có một là tùy chọn để vẽ kết quả, nhưng 1000 cuộn sẽ không tạo ra một cốt truyện có ý nghĩa). Điều này có vẻ hơi ít vũ lực, nhưng nội bộ nó tương tự như những gì đã được đăng. Bạn có thể sử dụng các hàm áp dụng hoặc có liên quan để thực hiện những việc như tổng hợp trên các cột của sự trở lại.