Nothing

```
library("model4you")
set.seed(2018)
n <- 160
trt <- factor(rep(0:1, each = n/2))
y <- 4 + (trt == 1) + rnorm(n)
z <- matrix(rnorm(n * 2), ncol = 2)
dat <- data.frame(y, trt, z)
mod <- lm(y ~ trt, data = dat)
## Note that ntree should usually be higher
frst <- pmforest(mod, ntree = 20)
pmods <- pmodel(frst, fun = identity)
sum(objfun(pmods))
## Note that B should be at least 100
## The low B is just for demonstration
## purposes.
tst <- pmtest(forest = frst,
pmodels = pmods,
B = 5)
tst
```

**Any scripts or data that you put into this service are public.**

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.