anchor.diff:          [|55.; 66.; 49.; 64.; 70.; 68.;|] ;   
       |]
   
  -let contrastMatrixDsd = contrastMatrix
  -    //[| 
  -     //[|1.0;-1.0;0.0;|] ; 
  -     //[|1.0;0.0;-1.0;|] ; 
  -     //[|0.0;1.0;-1.0;|]   
  -    //|]
  -
  -dsd |> Array.map Array.average   
  -   
  +let contrastMatrixDsd = contrastMatrix //used in docs for "Anova"
  +
   Anova.oneWayAnova dsd
   
   tukeyHSD contrastMatrixDsd dsd 
  
positive.code: ##Bartlett
  <a name="PostHoc"></a>
  ##PostHoc
  This test uses the data shown for Anova.
  *)
  
  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]
  
negative.code: |> List.rev
  
  
  
  