Medicine

#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

data()

data(package = .packages(all.available = TRUE))


library(psych)
mydata2 = force(sat.act)

head(mydata2)

##       gender education age ACT SATV SATQ
## 29442      2         3  19  24  500  500
## 29457      2         3  23  35  600  500
## 29498      2         3  20  21  480  470
## 29503      1         4  27  26  550  520
## 29504      1         2  33  31  600  550
## 29518      1         5  26  28  640  640