FLComp.R

# FLComp - Basic VIRTUAL class for all FLQuant-based classes
# FLCore/R/FLComp.R

# Copyright 2003-2018 FLR Team. Distributed under the GPL 2 or later
# Maintainer: Iago Mosqueira, EC JRC

# summary		{{{
#' @rdname summary-methods
#' @aliases summary,FLComp-method
setMethod("summary", signature(object="FLComp"),
	function(object, ...){

		cat("An object of class \"", class(object), "\"\n\n", sep="")
		cat("Name:", object@name, "\n")
		
    desc <- ifelse(nchar(object@desc) > getOption("width"),
      paste(substr(object@desc, 1, getOption("width") - 6), "[...]"),
      object@desc) 
		cat("Description:", desc, "\n")

    # character slots
		cnames <- getSlotNamesClass(object, 'character')
    cnames <- cnames[!cnames%in%c('name', 'desc')]
    if(length(cnames) > 0)
    {
      for (s in cnames)
        cat(paste(toupper(substring(s, 1,1)), substring(s, 2), sep=""), ": ",
          slot(object, s), "\n")
    }

    # FLArray slots
		qnames <- getSlotNamesClass(object, 'FLArray')
    dms <- dims(object)
    # Quant
		cat("Quant:", dms$quant, "\n")
    # Dims
		cat("Dims: ", quant(slot(object, qnames[1])), "\tyear\tunit\tseason\tarea\titer\n")
		cat("", dms[[2]], dim(slot(object, qnames[1]))[-c(1,6)],
      dms$iter, "\n\n", sep="\t")
		# Range
    cat("Range: ", paste(sub('plusgroup', 'pgroup', names(object@range)),
      collapse="\t"), "\n")
		cat("", object@range, "\n\n", sep="\t")
	}
)	# }}}

# show {{{
setMethod('show', signature('FLComp'),
  function(object) {
    summary(object)
  }
) # }}}

# print {{{
setMethod('print', signature(x='FLComp'),
  function(x) {
    showDefault(x, FALSE)
  }
) # }}}

# window    {{{
setMethod("window", signature(x="FLComp"),
	  function(x, start=dims(x)$minyear, end=dims(x)$maxyear, extend=TRUE, 
      frequency=1) {
      
      x <- qapply(x, window, start=start, end=end, extend=extend,
        frequency=frequency)
  		
      x@range["minyear"] <- as.numeric(dims(x)$minyear)
	    x@range["maxyear"] <- as.numeric(dims(x)$maxyear)

		return(x)
	}
)	# }}}

# propagate {{{
setMethod("propagate", signature(object="FLComp"),
	function(object, iter, fill.iter=TRUE) {

		# GET object iters
		mit <- unlist(qapply(object, function(x) dim(x)[6]))

		# CHECK iter can only be 1 or dim(object)[6]
		if(sum(mit) > length(mit) & !iter %in% mit)
			stop("incompatible number of iters requested")

		# GET slots to extend
		idx <- mit[mit != iter]

		for(sl in names(idx)) {
			slot(object, sl) <- propagate(slot(object, sl), iter, fill.iter=fill.iter)
		}

    # DO for FLPar
		pnms <- getSlots(class(object))
		pnames <- names(pnms)[pnms == "FLPar"]
		for(i in pnames)
			slot(object, i) <- propagate(slot(object, i), iter)
		return(object)
	}
) # }}}

# iter {{{
setMethod("iter", signature(obj="FLComp"),
  function(obj, iter) {

  # copy the iterate into the new slots
  names. <- c(getSlotNamesClass(obj, 'FLArray'), getSlotNamesClass(obj, 'FLPar'))
  
  for(s. in names.) {
    if(dims(slot(obj, s.))$iter == 1) {
      slot(obj, s.) <- propagate(iter(slot(obj, s.), 1), length(iter))
      dimnames(slot(obj, s.))$iter <- ac(iter)
    }
    else {
      slot(obj, s.) <- iter(slot(obj, s.), iter)
    }
  }
  return(obj)
  }
) # }}}

