Sometimes (usually?) relationships between variables are non-linear.
`simstudy`

can already accommodate that. But, if we want to
explicitly generate data from a piece-wise polynomial function to
explore spline methods in particular, or non-linear relationships more
generally. There are three functions that facilitate this:
`viewBasis`

, `viewSplines`

, and
`genSpline`

. The first two functions are more exploratory in
nature, and just provide plots of the B-spline basis functions and the
splines, respectively. The third function actually generates data and
adds to an existing data.table.

The shape of a spline is determined by three factors: (1) the
cut-points or knots that define the piecewise structure of the function,
(2) the polynomial degree, such as linear, quadratic, cubic, etc., and
(3) the linear coefficients that combine the basis functions, which is
contained in a vector or matrix *theta*.

First, we can look at the basis functions, which depend only the knots and degree. The knots are specified as quantiles, between 0 and 1:

```
<- c(0.25, 0.50, 0.75)
knots viewBasis(knots, degree = 2)
```

```
<- c(0.20, 0.40, 0.60, 0.80)
knots viewBasis(knots, degree = 3)
```

The splines themselves are specified as linear combinations of each
of the basis functions. The coefficients of those combinations are
specified in *theta*. Each individual spline curve represents a
specific linear combination of a particular set of basis functions. In
exploring, we can look at a single curve or multiple curves, depending
on whether or not we specify theta as a vector (single) or matrix
(multiple).

```
<- c(0.25, 0.5, 0.75)
knots
# number of elements in theta: length(knots) + degree + 1
= c(0.1, 0.8, 0.4, 0.9, 0.2, 1.0)
theta1
viewSplines(knots, degree = 2, theta1)
```

```
= matrix(c(0.1, 0.2, 0.4, 0.9, 0.2, 0.3, 0.6,
theta2 0.1, 0.3, 0.3, 0.8, 1.0, 0.9, 0.4,
0.1, 0.9, 0.8, 0.2, 0.1, 0.6, 0.1),
ncol = 3)
theta2
```

```
## [,1] [,2] [,3]
## [1,] 0.1 0.1 0.1
## [2,] 0.2 0.3 0.9
## [3,] 0.4 0.3 0.8
## [4,] 0.9 0.8 0.2
## [5,] 0.2 1.0 0.1
## [6,] 0.3 0.9 0.6
## [7,] 0.6 0.4 0.1
```

`viewSplines(knots, degree = 3, theta2)`

We can generate data using a predictor in an existing data set by
specifying the *knots* (in terms of quantiles), a vector of
coefficients in *theta*, the degree of the polynomial, as well as
a range

```
<- defData(varname = "age", formula = "20;60", dist = "uniform")
ddef
= c(0.1, 0.8, 0.6, 0.4, 0.6, 0.9, 0.9)
theta1 <- c(0.25, 0.5, 0.75) knots
```

Here is the shape of the curve that we want to generate data from:

`viewSplines(knots = knots, theta = theta1, degree = 3)`

Now we specify the variables in the data set and generate the data:

```
set.seed(234)
<- genData(1000, ddef)
dt <- genSpline(dt = dt, newvar = "weight",
dt predictor = "age", theta = theta1,
knots = knots, degree = 3,
newrange = "90;160",
noise.var = 64)
```

Hereâ€™s a plot of the data with a smoothed line fit to the data:

```
ggplot(data = dt, aes(x=age, y=weight)) +
geom_point(color = "grey65", size = 0.75) +
geom_smooth(se=FALSE, color="red", size = 1, method = "auto") +
geom_vline(xintercept = quantile(dt$age, knots)) +
theme(panel.grid.minor = element_blank())
```

Finally, we will fit three different spline models to the data - a linear, a quadratic, and a cubic - and plot the predicted values:

```
# normalize age for best basis functions
:= (age - min(age))/(max(age) - min(age))]
dt[, nage
# fit a cubic spline
<- lm(weight ~ bs(x = nage, knots = knots, degree = 3, intercept = TRUE) - 1, data = dt)
lmfit3
# fit a quadtratic spline
<- lm(weight ~ bs(x = nage, knots = knots, degree = 2), data = dt)
lmfit2
# fit a linear spline
<- lm(weight ~ bs(x = nage, knots = knots, degree = 1), data = dt)
lmfit1
# add predicted values for plotting
.3deg := predict(lmfit3)]
dt[, pred.2deg := predict(lmfit2)]
dt[, pred.1deg := predict(lmfit1)]
dt[, pred
ggplot(data = dt, aes(x=age, y=weight)) +
geom_point(color = "grey65", size = 0.75) +
geom_line(aes(x=age, y = pred.3deg), color = "#1B9E77", size = 1) +
geom_line(aes(x=age, y = pred.2deg), color = "#D95F02", size = 1) +
geom_line(aes(x=age, y = pred.1deg), color = "#7570B3", size = 1) +
geom_vline(xintercept = quantile(dt$age, knots)) +
theme(panel.grid.minor = element_blank())
```