if(!requireNamespace('lattice', quietly = TRUE)){
   install.packages('lattice', dependencies=TRUE, repos='http://cran.rstudio.com/')
}
library(lattice)

## This didn't work due to missing dependencies:
##
## if(!requireNamespace('RVAideMemoire', quietly = TRUE)){
##     install.packages('RVAideMemoire', dependencies=TRUE, repos='http://cran.rstudio.com/')
## }
## require('RVAideMemoire')
##
## Download the source code and install by hand:
##
## https://cran.r-project.org/web/packages/RVAideMemoire/index.html
##
## Source the code you will actuall use:
##
if( file.exists('~/Software/RVAideMemoire/R/mood.medtest.R') ){
    source('~/Software/RVAideMemoire/R/mood.medtest.R')
}

source('./utils.R')

DF = read.table('../../Data/Plant_Growth.txt')
DF$group = as.factor(DF$group)

histogram(~ weight | group, layout=c(1, 3), data=DF)

## Use a nonparametric method:
##
res = mood.medtest(weight ~ group, data=DF)
print('Median Test:')
print(res)

## Use the Pearson chi-squared test:
##
res = gen_median_test(DF, 'weight', 'group')
print(sprintf('pooled_median= %f', res$pooled_median))
print(res$T)
print(res$cs)