# Executive Summary

# Maternal Preconception Study

## Fecal Testosterone

# Introduction

The results of the analysis are presented using the following structure:

1. Searchable High‑Level Summary: Provides a brief overview of the key findings from our analysis with a table of effects
2. Visualisation, Descriptive Statistics and Diagnostics by Variable: Exploratory data analysis and model diagnostics for each behavioral measure
3. Searchable ROPE Analysis Table: Model diagnostics to check the fit of the models

# Effects of Interest

- All groups except for CRV have higher fecal testosterone following resident intruder behavioural testing.
- CRA has higher fecal testosterone, particularly prior to resident intruder behavioural testing.

# Searchable High‑Level Summary

### Enter search terms in the column headers to filter the table below

Note: TD in the Higher column indicates a Trend Difference for the continuous variable.

# Visualisation, Descriptive Statistics and Diagnostics by Variable

### Bar Plot

### Density Plot

Fecal Testosterone by Grouping Variables
Diet | n | Mean | SD | Min | Max | % of Zeros | Zero Flag
--- | --- | --- | --- | --- | --- | --- | ---
Control | 22.0 | 28.6 | 15.6 | 0.5 | 59.2 | 0.0 |
CR25 | 22.0 | 33.1 | 27.6 | 5.7 | 99.7 | 0.0 |
CRA | 20.0 | 41.1 | 20.6 | 9.6 | 85.0 | 0.0 |
CRV | 18.0 | 32.9 | 14.9 | 15.8 | 58.8 | 0.0 |
Warning symbols: * = >40% zeros in this group, ! = >20% zeros in entire sample Either flag triggers the use of a binary choice model.

### Summary Table

Model Summary
Family: gamma
Parameter | Estimate | Std. Error | 2.5% CI | 97.5% CI | Type
--- | --- | --- | --- | --- | ---
Fixed Effects
Intercept | 2.987 | 0.184 | 2.624 | 3.355 | Fixed Effects
DietCR25 | −0.156 | 0.251 | −0.649 | 0.340 | Fixed Effects
DietCRA | 0.446 | 0.256 | −0.050 | 0.954 | Fixed Effects
DietCRV | 0.368 | 0.261 | −0.138 | 0.884 | Fixed Effects
Time | 0.603 | 0.214 | 0.177 | 1.011 | Fixed Effects
DietCR25:Time | 0.412 | 0.294 | −0.160 | 0.994 | Fixed Effects
DietCRA:Time | −0.117 | 0.309 | −0.710 | 0.501 | Fixed Effects
DietCRV:Time | −0.378 | 0.311 | −0.975 | 0.249 | Fixed Effects

### Posterior Check

$data # A tibble: 902 × 7 y_id y_name rep_id rep_label is_y is_y_label value 1 1 1 1 italic(y)[rep] ( 1 ) FALSE italic(y)[rep] 21.2 2 1 1 2 italic(y)[rep] ( 2 ) FALSE italic(y)[rep] 21.5 3 1 1 3 italic(y)[rep] ( 3 ) FALSE italic(y)[rep] 20.9 4 1 1 4 italic(y)[rep] ( 4 ) FALSE italic(y)[rep] 11.8 5 1 1 5 italic(y)[rep] ( 5 ) FALSE italic(y)[rep] 9.10 6 1 1 6 italic(y)[rep] ( 6 ) FALSE italic(y)[rep] 12.2 7 1 1 7 italic(y)[rep] ( 7 ) FALSE italic(y)[rep] 13.9 8 1 1 8 italic(y)[rep] ( 8 ) FALSE italic(y)[rep] 16.1 9 1 1 9 italic(y)[rep] ( 9 ) FALSE italic(y)[rep] 11.2 10 1 1 10 italic(y)[rep] ( 10 ) FALSE italic(y)[rep] 24.0 # ℹ 892 more rows

\(layers\)layers[[1]] mapping: group = ~.data$rep_id, colour = yrep geom_line: na.rm = FALSE, orientation = NA stat_density: bw = nrd0, adjust = 1, kernel = gaussian, n = 1024, trim = FALSE, na.rm = FALSE, bounds = c(-Inf, Inf), orientation = NA position_identity

