simulare a rotazione due dadi

Question

Sto provando a simulare il lancio di due dadi. Ho usato:

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

e ottenere il risultato previsto. Ho anche provato

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

e che funziona anche, ma questi mi sento un po ' "forza bruta". C'è un modo più semplice ed elegante per farlo?

In termini più generali, esiste una funzione che richiede 2 (o più) eventi indipendenti e fa operazioni su di essi (aggiunta, moltiplicazione)?

Answer 1

Se il problema è che non si può rotolare qualsiasi numero arbitrario di dadi, qualcosa come:

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

sarebbe più flessibile.

Answer 2

Sono d'accordo con David sul fatto che nulla sembra particolarmente sbagliato con la prima opzione. Un altro modo per andare potrebbe essere questo, se siete davvero solo dopo la somma dei due dadi:

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

Answer 3

c'è una funzione dice nel pacchetto TeachingDemos che simula il rotolamento dei dadi (e v'è anche un opzione per tracciare i risultati, ma 1000 rotoli non renderebbero una trama significativa). Questo può sembrare un po 'meno forza bruta, ma internamente fa simile a quello che è già stato pubblicato. È possibile utilizzare le funzioni apply o correlate per eseguire operazioni come la somma tra le colonne del reso.