Nothing
### R code from vignette source 'circmax.Rnw'
###################################################
### code chunk number 1: preliminaries
###################################################
options(width = 70, prompt = "R> ", continue = "+ ")
library("circmax")
###################################################
### code chunk number 2: circmax.Rnw:158-162
###################################################
library("circmax")
sdat <- circmax_simulate(n = 1000, beta = c(3, 5, 2), gamma = c(3, 3))
head(sdat)
###################################################
### code chunk number 3: circmax.Rnw:171-173
###################################################
m.circmax <- circmax(y ~ x1 + x2 | x3, data = sdat)
print(m.circmax)
###################################################
### code chunk number 4: circmax.Rnw:209-212
###################################################
sdat <- circtree_simulate(n = 1000, mu = c(0, 2, 1), kappa = c(3, 3, 1))
m.circtree <- circtree(y ~ x1 + x2, data = sdat)
plot(m.circtree)
###################################################
### code chunk number 5: circmax.Rnw:232-235
###################################################
library("circmax")
sdat <- circtree_simulate(n = 1000, mu = c(0, 2, 1), kappa = c(3, 3, 1))
head(sdat)
###################################################
### code chunk number 6: circmax.Rnw:243-245
###################################################
m.circtree <- circtree(y ~ x1 + x2, data = sdat)
logLik(m.circtree)
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.