mirror of
https://github.com/wassname/catalyst.git
synced 2026-08-13 12:00:16 +08:00
Eddie Hebert
7904773d00
The latest flake8 release in now 1.5, which pulls in pep8: 1.3.4a0 The upgrade pep8 has changes to what it picks up as lint. Making code base compatible, so that new devs can install pep8 from PyPI and not have friction over the version difference. Currently using these ignores in the config file: ``` [pep8] ignore = E124,E125,E126 ``` Ignoring these since they are difficult to squash while maintaining an 80 char line length, and appear spurious. Should address later. Updates Travis config, README, and pip requirements to reflect change.
116 lines
3.6 KiB
Python
116 lines
3.6 KiB
Python
#
|
|
# Copyright 2012 Quantopian, Inc.
|
|
#
|
|
# Licensed under the Apache License, Version 2.0 (the "License");
|
|
# you may not use this file except in compliance with the License.
|
|
# You may obtain a copy of the License at
|
|
#
|
|
# http://www.apache.org/licenses/LICENSE-2.0
|
|
#
|
|
# Unless required by applicable law or agreed to in writing, software
|
|
# distributed under the License is distributed on an "AS IS" BASIS,
|
|
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
|
|
# See the License for the specific language governing permissions and
|
|
# limitations under the License.
|
|
|
|
from numbers import Number
|
|
from collections import defaultdict
|
|
from math import sqrt
|
|
|
|
from zipline.transforms.utils import EventWindow, TransformMeta
|
|
|
|
|
|
class MovingStandardDev(object):
|
|
"""
|
|
Class that maintains a dicitonary from sids to
|
|
MovingStandardDevWindows. For each sid, we maintain a the
|
|
standard deviation of all events falling within the specified
|
|
window.
|
|
"""
|
|
__metaclass__ = TransformMeta
|
|
|
|
def __init__(self, market_aware, days=None, delta=None):
|
|
|
|
self.market_aware = market_aware
|
|
|
|
self.delta = delta
|
|
self.days = days
|
|
|
|
# Market-aware mode only works with full-day windows.
|
|
if self.market_aware:
|
|
assert self.days and not self.delta,\
|
|
"Market-aware mode only works with full-day windows."
|
|
|
|
# Non-market-aware mode requires a timedelta.
|
|
else:
|
|
assert self.delta and not self.days, \
|
|
"Non-market-aware mode requires a timedelta."
|
|
|
|
# No way to pass arguments to the defaultdict factory, so we
|
|
# need to define a method to generate the correct EventWindows.
|
|
self.sid_windows = defaultdict(self.create_window)
|
|
|
|
def create_window(self):
|
|
"""
|
|
Factory method for self.sid_windows.
|
|
"""
|
|
return MovingStandardDevWindow(
|
|
self.market_aware,
|
|
self.days,
|
|
self.delta
|
|
)
|
|
|
|
def update(self, event):
|
|
"""
|
|
Update the event window for this event's sid. Return an ndict
|
|
from tracked fields to moving averages.
|
|
"""
|
|
# This will create a new EventWindow if this is the first
|
|
# message for this sid.
|
|
window = self.sid_windows[event.sid]
|
|
window.update(event)
|
|
return window.get_stddev()
|
|
|
|
|
|
class MovingStandardDevWindow(EventWindow):
|
|
"""
|
|
Iteratively calculates standard deviation for a particular sid
|
|
over a given time window. The expected functionality of this
|
|
class is to be instantiated inside a MovingStandardDev.
|
|
"""
|
|
|
|
def __init__(self, market_aware, days, delta):
|
|
# Call the superclass constructor to set up base EventWindow
|
|
# infrastructure.
|
|
EventWindow.__init__(self, market_aware, days, delta)
|
|
|
|
self.sum = 0.0
|
|
self.sum_sqr = 0.0
|
|
|
|
def handle_add(self, event):
|
|
assert 'price' in event
|
|
assert isinstance(event.price, Number)
|
|
|
|
self.sum += event.price
|
|
self.sum_sqr += event.price ** 2
|
|
|
|
def handle_remove(self, event):
|
|
assert 'price' in event
|
|
assert isinstance(event.price, Number)
|
|
|
|
self.sum -= event.price
|
|
self.sum_sqr -= event.price ** 2
|
|
|
|
def get_stddev(self):
|
|
# Sample standard deviation is undefined for a single event or
|
|
# no events.
|
|
if len(self) <= 1:
|
|
return None
|
|
|
|
else:
|
|
average = self.sum / len(self)
|
|
s_squared = (self.sum_sqr - self.sum * average) \
|
|
/ (len(self) - 1)
|
|
stddev = sqrt(s_squared)
|
|
return stddev
|