$layers[[2]] mapping: colour = y geom_line: na.rm = FALSE, orientation = NA, lineend = round stat_density: bw = nrd0, adjust = 1, kernel = gaussian, n = 1024, trim = FALSE, na.rm = FALSE, bounds = c(-Inf, Inf), orientation = NA position_identity

$scales <ggproto object: Class ScalesList, gg> add: function add_defaults: function add_missing: function backtransform_df: function clone: function find: function get_scales: function has_scale: function input: function map_df: function n: function non_position_scales: function scales: list set_palettes: function train_df: function transform_df: function super: <ggproto object: Class ScalesList, gg>

$guides <Guides[0] ggproto object>

\(mapping\)x expr: ^.data$value env: 0x0000017ad7e16200

attr(,“class”) [1] “uneval”

\(theme\)line $colour [1] “black”

$linewidth [1] 0.5

$linetype [1] 1

$lineend [1] “butt”

$arrow [1] FALSE

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_line” “element”

\(rect\)fill [1] “white”

$colour [1] “black”

$linewidth [1] 0.5

$linetype [1] 1

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_rect” “element”

\(text\)family [1] “”

$face [1] “plain”

$colour [1] “black”

$size [1] 11

$hjust [1] 0.5

$vjust [1] 0.5

$angle [1] 0

$lineheight [1] 0.9

$margin [1] 0points 0points 0points 0points

$debug [1] FALSE

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

$title NULL

$aspect.ratio NULL

$axis.title NULL

\(axis.title.x\)family NULL

$face NULL

$colour NULL

$size NULL

$hjust NULL

$vjust [1] 1

$angle NULL

$lineheight NULL

$margin [1] 2.75points 0points 0points 0points

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

\(axis.title.x.top\)family NULL

$face NULL

$colour NULL

$size NULL

$hjust NULL

$vjust [1] 0

$angle NULL

$lineheight NULL

$margin [1] 0points 0points 2.75points 0points

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

$axis.title.x.bottom NULL

\(axis.title.y\)family NULL

$face NULL

$colour NULL

$size NULL

$hjust NULL

$vjust [1] 1

$angle [1] 90

$lineheight NULL

$margin [1] 0points 2.75points 0points 0points

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

$axis.title.y.left NULL

\(axis.title.y.right\)family NULL

$face NULL

$colour NULL

$size NULL

$hjust NULL

$vjust [1] 1

$angle [1] -90

$lineheight NULL

$margin [1] 0points 0points 0points 2.75points

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

\(axis.text\)family NULL

$face NULL

$colour [1] “grey30”

$size [1] 0.8 *

$hjust NULL

$vjust NULL

$angle NULL

$lineheight NULL

$margin NULL

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

\(axis.text.x\)family NULL

$face NULL

$colour NULL

$size NULL

$hjust NULL

$vjust [1] 1

$angle NULL

$lineheight NULL

$margin [1] 2.2points 0points 0points 0points

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

\(axis.text.x.top\)family NULL

$face NULL

$colour NULL

$size NULL

$hjust NULL

$vjust [1] 0

$angle NULL

$lineheight NULL

$margin [1] 0points 0points 2.2points 0points

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

$axis.text.x.bottom NULL

\(axis.text.y\)family NULL

$face NULL

$colour NULL

$size NULL

$hjust [1] 1

$vjust NULL

$angle NULL

$lineheight NULL

$margin [1] 0points 2.2points 0points 0points

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

$axis.text.y.left NULL

\(axis.text.y.right\)family NULL

$face NULL

$colour NULL

$size NULL

$hjust [1] 0

$vjust NULL

$angle NULL

$lineheight NULL

$margin [1] 0points 0points 0points 2.2points

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

$axis.text.theta NULL

\(axis.text.r\)family NULL

$face NULL

$colour NULL

$size NULL

$hjust [1] 0.5

$vjust NULL

$angle NULL

$lineheight NULL

$margin [1] 0points 2.2points 0points 2.2points

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

$axis.ticks list() attr(,“class”) [1] “element_blank” “element”

$axis.ticks.x NULL

$axis.ticks.x.top NULL

$axis.ticks.x.bottom NULL

$axis.ticks.y NULL

$axis.ticks.y.left NULL

$axis.ticks.y.right NULL

$axis.ticks.theta NULL

