Open ylb613 opened 6 months ago
Dear authors, My program has a error when running this step. What is the variable logOR, which was not mentioned in the previous step
K_IPe_ij <- apply(logOR[logOR.id,], 2, function(x){
fit <- lm(y~m, data = data.frame(y = x, m=rowMeans(logOR)[logOR.id] )) y.est <- predict(fit, newdata = data.frame(m = rowMeans(logOR))) return( y.est - rowMeans(logOR) ) })
Best wishes Yang
Dear authors, My program has a error when running this step. What is the variable logOR, which was not mentioned in the previous step
calculate the offset term for IP efficiency
K_IPe_ij <- apply(logOR[logOR.id,], 2, function(x){
fit <- lm(y~m, data = data.frame(y = x, m=rowMeans(logOR)[logOR.id] )) y.est <- predict(fit, newdata = data.frame(m = rowMeans(logOR))) return( y.est - rowMeans(logOR) ) })
Best wishes Yang