search for: retrievegenotyp

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

Did you mean: retrievegenotypes
2020 Oct 29
1
R: sim1000G
...ent 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==0)),45) beta.sign = rep(1,45) c.value = 0.402 beta.abs = c.value*abs(log10(freq[causal])) beta.val = beta.sign*beta.abs x.bar = apply(gt[,causal],2,mean) x.bar = as.matrix(x.bar) beta.val =...
2020 Nov 01
0
R: sim1000G
...enetic 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==0)),45) > > beta.sign = rep(1,45) > c.value = 0.402 > beta.abs = c.value*abs(log10(freq[causal])) > beta.val = beta.sign*beta.abs > x.bar = apply(gt[,causal],...