This is a pre-production deployment of Warehouse. Changes made here affect the production instance of PyPI (
Help us improve Python packaging - Donate today!

A least-squares optimal method to find the saccadic reaction time and the saccade duration from tracked gaze points

Project Description

A least-squares optimal offline method to find saccadic reaction time and saccade duration from tracked gaze points.

You have tracked the gaze points of the following event sequence:

  1. A person looks at point (A). An image appears at (B).
  2. The person reacts to the image and starts to move their gaze toward (B). The saccade starts.
  3. The gaze arrives to (B). The saccade ends.
  4. The person looks at point (B).

Now you want to determine:

  1. The time between events 1 and 2, called the saccadic reaction time or SRT.
  2. The time between events 2 and 3, called the saccade duration or SD.

The saccademodel algorithm computes the times for you by fitting an ideal gaze path to the data. The ideal gaze path has the following structure. From t=0 to t=saccade_start the ideal gaze is exactly at point (A). From t=saccade_start to t=saccade_end the ideal gaze moves from (A) to (B) with constant velocity. From t=saccade_end to t=n the gaze remains at (B). The algorithm finds such times saccade_start and saccade_end that minimize the mean squared error between the ideal gaze path and the given tracked gaze points. In other words, the algorithm splits the data to three segments: source fixation, saccade, and target fixation.

As the greatest advantage, when compared to velocity-based saccade recognition methods, data does not need to be filtered beforehand because the squared error approach does that by itself. Even though filtering would yield smaller total squared error, it does not affect the estimates of saccade_start and saccade_end. However, if the noise in the data is nonstationary, some special noise filtering methods might be needed.

As the greatest disadvantage, the saccademodel algorithm is suitable only for offline analysis and therefore cannot be used in realtime setups.

1. Install

With pip:

$ pip install saccademodel

2. Usage

The data structure pointlist is used thoroughly. It is a list of points, where each point is a list [x, y].

The usage is simple:

>>> import saccademodel
>>> rawdata = [
    [130.012, 404.231],
    [129.234, 403.478],
    [None, None],
    [133.983, 450.044],
>>> results =
>>> print(results)
    'source_points': [[344.682, 200.115], ...],
    'saccade_points': [[324.233, 202.393], ...],
    'target_points': [[556.423, 112.607], ...],
    'mean_squared_error': 0.000166802

Note that the lengths of the returned lists can be used to determine saccadic reaction time and duration. For example, given the points from the appearance of stimulus, the saccadic reaction time is captured in the length of source_points and the saccade duration in saccade_points. If the frame rate is known, you can convert the lengths to seconds by:

>>> framerate = 300.0  # samples per second
>>> saccadic_reaction_time = len(results.source_points) / framerate
>>> saccade_duration = len(results.saccade_points) / framerate

3. API



  • gazepointlist: a list of [x, y] points i.e. a list of lists.

Return dict with following keys:

  • source_points: the points before the saccade
  • saccade_points: the points in the saccade
  • target_points: the points after the saccade.
  • mean_squared_error: the average squared error from the model for a point.

3.2. saccademodel.version

The current version string:

>>> saccademodel.version

4. For developers

4.1. Virtualenv

Use virtualenv:

$ virtualenv -p python3.5 saccademodel-py
$ cd saccademodel-py
$ source bin/activate
$ deactivate

4.2. Jupyter Notebook


$ cd explore
$ jupyter notebook

Install requirements:

$ pip install --editable .[notebook]

4.3. Testing

Follow instructions to install pyenv and then either run quick tests:

$ python3.5 test

or comprehensive tests for multiple Python versions in tox.ini:

$ pyenv local 2.6.9 2.7.10 3.1.5 3.2.6 3.3.6 3.4.3 3.5.0
$ eval "$(pyenv init -)"
$ pyenv rehash
$ tox

Install new pyenv environments by:

$ pyenv install 3.4.5

Validate README.rst at

4.4. Publishing to PyPI

Follow python packaging instructions:

  1. Create an unpacked sdist: $ python sdist
  2. Create a universal wheel: $ python bdist_wheel --universal
  3. Go to PyPI and register the project by filling the package form by uploading saccademodel.egg-info/PKG_INFO file.
  4. Upload the package with twine:
    1. Sign the dist: $ gpg --detach-sign -a dist/saccademodel-1.2.3*
    2. Upload: twine upload dist/saccademodel-1.2.3* (will ask your PyPI password)
  5. Package published!

Updating the package takes same steps except the 3rd.

6. License

MIT License

Release History

Release History

This version
History Node


Download Files

Download Files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

File Name & Checksum SHA256 Checksum Help Version File Type Upload Date
saccademodel-0.1.0-py2.py3-none-any.whl (14.5 kB) Copy SHA256 Checksum SHA256 py2.py3 Wheel Dec 9, 2015
saccademodel-0.1.0.tar.gz (10.4 kB) Copy SHA256 Checksum SHA256 Source Dec 9, 2015

Supported By

WebFaction WebFaction Technical Writing Elastic Elastic Search Pingdom Pingdom Monitoring Dyn Dyn DNS Sentry Sentry Error Logging CloudAMQP CloudAMQP RabbitMQ Heroku Heroku PaaS Kabu Creative Kabu Creative UX & Design Fastly Fastly CDN DigiCert DigiCert EV Certificate Rackspace Rackspace Cloud Servers DreamHost DreamHost Log Hosting