var unx = 1*0; var xjz = " "; var zjk = ""; var zbz = 4*1; var zbj = Math.pow(10,zbz); function rnd2(x){return Math.round(x*100)/100} function rnd4(x){return Math.round(x*zbj)/zbj} function unw() { if(analysis==2) { alert("Unweighted means analysis is not available for correlated samples.") }else{ document.setup.elements[4].value = "unweighted-means analysis"; document.qxr.elements[0].value = "unweighted-means analysis"; unx = 1*1; for(var i=0; i<13; i++) { document.f.elements[i].value = ""; } document.tuk.elements[0].value = "This test will be performed only if K>2
and the analysis of variance yields a
significant F-ratio."; } } function wt() { document.setup.elements[4].value = "standard weighted-means analysis"; document.qxr.elements[0].value = "standard weighted-means analysis"; unx = 1*0; for(var i=0; i<13; i++) { document.f.elements[i].value = ""; } document.tuk.elements[0].value = "This test will be performed only if K>2
and the analysis of variance yields a
significant F-ratio."; } function calc() { var x1=0, x2=0, x3=0, x4=0, x5=0; var sx1=0, sx2=0, sx3=0, sx4=0, sx5=0; var sum1=0, sum2=0, sum3=0, sum4=0, sum5=0, sumt=0; var sumsq1=0, sumsq2=0, sumsq3=0, sumsq4=0, sumsq5=0, sumsqt=0; var ss1=0, ss2=0, ss3=0, ss4=0, ss5=0, sst=0; var n1=0, n2=0, n3=0, n4=0, n5=0, nt=0; var m1=0, m2=0, m3=0, m4=0, m5=0, mt=0; var sswg=0, ssbg=0, sserr=0; var dfbg=0, dfwg=0, dft=0, dfss=0, dferr=0; var subj = new Array(); var sumsub=0, sssub=0; var mswg=0, msbg=0, msss=0, mserr=0, ff=0; zjk = ""; var nps = 1*0; var d12=1*0, d13=1*0, d14=1*0, d15=1*0, d23=1*0, d24=1*0, d25=1*0, d34=1*0, d35=1*0, d45=1*0; zbz = 4*1; zbj = Math.pow(10,zbz); for(var i=0; i<30; i++) { document.show1.elements[i].value = ""; } for(var i=0; i<12; i++) { document.f.elements[i].value = ""; } document.f.elements[0].value = document.setup.elements[3].value; for(var i=0; i<30; i++) { document.show1.elements[i].value = ""; } if(document.setup.elements[3].value == "") { alert("Please perform the Setup operation before proceeding."); } var vals1 = document.a1.elements[0].value; if(vals1.indexOf("\r")>-1){xjz="\r"}else{xjz="\n"} x1 = document.a1.elements[0].value; sx1 = x1.split(xjz); for(var i=0; i0) {n1=n1+1;} } x2 = document.a1.elements[1].value; sx2 = x2.split(xjz); for(var i=0; i0) {n2=n2+1;} } x3 = document.a1.elements[2].value; sx3 = x3.split(xjz); for(var i=0; i0) {n3=n3+1;} } x4 = document.a1.elements[3].value; sx4 = x4.split(xjz); for(var i=0; i0) {n4=n4+1;} } x5 = document.a1.elements[4].value; sx5 = x5.split(xjz); for(var i=0; i0) {n5=n5+1;} } for(var i=0; i } if(k>2) { m3 = sum3/n3; ss3 = sumsq3-((sum3*sum3)/n3); sswg = sswg+ss3*1; nt = nt+n3*1; sumt = sumt+sum3*1; sumsqt = sumsqt+sumsq3*1; document.show1.elements[2].value = n3; document.show1.elements[8].value = rnd4(sum3); document.show1.elements[14].value = rnd4(m3); document.show1.elements[20].value = rnd4(sumsq3); document.show1.elements[26].value =rnd4(ss3/(n3-1)); document.show1.elements[32].value =rnd4(Math.sqrt(ss3/(n3-1))); document.show1.elements[38].value =rnd4(Math.sqrt((ss3/(n3-1))/n3)); if(analysis==2) { for(var i=0; i } } if(k>3) { m4 = sum4/n4; ss4 = sumsq4-((sum4*sum4)/n4); sswg = sswg+ss4*1; nt = nt+n4*1; sumt = sumt+sum4*1; sumsqt = sumsqt+sumsq4*1; document.show1.elements[3].value = n4; document.show1.elements[9].value = rnd4(sum4); document.show1.elements[15].value = rnd4(m4); document.show1.elements[21].value = rnd4(sumsq4); document.show1.elements[27].value =rnd4(ss4/(n4-1)); document.show1.elements[33].value =rnd4(Math.sqrt(ss4/(n4-1))); document.show1.elements[39].value =rnd4(Math.sqrt((ss4/(n4-1))/n4)); if(analysis==2) { for(var i=0; i } } if(k>4) { m5 = sum5/n5; ss5 = sumsq5-((sum5*sum5)/n5); sswg = sswg+ss5*1; nt = nt+n5*1; sumt = sumt+sum5*1; sumsqt = sumsqt+sumsq5*1; document.show1.elements[4].value = n5; document.show1.elements[10].value = rnd4(sum5); document.show1.elements[16].value = rnd4(m5); document.show1.elements[22].value = rnd4(sumsq5); document.show1.elements[28].value =rnd4(ss5/(n5-1)); document.show1.elements[34].value =rnd4(Math.sqrt(ss5/(n5-1))); document.show1.elements[40].value =rnd4(Math.sqrt((ss5/(n5-1))/n5)); if(analysis==2) { for(var i=0; i } } sst = sumsqt-((sumt*sumt)/nt); ssbg = sst-sswg; dft = nt-1; dfbg = k-1; dfwg = nt-k; mt = sumt/nt; document.show1.elements[5].value = nt; document.show1.elements[11].value = rnd4(sumt); document.show1.elements[17].value = rnd4(mt); document.show1.elements[23].value = rnd4(sumsqt); document.show1.elements[29].value =rnd4(sst/(nt-1)); document.show1.elements[35].value =rnd4(Math.sqrt(sst/(nt-1))); document.show1.elements[41].value =rnd4(Math.sqrt((sst/(nt-1))/nt)); mswg = sswg/dfwg; mserr = mswg; if(n1>0){nps+=1/n1} if(n2>0){nps+=1/n2} if(n3>0&&k>2){nps+=1/n3} if(n4>0&&k>3){nps+=1/n4} if(n5>0&&k>4){nps+=1/n5} nps = k/nps; var mg = (m1+m2+m3+m4+m5)/k; var smg = 1*0; if(analysis==1*1&&unx==1*1) { if(k>0){smg+=(m1-mg)*(m1-mg)} if(k>1){smg+=(m2-mg)*(m2-mg)} if(k>2){smg+=(m3-mg)*(m3-mg)} if(k>3){smg+=(m4-mg)*(m4-mg)} if(k>4){smg+=(m5-mg)*(m5-mg)} smg = smg*nps; ssbg = smg*1; } msbg = ssbg/dfbg; document.f.elements[1].value = rnd4(ssbg); document.f.elements[2].value = dfbg; document.f.elements[3].value = rnd4(msbg); document.f.elements[11].value = rnd4(sst); document.f.elements[12].value = dft; dferr = dfwg; sserr = sswg; var zzz = 0*1; if(analysis==2) { for(var i=0; i sssub = (sumsub/k)-((sumt*sumt)/nt); sssub = Math.round(sssub*10000)/10000; dfss = n1-1; dferr = dfwg-dfss; msss = Math.round((sssub/dfss)*10000)/10000; sserr = Math.round((sswg-sssub)*10000)/10000; mserr = Math.round((sserr/dferr)*10000)/10000; document.f.elements[9].value = rnd4(sssub); document.f.elements[10].value = dfss; document.f.elements[12].value = dft; if(n1!=n2){zzz++} if(k>=3&&n2!=n3){zzz++} if(k>=4&&n3!=n4){zzz++} if(k>=5&&n4!=n5){zzz++} } if(zzz>0){ alert("The correlated-samples design requires equal n's per sample."); } if(mserr>0) {ff=msbg/mserr} else {ff=1*0} document.f.elements[6].value = rnd4(sserr); document.f.elements[7].value = dferr; document.f.elements[8].value = rnd4(mserr); document.f.elements[4].value = rnd2(ff); if(document.setup.elements[3].value == "") { for(var i=0; i<30; i++) { document.show1.elements[i].value = ""; } } Fprob(); if(pj==1*1&&k>2) { var t = 1*1.6; var crit0 = 1*0; var t05 = 1*0; var t01 = 1*0; var q05 = 1*0; var q01 = 1*0; var hsd05 = 1*0; var hsd01 = 1*0; var df = dferr; var tukey = ""; for (var j=0; j<10000; j++) { t+=.005; crit0 = Tcall(buzz(eval(t),eval(df))); if(crit0<=.05){t05=t;break;} } for (var j=0; j<10000; j++) { t+=.005; crit0 = Tcall(buzz(eval(t),eval(df))); if(crit0<=.01){t01=t;break;} } hsd05 = Math.round(hsd05*100)/100; if(k==3) { q05 = 2.106973*t05-0.817944; q01 = 1.946821*t01-0.8996111; } if(k==4) { q05 = 2.59017*t05-1.443969; q01 = 2.303537*t01-1.534593; } if(k==5) { q05 = 2.972018*t05-1.96814; q01 = 2.587002*t01-2.064423; } hsd05 = q05*Math.sqrt(mserr/nps); hsd01 = q01*Math.sqrt(mserr/nps); hsd05 = Math.round(hsd05*100)/100; hsd01 = Math.round(hsd01*100)/100; tukey+="HSD[.05]="+hsd05+"; HSD[.01]="+hsd01+"\r"; d12 = Math.abs(m1-m2); d13 = Math.abs(m1-m3); d14 = Math.abs(m1-m4); d15 = Math.abs(m1-m5); d23 = Math.abs(m2-m3); d24 = Math.abs(m2-m4); d25 = Math.abs(m2-m5); d34 = Math.abs(m3-m4); d35 = Math.abs(m3-m5); d45 = Math.abs(m4-m5); var gap = ""; if(d12>=hsd05){gap="P<.05"}else{gap="nonsignificant"} if(d12>=hsd01){gap="P<.01"} tukey+="M1 vs M2 "+gap+"\r"; if(d13>=hsd05){gap="P<.05"}else{gap="nonsignificant"} if(d13>=hsd01){gap="P<.01"} tukey+="M1 vs M3 "+gap+"\r"; if(k>3) { if(d14>=hsd05){gap="P<.05"}else{gap="nonsignificant"} if(d14>=hsd01){gap="P<.01"} tukey+="M1 vs M4 "+gap+"\r"; } if(k>4) { if(d15>=hsd05){gap="P<.05"}else{gap="nonsignificant"} if(d15>=hsd01){gap="P<.01"} tukey+="M1 vs M5 "+gap+"\r"; } if(d23>=hsd05){gap="P<.05"}else{gap="nonsignificant"} if(d23>=hsd01){gap="P<.01"} tukey+="M2 vs M3 "+gap+"\r"; if(k>3) { if(d24>=hsd05){gap="P<.05"}else{gap="nonsignificant"} if(d24>=hsd01){gap="P<.01"} tukey+="M2 vs M4 "+gap+"\r"; } if(k>4) { if(d25>=hsd05){gap="P<.05"}else{gap="nonsignificant"} if(d25>=hsd01){gap="P<.01"} tukey+="M2 vs M5 "+gap+"\r"; } if(k>3) { if(d34>=hsd05){gap="P<.05"}else{gap="nonsignificant"} if(d34>=hsd01){gap="P<.01"} tukey+="M3 vs M4 "+gap+"\r"; } if(k>4) { if(d35>=hsd05){gap="P<.05"}else{gap="nonsignificant"} if(d35>=hsd01){gap="P<.01"} tukey+="M3 vs M5 "+gap+"\r"; if(d45>=hsd05){gap="P<.05"}else{gap="nonsignificant"} if(d45>=hsd01){gap="P<.01"} tukey+="M4 vs M5 "+gap+"\r"; } document.tuk.elements[0].value = tukey; } } function unwipe() { if(analysis==1) {zjk="Independent"}; if(analysis==2) {zjk="Correlated"}; var jsum = new Array(42); for(var i=0; i<42; i++){jsum[i] = document.show1.elements[i].value} var n1 = document.show1.elements[0].value; var n2 = document.show1.elements[1].value; var n3 = document.show1.elements[2].value; var n4 = document.show1.elements[3].value; var n5 = document.show1.elements[4].value; var nt = document.show1.elements[5].value; var sum1 = document.show1.elements[6].value; var sum2 = document.show1.elements[7].value; var sum3 = document.show1.elements[8].value; var sum4 = document.show1.elements[9].value; var sum5 = document.show1.elements[10].value; var sumt = document.show1.elements[11].value; var sumsq1 = document.show1.elements[12].value; var sumsq2 = document.show1.elements[13].value; var sumsq3 = document.show1.elements[14].value; var sumsq4 = document.show1.elements[15].value; var sumsq5 = document.show1.elements[16].value; var sumsqt = document.show1.elements[17].value; var ss1 = document.show1.elements[18].value; var ss2 = document.show1.elements[19].value; var ss3 = document.show1.elements[20].value; var ss4 = document.show1.elements[21].value; var ss5 = document.show1.elements[22].value; var sst = document.show1.elements[23].value; var m1 = document.show1.elements[24].value; var m2 = document.show1.elements[25].value; var m3 = document.show1.elements[26].value; var m4 = document.show1.elements[27].value; var m5 = document.show1.elements[28].value; var mt = document.show1.elements[29].value; var ssbg = document.f.elements[1].value; var dfbg = document.f.elements[2].value; var msbg = document.f.elements[3].value; var ff = document.f.elements[4].value; var pp = document.f.elements[5].value; var sserr = document.f.elements[6].value; var dferr = document.f.elements[7].value; var mserr = document.f.elements[8].value; var ssss = document.f.elements[9].value; var dfss = document.f.elements[10].value; var sst = document.f.elements[11].value; var dft = document.f.elements[12].value; var tuck1 = ""; var tuck2 = document.tuk.elements[0].value.split(xjz); for(var i=0; i";} var jx1 = ""; var jz1 = document.a1.elements[0].value.split(xjz); for(var i=0; i";} var jx2 = ""; var jz2 = document.a1.elements[1].value.split(xjz); for(var i=0; i";} var jx3 = ""; var jz3 = document.a1.elements[2].value.split(xjz); if(k>=3) {for(var i=0; i";} } else {jx3="---";} var jx4 = ""; var jz4 = document.a1.elements[3].value.split(xjz); if(k>=4) {for(var i=0; i";} } else {jx4="---";} var jx5 = ""; var jz5 = document.a1.elements[4].value.split(xjz); if(k==5) {for(var i=0; i";} } else {jx5="---";} var date = new Date(); var aa = 'VassarStats Printable Report
One-Way ANOVA for '+k+' '+zjk+' Samples
'+date+'
'; document.open(); document.write(''); document.write(''+aa); document.write('Values Entered
Sample 1
Sample 2
Sample 3
Sample 4
Sample 5
'+jx1+''+jx2+''+jx3+''+jx4+''+jx5+'
'); document.write('
Data Summary
Samples
1
2
3
4
5
Total
N  
'+jsum[0]+'
'+jsum[1]+'
'+jsum[2]+'
'+jsum[3]+'
'+jsum[4]+'
'+jsum[5]+'
-X  
'+jsum[6]+'
'+jsum[7]+'
'+jsum[8]+'
'+jsum[9]+'
'+jsum[10]+'
'+jsum[11]+'
-Mean  
'+jsum[12]+'
'+jsum[13]+'
'+jsum[14]+'
'+jsum[15]+'
'+jsum[16]+'
'+jsum[17]+'
-X2  
'+jsum[18]+'
'+jsum[19]+'
'+jsum[20]+'
'+jsum[21]+'
'+jsum[22]+'
'+jsum[23]+'
Variance  
'+jsum[24]+'
'+jsum[25]+'
'+jsum[26]+'
'+jsum[27]+'
'+jsum[28]+'
'+jsum[29]+'
Std.Dev.  
'+jsum[30]+'
'+jsum[31]+'
'+jsum[32]+'
'+jsum[33]+'
'+jsum[34]+'
'+jsum[35]+'
Std.Err.  
'+jsum[36]+'
'+jsum[37]+'
'+jsum[38]+'
'+jsum[39]+'
'+jsum[40]+'
'+jsum[41]+'
'); document.write('
ANOVA Summary
Source
SS
df
MS
F
P
Treatment
[between groups]
'+ssbg+'
'+dfbg+'
'+msbg+'
'+ff+'
'+pp+'
Error
'+sserr+'
'+dferr+'
'+mserr+'
Ss/Bl
'+ssss+'
'+dfss+'
Total
'+sst+'
'+dft+'
Ss/Bl = Subjects or Blocks depending on the design.
Applicable only to correlated-samples ANOVA.
'); document.write('Tukey HSD Test
'+tuck1+'M1 = mean of Sample 1
M2 = mean of Sample 2
and so forth.
HSD = the absolute [unsigned]
difference between any two
sample means required for
significance at the designated
level. HSD[.05] for the .05 level;
HSD[.01] for the .01 level.
'); document.write('
HomeClick this link only if you did not arrive here via the VassarStats main page.
'); document.close(); }