• 検索結果がありません。

一対比較データの解析にも役立つ定差方程式

N/A
N/A
Protected

Academic year: 2021

シェア "一対比較データの解析にも役立つ定差方程式"

Copied!
2
0
0

読み込み中.... (全文を見る)

全文

(1)

健一!)IQ

メモランダム

このコラムは, ORにかかわる概念,知識(手法,原理),それらの図解.よい教材や問題,実学 ORの実施経験.そこから得られた智恵やアドバイス,失敗談と教訓,新しい視札視座,フレー ムワーク.未だ解けていない問題.面白い研究テー?などを,“新鮮にしかも“コンパクトに" 表現し,提示していただくものです.ユニークなアイディア,フレッシュな見方.発想,だれかと 意見をたたかわせたい問題提起など,ムるってご投稿くだきい, (原稿は,刷り上がり,半ページ から 3 ベ ジに納まるようにお書き〈ださい.簡単に F 加筆訂正をお願いする場合があります)

一対比較データの解析にも役立つ定差方程式

牧野都治

-一巡 3 角形の数の平均と分散

一対比較データの解析法の 1 っとして,選好多角形 内に生じている一巡 3 角形の数を調べ, 帰無仮設Ho; r判定がランダムである」 を検定する方法があることは,ケンドールの本 (Rank

C

o

r

r

e

l

a

t

i

o

n

Methods) などによって,よく知られてい る.ここで,選好多角形は各頂点を矢線で結んだ多角 形であるが,たとえば 2 つの頂点A, B について, A→ B のように, B に矢が向けられていれば, A より Bの方が 好まれていることを表わす.また, A と B とに関して, 判定がランダムであるということは,矢がどちらに向 くかが確率 1/2 で定まることを意味している. (これら のことについては,本誌'93年 7 月号の拙著「一対比較 データの図表示に思う j を参照されたい) きて,このような検定を行なうために,選好 n角形内 に生ずる一巡 3 角形の数民の確率分布がわかってい ると都合がよい.少なくとも,れの平均・分散ぐらい は知っておきたい.これらを,ケンドールの本では, 次のように求めている.まず平均についてであるが, 頂点が U , k である 1 つの 3 角形をとり出したとき,そ れが一巡 3 角形であれば Pijk=l , そうでなければ Pijk =0 となる確率変数P

iik

を考えて,

Tn=

~

P

i

j

k

i..i,k とおけば,これが一巡 3 角形の級数になる.ただし, ~は相異なるすべての U , k についての和をとること を意味している.ここで 1 つの 3 角形をとり出した とき,それが一巡 3 角形となる確率を考えてみると, 左まわりと右まわりがあるので

(tYX2=1

+1

2

/"~

X2=

4

まきの とじ東京理科大学経営学部 〒 346 久喜市下清久

2

2

4

(52) となる.よって

1

,^__

3 1

E(P

iik

)=lX 一一 +OX 一一=一一

4

' V "

4

4 ・

これより E(Tn) は

E( 九)=E( 玄九) =+.(~)

.

7

.

k

r

i

i

k

J

-4"¥

3

}

=in(n-1)(n-2) … H・H ・-… H・H ・..…(1)

2

4

と求まる.また,分散 V(Tn) については 2 次モーメ ント E( T;) を計算しておいて , V( お )=E( T;)-E2 (九)と求めればよい.そこでケンドールでは,次のよ うな計算を行なう.はじめにE[ (~Piik)2] を求める.そ のために(玄P抽 )2の展開式を考え,その期待値をとっ て整頓すると,

E[(ふPiik)2]=告げ).{( ~)十3}

となり,分散が

V(Tn)=in(n1)(n-2)-H ・H・-…H ・ H・..…(2)

3

2

と求まる. これらはいずれもよく知られていることであるが, 本稿では以下, (1)式や(2)式が差分を用いて容易に求め られること,および一巡 3 角形でなく一巡 4 角形のと きにはどうなるか,などについての考えを述べたい.

-定差方程式による求め方

選好n角形内の一巡 3 角形の数 Tn について,次の式 が成り立つ.

Tn= Tn-

1

+

Rn

(n 二~4) …....・ H ・...・ H ・-…・・ (3). ただし , Rn は選好(n-1)角形内の 2 つの頂点と第 n 番目の点とを結んだときできる一巡 3 角形の数を表わ す.そこで, (3)式の両辺の期待値をとると

E(

Tn}= E(

Tn-l)+ E(Rn).

ところで, (n- l)角形でたとえば頂点 1 , 2 を結ぶ辺 について,図 1 のように矢が 2 へ向いていれば n→ 1 © 日本オペレーションズ・リサーチ学会. 無断複写・複製・転載を禁ず.

(2)

