% D3Mdemo.m 
% provides examples of 3D dynamical density delay maps
% creating videos for 2 RR time series:
% a sleep RR and an atrial fibrillation RR

% Developed by A Burykin (burykin@gmail.com), L Citi, 
% MD Costa and AL Goldberger
% and modified by S Mariani (sara.mariani@wyss.harvard.edu)
% Wyss Institute at Harvard
% when using this script, please reference: Burykin, Anton et al.
% "Dynamical density delay maps: simple, new method for visualising 
% the behaviour of complex systems." BMC medical informatics and 
% decision making 14.1 (2014): 6.
clc
close all
clear all
display('Dynamical density delay maps: sleep RR')
load Sleep_RR.dat
RR=Sleep_RR;
T=500/60;      % min time interval to display in one frame
Tshift=10/60;  % min time shift between frames 
videoname='Sleep3D';
foldername='Sleep3Dframes';
step=0.00015;
D3M_3Dfun(RR,videoname, foldername, T, Tshift, step);

clear all
display('Dynamical density delay maps: typical AF')
load Typical_AF.dat
RR=Typical_AF;
T=20;      % min time interval to display in one frame
Tshift=5;  % min time shift between frames 
videoname='Typical_AF3D';
foldername='Typical_AF3Dframes';
step=0.00015;
D3M_3Dfun(RR,videoname, foldername, T, Tshift, step);