| plot_weights {dearseq} | R Documentation |
Display the variability with respect to the level of expression and the associated smoothed estimation of precision weights accounting for heteroscedasticity.
plot_weights(x)
x |
a list (such as outputed by the functions
|
a ggplot object
Boris Hejblum
G <- 10000
n <- 12
p <- 2
y <- sapply(1:n, FUN = function(x){rnbinom(n = G, size = 0.07, mu = 200)})
x <- sapply(1:p, FUN = function(x){rnorm(n = n, mean = n, sd = 1)})
if(interactive()){
w <- sp_weights(y, x, use_phi=FALSE, na.rm = TRUE, gene_based=TRUE)
plot_weights(w)
vw <- voom_weights(y, x)
plot_weights(vw)
}