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

whys lambda is so large #212

Open
koujiaodahan opened this issue Dec 13, 2023 · 2 comments
Open

whys lambda is so large #212

koujiaodahan opened this issue Dec 13, 2023 · 2 comments

Comments

@koujiaodahan
Copy link

I try to run fitted_rss1 <- susie_rss(bhat = sumread$beta_alt, shat = sumread$V10, n = N, R = ldmatrix, L = 10,
estimate_residual_variance = FALSE,verbose=TRUE)
based on my gwas summary statistics and the LD matrix of same samples.
After that,I run kriging_rss for diagnose. but the lambda is always above 0.4
and the plot is that. why?

image
@koujiaodahan
Copy link
Author

Besides that , running the susie_rss function at a locus(chr6), it breaks with "Error in susie_suff_stat(XtX = XtX, Xty = Xty, n = n, yty = (n - 1) * :
The estimated prior variance is unreasonably large.
This is usually caused by mismatch between the summary statistics and the LD matrix.
Please check the input.
Calls: susie_rss -> susie_suff_stat
Execution halted"

@stephens999
Copy link
Contributor

The results suggest something wrong with the way you are computing your ld matrix or z scores

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

2 participants