$axis.ticks.r NULL

$axis.minor.ticks.x.top NULL

$axis.minor.ticks.x.bottom NULL

$axis.minor.ticks.y.left NULL

$axis.minor.ticks.y.right NULL

$axis.minor.ticks.theta NULL

$axis.minor.ticks.r NULL

$axis.ticks.length [1] 2.75points

$axis.ticks.length.x NULL

$axis.ticks.length.x.top NULL

$axis.ticks.length.x.bottom NULL

$axis.ticks.length.y NULL

$axis.ticks.length.y.left NULL

$axis.ticks.length.y.right NULL

$axis.ticks.length.theta NULL

$axis.ticks.length.r NULL

$axis.minor.ticks.length [1] 0.75 *

$axis.minor.ticks.length.x NULL

$axis.minor.ticks.length.x.top NULL

$axis.minor.ticks.length.x.bottom NULL

$axis.minor.ticks.length.y NULL

$axis.minor.ticks.length.y.left NULL

$axis.minor.ticks.length.y.right NULL

$axis.minor.ticks.length.theta NULL

$axis.minor.ticks.length.r NULL

$axis.line list() attr(,“class”) [1] “element_blank” “element”

$axis.line.x NULL

$axis.line.x.top NULL

$axis.line.x.bottom NULL

$axis.line.y NULL

$axis.line.y.left NULL

$axis.line.y.right NULL

$axis.line.theta NULL

$axis.line.r NULL

$legend.background list() attr(,“class”) [1] “element_blank” “element”

$legend.margin [1] 5.5points 5.5points 5.5points 5.5points

$legend.spacing [1] 11points

$legend.spacing.x NULL

$legend.spacing.y NULL

$legend.key list() attr(,“class”) [1] “element_blank” “element”

$legend.key.size [1] 1.2lines

$legend.key.height NULL

$legend.key.width NULL

$legend.key.spacing [1] 5.5points

$legend.key.spacing.x NULL

$legend.key.spacing.y NULL

$legend.frame NULL

$legend.ticks NULL

$legend.ticks.length [1] 0.2 *

$legend.axis.line NULL

\(legend.text\)family NULL

$face NULL

$colour NULL

$size [1] 0.8 *

$hjust NULL

$vjust NULL

$angle NULL

$lineheight NULL

$margin NULL

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

$legend.text.position NULL

\(legend.title\)family NULL

$face NULL

$colour NULL

$size NULL

$hjust [1] 0

$vjust NULL

$angle NULL

$lineheight NULL

$margin NULL

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

$legend.title.position NULL

$legend.position [1] “right”

$legend.position.inside NULL

$legend.direction NULL

$legend.byrow NULL

$legend.justification [1] “center”

$legend.justification.top NULL

$legend.justification.bottom NULL

$legend.justification.left NULL

$legend.justification.right NULL

$legend.justification.inside NULL

$legend.location NULL

$legend.box NULL

$legend.box.just NULL

$legend.box.margin [1] 0cm 0cm 0cm 0cm

$legend.box.background list() attr(,“class”) [1] “element_blank” “element”

$legend.box.spacing [1] 11points

$panel.background list() attr(,“class”) [1] “element_blank” “element”

$panel.border list() attr(,“class”) [1] “element_blank” “element”

$panel.spacing [1] 5.5points

$panel.spacing.x NULL

$panel.spacing.y NULL

\(panel.grid\)colour [1] “grey92”

$linewidth NULL

$linetype NULL

$lineend NULL

$arrow [1] FALSE

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_line” “element”

$panel.grid.major NULL

\(panel.grid.minor\)colour NULL

$linewidth [1] 0.5 *

$linetype NULL

$lineend NULL

$arrow [1] FALSE

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_line” “element”

$panel.grid.major.x NULL

$panel.grid.major.y NULL

$panel.grid.minor.x NULL

$panel.grid.minor.y NULL

$panel.ontop [1] FALSE

$plot.background list() attr(,“class”) [1] “element_blank” “element”

\(plot.title\)family NULL

$face NULL

$colour NULL

$size [1] 1.2 *

$hjust [1] 0.5

$vjust [1] 1

$angle NULL

$lineheight NULL

$margin [1] 0points 0points 5.5points 0points

