anchor.diff:              let sStd = sStd |> toFloatList 
               let trainX = trnX |> toFloatList
               let trainT = trnT |> toFloatList
  +            let trainTUpper = trnT + trnSigma |> toFloatList
  +            let trainTLower = trnT - trnSigma |> toFloatList
               let upperStdev = List.map2 (fun m s-> m + s) sMean sStd |> List.rev
               let lowerStdev = List.map2 (fun m s-> m - s) sMean sStd
               let revX = List.rev sX
  
positive.code: 
  open SymTensor.Compiler.Cuda
  open RProvider.graphics
  open Models
  open ArrayNDNS
  open SymTensor
  
negative.code:         
              R.lock (fun () ->
                  namedParams [   
                       "x",       box sX
                       "y",       box sMean
                       "ylim",    box [minY; maxY]
                       "col",     box "blue"
                       "type",    box "n"
                       "xlab",    box "x"
                       "ylab",    box "y"]
                  |> R.plot |> ignore
                  namedParams [   
                       "x",       box (sX @ revX)
                       "y",       box (lowerStdev @ upperStdev)
                       "col",     box "beige"
                       "border" , box "NA"]
                  |> R.polygon |>ignore
                  namedParams [ 
                      "x",        box sX
                      "y",        box sMean
                      "col",      box "black"
                      "type",     box "l"
                      "size",     box 2]
                  |> R.lines |>ignore
                  namedParams [ 
                      "x",        box trainX
                      "y",        box trainT
                      "ui",       box trainTUpper
                      "li",       box trainTLower
                      "pch",      box "."
                      "col",      box "red"
                      "add",      box true]
                  |> R.plotCI |>ignore
              )
  
              ()
  