public StatSummary(IEnumerable<double> values)
{
var list = values.ToList();
N = list.Count;
if (N == 0)
throw new InvalidOperationException("StatSummary: Sequence contains no elements");
list.Sort();
if (N == 1)
Q1 = Median = Q3 = list[0];
else
{
Func<IList<double>, double> getMedian = x => x.Count % 2 == 0
? (x[x.Count / 2 - 1] + x[x.Count / 2]) / 2
: x[x.Count / 2];
Median = getMedian(list);
Q1 = getMedian(list.Take(N / 2).ToList());
Q3 = getMedian(list.Skip((N + 1) / 2).ToList());
}
Min = list.First();
Mean = list.Average();
Max = list.Last();
InterquartileRange = Q3 - Q1;
LowerFence = Q1 - 1.5 * InterquartileRange;
UpperFence = Q3 + 1.5 * InterquartileRange;
Outlier = list.Where(value => value < LowerFence || value > UpperFence).ToArray();
StandardDeviation = N == 1 ? 0 : Math.Sqrt(list.Sum(d => Math.Pow(d - Mean, 2)) / (N - 1));
StandardError = StandardDeviation / Math.Sqrt(N);
ConfidenceInterval = new ConfidenceInterval(Mean, StandardError);
}
public StatSummary(IEnumerable <double> values)
{
var list = values.ToList();
var n = list.Count;
if (n == 0)
{
throw new InvalidOperationException("Sequence contains no elements");
}
list.Sort();
if (n == 1)
{
Q1 = Median = Q3 = list[0];
}
else
{
Func <IList <double>, double> getMedian = x => x.Count % 2 == 0
? (x[x.Count / 2 - 1] + x[x.Count / 2]) / 2
: x[x.Count / 2];
Median = getMedian(list);
Q1 = getMedian(list.Take(n / 2).ToList());
Q3 = getMedian(list.Skip((n + 1) / 2).ToList());
}
Min = list.First();
Mean = list.Average();
Max = list.Last();
InterquartileRange = Q3 - Q1;
LowerFence = Q1 - 1.5 * InterquartileRange;
UpperFence = Q3 + 1.5 * InterquartileRange;
Outlier = list.Where(value => value <LowerFence || value> UpperFence).ToArray();
StandardDeviation = n == 1 ? 0 : Math.Sqrt(list.Sum(d => Math.Pow(d - Mean, 2)) / (n - 1));
StandardError = StandardDeviation / Math.Sqrt(n);
ConfidenceInterval = new ConfidenceInterval(Mean, 1.96 * StandardError);
}
