knitr::opts_chunk$set(echo = FALSE)
library(nlme)
library(tibble)
library(Dknitprintr)
library(knitr)
# Use characters as factor levels
set.seed(1)
df = as_tibble(expand.grid(pgroup = letters[1:2], app = LETTERS[1:3], dayc = seq(0,100,by = 10), Short = as.character(1:5), stringsAsFactors = FALSE ))
df$Hours = rnorm(nrow(df), 5,2)

# Use factors
set.seed(1010)
d = as_tibble(expand.grid(pgroup = letters[1:2], app = LETTERS[1:3], dayc = seq(0,100,by = 10), Short = as.character(1:5), stringsAsFactors = TRUE ))
d$Hours = rnorm(nrow(df), 5,2)

Simple output, first uses strings as factors

summary(lme(Hours~pgroup, random = ~1|Short, data = d))
summary(lme(Hours~pgroup, random = ~1|Short, data = df))

All categorical variables, first uses strings as factors

p = summary(lme(Hours~pgroup*app, random = ~1|Short, data = d))
p2 = summary(lme(Hours~pgroup*app, random = ~1|Short, data = df))
knit_print(p, parameter = "years", bold_p = 0.3)
knit_print(p2, parameter = "years", bold_p = 0.3)

All factors, explicit caption

knit_print(p, caption = "My short caption with strings as factors", moredec = 3)
knit_print(p2, caption = "My short caption with factors", moredec = 3)

With numeric covariable

summary(lme(Hours~pgroup*app+dayc, random = ~1|Short, data = d))
summary(lme(Hours~pgroup*app+dayc, random = ~1|Short, data = df))

With numeric covariable and all interactions

summary(lme(Hours~pgroup*app*dayc, random = ~1|Short, data = d))
summary(lme(Hours~pgroup*app*dayc, random = ~1|Short, data = df))

Extended linear model

summary(gls(Hours~pgroup*app*dayc, data = d))
summary(gls(Hours~pgroup*app*dayc, data = df))


dmenne/dknitprintr documentation built on Aug. 3, 2019, 4:14 p.m.