1

ある状態とは大きく違うデータが与えられた時に
判別するプログラムを書こうとしています。
1つの判別方法として、A_B間の標準偏差とA_C間の標準偏差を求めて、
Aの状態と似ている/異なる判断ができるのではないかと考えました。

しかしながら、現在のコードでは想定通りにはならず、
判別を行う閾値の設定方法もどうしたら良いかわからない状態です。

想定しているデータ

#Aは標準のデータとする
A =[9.43,5.04,8.86,2.48,5.27,5.21,7.95,1.96,1.53,3.13]
#Aの状態とほとんど同じ
B = [8.99,5.00,8,2.30,5.20,4.89,8.14,2.00,1.20,2.86]
#Aの状態とは大きく違う
C = [5.20,3.14,7.71,4.33,3.47,3.10,3.88,4.57,2.48,6.41]

実行コード

import numpy as np
#Aは標準のデータとする
A =[9.43,5.04,8.86,2.48,5.27,5.21,7.95,1.96,1.53,3.13]
#Aの状態とほとんど同じ
B = [8.99,5.00,8,12,2.30,5.20,4.89,8.14,2.00,1.20,2.86]
#Aの状態とは大きく違う
C = [5.20,3.14,7.71,4.33,3.47,3.10,3.88,4.57,2.48,6.41]

marginAB = 0
marginAC = 0
#A_B間の標準偏差とA_C間の標準偏差を求めて、Aの状態と似ている/異なるの判断をしようとしています
for i in range(len(A)):
        marginAB += (A[i] - B[i])**2
        marginAC += (A[i] - C[i])**2

#sd:標準偏差(standard deviation)
sdAB = marginAB / len(A) *2
sdAC = marginAC / len(A)*2

#sdACの方が同じ配列番号間の差は大きくなるからsdABよりも標準偏差が大きいと考えられる
#似ているか/異なるかの判断
#sdABにかけている値はどのように決定すべきかわからないため適当に想定
if sdAC > 2*sdAB:
    print('異なる')
else:
    print('似ている')

出力

似ている→異なる(修正後)

確認したこと
判別がうまくいかなかったことからsdABとsdACを確認しましたが、予想に反してsdACの値の方がsdABよりも小さかったことも理解できませんでした。(修正後解決しました)

sdAB,sdACの値(今回の例の場合)

0.2587999999999998
13.795580000000001
2
  • Bの要素の個数が他より一つ多いのですが、データの間違いではないですか。
    – Hideki
    Commented 2019年6月6日 13:36
  • ご指摘いただきましてありがとうございます。 解決しました。質問で「sdABにかけている値はどのように決定すべきかわからないため適当に想定」とあるのですが、今回は仮に2倍としましたが、この閾値はベースラインとなるsdABと比較対象のsdACを実際に求めてみて決めるべきでしょうか。
    – user32568
    Commented 2019年6月7日 0:15

1 件の回答 1

0

Bの4番目に12を入れているので、それ以降はABの値が全く違うものになっています。

A =[    B = [
9.43    8.99
5.04    5
8.86    8
2.48    12
5.27    2.3
5.21    5.2
7.95    4.89
1.96    8.14
1.53    2
3.13]   1.2
        2.86]
3
  • 8.128,12とtypoしたってことですかね。 Commented 2019年6月6日 14:01
  • ご指摘いただきましてありがとうございます。 解決しました。質問で「sdABにかけている値はどのように決定すべきかわからないため適当に想定」とあるのですが、今回は仮に2倍としましたが、この閾値はベースラインとなるsdABと比較対象のsdACを実際に求めてみて決めるべきでしょうか。
    – user32568
    Commented 2019年6月6日 22:17
  • 「似ている」かかどうかの判断基準は人により異なるので、決めようがない(逆に言えば主観で決めればいい)と思いますが。 Commented 2019年6月7日 0:57

この質問に回答するには、ログインする必要があります。