> descr(dat[, c("y", "P")], "P") %>% print()
  Variables     `1`        `2`        `3`        Total     p      Test          
1  y                                                                            
2    N           3          3          3          9         0.22~  F-test (ANOV~
3    mean        33         29         26         29                            
4    sd          5.3        4.9        2.6        4.9                           
5    median      34         30         26         28                            
6    Q1 - Q3     27 -- 37   24 -- 34   23 -- 28   26 -- 3~                      
7    min - max   27 -- 37   24 -- 34   23 -- 28   23 -- 3~                      

