![[BACK]](/cvsweb/icons/back.gif){kind=icon}
Return to anovaws CVS log ![[TXT]](/cvsweb/icons/text.gif){kind=icon}
![[DIR]](/cvsweb/icons/dir.gif){kind=icon}
Up to [CSRG BSD Unix] / 43BSD / contrib / apl / lib|
Return to anovaws CVS log | Up to [CSRG BSD Unix] / 43BSD / contrib / apl / lib |
1.1 root 1: n� �*.KB�& H ss r { ss y
2: r { (+/,y*2)%X/Ry
3:
4: sigma r { v sigma y
5: r { BJ((2XRv)R'+/'),'(|H|H(IRRy)Ev)\Oy'
6:
7: delta r { v delta y;p;q
8: r { y
9: nxt:}0 if 0=Rv
10: r { (rX(Rr)[1Yv])-(1Yv) sumx r
11: v { 1Uv
12: }nxt
13:
14: if r { a if b
15: r { b/a
16:
17: rm vc { v rm y
18: vc { (~(IRRy)Ev)/IRRy
19:
20: anova m anova y;a;b;n;s;t;aF;ym;yv;ye;dfs;nh;nv;um
21: start
22: aov1
23: table
24: start start
25: n { RRy
26: um { ~'l'E(`1+mI';')Um
27: yv { (y#0)V~';'Em
28: ym { (+/y)%+/yv
29: ye { yvXy-O\(O|Ry)RO\ym
30: nh { (X/`1URy)%+/,%+/yv
31: dfs {`1+(`1URy),(+/,yv)%X/`1URy
32: a{O|nY'abcdef'
33: aF{|OnY'aFbFcFdFeFfF'
34: m{1U(mEa,aF,'+')/m
35: b { (nR2)N`1+I2*n
36: aov1 aov1;x;e;et;er
37: s { 0 2 R0
38: t { 0 10 R' '
39: e {(n,0)R0
40: loop: t{t,[1] 10 Y x{n Y (`1+mI'+')Ym
41: s{s,[1] aov2 et {1 1 2 O\(aEx)J.V((aFEx)^.&b)/b
42: e { e,et
43: }loop if 0<Rm{(mI'+')Um
44: } 0 if 0=Rer{(~(I`1+2*n)E2Be)/I`1+2*n
45: t{t,[1] 'error '
46: s{s,[1] aov2 (nR2)Ner
47: aov2 r{ aov2 x;df;s;v;vc
48: df { s { 0
49: nxt:v {x[;1]/O|In
50: vc { (~x[;1])/O|In
51: df { df+X/dfs[v]
52: }error if (RRy)Ev
53: }umean if um
54: leastsq: s { s + ss2 vc sigma2 v delta2 y
55: }test
56: umean: s { s + nhX(ss (1Uvc) sigma v delta ym)%X/Rym
57: }test
58: error: s { s +(`1YRy)X (ss vc sigma (1Uv)delta ye)%X/`1URy
59: test: }nxt if 0<`1YRx{0 1Ux
60: r{ df,s
61: table table;dfe;dftt;sse;sst
62: ' df ss ms'
63: t, 3 0 16 5 16 5 NJ s,s[;2]%s[;1]
64:
65: means r { m means y;i;v;vc;x
66: v{,'abcdef'Im
67: vc{(~(IRRy)Ev)/IRRy
68: x{,(H|v)O\(vc sigma y)%X/(Ry)[vc]
69: i{O\1+(Ry)[v]N`1+IX/(Ry)[v]
70: r {(4 0 NJ i),12 3 NJ((Rx),1)Rx
71:
72: convert r { n convert data ;x
73: x{Rdata
74: r{2 1 3O\(x[2],(x[1]%n),n)RO\data
75: & m3 # & y=m+a+b+ab+c�+E newdata newdata
76: data { ((pXn),q)R0
77: 'enter data for each subject on a new line'
78: i { 1
79: nxt: 'subject: ';i
80: t { L
81: }error XI (R,t)#1URdata
82: data[i;]{t
83: }nxt XI (1YRdata)& i{ i+1
84: 'newdata complete'
85: }0
86: error:'wrong amount of data'
87: 1URdata;' elements required'
88: }nxt
89: & dat3 �& & - @A �A B �A B pB B `B �B B @A �@ @B B @B �A @A �B B �B B B �B �B B @A �A �A �B pB �B �B �B B B @B
90: sum r { v sum y;p;q
91: r { y
92: v { (~(IRRy)Ev)/IRRy
93: nxt:}0 if 0=Rv
94: p { |Hq {|H~(IRRr)=1Yv
95: r { q\O(Rr)[q]R (+/[1] p\Or)
96: v { 1Uv
97: }nxt
98:
99: sigma2 r { v sigma2 y
100: r { v sigma y
101: nv { v sigma nv
102:
103: delta2 r { v delta2 y;vc
104: r { y
105: nv { y#0
106: nxt:}0 if 0=Rv
107: vc { ((1Yv)#IRRr)/IRRr
108: r { r-((1Yv) sumx r)X(vc sumx nv)%+/,nv
109: v { 1Uv
110: }nxt
111: ss2 r { ss2 y
112: r { +/,(y*2)%nv
113: & m3u ( & y=m+a+b+ab+c�+E ; um& m3l ( & y=m+a+b+ab+c�+E ; ls testd testd
114: dat3
115: m3
116: m3 anova dat3
117: ''
118: m3u
119: m3u anova dat3
120: ''
121: m3l
122: m3l anova dat3
123:
124: sumx r { v sumx y;p;q
125: r { y
126: nxt:}0 if 0=R,v
127: p { |Hq {|H~(IRRr)=1Yv
128: r { q\O(Rr)[q]R (+/[1] p\Or)
129: v { 1Uv
130: }nxt
131: & Llx 3 & 'type ''describe'' for details' edtxt edtxt;fd;name;vr;r
132: L'{'variable name? '
133: name { L'
134: fd { 390 Lcreat 'apltxt.tmp'
135: }(fd < 0)/error
136: r { fd Lwrite BJ name
137: fd { Lclose fd
138: Lrun 'xed apltxt.tmp'
139: fd { 260 Lopen 'apltxt.tmp'
140: }(fd<0)/error
141: vr { ''
142: more: vr { vr,r { fd Lread 256
143: } (0#Rr)/more
144: fd { Lclose fd
145: r { BJ name,' { vr'
146: Lrun 'rm apltxt.tmp'
147: }0
148: error: 'Can''t open file: apltxt.tmp'
149: & describe � & d� �
150: these functions are described in:
151: algorithms for anova with unequal cell frequencies
152: by a.p. reeves and j. a. sniezek
153: in apl quote quad vol.10, no2, pp32-34, december 1979
154:
155: for an example usage type testd
156:
157: umeans r { m umeans y;v;vc;x
158: y { (+/y)%(`1YRy)-+/y=0
159: v {,'abcdef'Im
160: vc{(~(IRRy)Ev)/IRRy
161: x { (vc sigma y)%X/(Ry)[vc]
162: r { 12 3 NJ x
163:
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.