relist                 package:utils                 R Documentation

_A_l_l_o_w _R_e-_L_i_s_t_i_n_g _a_n _u_n_l_i_s_t()_e_d _O_b_j_e_c_t

_D_e_s_c_r_i_p_t_i_o_n:

     'relist()' is an S3 generic function with a few methods in order
     to allow easy inversion of 'unlist(obj)' when that is used with an
     object 'obj' of (S3) class '"relistable"'.

_U_s_a_g_e:

     relist(flesh, skeleton)
     ## Default S3 method:
     relist(flesh, skeleton = attr(flesh, "skeleton"))
     ## S3 method for class 'factor':
     relist(flesh, skeleton = attr(flesh, "skeleton"))
     ## S3 method for class 'list':
     relist(flesh, skeleton = attr(flesh, "skeleton"))
     ## S3 method for class 'matrix':
     relist(flesh, skeleton = attr(flesh, "skeleton"))

     as.relistable(x)
     is.relistable(x)

     ## S3 method for class 'relistable':
     unlist(x, recursive = TRUE, use.names = TRUE)

_A_r_g_u_m_e_n_t_s:

   flesh: a vector to be relisted

skeleton: a list, the structure of which determines the structure of
          the result

       x: an R object, typically a list (or vector).

recursive: logical.  Should unlisting be applied to list components of
          'x'?

use.names: logical.  Should names be preserved?

_D_e_t_a_i_l_s:

     Some functions need many parameters, which are most easily
     represented in complex structures, e.g., nested lists. 
     Unfortunately, many mathematical functions in R, including 'optim'
     and 'nlm' can only operate on functions whose domain is a vector. 
     R has 'unlist()' to convert nested list objects into a vector
     representation.  'relist()', it's methods and the functionality
     mentioned here provide the inverse operation to convert vectors
     back to the convenient structural representation. This allows
     structured functions (such as 'optim()') to have simple
     mathematical interfaces.

     For example, a likelihood function for a multivariate normal model
     needs a variance-covariance matrix and a mean vector.  It would be
     most convenient to represent it as a list containing a vector and
     a matrix.  A typical parameter might look like


           list(mean=c(0, 1), vcov=cbind(c(1, 1), c(1, 0))).

     However, 'optim' cannot operate on functions that take lists as
     input; it only likes numeric vectors.  The solution is conversion.
     Given a function 'mvdnorm(x, mean, vcov, log=FALSE)' which
     computes the required probability density, then


             ipar <- list(mean=c(0, 1), vcov=cbind(c(1, 1), c(1, 0)))
             initial.param <- as.relistable(ipar)

             ll <- function(param.vector)
             {
                param <- relist(param.vector, skeleton=ipar))
                -sum(mvdnorm(x, mean = param$mean, vcov = param$vcov,
                             log = TRUE))
             }

             optim(unlist(initial.param), ll)

     'relist' takes two parameters: skeleton and flesh.  Skeleton is a
     sample object that has the right 'shape' but the wrong content. 
     'flesh' is a vector with the right content but the wrong shape. 
     Invoking


         relist(flesh, skeleton)

     will put the content of flesh on the skeleton.  You don't need to
     specify skeleton explicitly if the skeleton is stored as an
     attribute inside flesh. In particular, if flesh was created from
     some object obj with 'unlist(as.relistable(obj))' then the
     skeleton attribute is automatically set.  (Note that this does not
     apply to the example here, as 'optim' is creating a new vector to
     pass to 'll' and not its 'par' argument.)

     As long as 'skeleton' has the right shape, it should be a precise
     inverse of 'unlist'.  These equalities hold:


        relist(unlist(x), x) == x
        unlist(relist(y, skeleton)) == y

        x <- as.relistable(x)
        relist(unlist(x)) == x


_V_a_l_u_e:

     an object of (S3) class '"relistable"' (and '"list"').

_A_u_t_h_o_r(_s):

     R Core, based on a code proposal by Andrew Clausen.

_S_e_e _A_l_s_o:

     'unlist'

_E_x_a_m_p_l_e_s:

      ipar <- list(mean=c(0, 1), vcov=cbind(c(1, 1), c(1, 0)))
      initial.param <- as.relistable(ipar)
      ul <- unlist(initial.param)
      relist(ul)
      stopifnot(identical(relist(ul), initial.param))

