blob: 0a325f5403e9661f690243320eab1ee778cce1fe (
plain)
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
|
//Caption:One-tailed Tests Concerning Difference between Two Proportions
//Test 1: Ho: p1-p2<=0 and H1: p1-p2>0
//Example8.26
//Page277
clc;
n1= 120; //sample size of the vendor-1
n2 = 100; //sample size of the vendor-2
N1 = 15;//number of defective pieces in the sample-1
N2 = 11;//number of defective pieces in the sample-2
p1 = N1/n1; // proportion w.r.to vendor-1
p2 = N2/n2; // proportion w.r.to vendor-2
Stdp1_p2 = sqrt(((p1*(1-p1))/n1)+((p2*(1-p2))/n2))
alpha = 0.05;
Test = input('Enter the type of test=')
Zp1_p2 = ((p1-p2)-(0))/Stdp1_p2; //calculated normal statistic for (p1-p2)
disp(Zp1_p2,'The calculated normal statistic for (p1-p2)=')
z_Stand = standard_normal_zstat(alpha);//standard normal statistic for (p1-p2)
disp(z_Stand,'The Standard normal statistic for (p1-p2)=')
if (Test==1) then
if(Zp1_p2 < z_Stand) then
disp('It falls in the Acceptance Region')
disp('Then Null Hypothesis Ho should be Accepted')
elseif(Zp1_p2 > z_Stand) then
disp('It falls in the Rejection Region')
disp('Then Null Hypothesis Ho should be Rejected')
end
elseif (Test==2) then
if(Zp1_p2 >-z_Stand)
disp('It falls in the Acceptance Region')
disp('Then Null Hypothesis Ho should be Accepted')
elseif(Zp1_p2<-z_Stand)
disp('It falls in the Rejection Region')
disp('Then Null Hypothesis Ho should be Rejected')
end
end
//Result
//Enter the type of test=1
//
// The calculated normal statistic for (p1-p2)=
//
// 0.3449910
//
// The Standard normal statistic for (p1-p2)=
//
// 1.64
//
// It falls in the Acceptance Region
//
// Then Null Hypothesis Ho should be Accepted
//
//-->p1
// p1 =
//
// 0.125
//
//-->p2
// p2 =
//
// 0.11
//
//-->Stdp1_p2
// Stdp1_p2 =
//
// 0.0434794
|