# iter<-  {{{
setMethod("iter<-", signature(object="FLComp", value="FLComp"),
	function(object, iter, value)
	{
		object[,,,,,iter] <- value
		return(object)
	}
)   # }}}

# transform {{{
setMethod("transform", signature(`_data`="FLComp"),
function(`_data`, ...)
  {

	env <- new.env(parent=parent.frame())

	for (i in slotNames(`_data`)) {
		assign(i, slot(`_data`, i), envir=env)
	}

	args <- eval(substitute(list(...)), env)

	# IF ... is only FLQuants
	if(is.null(names(args)) & is(args[[1]], 'FLQuants')) {
		args <- unlist(args, recursive=FALSE)
	}

	# IF both FLQuants and promises
	if(any(names(args) == "")) {
		fqs <- unlist(lapply(args, is, 'FLQuants'))
		args <- c(args[!fqs], unlist(args[fqs], recursive=FALSE))
	}

	for (i in 1:length(args)) {
		slot(`_data`, names(args)[i]) <- args[[i]]
	}

	if(validObject(`_data`))
		return(`_data`)
	stop('Attempt to modify object incorrectly: check input dimensions')
	}
)	# }}}

# qapply		{{{
setMethod('qapply', signature(X='FLComp', FUN='function'),
	function(X, FUN, ..., exclude=missing, simplify=FALSE) {
		
    FUN <- match.fun(FUN)
		slots <- getSlotNamesClass(X, 'FLArray')

		if(!missing(exclude))
      slots <- slots[!slots %in% exclude]

		if(is(do.call(FUN, list(slot(X, slots[1]), ...)), 'FLArray')) {
			res <- X
			for (i in slots)
				slot(res, i) <- do.call(FUN, list(slot(X,i), ...))
		}
		else {
			res  <- vector('list', 0)
			for (i in slots)
				res[[i]] <- do.call(FUN, list(slot(X,i), ...))
		}

    if(simplify)
      res <- unlist(res)

		return(res)
	}
)   # }}}

# trim     {{{

#' @rdname trim
#' @aliases trim,FLComp-method

setMethod("trim", signature("FLComp"),
	function(x, ...)
	{
	  args <- list(...)

    names <- getSlotNamesClass(x, 'FLArray')

    c1 <- args[[quant(slot(x, names[1]))]]
	  c2 <- args[["year"]]

    # FLQuants with quant
    x <- qapply(x, trim, ...)
    
    # range
  	if (length(c1) > 0)
    {
    	x@range["min"] <- as.numeric(c1[1])
	    x@range["max"] <- as.numeric(c1[length(c1)])
	  }
  	if (length(c2)>0 )
    {
    	x@range["minyear"] <- as.numeric(c2[1])
	    x@range["maxyear"] <- as.numeric(c2[length(c2)])
  	}
	  return(x)
	}
) # }}}

# units	    {{{

#' @rdname units-FLCore
#' @details The complete set of *units* for a complex object can be obtained
#' as a named *list*.

setMethod("units", signature(x="FLComp"), function(x)
	qapply(x, units)
)
#}}}

# units<-      {{{

#' @rdname units-FLCore
#' @details Assignment of *units* to the *FLQuant* slots of a complex object
#' can be carried out passing a named *list* or *character* vector containing
#' the units for the slots to be modified.
#' @examples
#' data(ple4)
#' units(ple4) <- list(harvest="hr")
#' units(ple4) <- c(harvest="hr")

setMethod("units<-", signature(x="FLComp", value="list"),
    function(x, value) {
      for(i in seq(along=value))
        if(is.character(value[[i]]))
          units(slot(x, names(value[i]))) <- value[[i]]
      return(x)
	}
)

#' @rdname units-FLCore
setMethod("units<-", signature(x="FLComp", value="character"),
    function(x, value) {
      value <- as.list(value)
      units(x) <- value
      return(x)
	}
)

