simulating rolling two dice

Question

I am trying to simulate rolling two dice. I used:

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

and get the expe

Answer 1

I agree with David that nothing seems particularly wrong with your first option. Another way to go might be this, if you're really just after the sum of the two dice:

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