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],...