#' @rdname units-FLCore
setMethod("units<-", signature(x="FLComp", value="function"),
  function(x, value) {
    value <- do.call(value, list(x))
    units(x) <- value
    return(x)
	}
) # }}}

# '['       {{{
#' @rdname Extract
#' @aliases [,FLComp,ANY,ANY,ANY-method
setMethod('[', signature(x='FLComp'),
	function(x, i, j, k, l, m, n, ..., drop=FALSE) {

    # FIND FLQs, FLAs or FLChs
		qnames <- names(getSlots(class(x))[
      getSlots(class(x))=="FLQuant" |
      getSlots(class(x))=="FLArray" |
      getSlots(class(x))=="FLCohort"])

		dx <- dim(slot(x, qnames[1]))
    args <- list(drop=FALSE)

		if (!missing(i)) {
      args <- c(args, list(i=i))
		}
		if (!missing(j))
      args <- c(args, list(j=j))
		if (!missing(k))
      args <- c(args, list(k=k))
		if (!missing(l))
      args <- c(args, list(l=l))
		if (!missing(m))
      args <- c(args, list(m=m))
		if (!missing(n))
      args <- c(args, list(n=n))

    for(q in qnames)
      slot(x, q) <- do.call('[', c(list(x=slot(x,q)), args))
    
    # range
		if (!missing(i)) {
    	x@range['min'] <- dims(slot(x, qnames[1]))$min
    	x@range['max'] <- dims(slot(x, qnames[1]))$max
		}
		if (!missing(j)) {
	    x@range['minyear'] <- dims(slot(x, qnames[1]))$minyear
  	  x@range['maxyear'] <- dims(slot(x, qnames[1]))$maxyear
		}
    return(x)
    }
)   # }}}

# '[<-'            {{{
#' @rdname Extract
#' @aliases [<-,FLComp,ANY,ANY,ANY-method
setMethod("[<-", signature(x="FLComp"),
	function(x, i, j, k, l, m, n, ..., value) {

    # SLOTS to work on
		qnames <- getSlotNamesClass(x, 'FLArray')

    # dims
		dx <- dim(slot(x, qnames[1]))

    # quant and iter by slot
    di <- qapply(x, function(y) seq(1, dim(y)[1]))
    dj <- qapply(x, function(y) seq(1, dim(y)[2]))
    dk <- qapply(x, function(y) seq(1, dim(y)[3]))
    dl <- qapply(x, function(y) seq(1, dim(y)[4]))
    dm <- qapply(x, function(y) seq(1, dim(y)[5]))
    dn <- qapply(x, function(y) seq(1, dim(y)[6]))
		
    if (!missing(i))
      di <- lapply(di, function(x) x <- i)
    if (!missing(j))
      dj <- lapply(dj, function(x) x <- j)
    if (!missing(k))
      dk <- lapply(dk, function(x) x <- k)
    if (!missing(l))
      dl <- lapply(dl, function(x) x <- l)
    if (!missing(m))
      dm <- lapply(dm, function(x) x <- m)
		if (!missing(n))
      dn <- lapply(dn, function(x) x <- n)
    
    for(q in qnames) {
      slot(x, q)[di[[q]],dj[[q]],dk[[q]],dl[[q]],dm[[q]],dn[[q]]] <- slot(value, q)
    }

   	return(x)
	}
)   # }}}

# as.data.frame        {{{
setMethod("as.data.frame", signature(x="FLComp", row.names="missing", optional="missing"),
	function(x, row.names, optional, drop=FALSE, cohort=FALSE, units=FALSE)
	{
    qnames <- getSlotNamesClass(x, 'FLArray')
    quant <- quant(slot(x, qnames[1]))
	  df   <- data.frame()
    for(s in qnames)
		{
      sdf <- as.data.frame(slot(x, s), cohort=cohort, units=units)
      sdf[[quant]] <- as.character(sdf[[quant]])
      dfq <- cbind(slot=s, sdf)

			df  <- rbind(df, dfq)
	  }
    # drop
    if(drop) {
      idx <- apply(df, 2, function(x) length(unique(x))) == 1
      df <- df[, !idx]
    }

		# add attributes
		attributes(df)$desc <- x@desc
		attributes(df)$name <- x@name
		attributes(df)$range <- x@range

		return(df)
	}
)   # }}}

