inst/doc/circtree.R

In circtree: Regression Trees and Forests for Circular Responses

### 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)