Probabilidad

Es el estudio de azar o incertidumbre en cualquier situación 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 modeda al aire.

Experimento: lanzar un dado.

Librerias

library(dice)
library(MASS)
library(purrr)

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

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

¿Cuál es la probabilidad de sumar 5 al lanzar dos dados?

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 cinco en dos lanzamientos 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

¿Qué número es más probable de 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 = "Probabildiad", xlab = "Suma de 2 dados", col = "Tomato")

Experimento Mano de Poker

Crear baraja inglesa

palo<-c("Picas", "Corazones", "Diamantes", "Treboles")
palos<-rep(palo,13)
numero<-c( "A",2, 3, 4, 5, 6, 7, 8, 9, "D", "J", "Q", "K")
numeros<-rep(numero,4)
baraja<-cbind(palos,numeros)
baraja
##       palos       numeros
##  [1,] "Picas"     "A"    
##  [2,] "Corazones" "2"    
##  [3,] "Diamantes" "3"    
##  [4,] "Treboles"  "4"    
##  [5,] "Picas"     "5"    
##  [6,] "Corazones" "6"    
##  [7,] "Diamantes" "7"    
##  [8,] "Treboles"  "8"    
##  [9,] "Picas"     "9"    
## [10,] "Corazones" "D"    
## [11,] "Diamantes" "J"    
## [12,] "Treboles"  "Q"    
## [13,] "Picas"     "K"    
## [14,] "Corazones" "A"    
## [15,] "Diamantes" "2"    
## [16,] "Treboles"  "3"    
## [17,] "Picas"     "4"    
## [18,] "Corazones" "5"    
## [19,] "Diamantes" "6"    
## [20,] "Treboles"  "7"    
## [21,] "Picas"     "8"    
## [22,] "Corazones" "9"    
## [23,] "Diamantes" "D"    
## [24,] "Treboles"  "J"    
## [25,] "Picas"     "Q"    
## [26,] "Corazones" "K"    
## [27,] "Diamantes" "A"    
## [28,] "Treboles"  "2"    
## [29,] "Picas"     "3"    
## [30,] "Corazones" "4"    
## [31,] "Diamantes" "5"    
## [32,] "Treboles"  "6"    
## [33,] "Picas"     "7"    
## [34,] "Corazones" "8"    
## [35,] "Diamantes" "9"    
## [36,] "Treboles"  "D"    
## [37,] "Picas"     "J"    
## [38,] "Corazones" "Q"    
## [39,] "Diamantes" "K"    
## [40,] "Treboles"  "A"    
## [41,] "Picas"     "2"    
## [42,] "Corazones" "3"    
## [43,] "Diamantes" "4"    
## [44,] "Treboles"  "5"    
## [45,] "Picas"     "6"    
## [46,] "Corazones" "7"    
## [47,] "Diamantes" "8"    
## [48,] "Treboles"  "9"    
## [49,] "Picas"     "D"    
## [50,] "Corazones" "J"    
## [51,] "Diamantes" "Q"    
## [52,] "Treboles"  "K"

Crear el mazo de barajas

mazo<-apply(format(baraja),1,paste,collapse="")
mazo
##  [1] "Picas    A        " "Corazones2        " "Diamantes3        "
##  [4] "Treboles 4        " "Picas    5        " "Corazones6        "
##  [7] "Diamantes7        " "Treboles 8        " "Picas    9        "
## [10] "CorazonesD        " "DiamantesJ        " "Treboles Q        "
## [13] "Picas    K        " "CorazonesA        " "Diamantes2        "
## [16] "Treboles 3        " "Picas    4        " "Corazones5        "
## [19] "Diamantes6        " "Treboles 7        " "Picas    8        "
## [22] "Corazones9        " "DiamantesD        " "Treboles J        "
## [25] "Picas    Q        " "CorazonesK        " "DiamantesA        "
## [28] "Treboles 2        " "Picas    3        " "Corazones4        "
## [31] "Diamantes5        " "Treboles 6        " "Picas    7        "
## [34] "Corazones8        " "Diamantes9        " "Treboles D        "
## [37] "Picas    J        " "CorazonesQ        " "DiamantesK        "
## [40] "Treboles A        " "Picas    2        " "Corazones3        "
## [43] "Diamantes4        " "Treboles 5        " "Picas    6        "
## [46] "Corazones7        " "Diamantes8        " "Treboles 9        "
## [49] "Picas    D        " "CorazonesJ        " "DiamantesQ        "
## [52] "Treboles K        "

Crear mano de cartas

