inst/doc/circmax.R

In circmax: Circular Regression with Maximum Likelihood Estimation and Regression Trees

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