Variables contained in part5 person summary output

[Kelsey Sewell]

Dear Vincent and Jairo,

We are using GGIR Version 2.10-1. 

We are wondering if GGIR does not produce mean-level columns in the part5 person summary output (e.g., ACC_day_VIG_unbt_mg_pla) if there are not enough IDs with vigorous level accelerations across days? This variable and a few other ACC_day_VIG variables are not being generated in our part5 person summary output despite at least two people having multiple days (including weekday and weekend) worth of vigorous acceleration data in part5 day summary output. We have all the behavioral level variables available (e.g., dur_day_VIG_unbt_min_pla) even though there is no column called (ACC_day_VIG_unbt_mg_pla). 

If this is a bug that has been addressed in version updates please let us know. We are unable to update our version due to ongoing longitudinal data collection but it will be good to have this information for our records moving forward.

Thank you so much.

Kind Regards,

Kelsey & Audrey 


See below our scripts for reference. 

GGIR(
  mode=c(1:5),
#=====================
# General + Part 1
#=====================
datadir="/Volumes/Data1/FLAME/ACTIGRAPH/GT3X/PITT/PITT_Data_Dir_Pitt",
outputdir="/Volumes/Data1/FLAME/ACTIGRAPH/GT3X/PITT/PITT_Output_Dir_Pitt",
desiredtz = "America/New_York",
overwrite = FALSE,
print.filename = TRUE,
do.parallel = TRUE,
#=====================
# Part 2
#=====================
idloc = 2,
strategy = 1,
maxdur = 0,
#maxdur = 7
#excludefirstlast = TRUE,
ilevels = seq(0,600,by=25),
iglevels = c(seq(0,4000,by=25),8000),
qlevels = c(960/1440, 1320/1440, 1380/1440, 1410/1440, 1430/1440, 1435/1440, 1438/1440),
#mvpathreshold =c(60, 100, 100.6),
mvpathreshold =c(100),
printsummary = TRUE,
do.part2.pdf = TRUE,
epochvalues2csv = TRUE,
winhr = c(5,10),
cosinor = TRUE,
#=====================
# Part 3+4
#=====================
do.part3.pdf = TRUE,
outliers.only = FALSE ,
#=====================
# Part 5
#=====================
#threshold.lig = c(18, 35, 44.8), threshold.mod = c(60, 100, 100.6), threshold.vig = c(430, 428.8),
threshold.lig = c(35), threshold.mod = c(100), threshold.vig = c(430),
boutcriter = 0.8, boutcriter.in = 0.9, boutcriter.lig = 0.8,
boutcriter.mvpa = 0.8, boutdur.in = c(10,20,30,60), boutdur.lig = c(1,5,10),
boutdur.mvpa = c(5,10),
includedaycrit.part5 = 2/3,
frag.metrics="all",
part5_agg2_60seconds = TRUE,
week_weekend_aggregate.part5 = TRUE,
#=====================
# Visual Report
#=====================
timewindow = c("MM", "WW"),
do.report=c(2,4,5),
visualreport = TRUE
)
quit(save = "no", status = 0, runLast = TRUE)

[Vincent van Hees]

Hi Kelsey and Audrey,

The average acceleration cannot be calculated if no time is spent in a behavioural class. The missing value is not a bug but intentional.

See also my response to the same question from last year: https://groups.google.com/g/rpackageggir/c/4-9_SiMtwcU/m/M61FsMRiAQAJ

The essence is that working with units of time requires a different interpretation than working with units of acceleration. Zero time spent in a behaviour always has a meaning, even if a person never spends time in it, while zero acceleration only has meaning if a person actually participants in a behaviour.

Best,

Vincent

--
You received this message because you are subscribed to the Google Groups "R package GGIR" group.
To unsubscribe from this group and stop receiving emails from it, send an email to RpackageGGIR...@googlegroups.com.
To view this discussion, visit https://groups.google.com/d/msgid/RpackageGGIR/b6d718d9-a5da-4d59-b1f5-bc1eaecd25a2n%40googlegroups.com.