1reset data
2reset plot control
3reset i/o
4dimension 500 30
5label size 3
6read mpc262.dat run wafer probe mo day op hum y sr
7let df = sr - sr + 5.
8let m3 = y - y
9let ybar = average y subset run 1 subset wafer 138
10let m3 = m3 + ybar subset run 1 subset wafer 138
11let ybar = average y subset run 1 subset wafer 139
12let m3 = m3 + ybar subset run 1 subset wafer 139
13let ybar = average y subset run 1 subset wafer 140
14let m3 = m3 + ybar subset run 1 subset wafer 140
15let ybar = average y subset run 1 subset wafer 141
16let m3 = m3 + ybar subset run 1 subset wafer 141
17let ybar = average y subset run 1 subset wafer 142
18let m3 = m3 + ybar subset run 1 subset wafer 142
19let ybar = average y subset run 2 subset wafer 138
20let m3 = m3 + ybar subset run 2 subset wafer 138
21let ybar = average y subset run 2 subset wafer 139
22let m3 = m3 + ybar subset run 2 subset wafer 139
23let ybar = average y subset run 2 subset wafer 140
24let m3 = m3 + ybar subset run 2 subset wafer 140
25let ybar = average y subset run 2 subset wafer 141
26let m3 = m3 + ybar subset run 2 subset wafer 141
27let ybar = average y subset run 2 subset wafer 142
28let m3 = m3 + ybar subset run 2 subset wafer 142
29let d = y - m3
30let bias1 = average d subset probe 2362 subset run 1
31let bias2 = average d subset probe 2362 subset run 2
32let b11=average d subset probe 2362 subset run 1 subset wafer 138
33let b21=average d subset probe 2362 subset run 1 subset wafer 139
34let b31=average d subset probe 2362 subset run 1 subset wafer 140
35let b41=average d subset probe 2362 subset run 1 subset wafer 141
36let b51=average d subset probe 2362 subset run 1 subset wafer 142
37let b12=average d subset probe 2362 subset run 2 subset wafer 138
38let b22=average d subset probe 2362 subset run 2 subset wafer 139
39let b32=average d subset probe 2362 subset run 2 subset wafer 140
40let b42=average d subset probe 2362 subset run 2 subset wafer 141
41let b52=average d subset probe 2362 subset run 2 subset wafer 142
42let z = 0 for i = 1 1 10
43let z = z +b11 for i = 1 1 1
44let z = z +b21 for i = 2 1 2
45let z = z +b31 for i = 3 1 3
46let z = z +b41 for i = 4 1 4
47let z = z +b51 for i = 5 1 5
48let z = z +b12 for i = 6 1 6
49let z = z +b22 for i = 7 1 7
50let z = z +b32 for i = 8 1 8
51let z = z +b42 for i = 9 1 9
52let z = z +b52 for i = 10 1 10
53let sd= standard deviation z
54let sdcorr=sd/(10**(1/2))
55let correction=-(bias1+bias2)/2.
56. correction for probe #2362, standard dev, and standard dev of corr
57print correction sd sdcorr
58. end of calculations
59
60