Reproducibility in Rmd

reproducible.Rmd

---
title: "reproducible rmd"
output: pdf_document
---

```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE)
```


```{r}
library(mrgsolve)

set.seed(11223)

data <- expand.ev(ID = 1:1000, amt = 100, ii = 24, addl = 27)

mod <- mrgsolve:::house(
  end = 24*28, 
  delta = 1, 
  omega = dmat(1,1,1,1),
  sigma = dmat(1)
)

sum(mrgsim(mod,data,output="matrix"))
```