C:/Documents and Settings/Stephen/Desktop/Biostatistics Library and Program/Biostatistics/ANOVAGroupsCollection.cs

Go to the documentation of this file.

using System;
using System.Collections.Generic;

namespace StephenAshley.Biostatistics
{
   using NUMBER = Decimal;

   public class ANOVAGroupsCollection : GroupsCollection
   {
      #region *** constructors ***
      public ANOVAGroupsCollection() : base()
      {
         elementsANOVAGroup = new Group(0);
      }

      // <summary>
      // Constructs an <see cref="ANOVAGroupsCollection"/> with a given capacity.
      // </summary>
      // <param name="capacity">The capacity of the <see cref="ANOVAGroupsCollection"/>.</param>
      //public ANOVAGroupsCollection(Int32 capacity)
      //   : base(capacity)
      //{
      //   elementsANOVAGroup = new Group(0);
      //}

      public ANOVAGroupsCollection(List<Group> groups)
         : base(groups)
      {
         elementsANOVAGroup = new Group(0);
         foreach (Group group in groups)
            elementsANOVAGroup.AddRange(group.ToArray());
      }
      #endregion
      #region *** fields ***
      protected Group elementsANOVAGroup;
      #endregion
      #region *** methods ***

      public override void Add(Group group)
      {
         base.Add(group);

         elementsANOVAGroup.AddRange(group.ToArray());
      }

      public override void AddRange(List<Group> lstGroups)
      {
         base.AddRange(lstGroups);
         foreach (Group grp in lstGroups)
            elementsANOVAGroup.AddRange(grp.ToArray());
      }

      public new Group this[int index]
      {
         get 
         {
            Group group;
            try
            {
               group = groupsArray[index]; 
            }
            catch (Exception e)
            {
               throw new BiostatisticsException(e.Source + ": " + e.Message, e);
            }
            return group;
         }
         set
         {
            try
            {
            groupsArray[index] = value;
            }
            catch (Exception e)
            {
               throw new BiostatisticsException(e.Source + ": " + e.Message, e);
            }
            elementsANOVAGroup = new Group(0);
            foreach (Group group in elementsANOVAGroup)
               elementsANOVAGroup.AddRange(group.ToArray());
         }
      }

      public Int32 DegreesOfFreedomBetweenGroups()
      {
         if (k() < 2)
            throw new BiostatisticsException(
         "ANOVAGroupsCollection DegreesOfFreedomBetweenGroups requires more than 1 Group.");
         return k() - 1;
      }

      public Int32 DegreesOfFreedomTotal()
      {
         if (k() < 2)
            throw new BiostatisticsException(
         "ANOVAGroupsCollection DegreesOfFreedomTotal requires more than 1 Group.");
         return nTotal() - 1;
      }

      public Int32 DegreesOfFreedomWithinGroups()
      {return nTotal() - k();}

      public NUMBER F()
      {
         NUMBER decF;
         try
         {
            decF = MeanSquareBetweenGroups() / MeanSquareWithinGroups();
         }
         catch ( Exception e)
         {
            throw new BiostatisticsException(e.Source + ": " + e.Message, e);
         }
         return decF;
      }

      public NUMBER MeanSquareBetweenGroups()
      {
         NUMBER decMeanSquareBetweenGroups;
         try
         {
            decMeanSquareBetweenGroups = SumOfSquaredDeviatesBetweenGroups() /
               DegreesOfFreedomBetweenGroups();
         }
         catch (Exception e)
         {
            throw new BiostatisticsException(e.Source + ": " + e.Message, e);
         }
         return decMeanSquareBetweenGroups;
      }
      public NUMBER MeanSquareWithinGroups()
      {
         NUMBER decMeanSquareWithinGroups;
         try
         {
            decMeanSquareWithinGroups = SumOfSquaredDeviatesWithinGroups() /
               DegreesOfFreedomWithinGroups();
         }
         catch (Exception e)
         {
            throw new BiostatisticsException(e.Source + ": " + e.Message, e);
         }
         return decMeanSquareWithinGroups;
      }

      public NUMBER MeanTotalElements()
      {return elementsANOVAGroup.Mean();}

      public Int32 nTotalElements()
      {return elementsANOVAGroup.n();}

      public virtual NUMBER p()
      {
         return Distributions.ProbabilityF(F(), DegreesOfFreedomBetweenGroups(), 
               DegreesOfFreedomWithinGroups());
      }

      public NUMBER SampleStandardDeviation()
      {return elementsANOVAGroup.SampleStandardDeviation();}

      public NUMBER StandardErrorOfMean()
      {return elementsANOVAGroup.StandardErrorOfMean();}

      public NUMBER SumTotalElements()
      {return elementsANOVAGroup.Sum();}

      public NUMBER SumOfSquaresTotal()
      {return elementsANOVAGroup.SumOfSquares();}

      public NUMBER SumOfSquaredDeviatesTotal()
      {return elementsANOVAGroup.SumOfSquaredDeviates();}

      public NUMBER SumOfSquaredDeviatesBetweenGroups()
      {
         return SumOfSquaredDeviatesTotal() - SumOfSquaredDeviatesWithinGroups();
      }

      public NUMBER SumOfSquaredDeviatesWithinGroups()
      {
         NUMBER decResult = 0m;
            foreach (Group group in groupsArray)
               decResult += group.SumOfSquaredDeviates();
            return decResult;
      }
      #endregion

   }
}

---