Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Reading PredictDB Weights | Don't hardcode rsid as the SNP identifier #7

Open
shugamoe opened this issue Oct 5, 2023 · 0 comments
Open

Comments

@shugamoe
Copy link
Contributor

shugamoe commented Oct 5, 2023

rownames(wgt.matrix) <- wgt$rsid

When reading in weights from a .db it can be inconvenient to have rsid hard coded as the SNP identifier. Especially since there is another column, varID, which could be used. While working with predictdb mashr gtex v8 eQTL and sQTL models, I've used that identifier instead. Some of the SNPs in these models, if rsid is missing for whatever reason, use varID as a backup in the rsid row anyway as well.

Here's an example of a db_id parameter added to allow the user to specify the SNP identifier from a .db weight source.
shugamoe@bd42fc0

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant