Definiciones

Probabilidad: Es el estudio de azar y la incertidumbre en cualquier situacion en la cual varios posibles sucesos pueden ocurrir. Es un valor entre 0 (imposible) y 1 (seguro) Ejemplo: la probabilidad de que llueva hoy es de 0.70 (70%).

Experimento: Cualquier acción cuyo resultado está sujeto a la incertidumbre. Ejemplo: Lanzar una moneda al aire.

Instalar libreria y llamar el paquete para el experimento de lanzar el dado.

#install.packages("dice")
library(dice)  
## Loading required package: gtools

Instalar libreria y llamar paquete para poner resultados en fracción.

#install.packages("MASS")
library(MASS) 

Ejemplos

¿Cuál es la probabilidad de obtener un 6 al lanzar un dado?

un_seis <- getEventProb(nrolls = 1, ndicePerRoll = 1, nsidesPerDie = 6, eventList = list(6))
un_seis  
## [1] 0.1666667
fractions(un_seis)  
## [1] 1/6

¿Cuál es la probabilidad de obtener un 5 al lanzar un dado?

un_cinco <- getEventProb(nrolls = 1, ndicePerRoll = 2, nsidesPerDie = 6, eventList = list(5))
un_cinco  
## [1] 0.1111111
fractions(un_cinco) 
## [1] 1/9

¿Cuál es la probabilidad de obtener un 5 en dos lanzammientos de dados consecutivos?

dos_cincos <- getEventProb(nrolls = 2, ndicePerRoll = 1, nsidesPerDie = 6, eventList = list(5,5))
dos_cincos 
## [1] 0.02777778
fractions(dos_cincos) 
## [1] 1/36

¿Que número es más probable alcanzar al lanzar dos dados?

sumar_dos <- getEventProb(nrolls = 1, ndicePerRoll = 2, nsidesPerDie = 6, eventList = list(2))
sumar_tres <- getEventProb(nrolls = 1, ndicePerRoll = 2, nsidesPerDie = 6, eventList = list(3))
sumar_cuatro <- getEventProb(nrolls = 1, ndicePerRoll = 2, nsidesPerDie = 6, eventList = list(4))
sumar_cinco <- getEventProb(nrolls = 1, ndicePerRoll = 2, nsidesPerDie = 6, eventList = list(5))
sumar_seis <- getEventProb(nrolls = 1, ndicePerRoll = 2, nsidesPerDie = 6, eventList = list(6))
sumar_siete <- getEventProb(nrolls = 1, ndicePerRoll = 2, nsidesPerDie = 6, eventList = list(7))
sumar_ocho <- getEventProb(nrolls = 1, ndicePerRoll = 2, nsidesPerDie = 6, eventList = list(8))
sumar_nueve <- getEventProb(nrolls = 1, ndicePerRoll = 2, nsidesPerDie = 6, eventList = list(9))
sumar_diez <- getEventProb(nrolls = 1, ndicePerRoll = 2, nsidesPerDie = 6, eventList = list(10))
sumar_once <- getEventProb(nrolls = 1, ndicePerRoll = 2, nsidesPerDie = 6, eventList = list(11))
sumar_doce <- getEventProb(nrolls = 1, ndicePerRoll = 2, nsidesPerDie = 6, eventList = list(12))
suma <- c(2,3,4,5,6,7,8,9,10,11,12) 
probabilidad <- c(sumar_dos, sumar_tres, sumar_cuatro, sumar_cinco, sumar_seis, sumar_siete, sumar_ocho, sumar_nueve, sumar_diez, sumar_once, sumar_doce)
tabla <- cbind(suma,probabilidad)  
barplot(probabilidad,names.arg = suma, main ="Probabilidad", xlab= "Suma de 2 Dados", col= "Tomato")

Experimento: Mano de Poker

Instalar libreria y llamar paquete para manejo de funciones y vectores

#install.packages("tidyverse")
library(purrr)  

Crear baraja inglesa T= 10

numero <- c(2,3,4,5,6,7,8,9,"T", "J", "Q", "K", "A")
numeros <- rep(numero,4) 
palo <- c("T", "C", "P", "D")
palos <- rep(palo,13)
baraja <- data.frame(numeros, palos)  

Crear el mazo de cartas

mazo <- apply(format(baraja), 1, paste,collapse="")
mazo
##    1    2    3    4    5    6    7    8    9   10   11   12   13   14   15   16 
## "2T" "3C" "4P" "5D" "6T" "7C" "8P" "9D" "TT" "JC" "QP" "KD" "AT" "2C" "3P" "4D" 
##   17   18   19   20   21   22   23   24   25   26   27   28   29   30   31   32 
## "5T" "6C" "7P" "8D" "9T" "TC" "JP" "QD" "KT" "AC" "2P" "3D" "4T" "5C" "6P" "7D" 
##   33   34   35   36   37   38   39   40   41   42   43   44   45   46   47   48 
## "8T" "9C" "TP" "JD" "QT" "KC" "AP" "2D" "3T" "4C" "5P" "6D" "7T" "8C" "9P" "TD" 
##   49   50   51   52 
## "JT" "QC" "KP" "AD"

Crear la mano de cartas

mano <- function(n) sample(mazo, n, rep=FALSE)
mi_mano <- mano(5)  
mi_mano 
##    1   36   38   12   35 
## "2T" "JD" "KC" "KD" "TP"

Conclusiones

A partir del este código observamos la utilidad que tiene R para diferentes aspectos de la estadística, como lo es la probabilidad. Inicialmente, observamos las definiciones de probabilidad y de experimenta, las cuales ayudarán próximamente en el código. Siguiendo con cortos ejemplos que ilustran para lo que sirve una probabilidad (como en Las Vegas) al momento de lanzar un dado; la gráfica expuesta demuestra la probabilidad de obtener cierta suma en dos dados.

Ahora, al ver un ejemplo logramos entender de mejor manera la utilidad de la probabilidad. En este caso, se toma en cuenta la baraja inglesa, en donde se quisiera ver cuál sería el resultado de una “mano” de poker. Para esto, se hizo el mazo de cartas, es decir, a cada valor se le asignó un signo para tener las 52 cartas totales.

Este tipo de herramientas y librerías podrán ser de gran utilidad para futuros análisis estadísticos.