# mcf	{{{
setMethod('mcf', signature(object='FLComp'),
	function(object, second) {

	qdims <- unlist(qapply(object, function(x) dim(x)[1]))
	qdims <- names(qdims[qdims==max(qdims)][1])

	dimnames <- list()
	dob <- dimnames(slot(object, qdims))
	dse <- dimnames(slot(second, qdims))

	for(i in names(dob))
		dimnames[[i]] <- unique(c(dob[[i]], dse[[i]]))

	foo <- function(x, dimnames) {
		if(all(dimnames(x)[[1]] == 'all'))
			return(FLQuant(x, dimnames=dimnames[-1]))
		else
			return(FLQuant(x))
	}

	res <- new('FLlst')
	res[[1]] <- qapply(object, foo, dimnames=dimnames)
	res[[2]] <- qapply(second, foo, dimnames=dimnames)

	return(res)
	}
)	# }}}

# dims {{{
setMethod("dims", signature(obj="FLComp"),
  # Returns a list with different parameters
  function(obj, ...) {
    qnames <- getSlotNamesClass(obj, 'FLArray')

    range <- as.list(range(obj))
    dimsl <- qapply(obj, dim)
    dnames <- qapply(obj, dimnames)
    dimsl <- dimsl[!names(dnames) %in% 'fbar']
    dims <- matrix(unlist(dimsl), ncol=6, byrow=TRUE)

    iter <- max(dims[,6])
    pnames <- getSlotNamesClass(obj, 'FLPar')
		if(length(pnames) > 0)
			for(p in pnames)
				iter <- max(iter, length(dimnames(slot(obj, p))$iter))

    # Hack for FLBRP, drop fbar slot
    dnames <- dnames[!names(dnames) %in% 'fbar']
    if(is(obj, "FLFishery"))
      quants <- c(hperiod=1:2)
    else
      quants <- lapply(dnames, function(x) x[[1]])[unlist(lapply(dimsl,
        function(x) x[1] == max(dims[,1])))][[2]]
    res <- list(
      quant = quant(slot(obj, qnames[1])),
      quants = max(dims[,1]),
      min = ifelse(!is.na(suppressWarnings(as.numeric(quants[1]))),
        suppressWarnings(as.numeric(quants[1])), as.numeric(NA)),
      max = ifelse(!is.na(suppressWarnings(as.numeric(quants[max(dims[,1])]))),
        suppressWarnings(as.numeric(quants[max(dims[,1])])), as.numeric(NA)),
      year = max(dims[,2]),
      minyear = as.numeric(dnames[[1]]$year[1]),
      maxyear = as.numeric(dnames[[1]]$year[max(dims[,2])]),
      plusgroup = ifelse('plusgroup' %in% names(range), range$plusgroup, NA),
      unit = max(dims[,3]),
      season = max(dims[,4]),
      area = max(dims[,5]),
      iter = iter)
    res <- lapply(res, function(x) if(is.null(x)) return(as.numeric(NA)) else return(x))
    names(res)[2] <- res$quant

    return(res)
  }
) # }}}

# model.frame {{{
setMethod('model.frame', signature(formula='FLComp'),
	function(formula, mcf=TRUE, ...)
  {
    lst <- FLQuants()
    names <- getSlotNamesClass(formula, 'FLQuant')
    for(i in names)
      lst[[i]] <- slot(formula, i)
    names(lst) <- names
    if(mcf)
      lst <- mcf(lst)
    return(model.frame(lst, ...))
  }
)
# }}}

# range {{{
setMethod("range", "FLComp",
  function(x, i='missing', ..., na.rm = FALSE) {
    if(missing(i))
      slot(x, 'range')
    else
      slot(x, 'range')[i]
  }
)

