search for: prvalance

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

2020 Oct 29
1
R: sim1000G
...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 = t(as.matrix(beta.val)) #disease prvalance = 1% #beta0 = -log(99)-beta.val %*% x.bar #disease prvalance = 1.5% beta0 = 0-beta.val %*% x.bar eta = beta.val %*% t(gt[,causal]) eta = as.vector(eta) + rep(beta0,nrow(gt)) prob = exp(eta)/(1+exp(eta)) genocase = rep(NA, sample.size) set.seed(seed.num) for(i in 1:sample.size){...
2020 Nov 01
0
R: sim1000G
...(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 = t(as.matrix(beta.val)) > #disease prvalance = 1% > #beta0 = -log(99)-beta.val %*% x.bar > #disease prvalance = 1.5% > beta0 = 0-beta.val %*% x.bar > > eta = beta.val %*% t(gt[,causal]) > eta = as.vector(eta) + rep(beta0,nrow(gt)) > prob = exp(eta)/(1+exp(eta)) > > genocase = rep(NA, sample.size) >...