From ffe2698259933f64223e90b705ee78bf09d5910f Mon Sep 17 00:00:00 2001 From: Peter Carbonetto Date: Fri, 16 Dec 2022 09:03:11 -0600 Subject: [PATCH] Fixed another bug in kriging_rss. --- DESCRIPTION | 4 ++-- R/susie_utils.R | 4 ++-- 2 files changed, 4 insertions(+), 4 deletions(-) diff --git a/DESCRIPTION b/DESCRIPTION index 32d4cd5..922584a 100644 --- a/DESCRIPTION +++ b/DESCRIPTION @@ -15,8 +15,8 @@ Description: Implements methods for variable selection in linear called "Iterative Bayesian Stepwise Selection" (IBSS), is simple and fast, allowing the SuSiE model be fit to large data sets (thousands of samples and hundreds of thousands of variables). -Date: 2022-12-15 -Version: 0.12.30 +Date: 2022-12-16 +Version: 0.12.31 Authors@R: c(person("Gao","Wang",role="aut",email="wang.gao@columbia.edu"), person("Yuxin","Zou",role="aut"), person("Kaiqian","Zhang",role="aut"), diff --git a/R/susie_utils.R b/R/susie_utils.R index a75383e..eaaddfd 100644 --- a/R/susie_utils.R +++ b/R/susie_utils.R @@ -864,14 +864,14 @@ kriging_rss = function (z, R, n, r_tol = 1e-08, z_std_diff = z_std_diff, logLR = logLRmix) p = ggplot(res,aes_string(y = "z",x = "condmean")) + - geom_point(aes(y = z, x = condmean)) + + geom_point() + labs(y = "Observed z scores", x = "Expected value") + geom_abline(intercept = 0, slope = 1) + theme_bw() idx = which(logLRmix > 2 & abs(z) > 2) if (length(idx) > 0) p = p + geom_point(data = res[idx,], - aes(y = z, x = condmean),col = "red") + aes_string(y = "z", x = "condmean"),col = "red") return(list(plot = p,conditional_dist = res)) }