Annotation of 43BSD/contrib/apl/lib/anovaws, revision 1.1

1.1     ! root        1: n��*.KB�&       H ssr { ss y
        !             2: r { (+/,y*2)%X/Ry
        !             3: 
        !             4: sigmar { v sigma y
        !             5: r { BJ((2XRv)R'+/'),'(|H|H(IRRy)Ev)\Oy'
        !             6: 
        !             7: deltar { 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: ifr { a if b
        !            15: r { b/a
        !            16: 
        !            17: rmvc { v rm y
        !            18: vc { (~(IRRy)Ev)/IRRy
        !            19: 
        !            20: anovam anova y;a;b;n;s;t;aF;ym;yv;ye;dfs;nh;nv;um
        !            21: start
        !            22: aov1
        !            23: table
        !            24: startstart
        !            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: aov1aov1;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:       aov2r{ 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: tabletable;dfe;dftt;sse;sst
        !            62: '           df          ss              ms'
        !            63: t, 3 0 16 5 16 5 NJ s,s[;2]%s[;1]
        !            64: 
        !            65: meansr { 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: convertr { 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�+Enewdatanewdata
        !            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 BpBB`B�BB@A�@@BB@B�A@A�B B�BBB�B�B B@A�A�A�BpB�B�B�BB B@B
        !            90: sumr { 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: sigma2r { v sigma2 y
        !           100: r { v sigma y
        !           101: nv { v sigma nv
        !           102: 
        !           103: delta2r { 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:       ss2r { 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 ; lstestdtestd
        !           114: dat3
        !           115: m3
        !           116: m3 anova dat3
        !           117: ''
        !           118: m3u
        !           119: m3u anova dat3
        !           120: ''
        !           121: m3l
        !           122: m3l anova dat3
        !           123: 
        !           124: sumxr { 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: &Llx3&  '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: umeansr { 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: 

unix.superglobalmegacorp.com

This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.