[R] Symulacja Monte Carlo

[secter]

Dzień dobry!
Mam taki problem. Za pomocą symulacji Monte Carlo mam obliczyć prawdopodobieństwo tego, że suma oczek wyrzuconych dwóch rzetelnych kostek wynosi co najmniej \(\displaystyle{ 5}\) i jest mniejsza od \(\displaystyle{ 10}\). Potrafię wykonać taką symulację tylko z jednym warunkiem np:

Kod: Zaznacz cały

zdarzsprz = function(k){ apply(matrix(sample(1:6, 2*k, replace=TRUE), nrow=2), 2, sum)>=5}

Jak do tego dodać warunek, że suma ma być również mniejsza od \(\displaystyle{ 10}\)?

[Igor V]

Np tak:

Kod: Zaznacz cały

zdarzsprz = function(k){ 
    x <- matrix(sample(1:6, 2*k, replace=TRUE), nrow=2)
    apply(x, 2, sum) >=5 & apply(x, 2, sum) < 10
}