knitr::opts_chunk$set( collapse = TRUE, comment = "#>" )
This is the vignette for
ActFrag. This package extracts commonly used fragmention features from minute level actigraphy data. Recent studies haven shown that, on top of total daily active/sedentary volumes, the time accumulation strategies provide more sensitive information. This package provides functions to extract commonly used fragmentation metrics to quantidy such time accumulation strategies based on minute level actigraphy-measured activity counts data.
To download the package from Github
And to load in the package into the R environment
The expected data should consider at least one data frame of minute level activity counts, stored in a format of
data.frame of dimension $(\sum_i d_i) \times 1442$, where $d_i$ is the number of available days for subject i.And the order of the 1442 columns (and corresponding column names) should be "ID","Day","MIN1",..."MIN1440".
It is preferrable that user can also provide a
data.frame of wear/non-wear flag as same dimension of the activity counts. This flag data can serve as the following purposes:
Thbe wear/nonwear flag data should only consist of entries 0/1 representing nonwear/wear, sleep/wake, regions of interest/regions of interests, respectively. This is especially crucial for calculating features like total sedentary time, fragmentation etc. Because, we are not supposed to mix sedentary with sleep.
If no wear/nonwear flag data is supplied, users can create one using the
wear_flag function providing the time region:
data(example_activity_data) count = example_activity_data$count weartime = wear_flag(count.data = count, start = "06:00", end = "23:00")
In this version, we only incorporate type 3 . For 1 and 2 there are more appropriate softwares to look for.
Frequently, studies extract total time spent in sedentary behavior (e.g., total sedentary minutes per day) or proportion of waking hours spent sedentary. More advanced techniques have examined the effect of replacing sedentary time with active time spent either in light or moderate-vigorous intensity. For example, isotemporal substitution modeling) and compositional data analysis examine the combined effects of time spent sedentary, light and moderate-vigorous activity, and sleep while taking into account the codependence due to finite time during a day. Yet, most studies typically only use total volume of sedentary time while ignoring the potential importance of accumulation patterns. Recent studies have suggested that such patterns (known as prolonged sedentary bouts) may provide additional sensitivity to predict health outcomes and provide additional information beyond total volume of sedentary time.
Fragmentation metrics study the accumulation pattern of TST and TAT by quantifying the alternating these sequences via summaries of duration of and fre- quency of switching between sedentary and active bouts. They provide unique translatable insights into accumulation patterns for sedentary and active time and lead to additional findings of associations between those patterns and health outcomes on top of total sedentary volume Here is the list of available fragmentation metrics:
We can calculate the above mentioned metrics for both sedentary and active bout. Details about fragmentations can be found here.
fragmentation_long calcualte fragmentation features, (for a single vector amd whole dataset respectively). The argument
metrics, which consists of "mean_bout", "TP", "Gini", "hazard", "power", and "all" decides which metrics to calcualte. "all" will lead to all metrics.
For a single day of count (a vector):
data(example_activity_data) count1 = c(t(example_activity_data$count[1,-c(1,2)])) wear1 = c(t(example_activity_data$wear[1,-c(1,2)])) mb = fragmentation(x = count1, w = wear1, thresh = 100, metrics = "mean_bout",bout.length = 1) tp = fragmentation(x = count1, w = wear1, thresh = 100, metrics = "TP",bout.length = 1)
Given all the activity and wear/nonwear flag data for the whole dataset, user can choose to calcualte framentation at daily level, or aggregate bouts across all available days by choosing from either "subject" and "day" for the argument
data(example_activity_data) count = example_activity_data$count wear = example_activity_data$wear frag_by_subject = fragmentation_long(count.data = count, weartime = wear,thresh = 100, metrics = "all",by = "subject",bout.length = 1) frag_by_day = fragmentation_long(count.data = count, weartime = wear,thresh = 100, metrics = "all",by = "day",bout.length = 1)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.