anchor.diff:  ##Bartlett
   <a name="PostHoc"></a>
   ##PostHoc
  +
   This test uses the data shown for Anova.
   *)
   
  +
  +open PostHoc
  +
   Testing.PostHoc.hays contrastMatrix dataOneWay 
  -
  -
  -open PostHoc
  -        
  -//let sumOfSquareContrast (sampleSizes:int[]) (sampleMeans:float[]) (contrast:float[]) =        
  -//    let l           =  Array.fold2 (fun state mi ai -> state + (mi * ai)) 0.0 sampleMeans contrast 
  -//    let denominator = (Array.map2 (fun a n -> (abs a) / (float n)) contrast sampleSizes) |> Array.sum
  -//    (l * l / denominator,l)
  -
  -hays contrastMatrix dataOneWay 
   
   
   //let m1 = Seq.mean dataOneWay.[0]
  
positive.code: |> List.rev
  
  
  (**
  ##Rest; needs sorting
  *)
  
  (*** hide ***)
  // Is this part from the PosHoc module? and why is it here? 
  
  
  
negative.code:         [|55.; 66.; 49.; 64.; 70.; 68.;|] ;   
      |]
  
  let contrastMatrixDsd = contrastMatrix //used in docs for "Anova"
  
  Anova.oneWayAnova dsd
  
  tukeyHSD contrastMatrixDsd dsd 
  