anchor.diff:  
   average_above_test_() ->
       Series = [pseudo_series([5.0, 7.0])],
  -    [?_assertEqual(Series, evaluate_call(<<"averageAbove">>, [Series, 6.0], 0, 0, 0)),
  +    Expected = [with_function_name(hd(Series), "averageAbove")],
  +    [?_assertEqual(Expected, evaluate_call(<<"averageAbove">>, [Series, 6.0], 0, 0, 0)),
        ?_assertEqual([], evaluate_call(<<"averageAbove">>, [Series, 6.5], 0, 0, 0))
       ].
   
   average_below_test_() ->
       Series = [pseudo_series([5.0, 7.0])],
  +    Expected = [with_function_name(hd(Series), "averageBelow")],
       [?_assertEqual([], evaluate_call(<<"averageBelow">>, [Series, 5.5], 0, 0, 0)),
  -     ?_assertEqual(Series, evaluate_call(<<"averageBelow">>, [Series, 6.0], 0, 0, 0)),
  -     ?_assertEqual(Series, evaluate_call(<<"averageBelow">>, [Series, 10], 0, 0, 0))
  +     ?_assertEqual(Expected, evaluate_call(<<"averageBelow">>, [Series, 6.0], 0, 0, 0)),
  +     ?_assertEqual(Expected, evaluate_call(<<"averageBelow">>, [Series, 10], 0, 0, 0))
       ].
   
   npercentile_test_() ->
  
positive.code:     S2 = [pseudo_series([3.0, 9.0, 6.0])], % avg 6.0
      S3 = [pseudo_series([3.0, 12.0, 9.0])], % avg 8.0
      Series = S1 ++ S2 ++ S3,
      [?_assertEqual(S1 ++ S3, evaluate_call(<<"averageOutsidePercentile">>, [Series, 80], 0, 0, 0)),
       ?_assertEqual(Series, evaluate_call(<<"averageOutsidePercentile">>, [Series, 50], 0, 0, 0))
      ].
  
  most_deviant_test_() ->
  
negative.code:     Avgs = lists:map(fun(#series{values=Values}) -> statser_calc:safe_average(Values) end, Series),
      LowPerc = percentile(Avgs, 100 - N),
      HighPerc = percentile(Avgs, N),
      filter_named(fun(#series{values=Values}) ->
                           Avg = statser_calc:safe_average(Values),
                           Avg =< LowPerc orelse Avg >= HighPerc
                   end, Series, "averageOutsidePercentile");
  
  % derivative
  evaluate_call(<<"derivative">>, [Series], _From, _Until, _Now) ->
  