$debug NULL

$inherit.blank [1] FALSE

attr(,“class”) [1] “element_text” “element”

$plot.title.position [1] “panel”

\(plot.subtitle\)family NULL

$face NULL

$colour NULL

$size NULL

$hjust [1] 0

$vjust [1] 1

$angle NULL

$lineheight NULL

$margin [1] 0points 0points 5.5points 0points

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

\(plot.caption\)family NULL

$face NULL

$colour NULL

$size [1] 0.8 *

$hjust [1] 1

$vjust [1] 1

$angle NULL

$lineheight NULL

$margin [1] 5.5points 0points 0points 0points

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

$plot.caption.position [1] “panel”

\(plot.tag\)family NULL

$face NULL

$colour NULL

$size [1] 1.2 *

$hjust [1] 0.5

$vjust [1] 0.5

$angle NULL

$lineheight NULL

$margin NULL

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

$plot.tag.position [1] “topleft”

$plot.tag.location NULL

$plot.margin [1] 5.5points 5.5points 5.5points 5.5points

$strip.background list() attr(,“class”) [1] “element_blank” “element”

$strip.background.x NULL

$strip.background.y NULL

$strip.clip [1] “inherit”

$strip.placement [1] “inside”

\(strip.text\)family NULL

$face NULL

$colour [1] “grey10”

$size [1] 0.8 *

$hjust NULL

$vjust NULL

$angle NULL

$lineheight NULL

$margin [1] 4.4points 4.4points 4.4points 4.4points

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

$strip.text.x NULL

$strip.text.x.bottom NULL

$strip.text.x.top NULL

\(strip.text.y\)family NULL

$face NULL

$colour NULL

$size NULL

$hjust NULL

$vjust NULL

$angle [1] -90

$lineheight NULL

$margin NULL

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

\(strip.text.y.left\)family NULL

$face NULL

$colour NULL

$size NULL

$hjust NULL

$vjust NULL

$angle [1] 90

$lineheight NULL

$margin NULL

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

$strip.text.y.right NULL

$strip.switch.pad.grid [1] 2.75points

$strip.switch.pad.wrap [1] 2.75points

attr(,“class”) [1] “theme” “gg”
attr(,“complete”) [1] TRUE attr(,“validate”) [1] TRUE

$coordinates <ggproto object: Class CoordCartesian, Coord, gg> aspect: function backtransform_range: function clip: on default: FALSE distance: function draw_panel: function expand: FALSE is_free: function is_linear: function labels: function limits: list modify_scales: function range: function render_axis_h: function render_axis_v: function render_bg: function render_fg: function reverse: none setup_data: function setup_layout: function setup_panel_guides: function setup_panel_params: function setup_params: function train_panel_guides: function transform: function super: <ggproto object: Class CoordCartesian, Coord, gg>

$facet <ggproto object: Class FacetNull, Facet, gg> attach_axes: function attach_strips: function compute_layout: function draw_back: function draw_front: function draw_labels: function draw_panel_content: function draw_panels: function finish_data: function format_strip_labels: function init_gtable: function init_scales: function map_data: function params: list set_panel_size: function setup_data: function setup_panel_params: function setup_params: function shrink: TRUE train_scales: function vars: function super: <ggproto object: Class FacetNull, Facet, gg>

$plot_env <environment: 0x0000017ad7e16200>

$layout <ggproto object: Class Layout, gg> coord: NULL coord_params: list facet: NULL facet_params: list finish_data: function get_scales: function layout: NULL map_position: function panel_params: NULL panel_scales_x: NULL panel_scales_y: NULL render: function render_labels: function reset_scales: function resolve_label: function setup: function setup_panel_guides: function setup_panel_params: function train_position: function super: <ggproto object: Class Layout, gg>

\(labels\)labels$title [1] “Posterior Predictive Check”

\(labels\)x NULL

\(labels\)y NULL

\(labels\)group [1] “rep_id”

\(labels\)colour [1] “colour”

\(labels\)fill [1] “fill” attr(,“fallback”) [1] TRUE

\(labels\)weight [1] “weight” attr(,“fallback”) [1] TRUE

attr(,“class”) [1] “gg” “ggplot”

### Convergence