setReplaceMethod("range", signature(x="FLComp", i="missing", value="numeric"),
  function(x, value) {
    slot(x, 'range')[names(value)] <- value
      return(x)
  }
)

setReplaceMethod("range", signature(x="FLComp", i="ANY", value="numeric"),
  function(x, i, value) {

    # HACK for "pgroup", as displayed by summary()
    if(any(i == "pgroup"))
      i[i == "pgroup"] <- "plusgroup"

    slot(x, 'range')[i] <- value
    if(validObject(x))
      return(x)
    else
      stop("Object not valid")
  }
) # }}}

# expand  {{{
setMethod('expand', signature(x='FLComp'),
  function(x, ...)
  {
    x <- qapply(x, expand, ...)

    # range
    range <- qapply(x, function(x) dimnames(x)[[1]])
    slot <- names(which.max(lapply(range, length)))
    dnames <- dimnames(slot(x, slot))
    range(x, c('min', 'max', 'minyear', 'maxyear')) <- c(as.numeric(dnames[[1]][1]),
      as.numeric(dnames[[1]][length(dnames[[1]])]), as.numeric(dnames[[2]][1]),
      as.numeric(dnames[[2]][length(dnames[[2]])]))

    return(x)
  }
) # }}}

# slots  {{{
setMethod(slots, signature(object='FLComp', name='character'),
  function(object, name, ...) {

		args <- list(...)

    # args
    if(length(args) > 0)
      if(all(unlist(lapply(args, function(x) is(x, 'character')))))
        name <- c(name, unlist(args))
      else
        stop(paste('Only character vectors for slot names allowed:',
          unlist(args[!unlist(lapply(args, function(x) is(x, 'character')))])))

		res <- vector(mode='list', length=length(name))
    names(res) <- name

    for (i in name)
      res[[i]] <- slot(object, i)


    return(new(getPlural(res[[1]]), res))
  }
) # }}}

# age/year vectors {{{
setMethod("rngyear", "FLComp", function(object){
	object@range[c("minyear","maxyear")]
})

setReplaceMethod("rngyear", "FLComp", function(object, value){
	object@range[c("minyear","maxyear")] <- value
	object
})

setMethod("rngage", "FLComp", function(object){
	object@range[c("min","max")]
})

setReplaceMethod("rngage", "FLComp", function(object, value){
	object@range[c("min","max")] <- value
	object
})

setMethod("vecyear", "FLComp", function(object){
	rng <- object@range[c("minyear","maxyear")]
	rng[1]:rng[2]
})

setMethod("vecage", "FLComp", function(object){
	rng <- object@range[c("min","max")]
	rng[1]:rng[2]
}) # }}}

# metrics {{{

#' @rdname metrics
#' @examples
#' data(ple4)
#' # missing
#' metrics(ple4)
#' # metrics = function
#' metrics(ple4, metrics=function(x) FLQuants(SSB=ssb(x), REC=rec(x),
#'   F=fbar(x), SSBREC=ssb(x) / rec(x)))
#' # metrics = formula
#' metrics(ple4, metrics=~ssb)
#' metrics(ple4, metrics=list(SSB=~ssb))
#' metrics(ple4, metrics=list(SBMSY=~ssb/SBMSY), FLPar(SBMSY=3.87e4))
#' # metrics = list
#' metrics(ple4, metrics=list(SSB=ssb, REC=rec, F=fbar))
#' metrics(ple4, metrics=list(SSB=~ssb, REC=rec, F=fbar))

setMethod("metrics", signature(object="FLComp", metrics="list"),
  function(object, metrics, ...) {
    return(FLQuants(lapply(metrics, function(x) {
      # CALL each function
      do.call("metrics", list(object=object, metrics=x, ...))
    })))
  }
)

setMethod("metrics", signature(object="FLComp", metrics="function"),
  function(object, metrics) {
    # CALL function
    return(do.call(metrics, list(object)))
  }
)

