Editing Finding relations between variables in time series

{{Topic Infobox}}
Most personal science projects require finding relationships between different variables of the type 'time series'<ref>Core-Guide_Longitudinal-Data-Analysis_10-05-17.pdf (duke.edu)</ref>. An example could be the question "does my daily chocolate consumption correlate with my daily focus score?".   

You could do experiments if you control everything rigidly or if the effects are strong and quick, like less than a week. Old data may be useable as Baseline and a baseline may rule out some issues. If both block (like 2 weeks) and daily mixed (random intervention every day)  produce the same results then issues of time series are probably not in your experiment.    

Finding more complicated relationships require better statistical tests and algorithms and data science skills. Apps that would do this automatically or at least easily are not yet available. See below. Most internet resources treat time series as (regular cyclical) series, which is not useful as most of the tracked variables have irregular patterns and don't even have a regularly cyclical component.  

To do anything mentioned above you need to have your data parsed, cleaned, all in one place, and ideally even visualized.   

== List of less technical tools ==
There are also a number of tools or apps that can semi-automatically perform these correlations and help in doing the analyses.

==== [[Open Humans]] and their Personal Analysis notebooks ====
Open Humans provides a library of notebooks that can be used to visualize data across data sources and find relations between different variables. It also supports the upload of generic data files through the [https://blog.openhumans.org/2021/05/17/file-upload-data-types/ File Uploader].

==== [[Zenobase]] ====
[https://blog.zenobase.com/post/81497604762 Zenobase] can test correlations based on user-specified questions. User must configure lag, regression method and aggregation method using a UI. Powerful filtering tools too.  

==== Curedao ====
[https://github.com/curedao/decentralized-fda Curedao]. Correlation over bins and lags selecting the biggest effect.

==== Data Flexor ====
[https://www.dataflexor.com/ DataFlexor] Lots of pretty pictures. Not super advanced statistics yet.

==== [[HALE.life]] ====
"Baysian nodes, 'do' semantics, AI and experts"<ref>https://www.hale.life/</ref> Right now only for sports teams.  

==== [[Exist.io]] <ref>https://github.com/ejain/n-of-1-ml</ref> ====
From the [https://exist.io/ Exist.io] main site :"Which habits go together? Correlations are the most powerful part of Exist. By combining your data, we can answer questions like: “What makes me happiest?”, “What can I do to be more active?”, “When am I most productive?”"

==== [[Habitdash]] ====

[https://habitdash.com/  Habitdash]'s ''Automatic data analysis searches for hidden patterns'' to find relationships between activity, sleep, weight and other habits.

==== [[Optimized app]] ====

[http://optimized-app.com/ Optimized] claims to do ''"automatic correlation mining"''

==== [[Lytiko]] ====
lytiko.com promises correlations connections deep insights and visualizations. 
==== [[Vital]] ====
Vital (https://tryvital.io). API for health and fitness data. Free to use API for collecting wearables and health data and standardising them into one API. You can also use Vital's API for delivering at-home test kits.  

==== young.ai and [http://www.aging.ai/ aging.ai] ====
Deep learning predictor of age based on human blood tests, young.ai also makes recommendations.

==== Sonar [https://www.sonarhealth.co sonarhealth.co] ====
Customizable aggregation and syncing like weigh fitbit twice as much as apple watch or average steps instead of sum.

====== tunum.health ======
pearson correlation, trend analysis and manual dichotomization

[[Gyroscope]]

[[Welltory]]

[[Bearable App]]

[https://www.realize.me/#sign-up Realize Me]

[https://headsuphealth.com/ Heads Up]

[[Inside Tracker]]

Wellness FX

Export from Apple Health<ref>github.com/Lybron/health-auto-export</ref> (no analysis)

ConnectorDB DIY OS no analysis

Heedy DIY OS no analysis

Zapier, Integromat, IFTTT, DIY no analysis

== List of very technical tools ==

Some people do all the data science by themselves, by using programming languages such as R and Python in notebooks or apps. Coding platforms such as the notebooks on [[Open Humans]], Kaggle or GitHub can help. 

Programming languages for statistics; Matlab, R, Python, Julia. 

Try [https://orange3-timeseries.readthedocs.io/en/latest/ Python GUI time series analysis] . 

==== [[List of Interesting Self-Tracking Results#Observational.2C%20Many%20variables|DIY Individuals]] ====
== Reasons time series analysis especially as applied to QS is hard ==
[https://forum.quantifiedself.com/t/my-baseline-network-physiology-10-days-of-eeg-egg-ekg-cgm-temperature-activity-and-food-logs/5671/19 Wavelet coherence] is one potential solution.

Really strong relationships will be detected even through most of these problems. 

[http://www.tylervigen.com/spurious-correlations Spurious Correlations] mostly shows that if two things are trending in one direction and are checked for correlation they will show a very significant correlation. Practice effect is a subset. Another is one instance of an event type A increases the chances of the same event type happening soon after. Economists suggest unit root.

Effects on target variable from outside known variables. In non time series this is compensated for with RCT but in time series such an effect may last a while and coincide with an intervention causing very false results. This problem makes baseline data gathering more difficult and also necessary. Sometimes a baseline will show that this issue does not occur for a particular target variable. Alternatively experimenter could compensate by strictly controlling all possible sources of variance.  

Lag. What if eating pizza on one day causes heartburn the next? 

Build up. What if it takes two days of eating pizza to cause heartburn? 

Rate of change. Trend. Opposite of build up; derivative instead of integral. Stopping or starting an all pizza diet causes heartburn. 

Bin. Window. Smooth. Variables only make domain sense as aggregate over some time. Variables have a really high sampling rate.  

Interpolate. Variables have different sampling rates so need to be interpolated to be compared.  

Types of data. [Exercised] is an event with specific occurrence moment and length while [tired] is a vaguer value user could use to try to describe feelings past 4 hours. 

All the [[Issues with Self Report]] .      

Few positive instances but they are important. Went to a specific restaurant twice got sick soon after twice. Only ever got sick with similar symptoms five times. Or. Two large rare humps happen almost one after the other, similar to previous example if treated as events, adding the fact that lots of samples showing their similarity in shape too.   

Since removing real effects of other variables on target variable makes the variable of interest's effect stand out, 'machine learning' needs to be used. Basic approach would be to bin predictor variables multiple ways based on time from effect being checked, mean or other aggregator method and window of the aggregator. 

Machine learning also has limits on the kind of patterns it can detect. 

=== What to expect from the complete analysis tool ===
User without experience in statistical analysis will not be able to tell the difference between correctly computed correlations and poorly computed ones. However, a genuinely complete analysis produces plots which should include at least some of the following:  

Interpolation for irregular time series.  

Change point or breakpoint detection.  

Outlier detection. Smoothing.  

Removal of effect of variables found to correlate with this one to show residuals. 

Cycles decomposition using a model like ARIMA. Ex. kayak season is in the summer or lunch is at exactly 1pm.  

Detection of repeated shapes implying similar events that are not cyclical; like dinner is anywhere between 4pm and 10pm and causes a particular 2 hour spike in glucose.

== References ==
[[Category:Data analysis]]