inst/doc/Introduction_FS.R

In tpbilton/GUSMap: Genotyping Uncertainty with Sequencing Data and Linkage Mapping (GUSMap)

## ----setup, include = FALSE---------------------------------------------------
knitr::opts_chunk$set( collapse = TRUE, comment = "#>", eval=TRUE)
library(GUSMap)

## ----simDS, eval = T----------------------------------------------------------
vcffile <- simDS() # extract filename
vcffile            # filename stored in object vcffile

## ----loadData, eval=T---------------------------------------------------------
# convert VCF file to an RA file
rafile <- VCFtoRA(infilename = vcffile$vcf, direct = "./", makePed=TRUE)

## ----rafile, eval=T-----------------------------------------------------------
rafile # file path of RA file

## ----readRA, eval=T-----------------------------------------------------------
RAdata <- readRA(rafile = rafile, sampthres = 0.01, excsamp = NULL)

## ----class_RAdata, eval=T-----------------------------------------------------
class(RAdata)

## ----RAdata, eval=T-----------------------------------------------------------
RAdata

## ----makeFS, eval=T-----------------------------------------------------------
mySpecies <- makeFS(RAobj = RAdata, pedfile = vcffile$ped, inferSNPs = FALSE,
                      filter = list(MAF = 0.05, MISS = 0.5,
                        BIN = 100, DEPTH = 5, PVALUE = 0.01, MAXDEPTH=500))

## ----help_makeFS, eval=F------------------------------------------------------
#  ?`$makeFS`

## ----2pt_rf-------------------------------------------------------------------
mySpecies$rf_2pt(nClust = 3)

## ----save, eval=FALSE---------------------------------------------------------
#  save(mySpecies, file="FSobject.Rdata")

## ----plotChr------------------------------------------------------------------
# plot 2-point rf matrix for BI and MI SNPs
mySpecies$plotChr(mat = "rf", parent = "maternal", lmai=0.5)

## ----plotChr2-----------------------------------------------------------------
# plot 2-point LOD matrix for BI, PI and MI SNPs
mySpecies$plotChr(mat = "LOD", parent = "both", lmai=0.5)

## ----createLG-----------------------------------------------------------------
mySpecies$createLG(parent = "both", LODthres = 5, nComp = 10, reset = FALSE)

## ----help_createLG, eval=F----------------------------------------------------
#  ?`$createLG`

## ----plotLG_mat---------------------------------------------------------------
mySpecies$plotLG(parent = "both", interactive = F, LG = NULL)

## ----createLG_pat-------------------------------------------------------------
mySpecies$createLG(parent = "paternal", LODthres = 10, nComp = 10)

## ----plotLG_pat2--------------------------------------------------------------
mySpecies$plotLG(parent = "both", interactive = F)

## ----mergeLG------------------------------------------------------------------
mySpecies$mergeLG(LG = c(3,7), mergeTo="paternal")
mySpecies$print(what = "LG-pts")

## ----removeLG, eval=FALSE-----------------------------------------------------
#  mySpecies$removeLG(LG = 3)

## ----removeSNP, eval=FALSE----------------------------------------------------
#  mySpecies$removeSNP(snps = c(526,527))

## ----addSNPs, eval=FALSE------------------------------------------------------
#  mySpecies$addSNPs(LODthres = 10, nComp = 10)

## ----addBIsnps----------------------------------------------------------------
mySpecies$addBIsnps(LODthres = 10, nComp = 10)

## ----orderLG, , fig.width=8, fig.height=6-------------------------------------
mySpecies$orderLG(mapfun = "haldane", weight = "LOD2", ndim = 30) 

## ----help_orderLG, eval=F-----------------------------------------------------
#  ?`$orderLG`

## ----removeSNP_2--------------------------------------------------------------
mySpecies$removeSNP(c(
  c(253,255,300,316,143,301,340,166,173,201,209,264,312,310,315,336,
    356,64,500,459,513,449,423,376,451,508,418,228),
  c(377,405,217,177,436,450,467,28,515,411,229,222,410,230,401,368,
    375,34,259,164,169,334,283,213,148,288,213),
  c(112,108,311,37,69,136,46,115,29,126,42,72,26,92,67,71,94))
)

## ----computeMap---------------------------------------------------------------
mySpecies$computeMap(chrom = NULL, nThreads = 2)

## ----plotLG-------------------------------------------------------------------
mySpecies$plotLM(LG = NULL, fun = "haldane", col=c("cyan","lightblue","blue"))

## ----plotSyn, fig.width=7, fig.height=7---------------------------------------
mySpecies$plotSyn()

## ----writeLG, eval=F----------------------------------------------------------
#  mySpecies$writeLM(file = "mySpecies")

## ----help_writeLM, eval=F-----------------------------------------------------
#  ?`$writeLM`