Displaying 2 results from an estimated 2 matches for "generateunrelatedindividu".
2020 Oct 29
1
R: sim1000G
...#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==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)...
2020 Nov 01
0
R: sim1000G
...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 = bet...