choishingwan / PRS-Tutorial

A tutorial on how to run basic polygenic risk score analysis
MIT License
71 stars 113 forks source link

'ind.col' and 'rows_along(A.col)' should have the same length. #30

Open Sandyyy123 opened 2 years ago

Sandyyy123 commented 2 years ago
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
  }
}`
Sandyyy123 commented 2 years ago

Tried beta_grid[ind.chr], but that also didnt work.