spatialFDR {cydar} | R Documentation

## Compute the spatial FDR

### Description

Computed adjusted p-values for all hyperspheres, using a density-weighted version of the Benjamini-Hochberg method.

### Usage

    spatialFDR(x, pvalues, neighbors=50, bandwidth=NULL)

### Arguments

`x` | A numeric matrix of hypersphere coordinates, containing the median intensity of each marker (column) in each hypersphere (row). Alternatively, a CyData object containing median intensities for groups of cells, such as that produced by `countCells`.
`pvalues` | A numeric vector of p-values for each hypersphere.
`neighbors` | An integer scalar specifying the number of neighbors with which to compute the bandwidth.
`bandwidth` | A numeric scalar specifying the bandwidth for density estimation.

### Details

Consider the set of significant hyperspheres, distributed in some manner across the M-dimensional space (for M markers). The aim is to control the FDR across the subspaces containing significant hyperspheres. This is subtly different from controlling the FDR across the hypersphere themselves, which will skew the results for densely occupied subspaces.

Control of the spatial FDR is achieved by weighting the hyperspheres inversely proportional to their local densities. This downweights hyperspheres in dense subspaces while upweighting hyperspheres in sparse subspaces. The computed weights are then used as frequency weights in the Benjamini-Hochberg method, to control the FDR across subspaces.

The local density is calculated using a tricube kernel and the specified `bandwidth`. If unspecified, `bandwidth` is set to the median of the distances to the `neighbors`-closest neighbor for all hyperspheres. This usually provides stable density estimates while maintaining sensitivity to fine-scale structure.

### Value

A numeric vector of adjusted p-values for all hyperspheres.

### Author(s)

Aaron Lun

### References

Lun ATL, Richard AC, Marioni JC (2017). Testing for differential abundance in mass cytometry data. Nat. Methods, 14, 7:707-709.

### Examples

    coords <- matrix(rgamma(10000, 2, 2), nrow=1000)
    pvalues <- rbeta(nrow(coords), 1, 2)
    out <- spatialFDR(coords, pvalues)

[Package cydar version 1.6.1 Index]