added example

Author: muschellij2

R/create_day_inclusion.R

@@ -94,6 +94,9 @@ add_day_inclusion = function(
     data,
     ...
 ) {
+  time = date = NULL
+  rm(list = c("time", "date"))
+
   day_data =   create_day_inclusion(data, ...)
   data = acti_standardize_data(data, subset_xyz = FALSE, check_xyz = FALSE)
   data = data %>%

example.R

@@ -5,13 +5,7 @@ sc = import("stepcount")
 sc$stepcount
 stepcount::stepcount_check()
 
-library(tidyverse)
-library(activerse)
-path = "~/Dropbox/Data/wearable_pilot/data/raw/ag_axivity_test/John/114890_0000000000.cwa"
-
-# path = "inst/extdata/TAS1H30182785_2019-09-17.gt3x.gz"
-data = acti_read_cwa(path)
-function(data, calibrate = FALSE) {
+acti_run = function(data, calibrate = FALSE) {
   if (calibrate) {
     data = acti_calibrate(data)
   }
@@ -25,19 +19,42 @@ function(data, calibrate = FALSE) {
   processed = processed %>%
     dplyr::full_join(steps)
 
-  res = create_day_inclusion(processed, min_required = 1368L)
+  inclusion = create_day_inclusion(processed, min_required = 1368L)
   processed = processed %>%
     mutate(date = as.Date(time)) %>%
-    left_join(res)
-  processed %>%
-    mutate(log10AC = log(counts)) %>%
+    left_join(inclusion)
+  na_sum = function(x) sum(x, na.rm = TRUE)
 
-    summarise
+  processed = processed %>%
+    mutate(log10AC = log(counts + 1))
+
+  daily = processed %>%
+    acti_create_date() %>%
+    group_by(date) %>%
+    summarise(
+      dplyr::across(any_of(
+        c("counts", "log10AC", "steps", "MIMS", "mims")),
+        na_sum
+      )
+    )
+  daily = daily %>%
+    full_join(inclusion)
+
+  list(
+    minute_level = processed,
+    inclusion = inclusion,
+    daily = daily
+  )
 }
 
+library(tidyverse)
+library(activerse)
+path = "~/Dropbox/Data/wearable_pilot/data/raw/ag_axivity_test/John/114890_0000000000.cwa"
 
+# path = "inst/extdata/TAS1H30182785_2019-09-17.gt3x.gz"
+data = acti_read_cwa(path)
 
-
-
-
-
+out = acti_run(data)
+out$minute_level
+out$inclusion
+out$daily