mano<-function(n) sample(mazo,n,rep=FALSE)
mi_mano<-mano(5)
mi_mano
## [1] "Corazones9        " "DiamantesQ        " "Diamantes6        "
## [4] "Treboles 8        " "CorazonesD        "
LS0tDQp0aXRsZTogIlByb2JhYmlsaWRhZCINCmF1dGhvcjogIkphdmllciBBeWFsYSINCmRhdGU6ICcyMDIyLTA5LTA1Jw0Kb3V0cHV0Og0KICBodG1sX2RvY3VtZW50Og0KICAgIHRvYzogdHJ1ZQ0KICAgIHRvY19mbG9hdDogdHJ1ZQ0KICAgIGNvZGVfZG93bmxvYWQ6IHRydWUNCi0tLQ0KDQojIyMgUHJvYmFiaWxpZGFkICAgDQpFcyBlbCBlc3R1ZGlvIGRlIGF6YXIgbyAgaW5jZXJ0aWR1bWJyZSBlbiBjdWFscXVpZXIgc2l0dWFjacOzbiBlbiBsYSBjdWFsIHZhcmlvcyBwb3NpYmxlcyBzdWNlc29zIHB1ZWRlbiBvY3VycmlyLiAgDQoNCkVzIHVuIHZhbG9yIGVudHJlIDAoaW1wb3NpYmxlKSB5IDEoc2VndXJvKS4gIA0KRWplbXBsbzogbGEgcHJvYmFiaWxpZGFkIGRlIHF1ZSBsbHVldmEgaG95IGVzIGRlIDAuNzAgKDcwJSkuICANCg0KRXhwZXJpbWVudG86IEN1YWxxdWllciBhY2Npw7NuIGN1eW8gcmVzdWx0YWRvIGVzdMOhIHN1amV0byBhIGxhIGluY2VydGlkdW1icmUuICANCkVqZW1wbG86IGxhbnphciB1bmEgbW9kZWRhIGFsIGFpcmUuICANCg0KIyMjIEV4cGVyaW1lbnRvOiBsYW56YXIgdW4gZGFkby4gIA0KDQojIyMjIExpYnJlcmlhcyAgDQpgYGB7ciBtZXNzYWdlPUZBTFNFLCB3YXJuaW5nPUZBTFNFfQ0KbGlicmFyeShkaWNlKQ0KbGlicmFyeShNQVNTKQ0KbGlicmFyeShwdXJycikNCmBgYA0KDQoNCiMjIyMgwr9DdcOhbCBlcyBsYSBwcm9iYWJpbGlkYWQgZGUgb2J0ZW5lciB1biA2IGFsIGxhbnphciB1biBkYWRvPyAgDQpgYGB7cn0NCnVuX3NlaXM8LWdldEV2ZW50UHJvYihucm9sbHMgPSAxLCBuZGljZVBlclJvbGwgPSAxLCBuc2lkZXNQZXJEaWUgPSA2LCBldmVudExpc3QgPSBsaXN0KDE6MikpDQp1bl9zZWlzDQpmcmFjdGlvbnModW5fc2VpcykNCmBgYA0KDQoNCiMjIyMgwr9DdcOhbCBlcyBsYSBwcm9iYWJpbGlkYWQgZGUgc3VtYXIgNSBhbCBsYW56YXIgZG9zIGRhZG9zPyAgDQpgYGB7cn0NCnVuX2NpbmNvPC1nZXRFdmVudFByb2IobnJvbGxzID0gMSwgbmRpY2VQZXJSb2xsID0gMiwgbnNpZGVzUGVyRGllID0gNiwgZXZlbnRMaXN0ID0gbGlzdCg1KSkNCnVuX2NpbmNvDQpmcmFjdGlvbnModW5fY2luY28pDQpgYGANCg0KDQojIyMjIMK/Q3XDoWwgZXMgbGEgcHJvYmFiaWxpZGFkIGRlIG9idGVuZXIgY2luY28gZW4gZG9zIGxhbnphbWllbnRvcyBkZSBkYWRvcyBjb25zZWN1dGl2b3M/ICANCmBgYHtyfQ0KZG9zX2NpbmNvczwtZ2V0RXZlbnRQcm9iKG5yb2xscyA9IDIsIG5kaWNlUGVyUm9sbCA9IDEsIG5zaWRlc1BlckRpZSA9IDYsIGV2ZW50TGlzdCA9IGxpc3QoNSw1KSkNCmRvc19jaW5jb3MNCmZyYWN0aW9ucyhkb3NfY2luY29zKQ0KYGBgDQoNCg0KIyMjIyDCv1F1w6kgbsO6bWVybyBlcyBtw6FzIHByb2JhYmxlIGRlIGFsY2FuemFyIGFsIGxhbnphciBkb3MgZGFkb3M/ICANCmBgYHtyfQ0Kc3VtYXJfZG9zPC1nZXRFdmVudFByb2IobnJvbGxzID0gMSwgbmRpY2VQZXJSb2xsID0gMiwgbnNpZGVzUGVyRGllID0gNiwgZXZlbnRMaXN0ID0gbGlzdCgyKSkNCnN1bWFyX3RyZXM8LWdldEV2ZW50UHJvYihucm9sbHMgPSAxLCBuZGljZVBlclJvbGwgPSAyLCBuc2lkZXNQZXJEaWUgPSA2LCBldmVudExpc3QgPSBsaXN0KDMpKQ0Kc3VtYXJfY3VhdHJvPC1nZXRFdmVudFByb2IobnJvbGxzID0gMSwgbmRpY2VQZXJSb2xsID0gMiwgbnNpZGVzUGVyRGllID0gNiwgZXZlbnRMaXN0ID0gbGlzdCg0KSkNCnN1bWFyX2NpbmNvPC1nZXRFdmVudFByb2IobnJvbGxzID0gMSwgbmRpY2VQZXJSb2xsID0gMiwgbnNpZGVzUGVyRGllID0gNiwgZXZlbnRMaXN0ID0gbGlzdCg1KSkNCnN1bWFyX3NlaXM8LWdldEV2ZW50UHJvYihucm9sbHMgPSAxLCBuZGljZVBlclJvbGwgPSAyLCBuc2lkZXNQZXJEaWUgPSA2LCBldmVudExpc3QgPSBsaXN0KDYpKQ0Kc3VtYXJfc2lldGU8LWdldEV2ZW50UHJvYihucm9sbHMgPSAxLCBuZGljZVBlclJvbGwgPSAyLCBuc2lkZXNQZXJEaWUgPSA2LCBldmVudExpc3QgPSBsaXN0KDcpKQ0Kc3VtYXJfb2NobzwtZ2V0RXZlbnRQcm9iKG5yb2xscyA9IDEsIG5kaWNlUGVyUm9sbCA9IDIsIG5zaWRlc1BlckRpZSA9IDYsIGV2ZW50TGlzdCA9IGxpc3QoOCkpDQpzdW1hcl9udWV2ZTwtZ2V0RXZlbnRQcm9iKG5yb2xscyA9IDEsIG5kaWNlUGVyUm9sbCA9IDIsIG5zaWRlc1BlckRpZSA9IDYsIGV2ZW50TGlzdCA9IGxpc3QoOSkpDQpzdW1hcl9kaWV6PC1nZXRFdmVudFByb2IobnJvbGxzID0gMSwgbmRpY2VQZXJSb2xsID0gMiwgbnNpZGVzUGVyRGllID0gNiwgZXZlbnRMaXN0ID0gbGlzdCgxMCkpDQpzdW1hcl9vbmNlPC1nZXRFdmVudFByb2IobnJvbGxzID0gMSwgbmRpY2VQZXJSb2xsID0gMiwgbnNpZGVzUGVyRGllID0gNiwgZXZlbnRMaXN0ID0gbGlzdCgxMSkpDQpzdW1hcl9kb2NlPC1nZXRFdmVudFByb2IobnJvbGxzID0gMSwgbmRpY2VQZXJSb2xsID0gMiwgbnNpZGVzUGVyRGllID0gNiwgZXZlbnRMaXN0ID0gbGlzdCgxMikpDQoNCnN1bWE8LWMoMiwzLDQsNSw2LDcsOCw5LDEwLDExLDEyKQ0KcHJvYmFiaWxpZGFkPC1jKHN1bWFyX2RvcyxzdW1hcl90cmVzLHN1bWFyX2N1YXRybyxzdW1hcl9jaW5jbyxzdW1hcl9zZWlzLHN1bWFyX3NpZXRlLHN1bWFyX29jaG8sc3VtYXJfbnVldmUsc3VtYXJfZGlleixzdW1hcl9vbmNlLHN1bWFyX2RvY2UpDQoNCnRhYmxhPC1jYmluZChzdW1hLHByb2JhYmlsaWRhZCkNCg0KYmFycGxvdChwcm9iYWJpbGlkYWQsIG5hbWVzLmFyZyA9IHN1bWEsIG1haW4gPSAiUHJvYmFiaWxkaWFkIiwgeGxhYiA9ICJTdW1hIGRlIDIgZGFkb3MiLCBjb2wgPSAiVG9tYXRvIikNCmBgYA0KDQoNCiMjIyBFeHBlcmltZW50byBNYW5vIGRlIFBva2VyICANCg0KIyMjIyBDcmVhciBiYXJhamEgaW5nbGVzYSAgDQpgYGB7cn0NCnBhbG88LWMoIlBpY2FzIiwgIkNvcmF6b25lcyIsICJEaWFtYW50ZXMiLCAiVHJlYm9sZXMiKQ0KcGFsb3M8LXJlcChwYWxvLDEzKQ0KbnVtZXJvPC1jKCAiQSIsMiwgMywgNCwgNSwgNiwgNywgOCwgOSwgIkQiLCAiSiIsICJRIiwgIksiKQ0KbnVtZXJvczwtcmVwKG51bWVybyw0KQ0KYmFyYWphPC1jYmluZChwYWxvcyxudW1lcm9zKQ0KYmFyYWphDQpgYGANCg0KDQojIyMjIENyZWFyIGVsIG1hem8gZGUgYmFyYWphcyAgDQpgYGB7cn0NCm1hem88LWFwcGx5KGZvcm1hdChiYXJhamEpLDEscGFzdGUsY29sbGFwc2U9IiIpDQptYXpvDQoNCmBgYA0KDQoNCiMjIyMgQ3JlYXIgbWFubyBkZSBjYXJ0YXMgIA0KYGBge3J9DQptYW5vPC1mdW5jdGlvbihuKSBzYW1wbGUobWF6byxuLHJlcD1GQUxTRSkNCm1pX21hbm88LW1hbm8oNSkNCm1pX21hbm8NCmBgYA0KDQo=