Error: Tested 13641 < 13641. Subscript out of bounds.

privefl / bigsnpr

R package for the analysis of massive SNP arrays.

https://privefl.github.io/bigsnpr/

183 stars 43 forks source link

Error: Tested 13641 < 13641. Subscript out of bounds. #463

Closed DUAN-GAO closed 2 months ago

DUAN-GAO commented 7 months ago

Hello, encountered following bug: Error: Tested 13641 < 13641. Subscript out of bounds. could u help, many thanks.

privefl commented 7 months ago

Which function? Can you provide a bit more detail?

DUAN-GAO commented 7 months ago

pred_auto <- big_prodVec(G_iPSCs_002, beta_auto * map_pgs2_iPSCs_002$beta,

ind.col = df_beta[["_NUMID"]],

ncores = 16)

Error: Tested 13641 < 13641. Subscript out of bounds. In addition: Warning message: In beta_auto * map_pgs2_iPSCs_002$beta : longer object length is not a multiple of shorter object length

It went all fine until this step.

DUAN-GAO commented 7 months ago

If change the ind.col parameter to the same object as beta_auto*

pred_auto <- big_prodVec(G_iPSCs_002, beta_auto * map_pgs2_iPSCs_002$beta,

ind.col = map_pgs2_iPSCs_002[["_NUMID"]],

ncores = 16)

then got this error instead

Error: Incompatibility between dimensions. 'y.col' and 'ind.col' should have the same length. In addition: Warning message: In beta_auto * map_pgs2_iPSCs_002$beta : longer object length is not a multiple of shorter object length

privefl commented 7 months ago

Two issues here I guess:

beta_auto and map_pgs2_iPSCs_002$beta do not have compatible dimensions
ind.col = df_beta[["_NUM_ID_"]] should be the column indices used in G_iPSCs_002; this means that max(df_beta[["_NUM_ID_"]]) should be lower or equal to ncol(G_iPSCs_002)

privefl commented 5 months ago

Any update on this?