simulateIBDsegmentsFabia
: R implementation of simulateIBDsegmentsFabia
.
Genotype data is simulated which contains rare variants and implanted IBD segments. Output is written for the bicluster algorithm fabia.
simulateIBDsegmentsFabia(fileprefix="dataSim", minruns=1,maxruns=1,snvs=1000,individualsN=100, avDistSnvs=100,avDistMinor=10,noImplanted=1, implanted=10,length=50,minors=30,mismatches=0, mismatchImplanted=0.5,overlap=50)
1/avDistMinor
is the average minor allele frequency (MAF).Data simulations for fabia focuses on rare variants but common variants are possible. Linkage disequilibrium and haplotype blocks are not simulated except by implanting IBD segments.
Simulated data is written to "...fabia_individuals.txt", "...fabia_annot.txt" "...fabia_mat.txt".
Implementation in R.
S. Hochreiter et al., FABIA: Factor Analysis for Bicluster Acquisition, Bioinformatics 26(12):1520-1527, 2010.
IBDsegment-class
,
IBDsegmentList-class
,
analyzeIBDsegments
,
compareIBDsegmentLists
,
extractIBDsegments
,
findDenseRegions
,
hapFabia
,
hapFabiaVersion
,
hapRes
,
chr1ASW1000G
,
IBDsegmentList2excel
,
identifyDuplicates
,
iterateIntervals
,
makePipelineFile
,
matrixPlot
,
mergeIBDsegmentLists
,
mergedIBDsegmentList
,
plotIBDsegment
,
res
,
setAnnotation
,
setStatistics
,
sim
,
simu
,
simulateIBDsegmentsFabia
,
simulateIBDsegments
,
split_sparse_matrix
,
toolsFactorizationClass
,
vcftoFABIA
old_dir <- getwd()
setwd(tempdir())
simulateIBDsegmentsFabia()
setwd(old_dir)
Run the code above in your browser using DataLab