createBindingFactor.layer_region: Create a binding factor object to match a given pattern of...

In davetgerrard/GenomicLayers: Simple, sequence-based simulation of epi-genomes.

Create a binding factor object to match a given pattern of layers, typically ignoring DNA

Description

Create a new binding factor based on a simple pattern of marks on others layers and may (when used) set marks on other layers.

Usage

createBindingFactor.layer_region(
  name,
  type = "layer_region",
  patternLength = 1,
  patternString = "N",
  mismatch.rate = 0,
  stateWidth = patternLength,
  profile.layers = "LAYER.1",
  profile.marks = 0,
  mod.layers = "LAYER.1",
  mod.marks = 1,
  offset = 0,
  offset.method = NULL,
  align = "centre",
  test.layer0.binding = FALSE,
  test.mismatch.rate = 0.1,
  max.pattern.tries = 1000,
  min.DM.length = 2,
  min.DR.length = 10,
  verbose = FALSE
)

Arguments

name	give the binding factor a name
type	"layer_region" to differentiate from other types
patternLength	[= 1] length of pattern
patternString	= NOT USED in this case
stateWidth	the width of pattern to recognise on other layers, default is same as patternLength
profile.layers	a vector of named layers to set as a match
profile.marks	a vector of 0/1 to match the layers in profile.layers
mod.layers	a vector of named layers to alter on a match
mod.marks	a vector of 0/1 to set on the mod.layers
offset	0 integer value to indicate relative distance from pattern to apply modifications
offset.method	NULL a function to apply to apply offset.
align	"centre"
test.layer0.binding	when creating, test if the DNA sequence has a match.
test.mismatch.rate	proportion of mismatches to tolerate when testing
max.pattern.tries	NA
min.DM.length	NA
min.DR.length	NA
verbose	set to TRUE for more output

Value

"hits"

Examples

bf.LR <- createBindingFactor.layer_region("layerBf", type="layer_region", 
         patternLength = 1, patternString = "N",
          stateWidth = patternLength, profile.layers = "LAYER.1",
          profile.marks = 0, mod.layers = "LAYER.1", mod.marks = 1)

davetgerrard/GenomicLayers documentation built on April 28, 2024, 2:53 p.m.