FS | R Documentation |
Class for storing RA data and associated functions for analysis of full-sib family populations.
## Create FS object
FSobj <- makeFS(RAobj, pedfile, family=NULL, inferSNPs = FALSE,
filter=list(MAF=0.05, MISS=0.2, BIN=100, DEPTH=5, PVALUE=0.01, MAXDEPTH=500))
## Functions (Methods) of an FS object
FSobj$addBIsnps(LODthres = 10, nComp = 10)
FSobj$addSNPs(LODthres = 10, nComp = 10)
FSobj$computeMap(chrom = NULL, init_r = 0.01, ep = 0.001, method = "optim", sexSpec = F, err = T, mapped = T, nThreads = 1)
FSobj$createLG(parent = "both", LODthres = 10, nComp = 10)
FSobj$maskSNP(snps)
FSobj$mergeLG(LG, where = NULL, mergeTo = NULL)
FSobj$orderLG(chrom = NULL, mapfun = "haldane", weight = "LOD2", ndim = 30, spar = NULL)
FSobj$plotChr(parent = "maternal", mat = "rf", filename = NULL, chrS = 2, lmai = 2)
FSobj$plotLG(parent = "maternal", LG = NULL, mat = "rf", filename = NULL, interactive = TRUE, what = NULL)
FSobj$plotLM(LG = NULL, fun = "haldane", col = "black")
FSobj$plotSyn()
FSobj$print(what = NULL, ...)
FSobj$removeLG(LG, where = NULL)
FSobj$removeSNP(snps, where = NULL)
FSobj$rf_2pt(nClust = 2, err = FALSE)
FSobj$unmaskSNP(snps)
FSobj$writeLM(file, direct = "./", LG = NULL, what = NULL, inferGeno=TRUE)
An FS object is created from the makeFS
function and contains RA data,
various statistics of the dataset that have been computed, and functions (methods)
for analyzing the data. Information in an FS object are specific to full-sib family populations.
$addBIsnps
Add Both-Informative (BI) snps to the maternal and/or paternal linkage groups and merge the parental linkage groups.
$addSNPs
Add MI, PI and/or SI SNPs to the existing linkage groups
$createLG
Create linkage group(s).
$computeMap
Compute linkage maps for a given marker order.
$maskSNP
Mask SNP(s) in the dataset.
$mergeLG
Merge linkage groups.
$orderLG
Order the SNPs in the linkage group(s).
$plotChr
Plot the heatmap of the 2-point recombination fraction estimates (or LOD scores) when SNPs are ordered according to the genome assembly.
$plotLG
Plot the heatmap of the 2-point recombination fraction estimates (or LOD scores) when SNPs are ordered according to their linkage groups.
$plotLM
Plot linkage maps.
$plotSyn
Produce a Synteny plot.
$print
Print summary information for FS object.
$removeLG
Remove linkage group(s).
$removeSNP
Remove SNP(s) from linkage group(s).
$rf_2pt
Compute the 2-point recombination fraction (and LOD score) between all SNP pairs.
$unmaskSNP
Unmask SNP(s) in the data set.
$writeLM
Write linkage mapping results to a file.
GUSbase::RA
-> FS
new()
FS$new(R6obj)
print()
FS$print(what = NULL, ...)
removeSNP()
FS$removeSNP(snps, where = NULL)
removeLG()
FS$removeLG(LG, where = NULL)
mergeLG()
FS$mergeLG(LG, where = NULL, mergeTo = NULL)
maskSNP()
FS$maskSNP(snps)
unmaskSNP()
FS$unmaskSNP(snps)
rf_2pt()
FS$rf_2pt(nClust = 2, err = FALSE)
createLG()
FS$createLG(parent = "both", LODthres = 10, nComp = 10, reset = FALSE)
addSNPs()
FS$addSNPs(LODthres = 10, nComp = 10)
addBIsnps()
FS$addBIsnps(LODthres = 10, nComp = 10, binsize = 0)
orderLG()
FS$orderLG( LG = NULL, mapfun = "haldane", weight = "LOD2", binsize = 0, ndim = 30, spar = NULL, filename = NULL )
plotLG()
FS$plotLG( parent = "maternal", LG = NULL, mat = "rf", filename = NULL, interactive = TRUE, what = NULL, ... )
plotChr()
FS$plotChr( chrom = NULL, parent = "maternal", mat = c("rf"), filename = NULL, chrS = 2, lmai = 2 )
plotLM()
FS$plotLM(LG = NULL, fun = "haldane", col = "black")
plotSyn()
FS$plotSyn()
computeMap()
FS$computeMap( chrom = NULL, init_r = 0.001, ep = 0.001, method = NULL, sexSpec = FALSE, err = TRUE, multiErr = FALSE, mapped = TRUE, nThreads = 1, inferOPGP = TRUE, rfthres = 0.1 )
writeLM()
FS$writeLM(file, direct = "./", LG = NULL, what = NULL, inferGeno = TRUE)
cometPlot()
FS$cometPlot( filename = NULL, cex = 1, maxdepth = 500, maxSNPs = 1e+05, res = 300, color = NULL, ind = FALSE, ncores = 1, ... )
rocketPlot()
FS$rocketPlot( ploid = 2, filename = NULL, cex = 1, maxdepth = 500, maxSNPs = 1e+05, res = 300, scaled = TRUE, ... )
RDDPlot()
FS$RDDPlot(filename = NULL, maxdepth = 500, maxSNPs = 1e+05, ...)
clone()
The objects of this class are cloneable with this method.
FS$clone(deep = FALSE)
deep
Whether to make a deep clone.
Timothy P. Bilton
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.