-
Notifications
You must be signed in to change notification settings - Fork 0
/
ex5.2c.r
43 lines (31 loc) · 844 Bytes
/
ex5.2c.r
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
#Remove todas as variáveis do workspace
rm(list=ls(all=TRUE))
#Carrega a biblioteca que permite "exponenciar" uma matriz quadrada
library(expm)
#Carrega as funções úteis
source("f.uteis.r")
#Exercício Propriamente Dito
exercicio<-c("Exercício 5.2c")
#Determinando a matriz de transição
P<-matrix(0,4,4)
P[1,]<-c(0.90, 0.05, 0.03, 0.02)
P[2,]<-c(0.00, 0.85, 0.09, 0.06)
P[3,]<-c(0.00, 0.00, 0.90, 0.10)
P[4,]<-c(1.00, 0.00, 0.00, 0.00)
#Obtendo o resultado por simulação
nsim<-10000
soma<-0
n<-3
for(isim in 1:nsim){
X<-f.gerasequencia(P,n,1)
if(sum(X==4)>0){
soma<-soma+1
}
}
prob.simulacao<-soma/nsim
#Obtendo o resultado teórico
P[4,]<-c(0.00, 0.00, 0.00, 1.00)
PN<-P%^%n
prob.teorica<-PN[1,4]
#Imprimindo os resultados na tela
f.resultado(exercicio,prob.teorica,prob.simulacao,nsim)