setMethod("metrics", signature(object="FLComp", metrics="character"),
  function(object, metrics) {
    # CALL function
    return(do.call(metrics, list(object)))
  }
)

setMethod("metrics", signature(object="FLComp", metrics="formula"),
  function(object, metrics, ...) {
    
    args <- list(...)

    # PARSE ...'s FLPar objects and vectors
    if(length(args) > 0) {
      
      names(args)[length(names(args)) == 0] <- "a"
      
      args <-  mapply(function(x, y) {
        if(is(x, "FLPar"))
          as(x, "list")
        else
          setNames(list(x), y)
      }, args, names(args), SIMPLIFY=FALSE)

      # COERCE into single-depth named list
      if(length(args) == 1)
        args <- args[[1]]
      else
        args <- do.call(append, unname(args))
    }

    # EVAL functions in metrics
    funs <- all.vars(metrics)[!all.vars(metrics) %in% names(args)]
    mets <- setNames(lapply(funs, do.call, list(object)), funs)

    # EVALUATE formula w/ object, mets and args
    res <- eval(metrics[[2]], c(object, mets, args)) 

    return(res)
  }
)
# }}}

# slim {{{

#' @rdname slim
#' @examples
#' data(ple4)
#' # Extend all of ple4 to 50 iters
#' ple4 <- propagate(ple4, 50)
#' # Add variability in catch.n
#' catch.n(ple4) <- rlnoise(50, log(catch.n(ple4)), log(catch.n(ple4))/10)
#' summary(ple4)
#' # slim object by dropping identical iters
#' sple4 <- slim(ple4)
#' summary(sple4)

setMethod("slim", signature(object="FLComp"),
  function(object, ...) {
    
    # FIND repeated iters
    res   <- qapply(object, function(x) {
      # CHECK sum(var) along iters == 0
      if(all(is.na(iterVars(x) / 0)))
        return(x[,,,,,1])
      else
        return(x)
      }, ...)

    return(res)
  }
) # }}}

# verify {{{

#' @rdname verify

setMethod("verify", signature(object="FLComp"),
  function(object, ..., report=TRUE) {
  
  rules <- list(...)
  
  # FIX names if missing
  nms <- names(rules)
  names(rules)[nms == ""] <- paste0("R", seq(length(rules)))[nms == ""]

  # PREPARE rules
  rules <- lapply(rules, function(x) {
    # formula, deparse
    if(is(x, 'formula')) return(list(rule=deparse(x[[length(x)]])))
    # character, create list
    else if(length(deparse(x))==1) return(list(rule=x))
    # list
    else if(is(x, "list"))
      # formula, deparse
      if(is(x[[1]], 'formula')) {
       x[[1]] <- deparse(x[[1]][[length(x)]])
        return(x)
      }
    else return(x)})
  
  # APPLY rules
  res <- lapply(rules, function(x) {
    expr <- parse(text=x$rule)
    args <- all.names(expr, functions=FALSE)

    env <- list2env(x["rule" != names(x)])

    # calls to methods
    inps <- lapply(args, do.call, list(object), envir=env)
    names(inps) <- args

    return(eval(expr, inps, enclos=env))
    })

  # CONVERT list to data.frame
  res <- data.frame(do.call(rbind, lapply(res, function(x) {
    # CREATE result summary vector
    c(items=length(x), passes=sum(x, na.rm=TRUE), fails=sum(!x, na.rm=TRUE),
      NAs=sum(is.na(x)))
    })),
    row.names=seq(1, length(res)))

  # data.frame: name, items, passes, fails, NAs, valid, rule
  res <- cbind(name=names(rules), res, valid=with(res, passes/items == 1),
    rule=unlist(lapply(rules, "[[", "rule")))

  if(report)
    return(res)
  else
    return(all(res$valid))
}) # }}}

# divide {{{
setMethod("divide", signature(object="FLComp"),
  function(object, names=dimnames(object)$iter) {

    its <- dims(object)$iter

    if(its == 1)
      return(object)

    res <- lapply(seq(its), iter, obj=object)

    names(res) <- names
 
    return(do.call(getPlural(res[[1]]), res))

  }
) # }}}

