search for: readgeneticmap

Displaying 2 results from an estimated 2 matches for "readgeneticmap".

2020 Oct 29
1
R: sim1000G
...vcf.gz" #nvariants = 442, ss=1000 vcf = readVCF( vcf_file, maxNumberOfVariants = 442 ,min_maf = 0.0005,max_maf = 0.01) #lowest MAF dim( vcf$gt1 ) #rows represent number of variants, columns represent number of individuals ## Download and use full chromosome genetic map downloadGeneticMap(4) readGeneticMap(4) sample.size=3000 startSimulation(vcf, totalNumberOfIndividuals = sample.size) data_sim = function(seed.num){ SIM$reset() id = generateUnrelatedIndividuals(sample.size) gt = retrieveGenotypes(id) freq = apply(gt,2,sum)/(2*nrow(gt)) causal = sample(setdiff(1:ncol(gt),which(freq==...
2020 Nov 01
0
R: sim1000G
...> vcf = readVCF( vcf_file, maxNumberOfVariants = 442 ,min_maf = > 0.0005,max_maf = 0.01) #lowest MAF > dim( vcf$gt1 ) #rows represent number of variants, columns represent > number of individuals > > ## Download and use full chromosome genetic map > downloadGeneticMap(4) > readGeneticMap(4) > > sample.size=3000 > > startSimulation(vcf, totalNumberOfIndividuals = sample.size) > > data_sim = function(seed.num){ > > SIM$reset() > > id = generateUnrelatedIndividuals(sample.size) > > gt = retrieveGenotypes(id) > > > freq = apply(gt,2,...