メイン　permtest4.html　　　Last modified: Sep 01, 2009
<html>
<head>
<meta http-equiv="Content-Type" content="text/html;CHARSET=EUC-JP">
<link rel="shortcut icon" href="../favicon.ico">
<title>JavaScript</title>
<script src="permutation.js">document.write("permutation.js ファイルが見つかりません??<br>")</script>
<script src="io.js">document.write("io.js ファイルが見つかりません??<br>")</script>

<script language="JavaScript">
<!--
//-----------　大域変数　-----------------------
// var array0, nr, nc, um, t, ntrue, denom, stat_val, nc1, nr1, array1, array2
// var EPSILON = 1e-14
// var MAXLOOP = 2000
// var CRITERION = 5000

function compareNumbers(a, b)
{
//  return a - b
  if (a > b) return 1
  else if (a == b) return 0
  else return -1
}

// median は直接関係しない
function median(value, n)
{
  var retv
  retv = 0
  value.sort(compareNumbers)
  if (Math.floor(n/2) != n/2) {
    return value[Math.floor(n/2)]
  }
  else { 
    return (value[Math.floor(n/2)]+value[Math.floor(n/2)-1])/2
  }
}


//-----------　検定に応じて定義　-----------------------
// stat 関数が引用される

function stat(array1, n1, array2, n2)
{
  var i, u
  u = 0
  for (i = 0; i < n1; i++) {
    u += array1[i]
  }
  return Math.abs((n1*n2+n1*(n1+1))/2-u)
}

//-----------　メイン　-----------------------
function calc(data1_string, data2_string)
{
  var n, n1, n2, data1, data2, i, j, array9, median1, median2, less, equal
  array1 = new Array()
  array2 = new Array()
  array0 = new Array()

  if ((data1 = getdata(data1_string, -1)) != false && (data2 = getdata(data2_string, -1)) != false) {
    n1 = decomp(data1, array1)
    n2 = decomp(data2, array2)
    n = n1+n2
    nr = 2
    nc = n

// 元データ
    array9 = new Array(n)

    for (i = 0; i < n1; i++) {
      array9[i] = array1[i]
    }
    for (i = 0; i < n2; i++) {
      array9[n1+i] = array2[i]
    }
    median1 = median(array1, n1)
    median2 = median(array2, n2)

// データではなくて順位を使う
    for (i = 0; i < n; i++) {
      equal = less = 0
      for (j = 0; j < n; j++) {
        if (array9[j] == array9[i]) {
          equal++
        }
        else if (array9[j] < array9[i]) {
          less++
        }
      }
      array0[i] = less+(equal+1)/2
    }
    for (i = 0; i < n1; i++) {
      array1[i] = array0[i]
    }
    for (i = 0; i < n2; i++) {
      array2[i] = array0[n1+i]
    }


// 以下の4行は検定に応じて追加・修正を

    printf("第一群の中央値 = %g (n1 = %i)\n", median1, n1)
    printf("第二群の中央値 = %g (n2 = %i)\n", median2, n2)
    stat_val = stat(array1, n1, array2, n2)
    printf("観察されたデータにおける U = %g\n", n1*n2/2-stat_val)

    if (com(n1+n2, n1) < CRITERION || confirm("計算時間がかかりそうです。正確な P 値を求めるときは「OK」，近似解を求めるときは「Cancel」をクリックします。")) {
      exact(n1, n2)
      printf("並べ替え検定による正確な P 値 = %.10g　( = %i / %i )\n", ntrue / denom, ntrue, denom)
    }
    else {
      approx(n1, n2)
      printf("並べ替え検定による近似的な P 値 = %.10g　( = %i / %i )\n", ntrue / denom, ntrue, denom)
    }
  }
  sep(40)
}
//-->
</script>
</head>

<body bgcolor="#ffffff">
<font size="+2"><b>並べ替え検定（U 統計量に基づく）</b></font>　<a  href="src/permtest4.html"><img src="png/src.png" width=35 height=11 alt="src" align=top></a>　　　　　Last modified: Jun 01, 2006<hr noshade><p>
<font color="#ff0000" size="+2">以下のプログラムのサポートは終了しました。自己責任でお使い下さい。</font>

<form name=Result>
<script language="JavaScript">
<!--
//-->JavaScript がサポートされていないブラウザですか？
</script>
注意：データの個数が多くなるにつれ，計算時間が長くなります。
<p>
<table><tr>
<td><input type="button" name="calcurate" value="計算開始" onClick="calc(this.form.data1.value,this.form.data2.value)">　　</td>
<td><input type="button" name="clear" value="入力欄クリア" onClick="this.form.data1.value=this.form.data2.value=''">　　</td>
<td><input type="button" name="clear" value="出力欄クリア" onClick="this.form.result.value=''">　　</td>
<td><img src="../gra/button3.png" width=9 height=9 alt="・"> <a href="exa/ttest2.html">使用法</a></td>
</tr></table>

<p>
<table>
<tr>
<td valign=top>第一群のデータ<br><textarea name="data1" rows=10 cols=20></textarea><p>
第二群のデータ<br><textarea name="data2" rows=10 cols=20></textarea></td>
<td valign=top rowspan=2>出力欄<br><textarea name="result" rows=24 cols=80></textarea></td>
</tr>
</table>
</form>

<p><hr noshade>
<img src="../gra/button3.png" width=9 height=9 alt="・"> <A HREF="javascript:history.go(-1)">直前のページへ戻る</A>　　<img src="../gra/button3.png" width=9 height=9 alt="・"> <a href="../mail.html">E-mail to Shigenobu AOKI</a>
<p><center><IMG SRC="../gra/ume5.png" width=121 height=37 ALT="Made with Macintosh"></center>
</body>
</html>

サブ　io.js　　　Last modified: Mar 25, 2004
サブ　permutation.js　　　Last modified: Mar 25, 2004
・ 直前のページへ戻る　　・ E-mail to Shigenobu AOKI

Made with Macintosh