gp

Purpose

Create a Gaussian Process.

Synopsis

net = gp(nin, covarfn)
net = gp(nin, covarfn, prior)

Description

net = gp(nin, covarfn) takes the number of inputs nin for a Gaussian Process model with a single output, together with a string covarfn which specifies the type of the covariance function, and returns a data structure net. The parameters are set to zero.

The fields in net are

  type = 'gp'
  nin = number of inputs
  nout = number of outputs: always 1
  nwts = total number of weights and covariance function parameters
  bias = logarithm of constant offset in covariance function
  noise = logarithm of output noise variance
  inweights = logarithm of inverse length scale for each input 
  covarfn = string describing the covariance function:
      'sqexp'
      'ratquad'
  fpar = covariance function specific parameters (1 for squared exponential,
   2 for rational quadratic)
  trin = training input data (initially empty)
  trtargets = training target data (initially empty)

net = gp(nin, covarfn, prior) sets a Gaussian prior on the parameters of the model. prior must contain the fields pr_mean and pr_variance. If pr_mean is a scalar, then the Gaussian is assumed to be isotropic and the additional fields net.pr_mean and pr_variance are set. Otherwise, the Gaussian prior has a mean defined by a column vector of parameters prior.pr_mean and covariance defined by a column vector of parameters prior.pr_variance. Each element of prmean corresponds to a separate group of parameters, which need not be mutually exclusive. The membership of the groups is defined by the matrix prior.index in which the columns correspond to the elements of prmean. Each column has one element for each weight in the matrix, in the order defined by the function gppak, and each element is 1 or 0 according to whether the parameter is a member of the corresponding group or not. The additional field net.index is set in this case.

See Also

gppak, gpunpak, gpfwd, gperr, gpcovar, gpgrad
Pages: Index

Copyright (c) Ian T Nabney (1996-9)