Nothing
### R code from vignette source 'circtree.Rnw'
###################################################
### code chunk number 1: preliminaries
###################################################
options(width = 70, prompt = "R> ", continue = "+ ")
library("circtree")
###################################################
### code chunk number 2: circtree.Rnw:158-162
###################################################
library("circtree")
sdat <- circmax_simulate(n = 1000, beta = c(3, 5, 2), gamma = c(3, 3))
head(sdat)
###################################################
### code chunk number 3: circtree.Rnw:171-173
###################################################
m.circmax <- circmax(y ~ x1 + x2 | x3, data = sdat)
print(m.circmax)
###################################################
### code chunk number 4: circtree.Rnw:209-213
###################################################
sdat <- circtree_simulate(n = 1000, mu = c(0, 2, 1), kappa = c(3, 3, 1))
m.circtree <- circtree(y ~ x1 + x2, data = sdat)
class(m.circtree) <- class(m.circtree)[class(m.circtree) != "circtree"] ## plot.circtree fails with cran version
plot(m.circtree)
###################################################
### code chunk number 5: circtree.Rnw:233-236
###################################################
library("circtree")
sdat <- circtree_simulate(n = 1000, mu = c(0, 2, 1), kappa = c(3, 3, 1))
head(sdat)
###################################################
### code chunk number 6: circtree.Rnw:244-246
###################################################
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.