{\rtf1\ansi\ansicpg1252\cocoartf1671\cocoasubrtf600 {\fonttbl\f0\fswiss\fcharset0 Helvetica;} {\colortbl;\red255\green255\blue255;} {\*\expandedcolortbl;;} \margl1440\margr1440\vieww10800\viewh8400\viewkind0 \pard\tx720\tx1440\tx2160\tx2880\tx3600\tx4320\tx5040\tx5760\tx6480\tx7200\tx7920\tx8640\pardirnatural\partightenfactor0 \f0\fs24 \cf0 > ## load libraries \ > library(brms)\ Loading required package: Rcpp\ Loading 'brms' package (version 2.14.4). Useful instructions\ can be found by typing help('brms'). A more detailed introduction\ to the package is available through vignette('brms_overview').\ \ Attaching package: \'91brms\'92\ \ The following object is masked from \'91package:stats\'92:\ \ ar\ \ > library(mice)\ \ Attaching package: \'91mice\'92\ \ The following object is masked from \'91package:stats\'92:\ \ filter\ \ The following objects are masked from \'91package:base\'92:\ \ cbind, rbind\ \ > \ > \ > # load data \ > load("/udd/nkabe/dissertation/aim2/imp1.Rda")\ > load("/udd/nkabe/dissertation/aim2/df_final.Rda")\ > \ > # normal model \ > m1 <- brm(purge01 ~ 1 + (1|s1), data=df_final)\ Compiling Stan program...\ Start sampling\ \ SAMPLING FOR MODEL '0c166cd5fcfd2d705d7f2b7d3a244928' NOW (CHAIN 1).\ Chain 1: \ Chain 1: Gradient evaluation took 0 seconds\ Chain 1: 1000 transitions using 10 leapfrog steps per transition would take 0 seconds.\ Chain 1: Adjust your expectations accordingly!\ Chain 1: \ Chain 1: \ Chain 1: Iteration: 1 / 2000 [ 0%] (Warmup)\ Chain 1: Iteration: 200 / 2000 [ 10%] (Warmup)\ Chain 1: Iteration: 400 / 2000 [ 20%] (Warmup)\ Chain 1: Iteration: 600 / 2000 [ 30%] (Warmup)\ Chain 1: Iteration: 800 / 2000 [ 40%] (Warmup)\ Chain 1: Iteration: 1000 / 2000 [ 50%] (Warmup)\ Chain 1: Iteration: 1001 / 2000 [ 50%] (Sampling)\ Chain 1: Iteration: 1200 / 2000 [ 60%] (Sampling)\ Chain 1: Iteration: 1400 / 2000 [ 70%] (Sampling)\ Chain 1: Iteration: 1600 / 2000 [ 80%] (Sampling)\ Chain 1: Iteration: 1800 / 2000 [ 90%] (Sampling)\ Chain 1: Iteration: 2000 / 2000 [100%] (Sampling)\ Chain 1: \ Chain 1: Elapsed Time: 31.98 seconds (Warm-up)\ Chain 1: 6.23 seconds (Sampling)\ Chain 1: 38.21 seconds (Total)\ Chain 1: \ \ SAMPLING FOR MODEL '0c166cd5fcfd2d705d7f2b7d3a244928' NOW (CHAIN 2).\ Chain 2: \ Chain 2: Gradient evaluation took 0 seconds\ Chain 2: 1000 transitions using 10 leapfrog steps per transition would take 0 seconds.\ Chain 2: Adjust your expectations accordingly!\ Chain 2: \ Chain 2: \ Chain 2: Iteration: 1 / 2000 [ 0%] (Warmup)\ Chain 2: Iteration: 200 / 2000 [ 10%] (Warmup)\ Chain 2: Iteration: 400 / 2000 [ 20%] (Warmup)\ Chain 2: Iteration: 600 / 2000 [ 30%] (Warmup)\ Chain 2: Iteration: 800 / 2000 [ 40%] (Warmup)\ Chain 2: Iteration: 1000 / 2000 [ 50%] (Warmup)\ Chain 2: Iteration: 1001 / 2000 [ 50%] (Sampling)\ Chain 2: Iteration: 1200 / 2000 [ 60%] (Sampling)\ Chain 2: Iteration: 1400 / 2000 [ 70%] (Sampling)\ Chain 2: Iteration: 1600 / 2000 [ 80%] (Sampling)\ Chain 2: Iteration: 1800 / 2000 [ 90%] (Sampling)\ Chain 2: Iteration: 2000 / 2000 [100%] (Sampling)\ Chain 2: \ Chain 2: Elapsed Time: 30.49 seconds (Warm-up)\ Chain 2: 6.12 seconds (Sampling)\ Chain 2: 36.61 seconds (Total)\ Chain 2: \ \ SAMPLING FOR MODEL '0c166cd5fcfd2d705d7f2b7d3a244928' NOW (CHAIN 3).\ Chain 3: \ Chain 3: Gradient evaluation took 0 seconds\ Chain 3: 1000 transitions using 10 leapfrog steps per transition would take 0 seconds.\ Chain 3: Adjust your expectations accordingly!\ Chain 3: \ Chain 3: \ Chain 3: Iteration: 1 / 2000 [ 0%] (Warmup)\ Chain 3: Iteration: 200 / 2000 [ 10%] (Warmup)\ Chain 3: Iteration: 400 / 2000 [ 20%] (Warmup)\ Chain 3: Iteration: 600 / 2000 [ 30%] (Warmup)\ Chain 3: Iteration: 800 / 2000 [ 40%] (Warmup)\ Chain 3: Iteration: 1000 / 2000 [ 50%] (Warmup)\ Chain 3: Iteration: 1001 / 2000 [ 50%] (Sampling)\ Chain 3: Iteration: 1200 / 2000 [ 60%] (Sampling)\ Chain 3: Iteration: 1400 / 2000 [ 70%] (Sampling)\ Chain 3: Iteration: 1600 / 2000 [ 80%] (Sampling)\ Chain 3: Iteration: 1800 / 2000 [ 90%] (Sampling)\ Chain 3: Iteration: 2000 / 2000 [100%] (Sampling)\ Chain 3: \ Chain 3: Elapsed Time: 35.78 seconds (Warm-up)\ Chain 3: 6.25 seconds (Sampling)\ Chain 3: 42.03 seconds (Total)\ Chain 3: \ \ SAMPLING FOR MODEL '0c166cd5fcfd2d705d7f2b7d3a244928' NOW (CHAIN 4).\ Chain 4: \ Chain 4: Gradient evaluation took 0 seconds\ Chain 4: 1000 transitions using 10 leapfrog steps per transition would take 0 seconds.\ Chain 4: Adjust your expectations accordingly!\ Chain 4: \ Chain 4: \ Chain 4: Iteration: 1 / 2000 [ 0%] (Warmup)\ Chain 4: Iteration: 200 / 2000 [ 10%] (Warmup)\ Chain 4: Iteration: 400 / 2000 [ 20%] (Warmup)\ Chain 4: Iteration: 600 / 2000 [ 30%] (Warmup)\ Chain 4: Iteration: 800 / 2000 [ 40%] (Warmup)\ Chain 4: Iteration: 1000 / 2000 [ 50%] (Warmup)\ Chain 4: Iteration: 1001 / 2000 [ 50%] (Sampling)\ Chain 4: Iteration: 1200 / 2000 [ 60%] (Sampling)\ Chain 4: Iteration: 1400 / 2000 [ 70%] (Sampling)\ Chain 4: Iteration: 1600 / 2000 [ 80%] (Sampling)\ Chain 4: Iteration: 1800 / 2000 [ 90%] (Sampling)\ Chain 4: Iteration: 2000 / 2000 [100%] (Sampling)\ Chain 4: \ Chain 4: Elapsed Time: 34.03 seconds (Warm-up)\ Chain 4: 6.2 seconds (Sampling)\ Chain 4: 40.23 seconds (Total)\ Chain 4: \ Warning message:\ Rows containing NAs were excluded from the model. \ > \ > # brms_multiple model \ > m2 <- brm_multiple(purge01 ~ 1 + (1|s1), data=imp1)\ Compiling the C++ model\ }