Set motif data — SetMotifData • Signac

Set motif matrix for given assay

SetMotifData(object, ...)

# S3 method for class 'Motif'
SetMotifData(object, slot, new.data, ...)

# S3 method for class 'ChromatinAssay'
SetMotifData(object, slot, new.data, ...)

# S3 method for class 'Seurat'
SetMotifData(object, assay = NULL, ...)

Arguments

object: A Seurat object
...: Arguments passed to other methods
slot: Name of slot to use
new.data: motif matrix to add. Should be matrix or sparse matrix class
assay: Name of assay whose data should be set

Value

Returns a Seurat object

Examples

motif.obj <- SeuratObject::GetAssayData(
  object = atac_small[['peaks']], slot = "motifs"
)
SetMotifData(object = motif.obj, slot = 'data', new.data = matrix(1:2))
#> A Motif object containing 1 motifs in 2 regions
new.data <- matrix(sample(c(0, 1), size = nrow(atac_small[["peaks"]]) * 10,
                   replace = TRUE), nrow = nrow(atac_small[["peaks"]]))
rownames(new.data) <- rownames(atac_small[["peaks"]])
SetMotifData(
  object = atac_small[['peaks']], slot = 'data', new.data = new.data
)
#> ChromatinAssay data with 323 features for 100 cells
#> Variable features: 323 
#> Genome: hg19 
#> Annotation present: TRUE 
#> Motifs present: TRUE 
#> Fragment files: 0 
motif.matrix <- GetMotifData(object = atac_small)
SetMotifData(
object = atac_small, assay = 'peaks', slot = 'data', new.data = motif.matrix
)
#> An object of class Seurat 
#> 1323 features across 100 samples within 3 assays 
#> Active assay: peaks (323 features, 323 variable features)
#>  2 layers present: counts, data
#>  2 other assays present: bins, RNA
#>  2 dimensional reductions calculated: lsi, umap