makemyprior
Science Score: 13.0%
This score indicates how likely this project is to be science-related based on various indicators:
-
○CITATION.cff file
-
✓codemeta.json file
Found codemeta.json file -
○.zenodo.json file
-
○DOI references
-
○Academic publication links
-
○Academic email domains
-
○Institutional organization owner
-
○JOSS paper metadata
-
○Scientific vocabulary similarity
Low similarity (8.0%) to scientific vocabulary
Repository
Basic Info
- Host: GitHub
- Owner: ingebogh
- Language: R
- Default Branch: master
- Size: 940 KB
Statistics
- Stars: 2
- Watchers: 1
- Forks: 0
- Open Issues: 0
- Releases: 0
Metadata Files
README.md
makemyprior
makemyprior is a tool for easy prior construction and visualization. It helps to formulate joint prior distributions for variance parameters in latent Gaussian models. The resulting prior is robust and can be created in an intuitive way. A graphical user interface (GUI) can be used to choose the joint prior, where the user can click through the model and select priors. An extensive guide is available in the GUI. The package allows for direct inference with the specified model and prior. Using a hierarchical variance decomposition, we formulate a joint variance prior that takes the whole model structure into account. In this way, existing knowledge can intuitively be incorporated at the level it applies to. Alternatively, one can use independent variance priors for each model components in the latent Gaussian model.
Installation
You can install the released version of makemyprior from CRAN with:
r
install.packages("makemyprior")
Example
This is an example showing how to implement a prior.
``` r
library(makemyprior)
set.seed(1) data <- list( a = rep(1:10, each = 10), b = rep(1:10, times = 10) ) data$y <- rnorm(10, 0, 0.4)[data$a] + rnorm(10, 0, 0.6)[data$b] + rnorm(100, 0, 1)
formula <- y ~ mc(a) + mc(b)
prior <- makeprior(formula, data, family = "gaussian", prior = list(tree = "s1 = (a, b); s2 = (s1, eps)", w = list(s2 = list(prior = "pc0", param = 0.25)), V = list(s2 = list(prior = "pc", param = c(3, 0.05)))), interceptprior = c(0, 1000))
summary(prior) # printing details plot(prior) # plotting prior
```
Owner
- Name: Ingeborg Hem
- Login: ingebogh
- Kind: user
- Repositories: 1
- Profile: https://github.com/ingebogh
GitHub Events
Total
- Watch event: 2
Last Year
- Watch event: 2
Packages
- Total packages: 1
-
Total downloads:
- cran 246 last-month
- Total dependent packages: 0
- Total dependent repositories: 0
- Total versions: 6
- Total maintainers: 1
cran.r-project.org: makemyprior
Intuitive Construction of Joint Priors for Variance Parameters
- Homepage: https://github.com/ingebogh/makemyprior
- Documentation: http://cran.r-project.org/web/packages/makemyprior/makemyprior.pdf
- License: GPL-2 | GPL-3 [expanded from: GPL (≥ 2)]
-
Latest release: 1.2.2
published almost 2 years ago