**TI89**mainRanova2ZO() Prgm ClrIO Disp "number of rows" Input r Disp "number of columns" Input s Disp "number per cell" Input c seq(i,i,1,r)list1 {sum(seq(a[1,j],j,1,s))/s}list2 For i,2,r augment(list2,{sum(seq(a[i,j],j,1,s))/s})list2 EndFor {sum(seq(list2[i],i,1,r))/r}list5 seq(j,j,1,s)list3 {sum(seq(a[i,1],i,1,r))/r}list4 For j,2,s augment(list4,{sum(seq(a[i,j],i,1,r))/r})list4 EndFor augment(list5,{c*s*sum(seq((list2[i]-list5[1])^2,i,1,r))})list5 augment(list5,{c*r*sum(seq((list4[j]-list5[1])^2,j,1,s))})list5 {sum(seq((a[1,j]+list5[1]-list2[1]-list4[j])^2,j,1,s))}list6 For i,2,r augment(list6,{sum(seq((a[i,j]+list5[1]-list2[i]-list4[j])^2,j,1,s))})list6 EndFor augment(list5,{c*sum(seq(list6[i],i,1,r))})list5 If c=1 Then (s-1)*list5[2]/list5[4]u (r-1)*list5[3]/list5[4]v Else {sum(seq((c-1)*b[1,j]^2,j,1,s))}list6 For i,2,r augment(list6,{sum(seq((c-1)*b[i,j]^2,j,1,s))})list6 EndFor augment(list5,{sum(seq(list6[i],i,1,r))})list5 (list5[2]/(r-1))/(list5[5]/(r*s*(c-1)))u (list5[3]/(s-1))/(list5[5]/(r*s*(c-1)))v (list5[4]/(r-1)/(S-1))/(list5[5]/(r*s*(c-1)))w EndIf NewData dist,list1,round(list2),list3,round(list4),round(list5) r-1m s-1o (r-1)*(s-1)z If c=1 Then (r-1)*(s-1)n Else r*s*(c-1)n EndIf 1-s\fcdff(0,u,m,n)t 1-s\fcdff(0,v,o,n)p If c1 Then 1-s\fcdff(0,w,z,n)q EndIf ClrIO Output 1,1,"row fstat" Output 1,65,round(u) Output 10,1,"row p val" Output 10,65,round(t) Output 26,1,"col fstat" Output 26,65,round(v) Output 36,1,"col p val" Output 36,65,round(p) If C1 Then Output 52,1,"int fstat" Output 52,65,round(w) Output 62,1,"int p val" Output 62,65,round(q) EndIf EndPrgm ܑ