def performEventDurationTest():
mahal_timeseries = readOutlierScores("results/outlier_scores.csv")
global_pace_timeseries = readGlobalPace("4year_features")
zscore_timeseries = readZScoresTimeseries("results/zscore.csv")
mahal_timeseries_fine = readOutlierScores("results/link_20_normalize_outlier_scores.csv")
threshold_vals = [.90,.91,.92,.93,.94,.95,.96,.97,.98,.99]
window_sizes = [1,2,3,4,6,8,12,24]
with open('results/threshold_experiment.csv', 'w') as f:
w = csv.writer(f)
w.writerow(['granularity', 'window','threshold', 'duration'])
for window_size in window_sizes:
for threshold in threshold_vals:
print (window_size, threshold)
events = detectWindowedEvents(mahal_timeseries, zscore_timeseries, global_pace_timeseries,
"results/events_windowed.csv", window_size=window_size, threshold_quant=threshold)
duration = getEventDuration(events, "2012-10-31")
w.writerow(["coarse", window_size, threshold, duration])
events= detectWindowedEvents(mahal_timeseries_fine, zscore_timeseries, global_pace_timeseries,
"results/link_20_normalize_events_windowed.csv", window_size=window_size,
threshold_quant=threshold)
duration = getEventDuration(events, "2012-10-31")
w.writerow(["fine", window_size, threshold, duration])
def performEventDurationTest():
mahal_timeseries = readOutlierScores("results/outlier_scores.csv")
global_pace_timeseries = readGlobalPace("4year_features")
zscore_timeseries = readZScoresTimeseries("results/zscore.csv")
mahal_timeseries_fine = readOutlierScores(
"results/link_20_normalize_outlier_scores.csv")
threshold_vals = [.90, .91, .92, .93, .94, .95, .96, .97, .98, .99]
window_sizes = [1, 2, 3, 4, 6, 8, 12, 24]
with open('results/threshold_experiment.csv', 'w') as f:
w = csv.writer(f)
w.writerow(['granularity', 'window', 'threshold', 'duration'])
for window_size in window_sizes:
for threshold in threshold_vals:
print(window_size, threshold)
events = detectWindowedEvents(mahal_timeseries,
zscore_timeseries,
global_pace_timeseries,
"results/events_windowed.csv",
window_size=window_size,
threshold_quant=threshold)
duration = getEventDuration(events, "2012-10-31")
w.writerow(["coarse", window_size, threshold, duration])
events = detectWindowedEvents(
mahal_timeseries_fine,
zscore_timeseries,
global_pace_timeseries,
"results/link_20_normalize_events_windowed.csv",
window_size=window_size,
threshold_quant=threshold)
duration = getEventDuration(events, "2012-10-31")
w.writerow(["fine", window_size, threshold, duration])
events= detectWindowedEvents(mahal_timeseries_fine, zscore_timeseries, global_pace_timeseries,
"results/link_20_normalize_events_windowed.csv", window_size=window_size,
threshold_quant=threshold)
duration = getEventDuration(events, "2012-10-31")
w.writerow(["fine", window_size, threshold, duration])
if(__name__=="__main__"):
#performEventDurationTest()
mahal_timeseries = readOutlierScores("results/outlier_scores.csv")
global_pace_timeseries = readGlobalPace("4year_features")
zscore_timeseries = readZScoresTimeseries("results/zscore.csv")
detectWindowedEvents(mahal_timeseries, zscore_timeseries, global_pace_timeseries,
"results/events_windowed.csv", window_size=8, threshold_quant=.95)
mahal_timeseries = readOutlierScores("results/link_20_normalize_outlier_scores.csv")
global_pace_timeseries = readGlobalPace("4year_features")
zscore_timeseries = readZScoresTimeseries("results/zscore.csv")
detectWindowedEvents(mahal_timeseries, zscore_timeseries, global_pace_timeseries,
"results/link_20_normalize_events_windowed.csv", window_size=8,
threshold_quant=.95)
logMsg("done")
events = detectWindowedEvents(
mahal_timeseries_fine,
zscore_timeseries,
global_pace_timeseries,
"results/link_20_normalize_events_windowed.csv",
window_size=window_size,
threshold_quant=threshold)
duration = getEventDuration(events, "2012-10-31")
w.writerow(["fine", window_size, threshold, duration])
if (__name__ == "__main__"):
#performEventDurationTest()
mahal_timeseries = readOutlierScores("results/outlier_scores.csv")
global_pace_timeseries = readGlobalPace("4year_features")
zscore_timeseries = readZScoresTimeseries("results/zscore.csv")
detectWindowedEvents(mahal_timeseries,
zscore_timeseries,
global_pace_timeseries,
"results/events_windowed.csv",
window_size=8,
threshold_quant=.95)
mahal_timeseries = readOutlierScores(
"results/link_20_normalize_outlier_scores.csv")
global_pace_timeseries = readGlobalPace("4year_features")
zscore_timeseries = readZScoresTimeseries("results/zscore.csv")
detectWindowedEvents(mahal_timeseries,
zscore_timeseries,