{

"nbformat": 4,
"nbformat_minor": 0,
"metadata": {
"colab": {
"name": "Untitled4.ipynb",
"provenance": [],
"collapsed_sections": []
},
"kernelspec": {
"name": "python3",
"display_name": "Python 3"
},
"language_info": {
"name": "python"
}
},
"cells": [
{
"cell_type": "code",
"execution_count": 64,
"metadata": {
"id": "70zCZb4VuYfk"
},
"outputs": [],
"source": [
"import scipy.stats as ss\n",
"import math\n",
"import numpy as np\n",
"\n"
]
},
{
"cell_type": "code",
"source": [
"#binomial\n",
"n=4\n",
"p=0.8078\n",
"mu=ss.binom.mean(n,p) #media o esperanza, numero esperado\n",
"sigma = ss.binom.std(n,p) # desviación tipica/desviacion estandar\n",
"var=ss.binom.var(n,p)\n",
"\n",
"#probabilidad en un punto\n",
"x=1\n",
"Pmf=ss.binom.pmf(x,n,p)\n",
"\n",
"#probabilidad acumulada a un punto P(X<=x)\n",
"Cdf = ss.binom.cdf(x,n,p)\n",
"\n",
"#probabilidad mayor a un punto P(X>x)\n",
"\n",
"Sf=ss.binom.sf(x,n,p)\n",
"\n",
"#Obtener x, para pertenecer al % inferior, ddado este %, el resultado
entrega X<x (debe ser menor al valor entregado para pertenecer)\n",
"\n",
"prob=0.117\n",
"Ppf=ss.binom(n,p).ppf(prob)\n",
"\n",
 "print(Ppf)\n",
"print(Pmf)\n",
"print(round(Cdf,4))\n",
"Cdf"
],
"metadata": {
"colab": {
"base_uri": "https://localhost:8080/"
},
"id": "c8X0BFISu2Vs",
"outputId": "17d2a863-be19-450a-ddd9-11ff15245989"
},
"execution_count": 78,
"outputs": [
{
"output_type": "stream",
"name": "stdout",
"text": [
"2.0\n",
"0.022941615152377616\n",
"0.0243\n"
]
},
{
"output_type": "execute_result",
"data": {
"text/plain": [
"0.024306240812283213"
]
},
"metadata": {},
"execution_count": 78
}
]
},
{
"cell_type": "code",
"source": [
"#POISSON\n",
"#Lambda = media = mu = esperanza\n",
"lamb=2\n",
"\n",
"P_cdf= round(ss.poisson.cdf(x,lamb),4)\n",
"\n",
"P_sf= round(ss.poisson.sf(1,lamb),4)\n",
"\n",
"P_pmf= ss.poisson.pmf(2,lamb)\n",
"\n",
"P_prob=0.117\n",
"Ppf=ss.poisson(lamb).ppf(P_prob)\n",
"\n",
"Ppf\n",
"ss.poisson.sf(1,lamb)\n",
"ss.poisson.pmf(2,lamb)\n",
"print(ss.poisson.cdf(2,lamb)-ss.poisson.cdf(1,lamb))\n",
"print(round(ss.poisson.pmf(2,lamb)/ss.poisson.sf(1,lamb),4))"
],
"metadata": {
"colab": {
 "base_uri": "https://localhost:8080/"
},
"id": "QWYuFuNA54ZR",
"outputId": "3d208784-38ee-43d2-e03a-ee8fbbe2eab5"
},
"execution_count": 74,
"outputs": [
{
"output_type": "stream",
"name": "stdout",
"text": [
"0.27067056647322546\n",
"0.4557\n"
]
}
]
},
{
"cell_type": "code",
"source": [
"#Normal\n",
"mu = 176\n",
"sigma = 4.5957\n",
"\n",
"y=180\n",
"\n",
"N_pdf = ss.norm(mu,sigma).pdf(y)\n",
"N_cdf = ss.norm(mu,sigma).cdf(y)\n",
"N_sf = ss.norm(mu,sigma).sf(y)\n",
"N_ppf=ss.norm(mu,sigma).ppf(0.117)\n",
"\n",
"print(\"Mayor a 36 =\", round(ss.norm(mu,sigma).sf(36),4))\n",
"print(\"Entre 21 y 27\", round(ss.norm(mu,sigma).cdf(27) -
ss.norm(mu,sigma).cdf(21),4))\n",
"print(\"Que tan alta debe ser la remuneración por hora para que el pago
este dentro del 10.2% superior\", ss.norm(mu,sigma).ppf(1-0.102))\n",
"\n",
"\n",
"1- ss.norm(mu,sigma).cdf(37.83499801405122)\n",
"\n",
"30*(ss.norm(mu,sigma).cdf(27) - ss.norm(mu,sigma).cdf(21))\n",
"\n",
"N_sf\n",
"\n"
],
"metadata": {
"colab": {
"base_uri": "https://localhost:8080/"
},
"id": "ReM1c3R5EjBv",
"outputId": "8f06530e-4184-4822-e81f-c7f6e5c01507"
},
"execution_count": 77,
"outputs": [
{
"output_type": "stream",
"name": "stdout",
"text": [
"Mayor a 36 = 1.0\n",
 "Entre 21 y 27 0.0\n",
"Que tan alta debe ser la remuneración por hora para que el pago este
dentro del 10.2% superior 181.83763104123219\n"
]
},
{
"output_type": "execute_result",
"data": {
"text/plain": [
"0.19204670532263335"
]
},
"metadata": {},
"execution_count": 77
}
]
},
{
"cell_type": "code",
"source": [
""
],
"metadata": {
"id": "xg2fzhCVEjJM"
},
"execution_count": null,
"outputs": []
}
]
}