Note: TD in the Higher column indicates a Trend Difference for the continuous variable.
Executive Summary
Maternal Preconception Study
Fecal Testosterone
Introduction
The results of the analysis are presented using the following structure:
- Searchable High‑Level Summary: Provides a brief overview of the key findings from our analysis with a table of effects
- Visualisation, Descriptive Statistics and Diagnostics by Variable: Exploratory data analysis and model diagnostics for each behavioral measure
- 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
Visualisation, Descriptive Statistics and Diagnostics by Variable
Bar Plot
Density Plot
| Fecal Testosterone by Grouping Variables | ||||||||
| Diet | Time | n | Mean | SD | Min | Max | % of Zeros | Zero Flag |
|---|---|---|---|---|---|---|---|---|
| Control | 0 | 11.0 | 20.5 | 13.7 | 0.5 | 41.9 | 0.0 | |
| Control | 1 | 11.0 | 36.7 | 13.4 | 14.0 | 59.2 | 0.0 | |
| CR25 | 0 | 11.0 | 17.7 | 14.6 | 5.7 | 55.7 | 0.0 | |
| CR25 | 1 | 11.0 | 48.4 | 29.5 | 17.4 | 99.7 | 0.0 | |
| CRA | 0 | 10.0 | 31.9 | 16.1 | 9.6 | 65.9 | 0.0 | |
| CRA | 1 | 10.0 | 50.2 | 21.3 | 19.5 | 85.0 | 0.0 | |
| CRV | 0 | 9.0 | 29.9 | 14.9 | 15.8 | 58.0 | 0.0 | |
| CRV | 1 | 9.0 | 35.8 | 15.2 | 15.9 | 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
\(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
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: 0x000001e4e4e1e200>
$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.45, 0.38] | 61.00 | 28.25 |
| Diet: Control vs CRA | ||||
| Diet: Control vs CRA | -0.38 | [-0.80, 0.03] | 96.90 | 4.85 |
| Diet: Control vs CRV | ||||
| Diet: Control vs CRV | -0.18 | [-0.64, 0.23] | 79.75 | 20.85 |
| Diet: CR25 vs CRA | ||||
| Diet: CR25 vs CRA | -0.34 | [-0.74, 0.10] | 94.50 | 8.85 |
| Diet: CR25 vs CRV | ||||
| Diet: CR25 vs CRV | -0.13 | [-0.57, 0.30] | 72.20 | 23.55 |
| Diet: CRA vs CRV | ||||
| Diet: CRA vs CRV | 0.21 | [-0.21, 0.67] | 82.55 | 19.65 |
| 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 × Diet: Control - CR25 | ||||
| Time × Diet: Control - CR25 0 | 0.15 | [-0.37, 0.60] | 72.60 | 22.55 |
| Time × Diet: Control - CR25 1 | -0.26 | [-0.81, 0.24] | 82.15 | 9.75 |
| Time × Diet: Control - CRA | ||||
| Time × Diet: Control - CRA 0 | -0.44 | [-0.91, 0.08] | 96.00 | 5.15 |
| Time × Diet: Control - CRA 1 | -0.33 | [-0.87, 0.16] | 90.70 | 5.30 |
| Time × Diet: Control - CRV | ||||
| Time × Diet: Control - CRV 0 | -0.38 | [-0.88, 0.13] | 93.25 | 8.65 |
| Time × Diet: Control - CRV 1 | 0.01 | [-0.54, 0.52] | 50.80 | 15.50 |
| Time × Diet: CR25 - CRA | ||||
| Time × Diet: CR25 - CRA 0 | -0.60 | [-1.10, -0.06] | 98.55 | 2.20 |
| Time × Diet: CR25 - CRA 1 | -0.08 | [-0.61, 0.44] | 61.05 | 14.00 |
| Time × Diet: CR25 - CRV | ||||
| Time × Diet: CR25 - CRV 0 | -0.53 | [-1.11, -0.06] | 97.60 | 3.80 |
| Time × Diet: CR25 - CRV 1 | 0.25 | [-0.28, 0.76] | 82.50 | 9.50 |
| Time × Diet: CRA - CRV | ||||
| Time × Diet: CRA - CRV 0 | 0.07 | [-0.49, 0.58] | 59.85 | 25.30 |
| Time × Diet: CRA - CRV 1 | 0.34 | [-0.21, 0.90] | 88.40 | 6.70 |
| Effect | Median | [95% HDI] | Direction Prob. | ROPE % |
|---|---|---|---|---|
| Time: Control | ||||
| Time: Control 0 - 1 | 0.60 | [0.18, 1.06] | 99.70 | 0.60 |
| Time: CR25 | ||||
| Time: CR25 0 - 1 | 1.01 | [0.55, 1.44] | 100.00 | 0.00 |
| Time: CRA | ||||
| Time: CRA 0 - 1 | 0.49 | [0.00, 0.90] | 98.35 | 2.30 |
| Time: CRV | ||||
| Time: CRV 0 - 1 | 0.22 | [-0.24, 0.70] | 82.00 | 10.45 |
| Effect | Median | [95% HDI] | Direction Prob. | ROPE % |
|---|---|---|---|---|
| Diet: Control vs CR25 | Time = 0 | ||||
| Diet: Control vs CR25 | Time = 0 0 | 0.15 | [-0.37, 0.60] | 72.60 | 22.55 |
| Diet: Control vs CRA | Time = 0 | ||||
| Diet: Control vs CRA | Time = 0 0 | -0.44 | [-0.91, 0.08] | 96.00 | 5.15 |
| Diet: Control vs CRV | Time = 0 | ||||
| Diet: Control vs CRV | Time = 0 0 | -0.38 | [-0.88, 0.13] | 93.25 | 8.65 |
| Diet: CR25 vs CRA | Time = 0 | ||||
| Diet: CR25 vs CRA | Time = 0 0 | -0.60 | [-1.10, -0.06] | 98.55 | 2.20 |
| Diet: CR25 vs CRV | Time = 0 | ||||
| Diet: CR25 vs CRV | Time = 0 0 | -0.53 | [-1.11, -0.06] | 97.60 | 3.80 |
| Diet: CRA vs CRV | Time = 0 | ||||
| Diet: CRA vs CRV | Time = 0 0 | 0.07 | [-0.49, 0.58] | 59.85 | 25.30 |
| Diet: Control vs CR25 | Time = 1 | ||||
| Diet: Control vs CR25 | Time = 1 1 | -0.26 | [-0.81, 0.24] | 82.15 | 9.75 |
| Diet: Control vs CRA | Time = 1 | ||||
| Diet: Control vs CRA | Time = 1 1 | -0.33 | [-0.87, 0.16] | 90.70 | 5.30 |
| Diet: Control vs CRV | Time = 1 | ||||
| Diet: Control vs CRV | Time = 1 1 | 0.01 | [-0.54, 0.52] | 50.80 | 15.50 |
| Diet: CR25 vs CRA | Time = 1 | ||||
| Diet: CR25 vs CRA | Time = 1 1 | -0.08 | [-0.61, 0.44] | 61.05 | 14.00 |
| Diet: CR25 vs CRV | Time = 1 | ||||
| Diet: CR25 vs CRV | Time = 1 1 | 0.25 | [-0.28, 0.76] | 82.50 | 9.50 |
| Diet: CRA vs CRV | Time = 1 | ||||
| Diet: CRA vs CRV | Time = 1 1 | 0.34 | [-0.21, 0.90] | 88.40 | 6.70 |