PermTestReg {pathwayPCA} | R Documentation

## AES-PCA permutation test of continuous response for pathway PCs

### Description

Given an `OmicsReg` object and a list of pathway PCs from the `ExtractAESPCs` function, test if each pathway with features recorded in the bio-assay design matrix is significantly related to the continuous response.

### Usage

    PermTestReg(
      OmicsReg,
      pathwayPCs_ls,
      numReps = 0L,
      parallel = FALSE,
      numCores = NULL,
      ...
    )

    ## S4 method for signature 'OmicsReg'
    PermTestReg(
      OmicsReg,
      pathwayPCs_ls,
      numReps = 0L,
      parallel = FALSE,
      numCores = NULL,
      ...
    )

### Arguments

`OmicsReg` | A data object of class `OmicsReg`, created by the `CreateOmics` function.
`pathwayPCs_ls` | A list of pathway PC matrices returned by the `ExtractAESPCs` function.
`numReps` | How many permutations to estimate the p-value? Defaults to 0 (that is, to estimate the p-value parametrically). If `numReps` > 0, then the non-parametric, permutation p-value will be returned based on the number of random samples specified.
`parallel` | Should the computation be completed in parallel? Defaults to `FALSE`.
`numCores` | If `parallel = TRUE`, how many cores should be used for computation? Internally defaults to the number of available cores minus 2.
`...` | Dots for additional internal arguments (currently unused).

### Details

This function takes in a list of the first principal components from each pathway and an object of class `OmicsReg`. This function will then calculate the AIC of a multivariate linear model (via the `lm` function) with the original observations as response and the pathway principal components as the predictor matrix. Note that the AIC and log-likelihood are proportional because the number of parameters in each pathway is constant.

Then, this function will create `numReps` permutations of the regression response, fit models to each of these permuted responses (holding the path predictor matrix fixed), and calculate the AIC of each model. This function will return a named vector of permutation p-values, where the value for each pathway is the proportion of models for which the AIC of the permuted response model is less than the AIC of the original model.

### Value

A named vector of pathway permutation p-values.

### See Also

`CreateOmics`; `ExtractAESPCs`; `lm`; `SampleReg`

### Examples

      # DO NOT CALL THIS FUNCTION DIRECTLY.
      # Use AESPCA_pVals() instead

    ## Not run:
      ###  Load the Example Data  ###
      data("colonSurv_df")
      data("colon_pathwayCollection")

      ###  Create an OmicsSurv Object  ###
      colon_Omics <- CreateOmics(
        assayData_df = colonSurv_df[, -(2:3)],
        pathwayCollection_ls = colon_pathwayCollection,
        response = colonSurv_df[, 1:2],
        respType = "reg"
      )

      ###  Extract Pathway PCs and Loadings  ###
      colonPCs_ls <- ExtractAESPCs(
        object = colon_Omics,
        parallel = TRUE,
        numCores = 2
      )

      ###  Pathway p-Values  ###
      PermTestReg(
        OmicsReg = colon_Omics,
        pathwayPCs_ls = colonPCs_ls$PCs,
        parallel = TRUE,
        numCores = 2
      )

    ## End(Not run)

[Package pathwayPCA version 1.6.3 Index]