# split {{{

#' @rdname split

setMethod("split", signature(x="FLComp", f="vector"),
  function(x, f) {

  return(do.call(getPlural(x), lapply(setNames(nm=unique(f)),
    function(i) {
      iter(x, f == i)
    }
  )))

  }
) # }}}

# combine {{{

setMethod('combine', signature(x='FLComp', y='FLComp'),
  function(x, y, ..., check=FALSE) {

    args <- c(list(x, y), list(...))

    # CHECK input classes match exactly
    if(length(unique(lapply(args, is))) > 1)
      stop("combine can only operate on objects of identical class")
    
    ds <- lapply(args, dims)

    # CHECK dimnames but iter
    if(check) {

		  idi <- names(ds[[1]])!="iter"
      
      # COMPARE dims(x)[-'iter')]
      if(length(unique(lapply(ds, "[", idi))) > 1)
        stop("Object dimensions but iter must match")
    }

    # CALCULATE iters
    its <- sum(unlist(lapply(ds, "[[", "iter")))

    # PROPAGATE object for output
    res <- propagate(iter(x, 1), its)

    # KEEP iter dimnames if unique
    itns <- unlist(lapply(args, function(x) dimnames(x)$iter))

    # CHECK iter dimnames are unique
    if(length(itns) > length(unique(itns)))
      itns <- ac(seq(1, its))

    # GET iter limits
    ite <- cumsum(unlist(lapply(ds, "[", "iter")))
    its <- ite - unlist(lapply(ds, "[", "iter")) + 1
    
    # ASSIGN
    for(i in seq(length(its)))
      iter(res, seq(its[i], ite[i])) <- args[[i]]
    
    return(res)
  }
)

# }}}

# dimnames {{{
setMethod("dimnames", signature(x="FLComp"),
  function(x) {
    
    # GET dimnames by FLQ/FLA slot
    dns <- qapply(x, dimnames)

    # FIND unique combinations
    udns <- unique(dns)

    # GET length by dim for each
    len <- lapply(udns, lapply, length)

    # TRAVERSE list by dim and return the longest element
    out <- lapply(names(udns[[1]]), function(x) {
        xlen <- unlist(lapply(len, "[[", x))
        idx <- which(xlen == max(xlen))[1]
        return(udns[[idx]][[x]])
    })

    # FIX missing names
    names(out) <- names(udns[[1]])

    return(out)
  }) # }}}

# dimnames<- {{{
setMethod('dimnames<-', signature(x='FLComp', value='list'),
  function(x, value)
  {
    res <- qapply(x, function(y) {
      dimnames(y) <- value
      return(y)
    })
  }
) # }}}

# update {{{
setMethod("update", signature(object="FLComp"),
  function(object, ...) {

    args <- list(...)
    
    # HANDLE list, e.g. FLQuants
    idx <- unlist(lapply(args, is, 'list'))
    args <- c(args[!idx], unlist(args[idx], recursive=FALSE))

    # CHECK all elements are named
    if(any(names(args) %in% c(NULL, character(1))))
      stop("Object to assign to slots must be named")

    # LOOP over value
    for(i in names(args)) {

     # GET dimnames to act on
      dms <- setNames(dimnames(args[[i]]), letters[9:14])

      # ASSIGN to subset of object the modified subset
      object[dms$i, dms$j, dms$k, dms$l, dms$m, dms$n] <- 
        do.call(paste0(i, "<-"), list(do.call("[", c(list(object), dms)),
          value=args[[i]]))
    }

    return(object)
  }
) # }}}

# dim {{{
setMethod("dim", signature(x="FLComp"),
  function(x) {

    # FIND FLQuant slots
    fqs <- which(getSlots(class(x)) == 'FLQuant')

    # TAKE first one
    slot <- names(fqs)[1]

    return(dim(slot(x, slot)))
  }
) # }}}