This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to find the SNPs used in final best.pgs in lassosum?

I'm working on PRS construction using lassosum package.

Once I obtain the best lambda value from the validation step, how can I see the final number of SNPs used in best.pgs in the following step?

out2 <- subset(out, s=v$best.s, lambda=v$best.lambda)
v2 <- validate(out2, covar=covar, test.bfile="Some_new_bfile")

My understanding is that we start with a certain number of SNPs from the summary stat based on p value cut off and then we select the SNPs with non-zero estimates.

TIA

lassoum lasso snps

0 answers

No answers yet.

Log in to answer this question.