VERSION 5.00
Object = "{F9043C88-F6F2-101A-A3C9-08002B2F49FB}#1.2#0"; "COMDLG32.OCX"
Begin VB.Form Form1
Caption = "Form1"
ClientHeight = 1830
ClientLeft = 60
ClientTop = 345
ClientWidth = 4530
LinkTopic = "Form1"
ScaleHeight = 1830
ScaleWidth = 4530
ShowInTaskbar = 0 'False
StartUpPosition = 3 'Windows の既定値
Begin VB.TextBox Text2
Height = 270
Left = 360
TabIndex = 2
Text = "４．ブロック化していない一元配置の分散分析"
Top = 480
Width = 3735
End
Begin MSComDlg.CommonDialog CommonDialog1
Left = 3840
Top = 1200
_ExtentX = 847
_ExtentY = 847
_Version = 393216
End
Begin VB.CommandButton Command1
Caption = "０入力；学生用ファイル(*.js)，その他の値；先生用ファイル(*.txt)作成 "
Height = 855
Left = 1320
TabIndex = 1
Top = 960
Width = 2055
End
Begin VB.TextBox Text1
Height = 495
Left = 240
TabIndex = 0
Top = 1200
Width = 735
End
End
Attribute VB_Name = "Form1"
Attribute VB_GlobalNameSpace = False
Attribute VB_Creatable = False
Attribute VB_PredeclaredId = True
Attribute VB_Exposed = False
Private Sub command1_click()
'分散分析，ブロック化していない一元配置．反復数が異なる．０入力は学生，それ以外は教師用ファイル作成'
Dim bango As Integer
bango = Val(Text1.Text)
CommonDialog1.ShowSave
Open CommonDialog1.FileName For Output As #1
Dim dUncan(15 To 30, 1 To 4) As Single
Dim hEikin(5) As Single
Dim dun(1 To 5) As Single
Dim iSeito(1 To 5, 1 To 7) As Integer
Dim iWa(5) As Integer
Dim sDun(6, 6) As String
Dim f1, f2, tValue, lSd, sA, du, sM, k, cF, hEihowa1, sS_shori, sS_gosa, mS_shori, mS_gosa, fValue, wOrk As Single
Dim yr, d, sOwa, iWa2, i, j, n1, n2, m, p, q, dF_total, dF_shori, dF_gosa As Integer
q = 0
dUncan(15, 1) = 3.014
dUncan(15, 2) = 3.16
dUncan(15, 3) = 3.25
dUncan(15, 4) = 3.312
dUncan(16, 1) = 2.998
dUncan(16, 2) = 3.144
dUncan(16, 3) = 3.235
dUncan(16, 4) = 3.298
dUncan(18, 1) = 2.971
dUncan(18, 2) = 3.118
dUncan(18, 3) = 3.21
dUncan(18, 4) = 3.274
dUncan(24, 1) = 2.919
dUncan(24, 2) = 3.066
dUncan(24, 3) = 3.16
dUncan(24, 4) = 3.226
dUncan(30, 1) = 2.888
dUncan(30, 2) = 3.035
dUncan(30, 3) = 3.131
dUncan(30, 4) = 3.199
Form1.AutoRedraw = True
yr = Val(Mid$(Date$, 3, 4))
Print yr
Randomize (yr + 14)
If bango = 0 Then
GoTo NOBANGO
Else
For m = 1 To 9
Print #1, ""; m; " ";
Next m
For m = 10 To 50
Print #1, ""; m; " ";
Next m
For m = 51 To 59
Print #1, ""; m; " ";
Next m
For m = 60 To 99
Print #1, ""; m; " ";
Next m
Print #1, ""
End If
NOBANGO:
For m = 1 To 50 '学生数'
If bango <> 0 Then Print #1, "-------------------------"
For p = 1 To 3 '繰り返し数'
'以下処理数ｎ1のサイクル変更'
d = m Mod 3
If p = 1 And d = 1 Then n1 = 5
If p = 1 And d = 2 Then n1 = 4
If p = 1 And d = 0 Then n1 = 3
If p = 2 And d = 1 Then n1 = 4
If p = 2 And d = 2 Then n1 = 3
If p = 2 And d = 0 Then n1 = 5
If p = 3 And d = 1 Then n1 = 3
If p = 3 And d = 2 Then n1 = 5
If p = 3 And d = 0 Then n1 = 4
'以下反復数n2'
If n1 = 3 Then n2 = 7
If n1 = 4 Then n2 = 7
If n1 = 5 Then n2 = 7
dF_total = n1 * n2 - 1
dF_shori = n1 - 1
dF_gosa = n1 * n2 - n1
f1 = 1.004 + 2.001 / dF_gosa + 41.79 / (dF_gosa ^ 2) + 2.318 / dF_shori - 0.957 / (dF_shori ^ 2) - 1.659 / (dF_shori * dF_gosa) + 2.395 / dF_gosa ^ (1 / 2) + 1.241 / dF_shori ^ (1 / 2)
f2 = -0.231 - 35.08 / dF_gosa + 214.2 / (dF_gosa ^ 2) + 4.645 / dF_shori - 1.552 / (dF_shori ^ 2) + 11.97 / (dF_shori * dF_gosa) + 16.6 / dF_gosa ^ (1 / 2) + 1.932 / dF_shori ^ (1 / 2)
SEITO:
For j = 1 To n1
iWa(j) = 0
Next j
sOwa = 0
iWa2 = 0
sS_shori = 0
For j = 1 To n1
For i = 1 To n2
If n1 <> 3 And j = 2 Then k = 0 Else k = 2
If n1 = 3 Then k = 2
iSeito(j, i) = 10 + Rnd * 10 + k * j 'データ作成'
Next i
Next j
For j = 1 To n1
For i = 1 To n2
iWa(j) = iWa(j) + iSeito(j, i)
iWa2 = iWa2 + iSeito(j, i) ^ 2
sOwa = sOwa + iSeito(j, i)
Next i
Next j
'同一平均値チェック'
For j = 1 To n1 - 1
For i = 1 To n1 - j
If iWa(j) = iWa(i + j) Then GoTo SEITO
Next i
Next j
For j = 1 To n1
hEikin(j) = iWa(j) / n2
sS_shori = sS_shori + iWa(j) ^ 2
Next j
cF = sOwa ^ 2 / (n1 * n2)
iWa2 = iWa2 - cF
sS_shori = sS_shori / n2 - cF
sS_gosa = iWa2 - sS_shori
mS_shori = sS_shori / dF_shori
mS_gosa = sS_gosa / dF_gosa
fValue = mS_shori / mS_gosa
sM = (mS_gosa / n2) ^ (1 / 2)
' Form1.Print m, p, n1, cF, fValue '
If Val(Len(cF)) < 10 Then GoTo SEITO
'n1=5は**'
If n1 = 5 And fValue < (f2 + 0.05) Then GoTo SEITO
'n1=4は**'
If n1 = 4 And fValue < (f2 + 0.05) Then GoTo SEITO
'n1=3は*'
If n1 = 3 And (fValue > (f2 - 0.05) Or fValue < (f1 + 0.05)) Then GoTo SEITO
If bango = 0 Then GoTo FILEFORSTUDENT
' FILEOUT FOR SENSEI '
If m < 10 And p = 1 Then Print #1, ""
If m > 9 And p = 1 Then Print #1, ""
Print #1, m; "-"; p
Print #1, m + 50; "-"; p
For j = 1 To n1
For i = 1 To n2
Print #1, iSeito(j, i);
Next i
Print #1,
Next j
Print #1, "平均値; ";
For j = 1 To n1
Print #1, Format(hEikin(j), "###.####"), " ";
Next j
Print #1,
Print #1, "総和 総2乗和 補正項 処理平均の標準偏差 ; ";
Print #1, sOwa, iWa2 + cF, Format(cF, "###.######"), Format(sM, "###.####")
Print #1, "Factor DF SS MS F-value (5%, 1%) "
Print #1, "Total ", dF_total, Format(iWa2, "000.000"); ""
Print #1, "Shori ", dF_shori, Format(sS_shori, "000.000"), Format(mS_shori, "000.000"), Format(fValue, "000.000"), Format(f1, "000.000"), Format(f2, "000.000")
Print #1, "Error ", dF_gosa, Format(sS_gosa, "000.000"), Format(mS_gosa, "000.000")
'mean values sorting
For i = 1 To n1 - 1
For j = i + 1 To n1
If hEikin(i) > hEikin(j) Then GoTo SORT
wOrk = hEikin(i): hEikin(i) = hEikin(j): hEikin(j) = wOrk
SORT:
Next j
Next i
Print #1, "Sorted means ";
For i = 1 To n1
Print #1, Format(hEikin(i), "###.####"), " ";
Next i
Print #1,
'Duncan '
For i = 1 To 4
dun(i) = sM * dUncan(dF_gosa, i)
Next i
lSd = dun(1)
tValue = lSd / (sM * 2 ^ (1 / 2))
Print #1, "LSD=t*SM*Root2= ";
Print #1, Format(tValue, "##.####");
Print #1, "*";
Print #1, Format(sM, "##.####");
Print #1, "* 1.4142 = ";
Print #1, Format(lSd, "##.####")
For j = 1 To n1 - 1
For i = 1 To n1 - j
Print #1, j; "-"; n1 - i + 1; "; ";
du = dun(n1 - i - j + 1)
Print #1, "dun("; n1 - i - j + 2; ")= ", Format(du, "###.####");
sA = hEikin(j) - hEikin(n1 - i + 1)
Print #1, " SA= ", Format(sA, "###.####");
If sA > du Then GoTo D2 Else GoTo D1
D2:
Print #1, " *"
sDun(n1 - i, j) = "* "
GoTo D3
D1:
Print #1, " NS"
sDun(n1 - i, j) = "NS"
D3:
Next i
Next j
For i = 1 To n1 - 1
For j = 1 To i
Print #1, sDun(i, j); " ";
Next j
Print #1,
Next i
GoTo TOBASU
' ====================== '
FILEFORSTUDENT:
If p > 1 Then GoTo TOBU3
Print #1, "jsDB[";
Print #1, q;
Print #1, "]='";
TOBU3:
Print #1, "No ";
Print #1, p;
Print #1, "
";
For j = 1 To n1
For i = 1 To n2
Print #1, iSeito(j, i);
Next i
Print #1, "
";
Next j
If p = 3 Then
Print #1, " No Answer
'"
GoTo TOBASU
End If
If p = 2 Then
GoTo FDAKE
End If
Print #1, "総和", " ", "総2乗和", " ", "CF ;";
Print #1, sOwa, " ", iWa2 + cF, " ", Format(cF, "###.#######"), "
";
Print #1, "------- ANOVA Table ----
";
Print #1, "Factor", " ", "DF ", " ", "SS", " ", "MS", " ", "F
";
Print #1, "Total ", " ", dF_total, " ", Format(iWa2, "###.###"), "
";
Print #1, "Shori ", " ", dF_shori, " ", Format(sS_shori, "###.####"), " ", Format(mS_shori, "###.####"), " ", Format(fValue, "###.####"), "
";
Print #1, "Error ", " ", dF_gosa, " ", Format(sS_gosa, "###.####"), " ", Format(mS_gosa, "###.####"), "
";
'mean values sorting
For i = 1 To n1 - 1
For j = i + 1 To n1
If hEikin(i) > hEikin(j) Then GoTo SORT2
wOrk = hEikin(i): hEikin(i) = hEikin(j): hEikin(j) = wOrk
SORT2:
Next j
Next i
'Duncan '
For i = 1 To 4
dun(i) = sM * dUncan(dF_gosa, i)
Next i
For j = 1 To n1 - 1
For i = 1 To n1 - j
du = dun(n1 - i - j + 1)
sA = hEikin(j) - hEikin(n1 - i + 1)
If sA > du Then GoTo D22 Else GoTo D12
D22:
sDun(n1 - i, j) = "* "
GoTo D32
D12:
sDun(n1 - i, j) = "N"
D32:
Next i
Next j
For i = 1 To n1 - 1
For j = 1 To i
Print #1, sDun(i, j); " ";
Next j
Print #1, "
";
Next i
GoTo TOBASU
FDAKE:
Print #1, "Ｆ値 ";
Print #1, Format(fValue, "##.####");
Print #1, "
";
TOBASU:
Next p
q = q + 1
Next m
Close #1
Form1.Print
Form1.Print "OK"
End Sub