generateObject: Load data

In rfael0cm/RTIGER: HMM-Based Model for Genotyping and Cross-Over Identification

Load data

Description

Load data

Usage

generateObject(experimentDesign = NULL,nstates = 3, rigidity=NULL,
seqlengths = NULL, verbose = TRUE)

Arguments

experimentDesign	a data Frame that contains minimum a column with the files direction (name of the column files) and another with a shorter name to be used inside the function.
nstates	the number of states to be fitted in the model. A standard setting would use 3 states (Homozygous1, Heterozygous, and Homozygous2).
rigidity	an integer number specifying the rigidity parameter to be used.
seqlengths	a named vector with the chromosome lenghts of the organism that the user is working with.
verbose	logical value. Whether to print info messages.

Value

RTIGER object

Examples

data("ATseqlengths")
path = system.file("extdata",  package = "RTIGER")
files = list.files(path, full.names = TRUE)
nam = sapply(list.files(path ), function(x) unlist(strsplit(x, split = "[.]"))[1])
expDesign = data.frame(files = files, name = nam)
names(ATseqlengths) = paste0("Chr", 1:5)
myres = generateObject(experimentDesign = expDesign,
              seqlengths = ATseqlengths,
              rigidity = 10
)

rfael0cm/RTIGER documentation built on Oct. 31, 2023, 1:25 p.m.