Runsigfit: Run sigfit extraction and attribution on a spectra catalog...

View source: R/Runsigfit.R

RunsigfitR Documentation

Run sigfit extraction and attribution on a spectra catalog file

Description

WARNING: sigfit can only do exposure attribution using SBS96 spectra catalog and signature catalog!

Usage

Runsigfit(
  input.catalog,
  out.dir,
  model = "nmf",
  CPU.cores = NULL,
  seedNumber = 1,
  K.exact = NULL,
  K.range = NULL,
  test.only = FALSE,
  overwrite = FALSE
)

Arguments

input.catalog

File containing input spectra catalog. Columns are samples (tumors), rows are mutation types.

out.dir

Directory that will be created for the output; abort if it already exits. Log files will be in paste0(out.dir, "/tmp").

model

Algorithm to be used to extract signatures and attribute exposures. Only "nmf" or "emu" is valid. Default: "nmf".

CPU.cores

Number of CPUs to use in running sigfit. For a server, 30 cores would be a good choice; while for a PC, you may only choose 2-4 cores. By default (CPU.cores = NULL), the CPU.cores would be equal to (parallel::detectCores())/2, total number of CPUs divided by 2.

seedNumber

Specify the pseudo-random seed number used to run sigfit. Setting seed can make the attribution of sigfit repeatable. Default: 1.

K.exact, K.range

K.exact is the exact value for the number of signatures active in spectra (K). Specify K.exact if you know exactly how many signatures are active in the input.catalog, which is the ICAMS-formatted spectra file.

K.range is A numeric vector (K.min,K.max) of length 2 which tell sigfit to search the best signature number active in spectra, K, in this range of Ks. Specify K.range if you don't know how many signatures are active in the input.catalog. K.max - K.min >= 3, otherwise an error will be thrown.

WARNING: You must specify only one of K.exact or K.range!

Default: NULL

test.only

If TRUE, only analyze the first 10 columns read in from input.catalog. Default: FALSE

overwrite

If TRUE, overwrite existing output. Default: FALSE

Details

Creates several files in out.dir. These are: TODO(Steve): list the files

TODO(Wuyang)

Value

The inferred exposure of sigfit, invisibly.


WuyangFF95/SynSigRun documentation built on Oct. 7, 2022, 1:16 p.m.