#install.packages("readxl")
library(readxl)
mydata <- read_xlsx("./Medicine.xlsx")
#install.packages("psych")
library(psych)
describeBy(mydata$Headache, group = mydata$Medicine)
##
## Descriptive statistics by group
## group: A
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 10 41.5 10.56 39.5 41 11.86 28 59 31 0.3 -1.45 3.34
## ------------------------------------------------------------
## group: B
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 10 37 10.37 40.5 37.75 5.93 18 50 32 -0.66 -1.1 3.28
## ------------------------------------------------------------
## group: C
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 10 26.8 8.84 24 26.75 11.12 15 39 24 0.14 -1.84 2.8