Skip to content

Commit

Permalink
Fixed another bug in kriging_rss.
Browse files Browse the repository at this point in the history
  • Loading branch information
pcarbo committed Dec 16, 2022
1 parent 459fa73 commit ffe2698
Show file tree
Hide file tree
Showing 2 changed files with 4 additions and 4 deletions.
4 changes: 2 additions & 2 deletions DESCRIPTION
Original file line number Diff line number Diff line change
Expand Up @@ -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="[email protected]"),
person("Yuxin","Zou",role="aut"),
person("Kaiqian","Zhang",role="aut"),
Expand Down
4 changes: 2 additions & 2 deletions R/susie_utils.R
Original file line number Diff line number Diff line change
Expand Up @@ -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))
}

Expand Down

0 comments on commit ffe2698

Please sign in to comment.