-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathexample-03-a.m
70 lines (49 loc) · 1.89 KB
/
example-03-a.m
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
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
%{
University: Universidad de Valladolid
Degree: Grado en Estadística
Subject: Muestreo Estadístico 1
Year: 2017/18
Teacher: Jesús Alberto Tapia García
Author: Sergio García Prado (garciparedes.me)
Name: simple-random-sampling/example-03-a.m
%}
clear;
addpath(genpath('./../utils/'));
% Tercer ejemplo de teoria
% Simulación para verificar la observación de m.a.s. donde se determina el
% tamaño de la población a partir del cual hay que encuestar a los mismos si
% mantenemos error de estimacion confianza y variabilidad
% Registros contables de un hospital
% Auditoria con el fin de estimar la cantidad total que falta por cobrar
% Tamaño poblacion registros por auditorar
N = 28000;
%Este tamaño lo iremos variando para comprobar la observación
% Simulamos datos poblacionales, cantidad falta por cobrar en la poblacion
% de N registros
Y = normrnd(1500, 100, N, 1);
% Ejercicio 1:
% Determinar el tamaño de muestra para estimar mu cometiendo un error de
% estimación B de 3 euros con una confianza del 95%
B = 3
alfa = 0.05;
z_value = norminv(1 - alfa / 2, 0,1);
n0 = z_value ^ 2 * var(Y) / B ^ 2;
n1 = ceil(n0 / ((n0 / N) +1))
% Ejercicio 2:
% Calcular el valor donde se alcanza la asintota segun la observación teórica N2
N2 = ceil(n0 * (n0 - 1))
Y2 = normrnd(3500, 100, N2, 1);
% Ejercicio 3:
% Determinar el tamaño de muestra para estimar mu en la poblacion e tamaño N2
% cometiendo el mismo error de estimacion B y confianza del 95%
n2 = calmumnmas(N2, B , 0.975, var(Y2))
% Ejercicio 4:
% Tomar poblaciones de tamaño superior a N2 y comprobar que el tamaño de muestra
% se mantiene
N3 = 2 * N2;
Y3 = normrnd(1500, 100, N3, 1);
n3 = calmumnmas(N3, B , 0.975, var(Y3))
% Ejercicio 5:
% Determinar el tamaño de muestra para estimar mu cometiendo un error de
% estimación B de 5000 euros con una confianza del 95%
n4 = calmumnmas(N3, 5000 , 0.975, var(Y3))