PRS-Tutorial
PRS-Tutorial copied to clipboard
'ind.col' and 'rows_along(A.col)' should have the same length.
pred_grid <- 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_prodMat( genotype,
beta_grid ,
ind.col = ind.chr)
if(is.null(pred_grid)){
pred_grid <- tmp
}else{
pred_grid <- pred_grid + tmp
}
}`