Last updated: 2018-12-14
workflowr checks: (Click a bullet for more information) ✔ R Markdown file: up-to-date
Great! Since the R Markdown file has been committed to the Git repository, you know the exact version of the code that produced these results.
✔ Environment: empty
Great job! The global environment was empty. Objects defined in the global environment can affect the analysis in your R Markdown file in unknown ways. For reproduciblity it’s best to always run the code in an empty environment.
✔ Seed:
set.seed(12345)
The command set.seed(12345) was run prior to running the code in the R Markdown file. Setting a seed ensures that any results that rely on randomness, e.g. subsampling or permutations, are reproducible.
✔ Session information: recorded
Great job! Recording the operating system, R version, and package versions is critical for reproducibility.
✔ Repository version: 73ed3b7
wflow_publish or wflow_git_commit). workflowr only checks the R Markdown file, but you know if there are other scripts or data files that it depends on. Below is the status of the Git repository when the results were generated:
Ignored files:
Ignored: .DS_Store
Ignored: .Rhistory
Ignored: .Rproj.user/
Ignored: analysis/.DS_Store
Ignored: analysis/BH_robustness_cache/
Ignored: analysis/FDR_Null_cache/
Ignored: analysis/FDR_null_betahat_cache/
Ignored: analysis/Rmosek_cache/
Ignored: analysis/StepDown_cache/
Ignored: analysis/alternative2_cache/
Ignored: analysis/alternative_cache/
Ignored: analysis/ash_gd_cache/
Ignored: analysis/average_cor_gtex_2_cache/
Ignored: analysis/average_cor_gtex_cache/
Ignored: analysis/brca_cache/
Ignored: analysis/cash_deconv_cache/
Ignored: analysis/cash_fdr_1_cache/
Ignored: analysis/cash_fdr_2_cache/
Ignored: analysis/cash_fdr_3_cache/
Ignored: analysis/cash_fdr_4_cache/
Ignored: analysis/cash_fdr_5_cache/
Ignored: analysis/cash_fdr_6_cache/
Ignored: analysis/cash_paper_fig1_cache/
Ignored: analysis/cash_paper_fig_avgecdf_cache/
Ignored: analysis/cash_paper_fig_deconv_cache/
Ignored: analysis/cash_paper_fig_leukemia_cache/
Ignored: analysis/cash_plots_2_cache/
Ignored: analysis/cash_plots_3_cache/
Ignored: analysis/cash_plots_4_cache/
Ignored: analysis/cash_plots_5_cache/
Ignored: analysis/cash_plots_cache/
Ignored: analysis/cash_sim_1_cache/
Ignored: analysis/cash_sim_2_cache/
Ignored: analysis/cash_sim_3_cache/
Ignored: analysis/cash_sim_4_cache/
Ignored: analysis/cash_sim_5_cache/
Ignored: analysis/cash_sim_6_cache/
Ignored: analysis/cash_sim_7_cache/
Ignored: analysis/correlated_z_2_cache/
Ignored: analysis/correlated_z_3_cache/
Ignored: analysis/correlated_z_cache/
Ignored: analysis/create_null_cache/
Ignored: analysis/cutoff_null_cache/
Ignored: analysis/design_matrix_2_cache/
Ignored: analysis/design_matrix_cache/
Ignored: analysis/diagnostic_ash_cache/
Ignored: analysis/diagnostic_correlated_z_2_cache/
Ignored: analysis/diagnostic_correlated_z_3_cache/
Ignored: analysis/diagnostic_correlated_z_cache/
Ignored: analysis/diagnostic_plot_2_cache/
Ignored: analysis/diagnostic_plot_cache/
Ignored: analysis/efron_leukemia_cache/
Ignored: analysis/fitting_normal_cache/
Ignored: analysis/gaussian_derivatives_2_cache/
Ignored: analysis/gaussian_derivatives_3_cache/
Ignored: analysis/gaussian_derivatives_4_cache/
Ignored: analysis/gaussian_derivatives_5_cache/
Ignored: analysis/gaussian_derivatives_cache/
Ignored: analysis/gd-ash_cache/
Ignored: analysis/gd_delta_cache/
Ignored: analysis/gd_lik_2_cache/
Ignored: analysis/gd_lik_cache/
Ignored: analysis/gd_w_cache/
Ignored: analysis/knockoff_10_cache/
Ignored: analysis/knockoff_2_cache/
Ignored: analysis/knockoff_3_cache/
Ignored: analysis/knockoff_4_cache/
Ignored: analysis/knockoff_5_cache/
Ignored: analysis/knockoff_6_cache/
Ignored: analysis/knockoff_7_cache/
Ignored: analysis/knockoff_8_cache/
Ignored: analysis/knockoff_9_cache/
Ignored: analysis/knockoff_cache/
Ignored: analysis/knockoff_var_cache/
Ignored: analysis/marginal_z_alternative_cache/
Ignored: analysis/marginal_z_cache/
Ignored: analysis/mosek_reg_2_cache/
Ignored: analysis/mosek_reg_4_cache/
Ignored: analysis/mosek_reg_5_cache/
Ignored: analysis/mosek_reg_6_cache/
Ignored: analysis/mosek_reg_cache/
Ignored: analysis/pihat0_null_cache/
Ignored: analysis/plot_diagnostic_cache/
Ignored: analysis/poster_obayes17_cache/
Ignored: analysis/real_data_simulation_2_cache/
Ignored: analysis/real_data_simulation_3_cache/
Ignored: analysis/real_data_simulation_4_cache/
Ignored: analysis/real_data_simulation_5_cache/
Ignored: analysis/real_data_simulation_cache/
Ignored: analysis/rmosek_primal_dual_2_cache/
Ignored: analysis/rmosek_primal_dual_cache/
Ignored: analysis/seqgendiff_cache/
Ignored: analysis/simulated_correlated_null_2_cache/
Ignored: analysis/simulated_correlated_null_3_cache/
Ignored: analysis/simulated_correlated_null_cache/
Ignored: analysis/simulation_real_se_2_cache/
Ignored: analysis/simulation_real_se_cache/
Ignored: analysis/smemo_2_cache/
Ignored: data/LSI/
Ignored: docs/.DS_Store
Ignored: docs/figure/.DS_Store
Ignored: output/fig/
Ignored: output/paper/
Unstaged changes:
Modified: analysis/cash_plots_2.rmd
Plotting the standard Gaussian and its standardized derivatives.
source("../code/gdash_lik.R")Loading required package: EQLLoading required package: ttutilsLoading required package: SQUAREMLoading required package: REBayesLoading required package: MatrixWarning: package 'Matrix' was built under R version 3.4.4Loading required package: CVXR
Attaching package: 'CVXR'The following object is masked from 'package:stats':
powerLoading required package: PolynomFLoading required package: RmosekLoading required package: ashrplotx <- seq(-10, 10, length = 1000)
gd.ord <- 10
hermite = Hermite(gd.ord)
gd0.std = dnorm(plotx)
matrix_lik_plot = cbind(gd0.std)
for (i in 1 : gd.ord) {
gd.std = (-1)^i * hermite[[i]](plotx) * gd0.std / sqrt(factorial(i))
matrix_lik_plot = cbind(matrix_lik_plot, gd.std)
}
method.col <- scales::hue_pal()(gd.ord)setEPS()
postscript("../output/paper/GD.eps", height = 4.5, width = 15)
#pdf("../output/paper/GD.pdf", height = 4.5, width = 15)
par(mfrow = c(1, 2))
par(mar = c(3, 6.5, 1, 1.5))
plot(plotx, gd0.std, type = "l", xlim = c(-5, 5), ylim = range(matrix_lik_plot), xlab = "", ylab = '', main = '', bty = 'n', lwd = 2)
title(ylab = latex2exp::TeX('$\\frac{1}{\\sqrt{\\textit{l}!}}\\varphi^{(\\textit{l})}(\\cdot)$'), line = 2.5, cex.lab = 1.5)
for (i in 1 : 4) {
lines(plotx, matrix_lik_plot[, (i + 1)], col = method.col[i], lty = i + 1, lwd = 2)
}
legend("topright", col = c(1, method.col[1 : 4]), lty = 1 : 5, lwd = 2, bty = "n", cex = 1.5,
c(latex2exp::TeX('$\\varphi$'),
latex2exp::TeX('$\\textit{l} = 1$'),
latex2exp::TeX('$\\textit{l} = 2$'),
latex2exp::TeX('$\\textit{l} = 3$'),
latex2exp::TeX('$\\textit{l} = 4$')
))
plot(plotx, gd0.std, type = "l", xlim = c(-5, 5), ylim = range(matrix_lik_plot), xlab = "", ylab = '', main = '', bty = 'n', lwd = 2)
title(ylab = latex2exp::TeX('$\\frac{1}{\\sqrt{\\textit{l}!}}\\varphi^{(\\textit{l})}(\\cdot)$'), line = 2.5, font = 1, cex.lab = 1.5)
for (i in 7 : 10) {
lines(plotx, matrix_lik_plot[, (i + 1)], col = method.col[i], lty = i + 1, lwd = 2)
}
legend("topright", col = c(1, method.col[7 : 10]), lty = c(1, 8 : 11), lwd = 2, bty = "n", cex = 1.5,
c(latex2exp::TeX('$\\varphi$'),
latex2exp::TeX('$\\textit{l} = 7$'),
latex2exp::TeX('$\\textit{l} = 8$'),
latex2exp::TeX('$\\textit{l} = 9$'),
latex2exp::TeX('$\\textit{l} = 10$')
))
dev.off()quartz_off_screen
2 sessionInfo()R version 3.4.3 (2017-11-30)
Platform: x86_64-apple-darwin15.6.0 (64-bit)
Running under: macOS 10.14.1
Matrix products: default
BLAS: /Library/Frameworks/R.framework/Versions/3.4/Resources/lib/libRblas.0.dylib
LAPACK: /Library/Frameworks/R.framework/Versions/3.4/Resources/lib/libRlapack.dylib
locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] ashr_2.2-7 Rmosek_8.0.69 PolynomF_1.0-2 CVXR_0.95
[5] REBayes_1.3 Matrix_1.2-14 SQUAREM_2017.10-1 EQL_1.0-0
[9] ttutils_1.0-1
loaded via a namespace (and not attached):
[1] latex2exp_0.4.0 gmp_0.5-13.1 Rcpp_0.12.16
[4] plyr_1.8.4 compiler_3.4.3 git2r_0.23.0
[7] workflowr_1.1.1 R.methodsS3_1.7.1 R.utils_2.6.0
[10] iterators_1.0.9 tools_3.4.3 digest_0.6.15
[13] bit_1.1-13 evaluate_0.10.1 lattice_0.20-35
[16] foreach_1.4.4 yaml_2.1.19 parallel_3.4.3
[19] Rmpfr_0.7-0 ECOSolveR_0.4 stringr_1.3.1
[22] knitr_1.20 rprojroot_1.3-2 bit64_0.9-7
[25] grid_3.4.3 R6_2.2.2 rmarkdown_1.9
[28] magrittr_1.5 whisker_0.3-2 scales_0.5.0
[31] MASS_7.3-50 backports_1.1.2 codetools_0.2-15
[34] htmltools_0.3.6 scs_1.1-1 colorspace_1.3-2
[37] stringi_1.2.2 munsell_0.4.3 pscl_1.5.2
[40] doParallel_1.0.11 truncnorm_1.0-8 R.oo_1.22.0 This reproducible R Markdown analysis was created with workflowr 1.1.1