しようとするには,判定がランダムであるとしたとき, 選好n角形内に生ずる一巡k角形の数れ(k) の平均・分 散を計算しておくのがよい.そこで,一巡 3 角形のと きと同様に,定差方程式をつくってみると,

Tn= T

n

-

1

+Rn

(n 注 h 注 4) ……...・ H ・...・ H ・ ..(6) となる. (6)式でのおや Rn は ,

Tn(k)

,

Rn(k) と書カ通れる べきものであるが,以下その意味での h を省略するこ とにする.きて,平均E(Tn) は E(Tn)=E( れ -l)+E(R

n

) ....・ H ・...・ H ・....・ H ・-…(7) から求めることができる.ただし,一巡h角形について n ー I は E(R

n

) を計算するのに,ループを数えるか数えない 図 n 番目の項点を結ぶ 3 角形

→ 2 を頂点とする 3 角形が一巡する確率は(+r 矢

が 2 から 1 へ向いていても n→ 2 → 1 を頂点とする 3 角形が一巡する(前とは逆順になるが)確率は

(すr

nを頂点とする,このような 3 角形が(γ)個でき

るので,

1

/n-

1

¥

E(R

n) =ー(

"

"

J. )

4 ¥ 2

J

よって,

1

fn-l¥

E( Tn

)

=

E( T

n

-

l

)

+

-

4

-

.

(

"

-

:

;

J. ) ..・ H ・...・ H ・.(4).

4 ¥ 2

J

ここで, E( 包)=十であることを考慮して(4)式を解

くと,結局(1)式と同じものが得られる. それでは分散はというと, (3)式の右辺で, Tn-1と Rn とは独立であると考えてよいので,

V(T

n

)= V(Tn-l)+ V(Rn)

(

n

~4). いまの場合, V(Rn) は

V(Rn)=(

n

-

:

;

1

)

X

l

(

1

2

xよ+02X+)

¥ 2

J" l\~"4 ,

v"4J

一(lX士+oxty}=会 (γ)

であるから

3 f

n-l

¥

V( Tn

)

=

V( T

n

-

l

)

+ ー(

"

-

:

;

J. )... ...・ H ・-…・・(5).

1

6

¥ 2

J

'こで, V(73)=iであることを考慮、して(5)式を解

1

6

くと, (2)式と同じ V(T

n

) の式が得られる.

-一巡h角形の数の平均と分散

一巡 3 角形に限らず,一般に一巡k角形 (k 孟 4) を活用 1994 年 4 月号 かを定めておかなければならない.この,ループを数 えるかどうかについては,前記「一対比較データの図 表示に思う J に詳述しであるが, r ループを数えない場 合j には

E(Rn)=(ザ X2.(~二 D

よって (7)式は

E( れ)=E(

T

n

-

l

)

+

(

-

*

-

-

r

1.(

~- ~)

¥2J ¥k-1J

となり,これを解いて

E民(仇川

η

刊)=ベ(す封r-

ト刊吋-→

1

が得られる. 一方, r ループを数える場合」には (8)式が,

E(Rn)=(

+

r

X{(k一俳句二 D

となるので, (7)式は

E( れ)=E( 九ーペ+r

X{(k- l) !}.(~-n

¥k-1J

となり,これを解いて E( お)の式が

E( 九)=(すrX{(k 俳句)

と求まる. それでは,分散 V(Tn) はどうかというと, (6)式右辺 Tn-1と Rn とが独立でないので, V( お)= V( お→)十 V(R

n

) とすることができない.そこで,この方法での解はあ きらめなければならない.この場合にはやはり,ケン ドールの本にあるような方法を拡張して,こつこつ解 くよりしかたないのかもしれない. (53)

2

2

5

© 日本オペレーションズ・リサーチ学会. 無断複写・複製・転載を禁ず.

参照

関連したドキュメント

[r]

この節では mKdV 方程式を興味の中心に据えて,mKdV 方程式によって統制されるような平面曲線の連 続朗変形,半離散 mKdV

Existence of weak solution for volume preserving mean curvature flow via phase field method. 13:55〜14:40 Norbert

 当社は取締役会において、取締役の個人別の報酬等の内容にかかる決定方針を決めておりま

In this paper, we consider the discrete deformation of the discrete space curves with constant torsion described by the discrete mKdV or the discrete sine‐Gordon equations, and

方式で 45 ~ 55 %、積上げ方式で 35 ~ 45% 又は純費用方式で 35 ~ 45 %)の選択制 (※一部例外を除く)

先に述べたように、このような実体の概念の 捉え方、および物体の持つ第一次性質、第二次

電子式の検知機を用い て、配管等から漏れるフ ロンを検知する方法。検 知機の精度によるが、他