> descr(dat[, c("y", "P")], "P") %>% print()
  Variables     `1`        `2`        `3`        Total      p        Test       
1  y                                                                            
2    N           24         24         24         72         >0.999   F-test (A~
3    mean        55         55         55         55                            
4    sd          6.2        6.2        6.2        6.1                           
5    median      56         56         56         56                            
6    Q1 - Q3     50 -- 59   50 -- 59   50 -- 59   50 -- 59                      
7    min - max   44 -- 70   44 -- 70   44 -- 70   44 -- 70                      

