Nothing
## ----include = FALSE----------------------------------------------------------
knitr::opts_chunk$set(
collapse = TRUE,
comment = "#>"
)
## ----setup--------------------------------------------------------------------
library(GLMMcosinor)
## ----message=F, warning=F, eval=F---------------------------------------------
# testdata <- simulate_cosinor(
# n = 200,
# mesor = 1,
# amp = 2,
# acro = 1.2,
# period = 12,
# n_period = 3,
# family = "poisson"
# )
#
# testdata
## ----message=F, warning=F, eval=F---------------------------------------------
# object <- cglmm(
# Y ~ amp_acro(times,
# n_components = 1,
# period = 12
# ),
# data = testdata,
# family = poisson()
# )
# summary(object)
# autoplot(object, superimpose.data = TRUE)
## ----message=F, warning=F, eval=F---------------------------------------------
# simulate_cosinor(
# n = 100,
# mesor = 1,
# amp = 2,
# acro = 1.2,
# period = 12,
# beta.group = TRUE,
# beta.mesor = 0.4,
# beta.amp = 0.5,
# beta.acro = 0.2,
# n_period = 3,
# n_components = 1,
# family = "poisson"
# )
## ----message=F, warning=F, echo=F---------------------------------------------
testdata <- simulate_cosinor(
n = 100,
mesor = 1,
amp = 2,
acro = 1.2,
period = 12,
beta.group = TRUE,
beta.mesor = 0.4,
beta.amp = 0.5,
beta.acro = 0.2,
n_components = 1,
n_period = 3,
family = "gaussian"
)
DT::datatable(testdata)
object <- cglmm(
Y ~ group + amp_acro(times,
n_components = 1,
period = 12,
group = "group"
),
data = testdata,
family = gaussian()
)
summary(object)
autoplot(object, superimpose.data = TRUE)
## ----message=F, warning=F, eval=F---------------------------------------------
# testdata <- simulate_cosinor(
# n = 200,
# mesor = 1,
# amp = c(0.2, 1),
# acro = c(1.2, 2),
# period = c(12, 6),
# n_components = 2,
# n_period = 2,
# family = "poisson"
# )
#
# testdata
# object <- cglmm(
# Y ~ amp_acro(times,
# n_components = 2,
# period = c(12, 6)
# ),
# data = testdata,
# family = poisson()
# )
# summary(object)
# autoplot(object, superimpose.data = TRUE)
## ----message=F, warning=F, fig.show = "hold"----------------------------------
testdata <- simulate_cosinor(100,
mesor = 7,
amp = c(0.1, 0.4, 0.5),
acro = c(1, 1.5, 0.1),
beta.mesor = 4.4,
beta.amp = c(2, 1, 0.4),
beta.acro = c(1, -1.5, -1),
family = "poisson",
period = c(12, 6, 8),
n_period = 2,
n_components = 3
)
object <- cglmm(Y ~ group + amp_acro(times,
n_components = 3,
period = c(12, 6, 8),
group = "group"
), data = testdata, family = poisson())
summary(object)
autoplot(object,
superimpose.data = TRUE,
x_str = "group",
predict.ribbon = FALSE
)
## ----message=F, warning=F, fig.show = "hold"----------------------------------
testdata <- simulate_cosinor(500,
mesor = 1,
amp = c(0.5, 0.5, 0.5),
acro = c(pi, pi / 2, pi),
alpha = 2,
beta.mesor = 2,
beta.amp = c(0.2, 0.2, 0.2),
beta.acro = c(pi / 2, pi, pi / 2),
beta.alpha = 3,
family = "gamma",
period = c(12, 6, 8),
n_period = 2,
n_components = 3
)
object <- cglmm(Y ~ group + amp_acro(times,
n_components = 3,
period = c(12, 6, 8),
group = "group"
), data = testdata, family = Gamma(link = "log"))
summary(object)
autoplot(object,
superimpose.data = TRUE,
x_str = "group",
predict.ribbon = FALSE
)
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.