LS0tCnRpdGxlOiA8c3BhbiBzdHlsZT0iY29sb3I6cmVkIj4qKlByb2JhYmlsaWRhZCoqCmF1dGhvcjogIkVsZW5hVmVsYV9BMDEyODM1MzUiCmRhdGU6ICIyMDIyLTEwLTA0IgpvdXRwdXQ6IAogIGh0bWxfZG9jdW1lbnQ6CiAgICB0b2M6IHRydWUKICAgIHRvY19mbG9hdDogdHJ1ZQogICAgdGhlbWU6IHVuaXRlZAogICAgaGlnaGxpZ2h0OiB0YW5nbwogICAgY29kZV9kb3dubG9hZDogdHJ1ZQotLS0KCjxpbWcgc3JjPSIvVXNlcnMvZWxlbmF2ZWxhL0Rlc2t0b3AvQ2FwdHVyYSBkZSBQYW50YWxsYSAyMDIyLTEwLTA0IGEgbGEocykgMTYuMDkuMjAucG5nIj4KCiMjIyAqKkRlZmluaWNpb25lcyoqCgoqKlByb2JhYmlsaWRhZDoqKiBFcyBlbCBlc3R1ZGlvIGRlIGF6YXIgeSBsYSBpbmNlcnRpZHVtYnJlIGVuIGN1YWxxdWllciBzaXR1YWNpb24gZW4gbGEgY3VhbCB2YXJpb3MgcG9zaWJsZXMgc3VjZXNvcyBwdWVkZW4gb2N1cnJpci4KKkVzIHVuIHZhbG9yIGVudHJlIDAgKGltcG9zaWJsZSkgeSAxIChzZWd1cm8pKgoqRWplbXBsbzoqIGxhIHByb2JhYmlsaWRhZCBkZSBxdWUgbGx1ZXZhIGhveSBlcyBkZSAwLjcwICg3MCUpLgoKKipFeHBlcmltZW50bzoqKiBDdWFscXVpZXIgYWNjacOzbiBjdXlvIHJlc3VsdGFkbyBlc3TDoSBzdWpldG8gYSBsYSBpbmNlcnRpZHVtYnJlLiAKKkVqZW1wbG86KiBMYW56YXIgdW5hIG1vbmVkYSBhbCBhaXJlLgoKSW5zdGFsYXIgbGlicmVyaWEgeSBsbGFtYXIgZWwgcGFxdWV0ZSBwYXJhIGVsIGV4cGVyaW1lbnRvIGRlIGxhbnphciBlbCBkYWRvLiAKYGBge3J9CiNpbnN0YWxsLnBhY2thZ2VzKCJkaWNlIikKbGlicmFyeShkaWNlKSAgCmBgYAoKSW5zdGFsYXIgbGlicmVyaWEgeSBsbGFtYXIgcGFxdWV0ZSBwYXJhIHBvbmVyIHJlc3VsdGFkb3MgZW4gZnJhY2Npw7NuLgpgYGB7cn0KI2luc3RhbGwucGFja2FnZXMoIk1BU1MiKQpsaWJyYXJ5KE1BU1MpIApgYGAKCiMjIyAqKkVqZW1wbG9zKioKCgrCv0N1w6FsIGVzIGxhIHByb2JhYmlsaWRhZCBkZSBvYnRlbmVyIHVuIDYgYWwgbGFuemFyIHVuIGRhZG8/IApgYGB7cn0KdW5fc2VpcyA8LSBnZXRFdmVudFByb2IobnJvbGxzID0gMSwgbmRpY2VQZXJSb2xsID0gMSwgbnNpZGVzUGVyRGllID0gNiwgZXZlbnRMaXN0ID0gbGlzdCg2KSkKdW5fc2VpcyAgCmZyYWN0aW9ucyh1bl9zZWlzKSAgCmBgYAoKCsK/Q3XDoWwgZXMgbGEgcHJvYmFiaWxpZGFkIGRlIG9idGVuZXIgdW4gNSBhbCBsYW56YXIgdW4gZGFkbz8gCmBgYHtyfQp1bl9jaW5jbyA8LSBnZXRFdmVudFByb2IobnJvbGxzID0gMSwgbmRpY2VQZXJSb2xsID0gMiwgbnNpZGVzUGVyRGllID0gNiwgZXZlbnRMaXN0ID0gbGlzdCg1KSkKdW5fY2luY28gIApmcmFjdGlvbnModW5fY2luY28pIApgYGAKCgrCv0N1w6FsIGVzIGxhIHByb2JhYmlsaWRhZCBkZSBvYnRlbmVyIHVuIDUgZW4gZG9zIGxhbnphbW1pZW50b3MgZGUgZGFkb3MgY29uc2VjdXRpdm9zPyAKYGBge3J9CmRvc19jaW5jb3MgPC0gZ2V0RXZlbnRQcm9iKG5yb2xscyA9IDIsIG5kaWNlUGVyUm9sbCA9IDEsIG5zaWRlc1BlckRpZSA9IDYsIGV2ZW50TGlzdCA9IGxpc3QoNSw1KSkKZG9zX2NpbmNvcyAKZnJhY3Rpb25zKGRvc19jaW5jb3MpIApgYGAKCsK/UXVlIG7Dum1lcm8gZXMgbcOhcyBwcm9iYWJsZSBhbGNhbnphciBhbCBsYW56YXIgZG9zIGRhZG9zPwpgYGB7cn0Kc3VtYXJfZG9zIDwtIGdldEV2ZW50UHJvYihucm9sbHMgPSAxLCBuZGljZVBlclJvbGwgPSAyLCBuc2lkZXNQZXJEaWUgPSA2LCBldmVudExpc3QgPSBsaXN0KDIpKQpzdW1hcl90cmVzIDwtIGdldEV2ZW50UHJvYihucm9sbHMgPSAxLCBuZGljZVBlclJvbGwgPSAyLCBuc2lkZXNQZXJEaWUgPSA2LCBldmVudExpc3QgPSBsaXN0KDMpKQpzdW1hcl9jdWF0cm8gPC0gZ2V0RXZlbnRQcm9iKG5yb2xscyA9IDEsIG5kaWNlUGVyUm9sbCA9IDIsIG5zaWRlc1BlckRpZSA9IDYsIGV2ZW50TGlzdCA9IGxpc3QoNCkpCnN1bWFyX2NpbmNvIDwtIGdldEV2ZW50UHJvYihucm9sbHMgPSAxLCBuZGljZVBlclJvbGwgPSAyLCBuc2lkZXNQZXJEaWUgPSA2LCBldmVudExpc3QgPSBsaXN0KDUpKQpzdW1hcl9zZWlzIDwtIGdldEV2ZW50UHJvYihucm9sbHMgPSAxLCBuZGljZVBlclJvbGwgPSAyLCBuc2lkZXNQZXJEaWUgPSA2LCBldmVudExpc3QgPSBsaXN0KDYpKQpzdW1hcl9zaWV0ZSA8LSBnZXRFdmVudFByb2IobnJvbGxzID0gMSwgbmRpY2VQZXJSb2xsID0gMiwgbnNpZGVzUGVyRGllID0gNiwgZXZlbnRMaXN0ID0gbGlzdCg3KSkKc3VtYXJfb2NobyA8LSBnZXRFdmVudFByb2IobnJvbGxzID0gMSwgbmRpY2VQZXJSb2xsID0gMiwgbnNpZGVzUGVyRGllID0gNiwgZXZlbnRMaXN0ID0gbGlzdCg4KSkKc3VtYXJfbnVldmUgPC0gZ2V0RXZlbnRQcm9iKG5yb2xscyA9IDEsIG5kaWNlUGVyUm9sbCA9IDIsIG5zaWRlc1BlckRpZSA9IDYsIGV2ZW50TGlzdCA9IGxpc3QoOSkpCnN1bWFyX2RpZXogPC0gZ2V0RXZlbnRQcm9iKG5yb2xscyA9IDEsIG5kaWNlUGVyUm9sbCA9IDIsIG5zaWRlc1BlckRpZSA9IDYsIGV2ZW50TGlzdCA9IGxpc3QoMTApKQpzdW1hcl9vbmNlIDwtIGdldEV2ZW50UHJvYihucm9sbHMgPSAxLCBuZGljZVBlclJvbGwgPSAyLCBuc2lkZXNQZXJEaWUgPSA2LCBldmVudExpc3QgPSBsaXN0KDExKSkKc3VtYXJfZG9jZSA8LSBnZXRFdmVudFByb2IobnJvbGxzID0gMSwgbmRpY2VQZXJSb2xsID0gMiwgbnNpZGVzUGVyRGllID0gNiwgZXZlbnRMaXN0ID0gbGlzdCgxMikpCnN1bWEgPC0gYygyLDMsNCw1LDYsNyw4LDksMTAsMTEsMTIpIApwcm9iYWJpbGlkYWQgPC0gYyhzdW1hcl9kb3MsIHN1bWFyX3RyZXMsIHN1bWFyX2N1YXRybywgc3VtYXJfY2luY28sIHN1bWFyX3NlaXMsIHN1bWFyX3NpZXRlLCBzdW1hcl9vY2hvLCBzdW1hcl9udWV2ZSwgc3VtYXJfZGlleiwgc3VtYXJfb25jZSwgc3VtYXJfZG9jZSkKdGFibGEgPC0gY2JpbmQoc3VtYSxwcm9iYWJpbGlkYWQpICAKYmFycGxvdChwcm9iYWJpbGlkYWQsbmFtZXMuYXJnID0gc3VtYSwgbWFpbiA9IlByb2JhYmlsaWRhZCIsIHhsYWI9ICJTdW1hIGRlIDIgRGFkb3MiLCBjb2w9ICJUb21hdG8iKQpgYGAKCiMjIyAqKkV4cGVyaW1lbnRvOiBNYW5vIGRlIFBva2VyKioKCkluc3RhbGFyIGxpYnJlcmlhIHkgbGxhbWFyIHBhcXVldGUgcGFyYSBtYW5lam8gZGUgZnVuY2lvbmVzIHkgdmVjdG9yZXMgCmBgYHtyfQojaW5zdGFsbC5wYWNrYWdlcygidGlkeXZlcnNlIikKbGlicmFyeShwdXJycikgIApgYGAKCgpDcmVhciBiYXJhamEgaW5nbGVzYSAKVD0gMTAKYGBge3J9Cm51bWVybyA8LSBjKDIsMyw0LDUsNiw3LDgsOSwiVCIsICJKIiwgIlEiLCAiSyIsICJBIikKbnVtZXJvcyA8LSByZXAobnVtZXJvLDQpIApwYWxvIDwtIGMoIlQiLCAiQyIsICJQIiwgIkQiKQpwYWxvcyA8LSByZXAocGFsbywxMykKYmFyYWphIDwtIGRhdGEuZnJhbWUobnVtZXJvcywgcGFsb3MpICAKYGBgCgpDcmVhciBlbCBtYXpvIGRlIGNhcnRhcyAKYGBge3J9Cm1hem8gPC0gYXBwbHkoZm9ybWF0KGJhcmFqYSksIDEsIHBhc3RlLGNvbGxhcHNlPSIiKQptYXpvCmBgYAoKICAgIApDcmVhciBsYSBtYW5vIGRlIGNhcnRhcyAKYGBge3J9Cm1hbm8gPC0gZnVuY3Rpb24obikgc2FtcGxlKG1hem8sIG4sIHJlcD1GQUxTRSkKbWlfbWFubyA8LSBtYW5vKDUpICAKbWlfbWFubyAKYGBgCgojIyMgKipDb25jbHVzaW9uZXMqKgpBIHBhcnRpciBkZWwgZXN0ZSBjw7NkaWdvIG9ic2VydmFtb3MgbGEgdXRpbGlkYWQgcXVlIHRpZW5lIFIgcGFyYSBkaWZlcmVudGVzIGFzcGVjdG9zIGRlIGxhIGVzdGFkw61zdGljYSwgY29tbyBsbyBlcyBsYSBwcm9iYWJpbGlkYWQuIEluaWNpYWxtZW50ZSwgb2JzZXJ2YW1vcyBsYXMgZGVmaW5pY2lvbmVzIGRlIHByb2JhYmlsaWRhZCB5IGRlIGV4cGVyaW1lbnRhLCBsYXMgY3VhbGVzIGF5dWRhcsOhbiBwcsOzeGltYW1lbnRlIGVuIGVsIGPDs2RpZ28uIFNpZ3VpZW5kbyBjb24gY29ydG9zIGVqZW1wbG9zIHF1ZSBpbHVzdHJhbiBwYXJhIGxvIHF1ZSBzaXJ2ZSB1bmEgcHJvYmFiaWxpZGFkIChjb21vIGVuIExhcyBWZWdhcykgYWwgbW9tZW50byBkZSBsYW56YXIgdW4gZGFkbzsgbGEgZ3LDoWZpY2EgZXhwdWVzdGEgZGVtdWVzdHJhIGxhIHByb2JhYmlsaWRhZCBkZSBvYnRlbmVyIGNpZXJ0YSBzdW1hIGVuIGRvcyBkYWRvcy4gIAoKQWhvcmEsIGFsIHZlciB1biBlamVtcGxvIGxvZ3JhbW9zIGVudGVuZGVyIGRlIG1lam9yIG1hbmVyYSBsYSB1dGlsaWRhZCBkZSBsYSBwcm9iYWJpbGlkYWQuIEVuICBlc3RlIGNhc28sIHNlIHRvbWEgZW4gY3VlbnRhIGxhIGJhcmFqYSBpbmdsZXNhLCBlbiBkb25kZSBzZSBxdWlzaWVyYSB2ZXIgY3XDoWwgc2Vyw61hIGVsIHJlc3VsdGFkbyBkZSB1bmEgIm1hbm8iIGRlIHBva2VyLiBQYXJhIGVzdG8sIHNlIGhpem8gZWwgbWF6byBkZSBjYXJ0YXMsIGVzIGRlY2lyLCBhIGNhZGEgdmFsb3Igc2UgbGUgYXNpZ27DsyB1biBzaWdubyBwYXJhIHRlbmVyIGxhcyA1MiBjYXJ0YXMgdG90YWxlcy4gIAoKRXN0ZSB0aXBvIGRlIGhlcnJhbWllbnRhcyB5IGxpYnJlcsOtYXMgcG9kcsOhbiBzZXIgZGUgZ3JhbiB1dGlsaWRhZCBwYXJhIGZ1dHVyb3MgYW7DoWxpc2lzIGVzdGFkw61zdGljb3MuIAogIAoKCgoKCgoKCgoKCgoKCgo=