Displaying 2 results from an estimated 2 matches for "prvalanc".
Did you mean:
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)
>...