# Create an initial D12 die roll

single_12_roll<-sample(1:12,1)

# Concatenate to form the next 19 dice rolls

single_12_roll<-c(single_12_roll, sample(1:12,1))
single_12_roll<-c(single_12_roll, sample(1:12,1))
single_12_roll<-c(single_12_roll, sample(1:12,1))
single_12_roll<-c(single_12_roll, sample(1:12,1))
single_12_roll<-c(single_12_roll, sample(1:12,1))
single_12_roll<-c(single_12_roll, sample(1:12,1))
single_12_roll<-c(single_12_roll, sample(1:12,1))
single_12_roll<-c(single_12_roll, sample(1:12,1))
single_12_roll<-c(single_12_roll, sample(1:12,1))
single_12_roll<-c(single_12_roll, sample(1:12,1))
single_12_roll<-c(single_12_roll, sample(1:12,1))
single_12_roll<-c(single_12_roll, sample(1:12,1))
single_12_roll<-c(single_12_roll, sample(1:12,1))
single_12_roll<-c(single_12_roll, sample(1:12,1))
single_12_roll<-c(single_12_roll, sample(1:12,1))
single_12_roll<-c(single_12_roll, sample(1:12,1))
single_12_roll<-c(single_12_roll, sample(1:12,1))
single_12_roll<-c(single_12_roll, sample(1:12,1))
single_12_roll<-c(single_12_roll, sample(1:12,1))
single_12_roll<-c(single_12_roll, sample(1:12,1))

# Display D12 data

single_12_roll

# Create an initial 2D6 die roll

double_6_roll<-sum(sample(1:6,2))

# Concatenate to form the next 19 dice rolls

double_6_roll<-c(double_6_roll, sum(sample(1:6,2)))
double_6_roll<-c(double_6_roll, sum(sample(1:6,2)))
double_6_roll<-c(double_6_roll, sum(sample(1:6,2)))
double_6_roll<-c(double_6_roll, sum(sample(1:6,2)))
double_6_roll<-c(double_6_roll, sum(sample(1:6,2)))
double_6_roll<-c(double_6_roll, sum(sample(1:6,2)))
double_6_roll<-c(double_6_roll, sum(sample(1:6,2)))
double_6_roll<-c(double_6_roll, sum(sample(1:6,2)))
double_6_roll<-c(double_6_roll, sum(sample(1:6,2)))
double_6_roll<-c(double_6_roll, sum(sample(1:6,2)))
double_6_roll<-c(double_6_roll, sum(sample(1:6,2)))
double_6_roll<-c(double_6_roll, sum(sample(1:6,2)))
double_6_roll<-c(double_6_roll, sum(sample(1:6,2)))
double_6_roll<-c(double_6_roll, sum(sample(1:6,2)))
double_6_roll<-c(double_6_roll, sum(sample(1:6,2)))
double_6_roll<-c(double_6_roll, sum(sample(1:6,2)))
double_6_roll<-c(double_6_roll, sum(sample(1:6,2)))
double_6_roll<-c(double_6_roll, sum(sample(1:6,2)))
double_6_roll<-c(double_6_roll, sum(sample(1:6,2)))
                                                                                                     
# Display 2D6 data

double_6_roll

# Calculate means

mean(single_12_roll)
mean(double_6_roll)


