#HW chpt 6

#6.1
x1 <- c(0,1,3,3,5,10,13,17,26) # control subjects
y1 <- c(0,6,7,9,11,13,20,20,24)      # TR
z1 <- c(0,5,8,9,11,13,16,17,20) # VTP
zz1 <- c(1,5,12,13,19,22,25,27,29) # RH

kruskal.test(list(x1, y1,z1,zz1))

#pvalue
pchisq(4.2646,3, ncp=0, lower.tail=F)

#OR
x <- c(x1, y1, z1, zz1)
g <- factor(rep(1:4, c(length(x1),length(y1),length(z1),length(zz1))))
kruskal.test(x, g)


#===========================================================

#6.4
x1 <- c(46,28,46,37,32,41,42,45,38,44) # site1
y1 <- c(42,60,32,42,45,58,27,51,42,52)      # site2
z1 <- c(38,33,26,25,28,28,26,27,27,27) # site3
zz1 <- c(31,30,27,29,30,25,25,24,27,30) # site4

kruskal.test(list(x1, y1,z1,zz1))

#OR
x <- c(x1, y1, z1, zz1)
g <- factor(rep(1:4, c(length(x1),length(y1),length(z1),length(zz1))))
kruskal.test(x, g)


#===========================================================

#6.8
a1 <- c(3500,3500,3500,4000,4000,4000,4300,4500,4900,5200,6000,6750,8000) # normal
b1 <- c(5710,6110,8060,11400)      # Anemia
c1 <- c(2930,3330,3880,4280,5120) # Lym
d1 <- c(6320,6860,11400,14000) # Mye
e1 <- c(3230,3880,7640,8280,16200,18250,29900) # Acute

kruskal.test(list(a1, b1,c1,d1,e1))

#OR
x <- c(a1, b1,c1,d1,e1)
g <- factor(rep(1:5, c(length(a1),length(b1),length(c1),length(d1),length(e1))))
kruskal.test(x, g)


#===========================================================

#6.12
a1 <- c(52.2,56.4,57.1,46.9,52.5,63,52,61.1,55.3,46.2,57.2) # scotch
b1 <- c(46.7,60.5,58.9,82.9,65.8,93.3,66.9,70.9,73.7,65.8,90.2,68.9) # Blue
c1 <- c(75.2,63.7,73.2,66.2,67.4,69.4,70.4,72.3,63.6,61.9,74.4,70.1) # White

kruskal.test(list(a1, b1,c1))

#OR
x <- c(a1, b1,c1)
g <- factor(rep(1:3, c(length(a1),length(b1),length(c1))))
kruskal.test(x, g)

#===========================================================

#6.42