activityTacsatAnalyse() by VE_REF help please

[TimWOceanSciences]

Hello

I want to use activityTacsatAnalyse() , but I want to do the analysis by vessel. I'm only interested in one gear type. I modified the function to use the SI_SPCA instead of SI_SP, but I can't get the analysis to work using VE_REF which from the documentation/function does seem possible. Could you possibly give an example of using activityTacsatAnalyse() but doing by vessel instead or with year? The below is from your wiki.

Thank you for any help.

Tim Whitton

https://github.com/nielshintzen/vmstools/wiki/Practicals3
tacsat$YEAR <- format(tacsat$SI_DATIM,format='%Y')
#In order to get this working, I cannot have NA's as gear type, so we can assign them as 'NO_GEAR' by:
tacsatp$LE_GEAR[which(is.na(tacsatp$LE_GEAR)==T)] <- "NO_GEAR"

#Analyse a selection of the data, a window will pop-up!
require(tcltk) #- install.packages("tcltk",repos=getOption("repos"))
X11()
storeScheme <- activityTacsatAnalyse(subset(tacsatp,LE_GEAR == "OTM" & format(SI_DATIM,format='%Y') == 1801),units="year",analyse.by="LE_GEAR")