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           =   2295
   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()
   '13.一元配置（反復数が異なる）．０入力は学生，それ以外は教師用ファイル作成'
 Dim bango As Integer
 bango = Val(Text1.Text)
     CommonDialog1.ShowSave
     Open CommonDialog1.FileName For Output As #1
  Dim dUncan(14 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 n0, 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, n3, n4, n5, m, p, q, dF_total, dF_shori, dF_gosa As Integer
  q = 0
  dUncan(14, 1) = 3.033
  dUncan(14, 2) = 3.178
  dUncan(14, 3) = 3.268
  dUncan(14, 4) = 3.329
  dUncan(18, 1) = 2.971
  dUncan(18, 2) = 3.118
  dUncan(18, 3) = 3.21
  dUncan(18, 4) = 3.274
  dUncan(22, 1) = 2.95
  dUncan(22, 2) = 3.097
  dUncan(22, 3) = 3.19
  dUncan(22, 4) = 3.255
  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, "<a href=#"; Format(m, "#"); ">"; m; "</a> ";
  Next m
  For m = 10 To 50
  Print #1, "<a href=#"; Format(m, "##"); ">"; m; "</a> ";
  Next m
  For m = 51 To 59
  Print #1, "<a href=#"; Format(m - 50, "#"); ">"; m; "</a> ";
  Next m
  For m = 60 To 99
  Print #1, "<a href=#"; Format(m - 50, "##"); ">"; m; "</a> ";
  Next m
  Print #1, "<pre>"
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'
n2 = 6
If n1 = 3 Then dF_total = 16
If n1 = 4 Then dF_total = 21
If n1 = 5 Then dF_total = 26
    dF_shori = n1 - 1
    dF_gosa = dF_total - dF_shori
    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:
'   Form1.Print "FAIL", p, m, n1  '
  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
    If i = 6 And j >= 3 Then GoTo W1
      iWa(j) = iWa(j) + iSeito(j, i)
      iWa2 = iWa2 + iSeito(j, i) ^ 2
      sOwa = sOwa + iSeito(j, i)
W1:
  Next i
  Next j
  n4 = 0
  For j = 1 To n1
    If j >= 3 Then
      n3 = n2 - 1
      Else
      n3 = n2
    End If
    hEikin(j) = iWa(j) / n3
    n4 = n4 + n3 ^ 2
    sS_shori = sS_shori + iWa(j) ^ 2 / n3
  Next j
   '同一平均値チェック'
  For j = 1 To n1 - 1
  For i = 1 To n1 - j
      If hEikin(j) = hEikin(i + j) Then GoTo SEITO
  Next i
    cF = sOwa ^ 2 / (dF_total + 1)
  Next j
    iWa2 = iWa2 - cF
    sS_shori = sS_shori - cF
    sS_gosa = iWa2 - sS_shori
    mS_shori = sS_shori / dF_shori
    mS_gosa = sS_gosa / dF_gosa
    fValue = mS_shori / mS_gosa
    n0 = (1 / (n1 - 1)) * ((dF_total + 1) - n4 / (dF_total + 1))
 '    Form1.Print m, p, n1, cF, s_shori, sS_gosa, n0, sM '
    sM = (mS_gosa / n0) ^ (1 / 2)
     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, "<a name=#"; Format(m, "#"); "></a>"
  If m > 9 And p = 1 Then Print #1, "<a name=#"; Format(m, "##"); "></a>"
   Print #1, m; "-"; p
   Print #1, m + 50; "-"; p
   For j = 1 To n1
   For i = 1 To n2
    If i = 6 And j >= 3 Then GoTo W2
    Print #1, iSeito(j, i);
W2:
   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, "<BR>";
   For j = 1 To n1
   For i = 1 To n2
   If i = 6 And j >= 3 Then GoTo W3
     Print #1, iSeito(j, i);
W3:
   Next i
    Print #1, "<BR>";
   Next j
  If p = 3 Then
   Print #1, " No Answer<BR>'"
   GoTo TOBASU
  End If
  If p = 2 Then
   GoTo FDAKE
  End If
   Print #1, "総和", "&nbsp;&nbsp;&nbsp", "総2乗和", "&nbsp;&nbsp;&nbsp", "CF  ;";
   Print #1, sOwa, "&nbsp;&nbsp;&nbsp", iWa2 + cF, "&nbsp;&nbsp;&nbsp", Format(cF, "###.#######"), "<BR>";
   Print #1, "-------  ANOVA Table  ---- <BR>";
   Print #1, "Factor", "&nbsp;&nbsp;&nbsp", "DF ", "&nbsp;&nbsp;&nbsp", "SS", "&nbsp;&nbsp;&nbsp", "MS", "&nbsp;&nbsp;&nbsp", "F<BR>";
   Print #1, "Total ", "&nbsp;&nbsp;&nbsp", dF_total, "&nbsp;&nbsp;&nbsp", Format(iWa2, "###.###"), "<BR>";
   Print #1, "Shori ", "&nbsp;&nbsp;&nbsp;&nbsp", dF_shori, "&nbsp;&nbsp;&nbsp", Format(sS_shori, "###.####"), "&nbsp;&nbsp;&nbsp", Format(mS_shori, "###.####"), "&nbsp;&nbsp;&nbsp", Format(fValue, "###.####"), "<BR>";
   Print #1, "Error ", "&nbsp;&nbsp;&nbsp", dF_gosa, "&nbsp;&nbsp;&nbsp", Format(sS_gosa, "###.####"), "&nbsp;&nbsp;&nbsp", Format(mS_gosa, "###.####"), "<BR>";
   '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, "<BR>";
  Next i
  GoTo TOBASU
FDAKE:
  Print #1, "Ｆ値  ";
  Print #1, Format(fValue, "##.####");
  Print #1, "<BR>";
TOBASU:
 Next p
  q = q + 1
 Next m
  Close #1
  Form1.Print
  Form1.Print "OK"
End Sub