Nothing
library(bayes4psy)
# set seed
seed <- 0
set.seed(0)
# set tolerance
tol <- 0.2
# priors for rgb
mu_prior <- b_prior(family="uniform", pars=c(0, 255))
sigma_prior <- b_prior(family="uniform", pars=c(0, 100))
# attach priors to relevant parameters
priors_rgb <- list(c("mu_r", mu_prior),
c("sigma_r", sigma_prior),
c("mu_g", mu_prior),
c("sigma_g", sigma_prior),
c("mu_b", mu_prior),
c("sigma_b", sigma_prior))
# generate data (rgb) and fit
r <- as.integer(rnorm(50, mean=250, sd=20))
r[r > 255] <- 255
r[r < 0] <- 0
g <- as.integer(rnorm(50, mean=20, sd=20))
g[g > 255] <- 255
g[g < 0] <- 0
b <- as.integer(rnorm(50, mean=40, sd=20))
b[b > 255] <- 255
b[b < 0] <- 0
colors <- data.frame(r=r, g=g, b=b)
fit1 <- b_color(colors=colors, priors=priors_rgb, chains=1, seed=seed, refresh=0)
# priors for hsv
h_prior <- b_prior(family="uniform", pars=c(0, 2*pi))
sv_prior <- b_prior(family="uniform", pars=c(0, 1))
kappa_prior <- b_prior(family="uniform", pars=c(0, 500))
sigma_prior <- b_prior(family="uniform", pars=c(0, 1))
# attach priors to relevant parameters
priors_hsv <- list(c("mu_h", h_prior),
c("kappa_h", kappa_prior),
c("mu_s", sv_prior),
c("sigma_s", sigma_prior),
c("mu_v", sv_prior),
c("sigma_v", sigma_prior))
# generate data (hsv) and fit
h <- rnorm(100, mean=2*pi/3, sd=0.5)
h[h > 2*pi] <- 2*pi
h[h < 0] <- 0
s <- rnorm(100, mean=0.9, sd=0.2)
s[s > 1] <- 1
s[s < 0] <- 0
v <- rnorm(100, mean=0.9, sd=0.2)
v[v > 1] <- 1
v[v < 0] <- 0
colors <- data.frame(h=h, s=s, v=v)
fit2 <- b_color(colors=colors, hsv=TRUE, priors=priors_hsv, chains=1, seed=seed, refresh=0)
# summary
test_that("color summary", {
expect_output(summary(fit1), regexp="mu_h")
expect_output(summary(fit1), regexp="kappa_h")
})
# print and show
test_that("color print and show", {
expect_output(print(fit1), regexp="mu_h")
expect_output(print(fit1), regexp="kappa_h")
})
# get_parameters
test_that("color get_parameters", {
parameters <- get_parameters(fit1)
expect_equal(mean(parameters$h), -0.06, tolerance=tol)
expect_equal(mean(parameters$s), 0.93, tolerance=tol)
expect_equal(mean(parameters$v), 0.96, tolerance=tol)
})
# compare_means two fits
test_that("color compare_means two fits", {
o <- capture.output(output <- compare_means(fit1, fit2=fit2))
r <- c(0, 1, NA)
g <- c(1, 0, NA)
b <- c(1, 0, NA)
h <- c(1, 0, NA)
s <- c(0, 1, NA)
v <- c(0, 1, NA)
compare <- rbind(r, g, b, h, s, v)
expect_equal(output, compare, tolerance=tol)
})
# compare_distributions two fits
test_that("color compare_distributions two fits", {
o <- capture.output(output <- compare_distributions(fit1, fit2=fit2))
r <- c(0.00, 1.00, NA)
g <- c(1.00, 0.00, NA)
b <- c(0.72, 0.28, NA)
h <- c(1.00, 0.00, NA)
s <- c(0.32, 0.68, NA)
v <- c(0.27, 0.73, NA)
compare <- rbind(r, g, b, h, s, v)
expect_equal(output, compare, tolerance=tol)
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.