def testRootMeanSquareError():
myFunc = atsFunc()
valList1 = []
valList2 = []
valList3 = []
for val in range(0, 10):
valList1.append(val)
valList2.append(val + 1)
valList3.append(val + 2)
print valList1
print myFunc.rootMeanSquareError(valList2, valList1, [1, 1, 0])
print valList2
print myFunc.rootMeanSquareError(valList2, valList2, [1, 1, 0])
print valList3
print myFunc.rootMeanSquareError(valList2, valList3, [1, 1, 0])
def testCorrcoef():
myFunc = atsFunc()
for testCount in range(0, 10):
valList1 = []
valList2 = []
random.seed(testCount)
for idx in range(0, 10):
valList1.append(random.random())
valList2.append(random.random())
print valList1
print valList2
corrcoefVal = myFunc.corrcoef(valList1, valList2)
print u"相関係数の値(numpy) \t: ", numpy.corrcoef(valList1, valList2)[0][1]
print u"相関係数の値(atsFunc)\t: ", corrcoefVal
valList1 = [10, 11, 12, 13, 14, 15, 14, 13, 12, 11, 10, 11]
valList2 = [10, 11, 12, 13, 14, 15, 14, 13, 12, 11, 10, 11]
corrcoefVal = myFunc.corrcoef(valList1, valList2)
print numpy.corrcoef(valList1, valList2)
print u"相関係数の値\t: ", corrcoefVal
def testGesdOutlierTest():
valueList = [
-0.25,
0.68,
0.94,
1.15,
1.20,
1.26,
1.26,
1.34,
1.38,
1.43,
1.49,
1.49,
1.55,
1.56,
1.58,
1.65,
1.69,
1.70,
1.76,
1.77,
1.81,
1.91,
1.94,
1.96,
1.99,
2.06,
2.09,
2.10,
2.14,
2.15,
2.23,
2.24,
2.26,
2.35,
2.37,
2.40,
2.47,
2.54,
2.62,
2.64,
2.90,
2.92,
2.92,
2.93,
3.21,
3.26,
3.30,
3.59,
3.68,
4.30,
4.64,
5.34,
5.42,
6.01,
]
correctRList = [0, 3.118, 2.942, 3.179, 2.810, 2.815, 2.848, 2.279, 2.310, 2.101, 2.067]
correctThresholdList = [0, 3.158, 3.151, 3.143, 3.136, 3.128, 3.120, 3.111, 3.103, 3.094, 3.085]
print "valueList ", valueList
print "len(valueList) ", len(valueList)
print "R should have these values : ", correctRList[1:]
print "critical values should have these values : ", correctThresholdList[1:]
testValueList = copy.deepcopy(valueList)
maxOutlierCount = 10
print u"-----^"
myFunc = atsFunc()
gesdResultList = myFunc.gesdOutlierTest(valueList, maxOutlierCount=maxOutlierCount, debugFlag=True)
print u"外れ値とする値の個数の最大 :\t", maxOutlierCount
print u"外れ値ではないと判断された値のリスト:\t", gesdResultList[0]
print u"外れ値ではないと判断された値の個数 :\t", len(gesdResultList[0])
print u"外れ値となった値のリスト :\t", gesdResultList[1]
print u"外れ値となった値の個数 :\t", len(gesdResultList[1])
