g.shell.GGIR(
               mode=c(1,2,3,4,5),
               datadir="C:/Users/psdmp/Reduction/test/data",
               outputdir="C:/Users/psdmp/Reduction/test/results",
               desiredtz ="Australia/NSW",
               do.report=c(2,4,5),
               #=====================
               # Part 1
               #=====================
               windowsizes = c(5,900,3600), #default: Epoch length, non-wear detection resolution, non-wear detection, evaluation window 
               do.cal=TRUE,   # Apply autocalibration? (recommended)
               do.enmo = TRUE,
               do.anglez=TRUE, #needed for sleep analyses with HDCZA
               #chunksize=1,   #size of data chunks to be read (value = 1 is maximum)
               printsummary=TRUE,
               #===================== 
               # Part 2
               #=====================
               strategy = 1,   #limit anlyses to times specified by hrs.del.start, hrs.del.end and maxdur
               hrs.del.start = 0,      
               hrs.del.end = 0,
               maxdur = 7,    # How many days of analyses do yuo expect?
               includedaycrit = 16,
               qlevels = c(c(1380/1440),c(1410/1440),c(1438/1440)),
               qwindow=c(0,6, 8.5, 15.5, 18.5),
               ilevels = c(seq(0,400,by=50),8000),
               iglevels = c(seq(0,4000,by=25),8000),
               mvpathreshold =c(201),
               bout.metric = 4,    # using sliding window across the, data to test bout criteria per window and do not allow for breaks of 1 minute or
                                   # longer. Requires the first and last epoch to require the threshold criteria.
               closedbout = FALSE, # If TRUE then count breaks in a bout towards the bout duration. If FALSE then
                                   # only count time spent above the threshold towards the bout duration.
               printsummary=TRUE,
               #=====================
               # Part 3
               #=====================
               anglethreshold = 5,
               timethreshold = 5,
               acc.metric="ENMO", 
               ignorenonwear=TRUE, # if TRUE non-wear is not detected as sleep (if FALSE then it will work with imputed data)
               overwrite=FALSE,
               constrain2range=TRUE,
               do.part3.pdf=TRUE, 
               do.parallel=FALSE,
               printsummary=TRUE,
               #=====================
               # Part 4
               #=====================
               excludefirstlast = FALSE,
               includenightcrit = 16, # number of minimum valid hours over 24 hr to attempt sleep analysis
               def.noc.sleep = 1,     # Apply HDCZA algorithm to sleep detection
               outliers.only = FALSE,
               criterror = 4,
               relyonsleeplog = FALSE,
               do.visual = TRUE,
               printsummary=TRUE,
               #=====================
               # Part 5
               #=====================
               excludefirstlast.part5 =TRUE,
               threshold.lig = c(30,40,50), #thresholds for inactivity 
               threshold.mod = c(101,201),  #thresholds for moderate 
               threshold.vig = c(429,707),  #thresholds for vigorous 
               boutcriter = 0.8,            
               boutcriter.in = 0.9,         
               boutcriter.lig = 0.8,1.0,
               boutcriter.mvpa = 0.8,       
               boutcriter.vig = 0.8,        
               boutdur.in = c(1,10,20,30), 
               boutdur.lig = c(1,10),
               boutdur.mvpa = c(1,10),
               printsummary=TRUE,
               #=====================
               # Visual report
               #=====================
               timewindow = c("WW","MM"),
               visualreport=TRUE)