# Generate example data
set.seed(1199)
grp1 <- rnorm(20, 45, 5)
grp2 <- c(rnorm(10, 45, 10), rnorm(10, 52, 20))
df <- data.frame(cbind(grp1=grp1, grp2=grp2))
head(df)
bivarplot(df)
# Reshape the data into long format
df <- wide2long(df)
head(df)
# Differences between the basic stats
calcstatdif(df)
Run the code above in your browser using DataLab