Closed Sandyyy123 closed 3 years ago
[pred_inf <- NULL for(chr in 1:22){ obj.bigSNP <- snp_attach(paste0("courage_eur.qc01_chr",chr,".rds")) genotype <- obj.bigSNP$genotypes ind.test <- 1:nrow(genotype) chr.idx <- which(info_snp$chr == chr) ind.chr <- info_snp$_NUM_ID_[chr.idx] tmp <- big_prodVec(genotype, beta_inf, ind.row = ind.test, ind.col = ind.chr) if(is.null(pred_inf)){ pred_inf <- tmp }else{ pred_inf <- pred_inf + tmp } } ]
_NUM_ID_
[pred_inf <- NULL for(chr in 1:22){ obj.bigSNP <- snp_attach(paste0("courage_eur.qc01_chr",chr,".rds")) genotype <- obj.bigSNP$genotypes ind.test <- 1:nrow(genotype) chr.idx <- which(info_snp$chr == chr) ind.chr <- info_snp$
_NUM_ID_[chr.idx] tmp <- big_prodVec(genotype, beta_inf, ind.row = ind.test, ind.col = ind.chr) if(is.null(pred_inf)){ pred_inf <- tmp }else{ pred_inf <- pred_inf + tmp } } ]