count_within_many: Count cells within a radius for multiple tissue categories...

Description Usage Arguments Details Value See Also Examples

Description

This is a wrapper around count_within() which supports counting multiple phenotype pairs and tissue categories within a single field. For each given tissue category, pair of 'from' phenotype and 'to' phenotype, and radius, it counts the number of 'from' cells having a 'to' cell within radius microns.

Usage

1
2
3
4
5
6
7
8
count_within_many(
  csd,
  pairs,
  radius,
  category = NA,
  phenotype_rules = NULL,
  verbose = TRUE
)

Arguments

csd

A cell seg data table.

pairs

A list of pairs of phenotypes. Each entry is a two-element vector. The result will contain values for each pair.

radius

The radius or radii to search within.

category

Optional tissue categories to restrict both from and to phenotypes.

phenotype_rules

(Optional) A named list. Item names are phenotype names and must match entries in pairs. Item values are selectors for select_rows().

verbose

If TRUE, display progress.

Details

The category parameter may be a single category or a list of categories.

See the tutorial Selecting cells within a cell segmentation table for more on the use of pairs and phenotype_rules.

Value

A tibble containing these columns:

slide_id

Slide ID from the data, if available.

source

Source field name.

field

Name of the individual field, if available.

category

Tissue category, if provided as a parameter, or "all".

from

From phenotype.

to

To phenotype.

radius, from_count, to_count, from_with, within_mean

Results from count_within for this data file and tissue category.

See Also

Other distance functions: compute_all_nearest_distance(), count_touching_cells(), count_within_batch(), count_within(), distance_matrix(), find_nearest_distance(), spatial_distribution_report(), subset_distance_matrix()

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
csd <- sample_cell_seg_data

# Count tumor cells near macrophages, and tumor cells near CD8 separately,
# in tumor and stroma tissue categories separately.
pairs <- list(c('CK+', 'CD68+'),
             c('CK+', 'CD8+'))
radius <- c(10, 25)
category <- list('Tumor', 'Stroma')
count_within_many(csd, pairs, radius, category)

# Count tumor cells near any T cell in all tissue categories.
# Use `phenotype_rules` to define the T cell phenotype
pairs <- c('CK+', 'T cell')
rules <- list(
'T cell'=c('CD8+', 'FoxP3+'))
count_within_many(csd, pairs, radius, phenotype_rules=rules)

akoyabio/phenoptr documentation built on Jan. 7, 2022, 5:37 p.m.