R: Smooth principal components analysis

sm.pca {sm} | R Documentation

## Smooth principal components analysis

### Description

This function calculates principal components in a manner which changes smoothly with a covariate. The smooth eigenvalues and eigenvector loadings can be plotted. A permutation test of equality of the components, both eigenvalues and eigenvectors, can be carried out.

### Usage

    sm.pca(x, Y, h, cor = TRUE, nperm = 100, pc = 1, ...)

### Arguments

`x` | either a vector of covariate values or a list object which is the output of a previous call to `sm.pca`. In the latter case, previously computed information is used to create plots and tests and the arguments `Y`, `h`, `cor` and `nperm` are not required.
`Y` | a matrix of responses whose rows correspond to the caovariate values.
`h` | the smoothing parameter which controls the smoothness of estimation with respect to the covariate `x`.
`cor` | a logical value indicating whether the correlation, rather than covariance, matrix should be used.
`nperm` | the number of permutations used in the permutation test and graphical reference band.
`pc` | an integar value indicating the component to be plotted against the covariate.
`...` | other optional parameters are passed to the `sm.options` function, through a mechanism which limits their effect only to this call of the function. Those relevant for this function are the following: `display` (here set to `"eigevalues"` or `"eigenvectors"`) `ngrid`, `xlab`; see the documentation of `sm.options` for their description.

### Details

Several further arguments may be set and these are passed to `sm.options`. Reelevant arguments for this function are `display` (`"eigenvalues"`, `"eigenvectors"`), `ngrid` and `df`. See `link{sm.options}` for further details.

The smoothing is performed by the local constant kernel method and the smoothing parameter corresponds to the standard deviation of a normal kernel function. If `h` is left unspecified then it is selected to correspond to the degrees of freedom set by the parameter `df`.

The reference band for a constant eigenvalue is constructed from the upper and lower pointwise 2.5 percentiles of the smooth eigenvalue curves from the data with permuted covariate values. The p-value compares the observed value of the difference between the smoothed and constant eigenvalues, summed over the covariate grid in `eval.points`, with the values generated from the permuted data.

In the eigenvector case, a reference band is computed from the percentiles of the curves from the permuted data, for each of the loadings. In order to plot all the loadings curves simultaneously, the locations where each curve lies inside its respective reference band are indicated by pale colour. The p-value compares the observed value of `1 - sum(e*e0)^2`, where `e` and `e0` are the eigenvectors under the smooth and constant scenarios (summed over the covariate grid), with the values generated from the permuted data. This test statistic differs from the one described in the Miller and Bowman (2012) reference below. It has been used as it conveniently handles the arbitrary sign of principal components.

When some components explain similar proportions of variance, the eigenvalues and eigenvectors can easily interchange, causing apparent sharp changes in the eigenvalue and eigenvector curves. It is difficult to track the components to avoid this.

### Value

a list with the following components:

xgrid

  a vector of values on the covariate scale at which the smooth estimates are constructed.

evals

  a matrix whose columns give the smooth eigenvalues for each component at the covariate values.

evecs

  a three-dimensional array whose third dimension corresponds to the covariate values and whose second dimension indexes the smooth components.

mhat

  a matrix whose columns give the estimated smooth means for each dimension of `Y` at the covariate values.

var.explained

  a matrix whose rows give the proportions of variance exaplined by each component at each covariate value.

xlab

  the label attached to the x-axis.

h

  the smoothing parameter used.

x

  the covariate values, after removal of missing cases.

Y

  the matrix of response values, after removal of missing cases.

cor

  a logical indicator of whether the correlation, rather than covariance, matrix is used in the construction of the eigenvalues and eigenvectors.

When a test or reference band is computed, the list has the additional components:

nperm

  the number of permutations used.

evals.perm

  the eigenvalues computed from the permuted data.

evecs.perm

  the eigenvectors computed from the permuted data.

When display contains `"eigenvalues"` or `"eigenvectors"`, the list has the additional components:

p.values

  the p-value for a test of constant eigenvalue for the component identified by `pc`.

p.vectors

  the p-value for a test of constant eigenvectors for the component identified by `pc`.

When display contains `"eigenvalues"`, the list has the additional component:

band

  a matrix whose two columns contain the boundaries of a reference band which indicates where the smooth eigenvalue curve should like if the hypothesis of no change in the eigenvalues with the ovariate is correct.

When display contains `"eigenvectors"`, the list has the additional components:

xgrid.plot

  a vector of values used for plotting the smooth eigenvectors.

evecs.plot

  a matrix whose rows contain the smooth eigenvectors at each value of `xgrid.plot`.

evecs.plot

  a matrix whose columns contain the colours for the line segments in each smooth eigenvector component.

### Side Effects

a plot on the current graphical device is produced, unless `display="none"`

### References

Miller, C. and Bowman, A.W. (2012). Smooth principal components for investigating changes in covarinces over time. Applied Statistics 61, 693–714.

### See Also

`sm.regression`, `sm.options`

### Examples

    ## Not run:
    Y    <- log(as.matrix(aircraft[ , -(1:2)]))
    year <- aircraft$Yr
    h    <- h.select(year, Y[ , 1], method = "df", df = 4)
    spca <- sm.pca(year, Y, h, display = "none")
    sm.pca(year, Y, h, display = "eigenvalues")
    sm.pca(year, Y, h, display = "eigenvectors", ylim = c(-1, 1))

    # The following code shows how the plots can be redrawn from the returned object

    spca <- sm.pca(year, Y, h, display = "eigenvalues")
    spca <- sm.pca(year, Y, h, display = "eigenvectors", ylim = c(-1, 1))

    with(spca, {
       ylim <- range(evals[ , 1], band)
       plot(xgrid, evals[ , 1], type = "n", ylab = "Variance", ylim = ylim)
       polygon(c(xgrid, rev(xgrid)), c(band[ , 1], rev(band[ , 2])),
               col = "lightgreen", border = NA)
       lines(xgrid, evals[ , 1], col = "red")
    })

    with(spca, {
       pc <- 1
       plot(range(xgrid.plot), range(evecs.plot), type = "n",
            xlab = "x", ylab = "PC loadings")
       for (i in 1:ncol(Y))
          segments(xgrid.plot[-length(xgrid.plot)],
                   evecs.plot[-nrow(evecs.plot), i],
                   xgrid.plot[-1], evecs.plot[-1, i],
                   col = col.plot[ , i], lty = i)
    })

    ## End(Not run)

[Package sm version 2.2-5.4 Index]