RStudio - frequenza relativa

[abaco90]

Ciao a tutti,

sto approcciando a Rstudio e ho un esercizio che mi chiede di creare un'urna con 7 palline: 4 bianche e 3 rosse.
Successivamente devo calcolare la probabilità, usando la frequenza relativa, di estrarre 2 palline entrambe rosse.

Il mio codice è

urna <- c('r', 'r', 'r', 'b', 'b', 'b', 'b')
estrazione_urna <- sample(urna, 2)
table(estrazione_urna)/7

Il problema è che eseguendo il codice, quando mi capita il caso con due palline rosse, mi dà probabilità:

estrazioni
        r 
0.2857143 

mentre dovrei ottenere: $ 3/7 * 2/6 = 0.1428571 $.

Dove ho sbagliato? Grazie

[dasalv12]

Un modo può essere il seguente (ci sono anche dei pacchetti che calcolano le permutazioni).

length(urna[which(urna == "r")]) / length(urna) * 
  (length(urna[which(urna == "r")]) -1) / (length(urna) -1)

Quello che sbagli nel codice precedente è che con "sample" estrai due elementi a caso 1) senza indicare che non vi reinserimento e 2) senza indicarne le caratteristiche volute. La funzione "table" ti indica 0.28 semplicemente perché ti conta 2/7, ovvero i due elementi estratti.