#
# Written by:
# -- 
# John L. Weatherwax                2009-04-21
# 
# email: wax@alum.mit.edu
# 
# Please send comments and especially bug reports to the
# above email address.
#
#-----

source('wilcoxon_signed_rank_interval.R')
source('wilcoxon_rank_sum_interval.R')

# Ex 15.17: EPage 618; Anwsers EPage 708 
#
data = c( 17.0, 21.4, 30.6, 5.0, 12.2, 11.8, 17.3, 18.8 ) 
c = 32
out = wilcoxon_signed_rank_interval( data, c ) 
print( out )

# Ex 15.18: EPage 618; Anwsers EPage 708 
#
DF = read.csv( "../../Data/CH15/ex15-03.txt", header=TRUE, quote="'" )
c = 93
out = wilcoxon_signed_rank_interval( DF$Observations, c ) 
print( out )

# Ex 15.19: EPage 618; Anwsers EPage 708; Example 15.2 EPage 604
#
data = c( -0.33, -0.41, -0.71, 0.19, -0.52, 0.20, -0.65, -0.14 )
c = 32
out = wilcoxon_signed_rank_interval( data, c ) 
print( out )

# Ex 15.20: EPage 618; Anwsers EPage 708 
#x
#data = c( 0.045, 0.117, 0.062, 0.072 ) 

# Ex 15.21: EPage 618; Anwsers EPage 708; Exercise 10 EPage 613
#
x_data = c( 229, 286, 245, 299, 250 )
y_data = c( 213, 179, 163, 247, 225 ) 
c = 21
out = wilcoxon_rank_sum_interval( x_data, y_data, c )
print( out )

# Ex 15.22: EPage 618; Anwsers EPage 708; Exercise 11 EPage 613
#
DF = read.csv( "../../Data/CH15/ex15-11.txt", header=TRUE, quote="'" )
x_data = DF[ DF$Wood == 'Pine', ]$Time # pick length(x_data)<length(y_data)
y_data = DF[ DF$Wood == 'Oak', ]$Time
c = 44
out = wilcoxon_rank_sum_interval( x_data, y_data, c )
print( out )