seqbias (version 1.20.0)

seqbias.load: Loading seqbias models

Description

Loads a seqbias model from a file written with 'seqbias.save'.

Usage

seqbias.load(ref_fn, model_fn)

Arguments

ref_fn

filename of a reference sequence against which the reads are aligned in FASTA format.

model_fn

filename of a saved seqbias model

Value

Details

A large seqbias model can take some time (several minutes) to fit. It is often preferable to do this just once. This function load the model from a file in YAML format, having been written with 'seqbias.save'.

See Also

seqbias.save

Examples

  reads_fn <- system.file( "extra/example.bam", package = "seqbias" )
  ref_fn <- system.file( "extra/example.fa", package = "seqbias" )

  sb <- seqbias.fit( ref_fn, reads_fn )

  seqbias.save( sb, "my_seqbias_model.yml" )

  # load sometime later
  sb <- seqbias.load( ref_fn, "my_seqbias_model.yml" )