> descr(dat[, c("Sex", "distance")], "Sex", test_options = list(paired = TRUE,
+   indices = dat$Subject)) %>% print()
  Variables     Male       Female    Divers    Total     p      Test            
1  distance                                                                     
2    N           64         44        64        172       0.00~  Mixed model AN~
3    mean        25         23        25        24                              
4    sd          2.9        2.4       2.9       3                               
5    median      25         23        25        24                              
6    Q1 - Q3     23 -- 26   21 -- 2~  23 -- 2~  22 -- 2~                        
7    min - max   17 -- 32   16 -- 2~  17 -- 3~  16 -- 3~                        

