alexhallam · June 28, 2021 11:20
diff --git a/monte_carlo_sum.R b/monte_carlo_sum.R
 # calculate the sum of a, b, b2
 n <- 1000
 m <- 10
 a <- rnorm(n, 10)
 b <- rnorm(n, 20)
 b2 <- .1*a + .1*b
 c <- sample(a, size = m*n, replace=T) + sample(b, size = m*n, replace=T) + sample(b2, size = m*n, replace=T)
	# calculate the sum of a, b, b2
	n <- 1000
	m <- 10
	a <- rnorm(n, 10)
	b <- rnorm(n, 20)
	b2 <- .1a + .1b
	c <- sample(a, size = mn, replace=T) + sample(b, size = mn, replace=T) + sample(b2, size = m*n, replace=T)