P2

Sea la cadena de Markov cuya matriz de transición en un paso viene dada por

La matriz de transición de compone de 5 estados de los diferentes tipos de productos:congelados, conservas, curados, harina de pescado y aceite respectivamente\(S=\{1,2,3,4,5\}\), los cuales se pueden modelar en R, gracias a la librería markovchain.

library(markovchain)

## Package:  markovchain
## Version:  0.8.5-2
## Date:     2020-09-07
## BugReport: https://github.com/spedygiorgio/markovchain/issues

library(igraph)

## 
## Attaching package: 'igraph'

## The following objects are masked from 'package:stats':
## 
##     decompose, spectrum

## The following object is masked from 'package:base':
## 
##     union

# Cargamos la matriz de transición en un paso una fila (row) a la vez
P1 = matrix (c(0.5,0.5,0,0,0,
               0.5,0.25,0.25,0,0,
               0,0,0.375,0.625,0,
               0,0,0,0.1429,0.8571,
               0,0,0.3333,0.3333,0.3334), nrow=5, byrow = TRUE)
# Creamos un objeto tipo markovchain llamado cm1 
# Describimos su conjunto de estados con el parámetro states
# Le damos el nombre: Cadena de Markov No. 1 con el parámetro name
mc1 = new("markovchain", transitionMatrix = P1, states = c("1","2","3","4","5"), name = "Cadena de Markov No. 1")
# Generamos el diagrama de transición en un paso para una mejor visualización de la cadena mc1
plot(mc1, layout=layout_with_gem, vertex.color = "white", vertex.size = 40, edge.color = "blue", vertex.frame.color = "red", vertex.label.cex = 2, edge.arrow.size = 0.5, edge.curved = 0, edge.label.cex = 1.5, edge.label.x = NA, edge.label.y = NA, edge.label.color = "black")

La probabilidad incondicional de que la cadena en 3 órdenes ordena 4(harina de pescado)

#Distribución inicial
alpha0<-c(0.13,0.2,0.12,0.37,0.18)
alpha3<-alpha0*mc1^3
alpha3

##             1         2         3         4         5
## [1,] 0.125625 0.0978125 0.2031042 0.2640768 0.3093815

De acá a 3 órdenes qué la fracción de los clientes estará ordenando harina de pescado(4) es 0.2640768