1 files changed, 25 insertions, 0 deletions
diff --git a/291/CH7/EX7.4b/eg7_4b.sce b/291/CH7/EX7.4b/eg7_4b.sce
new file mode 100755
index 000000000..f6eaba5fc
--- /dev/null
+++ b/291/CH7/EX7.4b/eg7_4b.sce
@@ -0,0 +1,25 @@
+tech1 = [140 136 138 150 152 144 132 142 150 154 136 142];
+tech2 = [144 132 136 140 128 150 130 134 130 146 128 131 137 135];
+num1= 12;
+num2= 14;
+mean1= mean(tech1);
+mean2= mean(tech2);
+//disp(mean1)
+//disp(Sp)
+alpha = 0.9;
+S1 = variance(tech1) //*num1/(num1-1);
+S2 = variance(tech2) // *num2/(num2-1);
+Sp = (((num1-1)*S1) + ((num2-1)*S2))/(num1+ num2 -2);
+Sp= sqrt(Sp);
+num= (1/num1)+(1/num2);
+betaa = (1-alpha)/2;
+tval = cdft("T", num1+num2-2, 1-betaa, betaa);
+upperlim = mean1-mean2 + (tval*Sp)*sqrt(num);
+lowerlim = mean1-mean2 - (tval*Sp)*sqrt(num);
+disp(upperlim, "to ",lowerlim,"The 90% confidence interval is "   )
+alpha = 0.95
+betaaa = 1-alpha;
+tval = cdft("T", num1+num2-2, 1-betaaa, betaaa);
+lowerlim = mean1-mean2 - (tval*Sp)*sqrt(num);
+disp("the upper confidence interval is")
+disp(" to infinity", lowerlim)
+\ No newline at end of file