Effect | Median | [95% HDI] | Direction Prob. | ROPE %
--- | --- | --- | --- | ---
Diet: Control vs CR25
Diet: Control vs CR25 | -0.05 | [-0.44, 0.35] | 58.65 | 30.05
Diet: Control vs CRA
Diet: Control vs CRA | -0.38 | [-0.80, 0.02] | 97.05 | 4.50
Diet: Control vs CRV
Diet: Control vs CRV | -0.18 | [-0.57, 0.26] | 81.05 | 19.55
Diet: CR25 vs CRA
Diet: CR25 vs CRA | -0.34 | [-0.73, 0.09] | 94.90 | 8.40
Diet: CR25 vs CRV
Diet: CR25 vs CRV | -0.14 | [-0.55, 0.30] | 74.15 | 24.35
Diet: CRA vs CRV
Diet: CRA vs CRV | 0.21 | [-0.24, 0.63] | 82.50 | 19.15

Effect | Median | [95% HDI] | Direction Prob. | ROPE %
--- | --- | --- | --- | ---
Time: Overall Change
Time: Overall Change 0 - 1 | 0.58 | [0.35, 0.82] | 100 | 0

Effect | Median | [95% HDI] | Direction Prob. | ROPE %
--- | --- | --- | --- | ---
Time: Control
Time: Control 0 | 2.99 | [2.67, 3.36] | 100.00 | 0.00
Time: Control 1 | 3.58 | [3.26, 3.95] | 100.00 | 0.00
Time: Control 0 - 1 | 0.60 | [0.20, 1.01] | 99.80 | 0.90
Time: CR25
Time: CR25 0 | 2.83 | [2.45, 3.20] | 100.00 | 0.00
Time: CR25 1 | 3.85 | [3.50, 4.23] | 100.00 | 0.00
Time: CR25 0 - 1 | 1.02 | [0.58, 1.45] | 100.00 | 0.00
Time: CRA
Time: CRA 0 | 3.42 | [3.07, 3.82] | 100.00 | 0.00
Time: CRA 1 | 3.91 | [3.54, 4.29] | 100.00 | 0.00
Time: CRA 0 - 1 | 0.49 | [0.03, 0.94] | 98.30 | 3.00
Time: CRV
Time: CRV 0 | 3.36 | [2.99, 3.77] | 100.00 | 0.00
Time: CRV 1 | 3.59 | [3.18, 4.02] | 100.00 | 0.00
Time: CRV 0 - 1 | 0.22 | [-0.24, 0.69] | 82.35 | 16.45

Effect | Median | [95% HDI] | Direction Prob. | ROPE %
--- | --- | --- | --- | ---
Time × Diet: Control - CR25
Time × Diet: Control - CR25 0 | 0.15 | [-0.33, 0.61] | 73.70 | 21.55
Time × Diet: Control - CR25 1 | -0.27 | [-0.74, 0.23] | 84.90 | 13.85
Time × Diet: Control - CRA
Time × Diet: Control - CRA 0 | -0.44 | [-0.94, 0.07] | 95.85 | 5.20
Time × Diet: Control - CRA 1 | -0.34 | [-0.82, 0.16] | 90.90 | 10.05
Time × Diet: Control - CRV
Time × Diet: Control - CRV 0 | -0.38 | [-0.91, 0.12] | 92.60 | 7.00
Time × Diet: Control - CRV 1 | 0.00 | [-0.52, 0.53] | 50.85 | 23.30
Time × Diet: CR25 - CRA
Time × Diet: CR25 - CRA 0 | -0.60 | [-1.15, -0.09] | 98.25 | 1.90
Time × Diet: CR25 - CRA 1 | -0.06 | [-0.56, 0.50] | 58.85 | 23.85
Time × Diet: CR25 - CRV
Time × Diet: CR25 - CRV 0 | -0.54 | [-1.08, 0.01] | 96.85 | 4.15
Time × Diet: CR25 - CRV 1 | 0.27 | [-0.27, 0.82] | 82.55 | 14.90
Time × Diet: CRA - CRV
Time × Diet: CRA - CRV 0 | 0.06 | [-0.46, 0.63] | 59.40 | 22.70
Time × Diet: CRA - CRV 1 | 0.33 | [-0.19, 0.92] | 87.80 | 11.80

# End of Document