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

Go to the documentation of this file.

using System;
using System.Collections.Generic;

namespace StephenAshley.Biostatistics
{
   using NUMBER = Decimal;

   public class RepeatedMeasuresANOVAGroupsCollection : ANOVAGroupsCollection
   {
      public RepeatedMeasuresANOVAGroupsCollection()
         : base()
      {
         N = 0;
      }

      // <summary>
      // Constructs a <see cref="RepeatedMeasuresANOVAGroupsCollection"/> 
      // with a specified capacity.
      // </summary>
      // <param name="capacity"></param>The capacity of the 
      // <see cref="RepeatedMeasuresANOVAGroupsCollection"/>.
      //public RepeatedMeasuresANOVAGroupsCollection(Int32 capacity)
      //   : base(capacity)
      //{
      //   N = 0;
      //}

      public RepeatedMeasuresANOVAGroupsCollection(List<Group> groups)
         : base(groups)
      {
         N = groups[0].n();
         Int32 iSize = N;
         for (Int32 i = 1; i < k(); i++)
            if (groups[i].n() != N)
               throw new BiostatisticsException(
                  "Groups in RepeatedMeasuresANOVAGroupsCollection have unequal sizes.");
      }

      private Int32 N;

      public override void Add(Group group)
      {
         if (group == null)
            throw new BiostatisticsException(
   "group parameter in Add method of RepeatedMeasuresGroupsCollection is null.");
         if (N == 0)
            N = group.n();
         else if (N != group.n())
         {
            throw new BiostatisticsException(
"Group to be added to RepeatedMeasuresGroupsCollection does not have the same number of elements as the other Groups in the collection.");
         }
         base.Add(group);
      }

      public override void AddRange(List<Group> lstGroups)
      {
         base.AddRange(lstGroups);
      }

      public Int32 DegreesOfFreedomError()
      {
         return DegreesOfFreedomWithinGroups()
               - DegreesOfFreedomSubjects();
      }

      public Int32 DegreesOfFreedomSubjects()
      { return N - 1; }

      public new NUMBER F()
      { return MeanSquareBetweenGroups() / MeanSquareError(); }

      public NUMBER MeanSquareError()
      {
         return SumOfSquaredDeviatesError()
               / DegreesOfFreedomError();
      }

      public override NUMBER p()
      {
         return Distributions.ProbabilityF(F(), DegreesOfFreedomBetweenGroups(),
               DegreesOfFreedomError());
      }
      public NUMBER SumOfSquaredDeviatesError()
      {
         return SumOfSquaredDeviatesWithinGroups() - SumOfSquaredDeviatesForSubjects();
      }

      public NUMBER SumOfSquaredDeviatesForSubjects()
      {
         Group grpSubjects = new Group();
         NUMBER decSumOfScores, decSumOfSquaredScores;
         for (Int32 i = 0; i < N; i++)
         {
            decSumOfScores = 0.0m;
            for (Int32 j = 0; j < k(); j++)
            {
               decSumOfScores += groupsArray[j][i];
            }
            grpSubjects.Add(decSumOfScores * decSumOfScores);
         }
         decSumOfSquaredScores = 0.0m;

         for (Int32 i = 0; i < N; i++)
            decSumOfSquaredScores += grpSubjects[i];
         NUMBER decSumOfSquaredDeviatesForSubjects;
         try
         {
         NUMBER dec1stTerm = decSumOfSquaredScores / k();

         NUMBER dec2ndTerm = (elementsANOVAGroup.Sum() * elementsANOVAGroup.Sum())
            / elementsANOVAGroup.n();

            decSumOfSquaredDeviatesForSubjects = dec1stTerm - dec2ndTerm;
         }
         catch (Exception e)
         {
            throw new BiostatisticsException(e.Source + ": " + e.Message, e);
         }
         return decSumOfSquaredDeviatesForSubjects;
      }
   }
}

---