/Epi Info 7/Epi.UnitTests/AnalysisTest.cs

# · C# · 501 lines · 490 code · 11 blank · 0 comment · 17 complexity · 65687cf1dd977772797534c034687f6e MD5 · raw file 

    

  1. using System;

    2. 

  2. using Microsoft.VisualStudio.TestTools.UnitTesting;

    3. 

  3. using System.IO;

    4. 

  4. using Epi.Windows;

    5. 

  5. using EpiDashboard;

    6. 

  6. using Epi.Core;

    7. 

  7. using StatisticsRepository;

    8. 

  8. using System.Collections.Generic;

    9. 

  9. using System.Data;

    10. 

  10. 

    11. 

  11. namespace Epi.UnitTests

    12. 

  12. {

    13. 

  13.     [TestClass]

    14. 

  14.     public class AnalysisTest : StatisticsTests

    15. 

  15.     {

    16. 

  16.         [TestMethod]

    17. 

  17.         public void MeansTestTwoValues()

    18. 

  18.         {

    19. 

  19.             int iterations = 100;

    20. 

  20. 

    21. 

  21.             var reader = new StreamReader(File.OpenRead(@"..\..\Data\MeansTwoValueResults.csv"));

    22. 

  22.             int linesinfile = iterations;

    23. 

  23.             double[] obs0s = new double[linesinfile];

    24. 

  24.             double[] obs1s = new double[linesinfile];

    25. 

  25.             double[] total0s = new double[linesinfile];

    26. 

  26.             double[] total1s = new double[linesinfile];

    27. 

  27.             double[] mean0s = new double[linesinfile];

    28. 

  28.             double[] mean1s = new double[linesinfile];

    29. 

  29.             double[] var0s = new double[linesinfile];

    30. 

  30.             double[] var1s = new double[linesinfile];

    31. 

  31.             double[] sd0s = new double[linesinfile];

    32. 

  32.             double[] sd1s = new double[linesinfile];

    33. 

  33.             double[] min0s = new double[linesinfile];

    34. 

  34.             double[] min1s = new double[linesinfile];

    35. 

  35.             double[] q10s = new double[linesinfile];

    36. 

  36.             double[] q11s = new double[linesinfile];

    37. 

  37.             double[] median0s = new double[linesinfile];

    38. 

  38.             double[] median1s = new double[linesinfile];

    39. 

  39.             double[] q30s = new double[linesinfile];

    40. 

  40.             double[] q31s = new double[linesinfile];

    41. 

  41.             double[] max0s = new double[linesinfile];

    42. 

  42.             double[] max1s = new double[linesinfile];

    43. 

  43.             double[] mode0s = new double[linesinfile];

    44. 

  44.             double[] mode1s = new double[linesinfile];

    45. 

  45.             double[] poolmeans = new double[linesinfile];

    46. 

  46.             double[] poollcls = new double[linesinfile];

    47. 

  47.             double[] poolucls = new double[linesinfile];

    48. 

  48.             double[] sattermeans = new double[linesinfile];

    49. 

  49.             double[] satterlcls = new double[linesinfile];

    50. 

  50.             double[] satterucls = new double[linesinfile];

    51. 

  51.             double[] sds = new double[linesinfile];

    52. 

  52.             double[] equaldfs = new double[linesinfile];

    53. 

  53.             double[] equalts = new double[linesinfile];

    54. 

  54.             double[] equaltps = new double[linesinfile];

    55. 

  55.             double[] unequaldfs = new double[linesinfile];

    56. 

  56.             double[] unequalts = new double[linesinfile];

    57. 

  57.             double[] unequaltps = new double[linesinfile];

    58. 

  58.             double[] betweensss = new double[linesinfile];

    59. 

  59.             double[] withinsss = new double[linesinfile];

    60. 

  60.             double[] totalsss = new double[linesinfile];

    61. 

  61.             double[] betweendfs = new double[linesinfile];

    62. 

  62.             double[] withindfs = new double[linesinfile];

    63. 

  63.             double[] totaldfs = new double[linesinfile];

    64. 

  64.             double[] betweenmss = new double[linesinfile];

    65. 

  65.             double[] withinmss = new double[linesinfile];

    66. 

  66.             double[] fanovas = new double[linesinfile];

    67. 

  67.             double[] pfanovas = new double[linesinfile];

    68. 

  68.             double[] dfbarts = new double[linesinfile];

    69. 

  69.             double[] chibarts = new double[linesinfile];

    70. 

  70.             double[] chipbarts = new double[linesinfile];

    71. 

  71.             double[] chikwts = new double[linesinfile];

    72. 

  72.             double[] dfkwts = new double[linesinfile];

    73. 

  73.             double[] pchikwts = new double[linesinfile];

    74. 

  74.             int i = 0;

    75. 

  75.             while (!reader.EndOfStream)

    76. 

  76.             {

    77. 

  77.                 var line = reader.ReadLine();

    78. 

  78.                 var values = line.Split(',');

    79. 

  79.                 bool rv = Double.TryParse(values[1], out obs0s[i]);

    80. 

  80.                 rv = Double.TryParse(values[2], out obs1s[i]);

    81. 

  81.                 rv = Double.TryParse(values[3], out total0s[i]);

    82. 

  82.                 rv = Double.TryParse(values[4], out total1s[i]);

    83. 

  83.                 rv = Double.TryParse(values[5], out mean0s[i]);

    84. 

  84.                 rv = Double.TryParse(values[6], out mean1s[i]);

    85. 

  85.                 rv = Double.TryParse(values[7], out var0s[i]);

    86. 

  86.                 rv = Double.TryParse(values[8], out var1s[i]);

    87. 

  87.                 rv = Double.TryParse(values[9], out sd0s[i]);

    88. 

  88.                 rv = Double.TryParse(values[10], out sd1s[i]);

    89. 

  89.                 rv = Double.TryParse(values[11], out min0s[i]);

    90. 

  90.                 rv = Double.TryParse(values[12], out min1s[i]);

    91. 

  91.                 rv = Double.TryParse(values[13], out q10s[i]);

    92. 

  92.                 rv = Double.TryParse(values[14], out q11s[i]);

    93. 

  93.                 rv = Double.TryParse(values[15], out median0s[i]);

    94. 

  94.                 rv = Double.TryParse(values[16], out median1s[i]);

    95. 

  95.                 rv = Double.TryParse(values[17], out q30s[i]);

    96. 

  96.                 rv = Double.TryParse(values[18], out q31s[i]);

    97. 

  97.                 rv = Double.TryParse(values[19], out max0s[i]);

    98. 

  98.                 rv = Double.TryParse(values[20], out max1s[i]);

    99. 

  99.                 rv = Double.TryParse(values[21], out mode0s[i]);

    100. 

  100.                 rv = Double.TryParse(values[22], out mode1s[i]);

    101. 

  101.                 rv = Double.TryParse(values[23], out poolmeans[i]);

    102. 

  102.                 rv = Double.TryParse(values[24], out poollcls[i]);

    103. 

  103.                 rv = Double.TryParse(values[25], out poolucls[i]);

    104. 

  104.                 rv = Double.TryParse(values[26], out sattermeans[i]);

    105. 

  105.                 rv = Double.TryParse(values[27], out satterlcls[i]);

    106. 

  106.                 rv = Double.TryParse(values[28], out satterucls[i]);

    107. 

  107.                 rv = Double.TryParse(values[29], out sds[i]);

    108. 

  108.                 rv = Double.TryParse(values[30], out equaldfs[i]);

    109. 

  109.                 rv = Double.TryParse(values[31], out equalts[i]);

    110. 

  110.                 rv = Double.TryParse(values[32], out equaltps[i]);

    111. 

  111.                 rv = Double.TryParse(values[33], out unequaldfs[i]);

    112. 

  112.                 rv = Double.TryParse(values[34], out unequalts[i]);

    113. 

  113.                 rv = Double.TryParse(values[35], out unequaltps[i]);

    114. 

  114.                 rv = Double.TryParse(values[36], out betweensss[i]);

    115. 

  115.                 rv = Double.TryParse(values[37], out withinsss[i]);

    116. 

  116.                 rv = Double.TryParse(values[38], out totalsss[i]);

    117. 

  117.                 rv = Double.TryParse(values[39], out betweendfs[i]);

    118. 

  118.                 rv = Double.TryParse(values[40], out withindfs[i]);

    119. 

  119.                 rv = Double.TryParse(values[41], out totaldfs[i]);

    120. 

  120.                 rv = Double.TryParse(values[42], out betweenmss[i]);

    121. 

  121.                 rv = Double.TryParse(values[43], out withinmss[i]);

    122. 

  122.                 rv = Double.TryParse(values[44], out fanovas[i]);

    123. 

  123.                 rv = Double.TryParse(values[45], out pfanovas[i]);

    124. 

  124.                 rv = Double.TryParse(values[46], out dfbarts[i]);

    125. 

  125.                 rv = Double.TryParse(values[47], out chibarts[i]);

    126. 

  126.                 rv = Double.TryParse(values[48], out chipbarts[i]);

    127. 

  127.                 rv = Double.TryParse(values[49], out chikwts[i]);

    128. 

  128.                 rv = Double.TryParse(values[50], out dfkwts[i]);

    129. 

  129.                 rv = Double.TryParse(values[51], out pchikwts[i]);

    130. 

  130.                 i++;

    131. 

  131.                 if (i == linesinfile)

    132. 

  132.                     break;

    133. 

  133.             }

    134. 

  134. 

    135. 

  135.             EpiDashboard.MeansParameters Parameters = new MeansParameters();

    136. 

  136.             Parameters.CrosstabVariableName = "crosstabvar2";

    137. 

  137.             Parameters.ColumnNames.Add("meanvar");

    138. 

  138.             Parameters.ColumnNames.Add("crosstabvar2");

    139. 

  139.             Parameters.IncludeFullSummaryStatistics = true;

    140. 

  140. 

    141. 

  141.             for (int j = 0; j < iterations; j++)

    142. 

  142.             {

    143. 

  143.                 getData("MeansTestData", "iteration = " + j);

    144. 

  144.                 Parameters.CustomFilter = "iteration = " + j;

    145. 

  145.                 Dictionary<DataTable, List<DescriptiveStatistics>> stratifiedFrequencyTables = DashboardHelper.GenerateFrequencyTable(Parameters);

    146. 

  146.                 foreach (object key in stratifiedFrequencyTables.Keys)

    147. 

  147.                 {

    148. 

  148.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].observations, obs0s[j]);

    149. 

  149.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][1].observations, obs1s[j]);

    150. 

  150.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].sum.Value, total0s[j], 0.0001);

    151. 

  151.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][1].sum.Value, total1s[j], 0.0001);

    152. 

  152.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].mean.Value, mean0s[j], 0.0001);

    153. 

  153.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][1].mean.Value, mean1s[j], 0.0001);

    154. 

  154.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].variance.Value, var0s[j], 0.0001);

    155. 

  155.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][1].variance.Value, var1s[j], 0.0001);

    156. 

  156.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].stdDev.Value, sd0s[j], 0.0001);

    157. 

  157.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][1].stdDev.Value, sd1s[j], 0.0001);

    158. 

  158.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].min.Value, min0s[j], 0.0001);

    159. 

  159.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][1].min.Value, min1s[j], 0.0001);

    160. 

  160.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].q1.Value, q10s[j], 0.0001);

    161. 

  161.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][1].q1.Value, q11s[j], 0.0001);

    162. 

  162.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].median.Value, median0s[j], 0.0001);

    163. 

  163.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][1].median.Value, median1s[j], 0.0001);

    164. 

  164.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].q3.Value, q30s[j], 0.0001);

    165. 

  165.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][1].q3.Value, q31s[j], 0.0001);

    166. 

  166.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].max.Value, max0s[j], 0.0001);

    167. 

  167.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][1].max.Value, max1s[j], 0.0001);

    168. 

  168.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].mode.Value, mode0s[j], 0.0001);

    169. 

  169.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][1].mode.Value, mode1s[j], 0.0001);

    170. 

  170.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].meansDiff, poolmeans[j], 0.0001);

    171. 

  171.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].equalLCLMean, poollcls[j], 0.0001);

    172. 

  172.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].equalUCLMean, poolucls[j], 0.0001);

    173. 

  173.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].meansDiff, sattermeans[j], 0.0001);

    174. 

  174.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].unequalLCLMean, satterlcls[j], 0.001);

    175. 

  175.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].unequalUCLMean, satterucls[j], 0.001);

    176. 

  176.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].stdDevDiff, sds[j], 0.0001);

    177. 

  177.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].df, equaldfs[j], 0.0001);

    178. 

  178.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].tStatistic, equalts[j], 0.0001);

    179. 

  179.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].pEqual, equaltps[j], 0.0001);

    180. 

  180.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].SatterthwaiteDF, unequaldfs[j], 0.0001);

    181. 

  181.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].tStatisticUnequal, unequalts[j], 0.0001);

    182. 

  182.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].pUneqal, unequaltps[j], 0.0001);

    183. 

  183.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].ssBetween.Value, betweensss[j], 0.0001);

    184. 

  184.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].ssWithin.Value, withinsss[j], 0.0001);

    185. 

  185.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].ssBetween.Value + stratifiedFrequencyTables[(System.Data.DataTable)key][0].ssWithin.Value, totalsss[j], 0.0001);

    186. 

  186.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].dfBetween.Value, betweendfs[j], 0.0001);

    187. 

  187.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].dfWithin.Value, withindfs[j], 0.0001);

    188. 

  188.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].dfBetween.Value + stratifiedFrequencyTables[(System.Data.DataTable)key][0].dfWithin.Value, totaldfs[j], 0.0001);

    189. 

  189.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].msBetween.Value, betweenmss[j], 0.0001);

    190. 

  190.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].msWithin.Value, withinmss[j], 0.0001);

    191. 

  191.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].fStatistic.Value, fanovas[j], 0.0001);

    192. 

  192.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].anovaPValue.Value, pfanovas[j], 0.0001);

    193. 

  193.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].dfBetween.Value, dfbarts[j], 0.0001);

    194. 

  194.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].chiSquare.Value, chibarts[j], 0.0001);

    195. 

  195.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].bartlettPValue.Value, chipbarts[j], 0.0001);

    196. 

  196.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].kruskalWallisH.Value, chikwts[j], 0.0001);

    197. 

  197.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].dfBetween.Value, dfkwts[j], 0.0001);

    198. 

  198.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].kruskalPValue.Value, pchikwts[j], 0.0001);

    199. 

  199.                 }

    200. 

  200.             }

    201. 

  201.         }

    202. 

  202. 

    203. 

  203.         [TestMethod]

    204. 

  204.         public void MeansTestFourValues()

    205. 

  205.         {

    206. 

  206.             int iterations = 100;

    207. 

  207. 

    208. 

  208.             var reader = new StreamReader(File.OpenRead(@"..\..\Data\MeansFourValueResults.csv"));

    209. 

  209.             int linesinfile = iterations;

    210. 

  210.             double[] obs0s = new double[linesinfile];

    211. 

  211.             double[] obs1s = new double[linesinfile];

    212. 

  212.             double[] obs2s = new double[linesinfile];

    213. 

  213.             double[] obs3s = new double[linesinfile];

    214. 

  214.             double[] total0s = new double[linesinfile];

    215. 

  215.             double[] total1s = new double[linesinfile];

    216. 

  216.             double[] total2s = new double[linesinfile];

    217. 

  217.             double[] total3s = new double[linesinfile];

    218. 

  218.             double[] mean0s = new double[linesinfile];

    219. 

  219.             double[] mean1s = new double[linesinfile];

    220. 

  220.             double[] mean2s = new double[linesinfile];

    221. 

  221.             double[] mean3s = new double[linesinfile];

    222. 

  222.             double[] var0s = new double[linesinfile];

    223. 

  223.             double[] var1s = new double[linesinfile];

    224. 

  224.             double[] var2s = new double[linesinfile];

    225. 

  225.             double[] var3s = new double[linesinfile];

    226. 

  226.             double[] sd0s = new double[linesinfile];

    227. 

  227.             double[] sd1s = new double[linesinfile];

    228. 

  228.             double[] sd2s = new double[linesinfile];

    229. 

  229.             double[] sd3s = new double[linesinfile];

    230. 

  230.             double[] min0s = new double[linesinfile];

    231. 

  231.             double[] min1s = new double[linesinfile];

    232. 

  232.             double[] min2s = new double[linesinfile];

    233. 

  233.             double[] min3s = new double[linesinfile];

    234. 

  234.             double[] q10s = new double[linesinfile];

    235. 

  235.             double[] q11s = new double[linesinfile];

    236. 

  236.             double[] q12s = new double[linesinfile];

    237. 

  237.             double[] q13s = new double[linesinfile];

    238. 

  238.             double[] median0s = new double[linesinfile];

    239. 

  239.             double[] median1s = new double[linesinfile];

    240. 

  240.             double[] median2s = new double[linesinfile];

    241. 

  241.             double[] median3s = new double[linesinfile];

    242. 

  242.             double[] q30s = new double[linesinfile];

    243. 

  243.             double[] q31s = new double[linesinfile];

    244. 

  244.             double[] q32s = new double[linesinfile];

    245. 

  245.             double[] q33s = new double[linesinfile];

    246. 

  246.             double[] max0s = new double[linesinfile];

    247. 

  247.             double[] max1s = new double[linesinfile];

    248. 

  248.             double[] max2s = new double[linesinfile];

    249. 

  249.             double[] max3s = new double[linesinfile];

    250. 

  250.             double[] mode0s = new double[linesinfile];

    251. 

  251.             double[] mode1s = new double[linesinfile];

    252. 

  252.             double[] mode2s = new double[linesinfile];

    253. 

  253.             double[] mode3s = new double[linesinfile];

    254. 

  254.             double[] betweensss = new double[linesinfile];

    255. 

  255.             double[] withinsss = new double[linesinfile];

    256. 

  256.             double[] totalsss = new double[linesinfile];

    257. 

  257.             double[] betweendfs = new double[linesinfile];

    258. 

  258.             double[] withindfs = new double[linesinfile];

    259. 

  259.             double[] totaldfs = new double[linesinfile];

    260. 

  260.             double[] betweenmss = new double[linesinfile];

    261. 

  261.             double[] withinmss = new double[linesinfile];

    262. 

  262.             double[] fanovas = new double[linesinfile];

    263. 

  263.             double[] pfanovas = new double[linesinfile];

    264. 

  264.             double[] dfbarts = new double[linesinfile];

    265. 

  265.             double[] chibarts = new double[linesinfile];

    266. 

  266.             double[] chipbarts = new double[linesinfile];

    267. 

  267.             double[] chikwts = new double[linesinfile];

    268. 

  268.             double[] dfkwts = new double[linesinfile];

    269. 

  269.             double[] pchikwts = new double[linesinfile];

    270. 

  270.             int i = 0;

    271. 

  271.             while (!reader.EndOfStream)

    272. 

  272.             {

    273. 

  273.                 var line = reader.ReadLine();

    274. 

  274.                 var values = line.Split(',');

    275. 

  275.                 bool rv = Double.TryParse(values[1], out obs0s[i]);

    276. 

  276.                 rv = Double.TryParse(values[2], out obs1s[i]);

    277. 

  277.                 rv = Double.TryParse(values[3], out obs2s[i]);

    278. 

  278.                 rv = Double.TryParse(values[4], out obs3s[i]);

    279. 

  279.                 rv = Double.TryParse(values[5], out total0s[i]);

    280. 

  280.                 rv = Double.TryParse(values[6], out total1s[i]);

    281. 

  281.                 rv = Double.TryParse(values[7], out total2s[i]);

    282. 

  282.                 rv = Double.TryParse(values[8], out total3s[i]);

    283. 

  283.                 rv = Double.TryParse(values[9], out mean0s[i]);

    284. 

  284.                 rv = Double.TryParse(values[10], out mean1s[i]);

    285. 

  285.                 rv = Double.TryParse(values[11], out mean2s[i]);

    286. 

  286.                 rv = Double.TryParse(values[12], out mean3s[i]);

    287. 

  287.                 rv = Double.TryParse(values[13], out var0s[i]);

    288. 

  288.                 rv = Double.TryParse(values[14], out var1s[i]);

    289. 

  289.                 rv = Double.TryParse(values[15], out var2s[i]);

    290. 

  290.                 rv = Double.TryParse(values[16], out var3s[i]);

    291. 

  291.                 rv = Double.TryParse(values[17], out sd0s[i]);

    292. 

  292.                 rv = Double.TryParse(values[18], out sd1s[i]);

    293. 

  293.                 rv = Double.TryParse(values[19], out sd2s[i]);

    294. 

  294.                 rv = Double.TryParse(values[20], out sd3s[i]);

    295. 

  295.                 rv = Double.TryParse(values[21], out min0s[i]);

    296. 

  296.                 rv = Double.TryParse(values[22], out min1s[i]);

    297. 

  297.                 rv = Double.TryParse(values[23], out min2s[i]);

    298. 

  298.                 rv = Double.TryParse(values[24], out min3s[i]);

    299. 

  299.                 rv = Double.TryParse(values[25], out q10s[i]);

    300. 

  300.                 rv = Double.TryParse(values[26], out q11s[i]);

    301. 

  301.                 rv = Double.TryParse(values[27], out q12s[i]);

    302. 

  302.                 rv = Double.TryParse(values[28], out q13s[i]);

    303. 

  303.                 rv = Double.TryParse(values[29], out median0s[i]);

    304. 

  304.                 rv = Double.TryParse(values[30], out median1s[i]);

    305. 

  305.                 rv = Double.TryParse(values[31], out median2s[i]);

    306. 

  306.                 rv = Double.TryParse(values[32], out median3s[i]);

    307. 

  307.                 rv = Double.TryParse(values[33], out q30s[i]);

    308. 

  308.                 rv = Double.TryParse(values[34], out q31s[i]);

    309. 

  309.                 rv = Double.TryParse(values[35], out q32s[i]);

    310. 

  310.                 rv = Double.TryParse(values[36], out q33s[i]);

    311. 

  311.                 rv = Double.TryParse(values[37], out max0s[i]);

    312. 

  312.                 rv = Double.TryParse(values[38], out max1s[i]);

    313. 

  313.                 rv = Double.TryParse(values[39], out max2s[i]);

    314. 

  314.                 rv = Double.TryParse(values[40], out max3s[i]);

    315. 

  315.                 rv = Double.TryParse(values[41], out mode0s[i]);

    316. 

  316.                 rv = Double.TryParse(values[42], out mode1s[i]);

    317. 

  317.                 rv = Double.TryParse(values[43], out mode2s[i]);

    318. 

  318.                 rv = Double.TryParse(values[44], out mode3s[i]);

    319. 

  319.                 rv = Double.TryParse(values[45], out betweensss[i]);

    320. 

  320.                 rv = Double.TryParse(values[46], out withinsss[i]);

    321. 

  321.                 rv = Double.TryParse(values[47], out totalsss[i]);

    322. 

  322.                 rv = Double.TryParse(values[48], out betweendfs[i]);

    323. 

  323.                 rv = Double.TryParse(values[49], out withindfs[i]);

    324. 

  324.                 rv = Double.TryParse(values[50], out totaldfs[i]);

    325. 

  325.                 rv = Double.TryParse(values[51], out betweenmss[i]);

    326. 

  326.                 rv = Double.TryParse(values[52], out withinmss[i]);

    327. 

  327.                 rv = Double.TryParse(values[53], out fanovas[i]);

    328. 

  328.                 rv = Double.TryParse(values[54], out pfanovas[i]);

    329. 

  329.                 rv = Double.TryParse(values[55], out dfbarts[i]);

    330. 

  330.                 rv = Double.TryParse(values[56], out chibarts[i]);

    331. 

  331.                 rv = Double.TryParse(values[57], out chipbarts[i]);

    332. 

  332.                 rv = Double.TryParse(values[58], out chikwts[i]);

    333. 

  333.                 rv = Double.TryParse(values[59], out dfkwts[i]);

    334. 

  334.                 rv = Double.TryParse(values[60], out pchikwts[i]);

    335. 

  335.                 i++;

    336. 

  336.                 if (i == linesinfile)

    337. 

  337.                     break;

    338. 

  338.             }

    339. 

  339. 

    340. 

  340.             EpiDashboard.MeansParameters Parameters = new MeansParameters();

    341. 

  341.             Parameters.CrosstabVariableName = "crosstabvar4";

    342. 

  342.             Parameters.ColumnNames.Add("meanvar");

    343. 

  343.             Parameters.ColumnNames.Add("crosstabvar4");

    344. 

  344.             Parameters.IncludeFullSummaryStatistics = true;

    345. 

  345. 

    346. 

  346.             for (int j = 0; j < iterations; j++)

    347. 

  347.             {

    348. 

  348.                 getData("MeansTestData", "iteration = " + j);

    349. 

  349.                 Parameters.CustomFilter = "iteration = " + j;

    350. 

  350.                 Dictionary<DataTable, List<DescriptiveStatistics>> stratifiedFrequencyTables = DashboardHelper.GenerateFrequencyTable(Parameters);

    351. 

  351.                 foreach (object key in stratifiedFrequencyTables.Keys)

    352. 

  352.                 {

    353. 

  353.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].observations, obs0s[j]);

    354. 

  354.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][1].observations, obs1s[j]);

    355. 

  355.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][2].observations, obs2s[j]);

    356. 

  356.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][3].observations, obs3s[j]);

    357. 

  357.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].sum.Value, total0s[j], 0.0001);

    358. 

  358.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][1].sum.Value, total1s[j], 0.0001);

    359. 

  359.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][2].sum.Value, total2s[j], 0.0001);

    360. 

  360.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][3].sum.Value, total3s[j], 0.0001);

    361. 

  361.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].mean.Value, mean0s[j], 0.0001);

    362. 

  362.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][1].mean.Value, mean1s[j], 0.0001);

    363. 

  363.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][2].mean.Value, mean2s[j], 0.0001);

    364. 

  364.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][3].mean.Value, mean3s[j], 0.0001);

    365. 

  365.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].variance.Value, var0s[j], 0.0001);

    366. 

  366.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][1].variance.Value, var1s[j], 0.0001);

    367. 

  367.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][2].variance.Value, var2s[j], 0.0001);

    368. 

  368.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][3].variance.Value, var3s[j], 0.0001);

    369. 

  369.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].stdDev.Value, sd0s[j], 0.0001);

    370. 

  370.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][1].stdDev.Value, sd1s[j], 0.0001);

    371. 

  371.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][2].stdDev.Value, sd2s[j], 0.0001);

    372. 

  372.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][3].stdDev.Value, sd3s[j], 0.0001);

    373. 

  373.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].min.Value, min0s[j], 0.0001);

    374. 

  374.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][1].min.Value, min1s[j], 0.0001);

    375. 

  375.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][2].min.Value, min2s[j], 0.0001);

    376. 

  376.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][3].min.Value, min3s[j], 0.0001);

    377. 

  377.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].q1.Value, q10s[j], 0.0001);

    378. 

  378.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][1].q1.Value, q11s[j], 0.0001);

    379. 

  379.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][2].q1.Value, q12s[j], 0.0001);

    380. 

  380.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][3].q1.Value, q13s[j], 0.0001);

    381. 

  381.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].median.Value, median0s[j], 0.0001);

    382. 

  382.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][1].median.Value, median1s[j], 0.0001);

    383. 

  383.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][2].median.Value, median2s[j], 0.0001);

    384. 

  384.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][3].median.Value, median3s[j], 0.0001);

    385. 

  385.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].q3.Value, q30s[j], 0.0001);

    386. 

  386.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][1].q3.Value, q31s[j], 0.0001);

    387. 

  387.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][2].q3.Value, q32s[j], 0.0001);

    388. 

  388.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][3].q3.Value, q33s[j], 0.0001);

    389. 

  389.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].max.Value, max0s[j], 0.0001);

    390. 

  390.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][1].max.Value, max1s[j], 0.0001);

    391. 

  391.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][2].max.Value, max2s[j], 0.0001);

    392. 

  392.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][3].max.Value, max3s[j], 0.0001);

    393. 

  393.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].mode.Value, mode0s[j], 0.0001);

    394. 

  394.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][1].mode.Value, mode1s[j], 0.0001);

    395. 

  395.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][2].mode.Value, mode2s[j], 0.0001);

    396. 

  396.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][3].mode.Value, mode3s[j], 0.0001);

    397. 

  397.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].ssBetween.Value, betweensss[j], 0.0001);

    398. 

  398.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].ssWithin.Value, withinsss[j], 0.0001);

    399. 

  399.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].ssBetween.Value + stratifiedFrequencyTables[(System.Data.DataTable)key][0].ssWithin.Value, totalsss[j], 0.0001);

    400. 

  400.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].dfBetween.Value, betweendfs[j], 0.0001);

    401. 

  401.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].dfWithin.Value, withindfs[j], 0.0001);

    402. 

  402.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].dfBetween.Value + stratifiedFrequencyTables[(System.Data.DataTable)key][0].dfWithin.Value, totaldfs[j], 0.0001);

    403. 

  403.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].msBetween.Value, betweenmss[j], 0.0001);

    404. 

  404.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].msWithin.Value, withinmss[j], 0.0001);

    405. 

  405.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].fStatistic.Value, fanovas[j], 0.0001);

    406. 

  406.                     if (pfanovas[j] < 0.001)

    407. 

  407.                         Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].anovaPValue.Value, pfanovas[j], 0.0001);

    408. 

  408.                     else if (pfanovas[j] < 0.01)

    409. 

  409.                         Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].anovaPValue.Value, pfanovas[j], 0.001);

    410. 

  410.                     else if (pfanovas[j] < 0.1)

    411. 

  411.                         Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].anovaPValue.Value, pfanovas[j], 0.01);

    412. 

  412.                     else

    413. 

  413.                         Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].anovaPValue.Value, pfanovas[j], 0.1);

    414. 

  414.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].dfBetween.Value, dfbarts[j], 0.0001);

    415. 

  415.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].chiSquare.Value, chibarts[j], 0.0001);

    416. 

  416.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].bartlettPValue.Value, chipbarts[j], 0.0001);

    417. 

  417.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].kruskalWallisH.Value, chikwts[j], 0.0001);

    418. 

  418.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].dfBetween.Value, dfkwts[j], 0.0001);

    419. 

  419.                     Assert.AreEqual(stratifiedFrequencyTables[(System.Data.DataTable)key][0].kruskalPValue.Value, pchikwts[j], 0.0001);

    420. 

  420.                 }

    421. 

  421.             }

    422. 

  422.         }

    423. 

  423. 

    424. 

  424.         [TestMethod]

    425. 

  425.         public void TwoByTwoTest()

    426. 

  426.         {

    427. 

  427.             int iterations = 100;

    428. 

  428.             StatisticsRepository.cTable cTable = new StatisticsRepository.cTable();

    429. 

  429. 

    430. 

  430.             var reader = new StreamReader(File.OpenRead(@"..\..\Data\TwoBy2Stats.csv"));

    431. 

  431.             int linesinfile = iterations;

    432. 

  432.             double[] chisqs = new double[linesinfile];

    433. 

  433.             double[] chisqps = new double[linesinfile];

    434. 

  434.             double[] mhchisqs = new double[linesinfile];

    435. 

  435.             double[] mhchisqps = new double[linesinfile];

    436. 

  436.             double[] yateschisqs = new double[linesinfile];

    437. 

  437.             double[] yateschisqps = new double[linesinfile];

    438. 

  438.             double[] ors = new double[linesinfile];

    439. 

  439.             double[] orlowers = new double[linesinfile];

    440. 

  440.             double[] oruppers = new double[linesinfile];

    441. 

  441.             double[] rrs = new double[linesinfile];

    442. 

  442.             double[] rrlowers = new double[linesinfile];

    443. 

  443.             double[] rruppers = new double[linesinfile];

    444. 

  444.             double[] fisherlowers = new double[linesinfile];

    445. 

  445.             double[] fisheruppers = new double[linesinfile];

    446. 

  446.             double[] fisher1s = new double[linesinfile];

    447. 

  447.             double[] fisher2s = new double[linesinfile];

    448. 

  448.             int i = 0;

    449. 

  449.             while (!reader.EndOfStream)

    450. 

  450.             {

    451. 

  451.                 var line = reader.ReadLine();

    452. 

  452.                 var values = line.Split(',');

    453. 

  453.                 bool rv = Double.TryParse(values[1], out chisqs[i]);

    454. 

  454.                 rv = Double.TryParse(values[2], out chisqps[i]);

    455. 

  455.                 rv = Double.TryParse(values[3], out mhchisqs[i]);

    456. 

  456.                 rv = Double.TryParse(values[4], out mhchisqps[i]);

    457. 

  457.                 rv = Double.TryParse(values[5], out yateschisqs[i]);

    458. 

  458.                 rv = Double.TryParse(values[6], out yateschisqps[i]);

    459. 

  459.                 rv = Double.TryParse(values[7], out ors[i]);

    460. 

  460.                 rv = Double.TryParse(values[8], out orlowers[i]);

    461. 

  461.                 rv = Double.TryParse(values[9], out oruppers[i]);

    462. 

  462.                 rv = Double.TryParse(values[10], out rrs[i]);

    463. 

  463.                 rv = Double.TryParse(values[11], out rrlowers[i]);

    464. 

  464.                 rv = Double.TryParse(values[12], out rruppers[i]);

    465. 

  465.                 rv = Double.TryParse(values[13], out fisherlowers[i]);

    466. 

  466.                 rv = Double.TryParse(values[14], out fisheruppers[i]);

    467. 

  467.                 rv = Double.TryParse(values[15], out fisher1s[i]);

    468. 

  468.                 rv = Double.TryParse(values[16], out fisher2s[i]);

    469. 

  469.                 i++;

    470. 

  470.                 if (i == linesinfile)

    471. 

  471.                     break;

    472. 

  472.             }

    473. 

  473. 

    474. 

  474.             for (int j = 0; j < iterations; j++)

    475. 

  475.             {

    476. 

  476.                 getData("LogisticTestData", "iteration = " + j);

    477. 

  477.                 int yy = mainTable.Select("depvar = 1 and indepvar1 = 1").Length;

    478. 

  478.                 int yn = mainTable.Select("depvar = 1 and indepvar1 = 0").Length;

    479. 

  479.                 int ny = mainTable.Select("depvar = 0 and indepvar1 = 1").Length;

    480. 

  480.                 int nn= mainTable.Select("depvar = 0 and indepvar1 = 0").Length;

    481. 

  481.                 StatisticsRepository.cTable.SingleTableResults singleTableResults = cTable.SigTable((double)yy, (double)yn, (double)ny, (double)nn, 0.95);

    482. 

  482.                 Assert.AreEqual(singleTableResults.ChiSquareUncorrectedVal, chisqs[j], 0.0001);

    483. 

  483.                 Assert.AreEqual(singleTableResults.ChiSquareUncorrected2P, chisqps[j], 0.0001);

    484. 

  484.                 Assert.AreEqual(singleTableResults.ChiSquareMantelVal, mhchisqs[j], 0.0001);

    485. 

  485.                 Assert.AreEqual(singleTableResults.ChiSquareMantel2P, mhchisqps[j], 0.0001);

    486. 

  486.                 Assert.AreEqual(singleTableResults.ChiSquareYatesVal, yateschisqs[j], 0.0001);

    487. 

  487.                 Assert.AreEqual(singleTableResults.ChiSquareYates2P, yateschisqps[j], 0.0001);

    488. 

  488.                 Assert.AreEqual(singleTableResults.OddsRatioEstimate.Value, ors[j], 0.0001);

    489. 

  489.                 Assert.AreEqual(singleTableResults.OddsRatioLower.Value, orlowers[j], 0.0001);

    490. 

  490.                 Assert.AreEqual(singleTableResults.OddsRatioUpper.Value, oruppers[j], 0.001);

    491. 

  491.                 Assert.AreEqual(singleTableResults.RiskRatioEstimate.Value, rrs[j], 0.0001);

    492. 

  492.                 Assert.AreEqual(singleTableResults.RiskRatioLower.Value, rrlowers[j], 0.0001);

    493. 

  493.                 Assert.AreEqual(singleTableResults.RiskRatioUpper.Value, rruppers[j], 0.0001);

    494. 

  494.                 Assert.AreEqual(singleTableResults.OddsRatioMLEFisherLower, fisherlowers[j], 0.0001);

    495. 

  495.                 Assert.AreEqual(singleTableResults.OddsRatioMLEFisherUpper, fisheruppers[j], 0.0001);

    496. 

  496.                 Assert.AreEqual(singleTableResults.FisherExactP, fisher1s[j], 0.0001);

    497. 

  497.                 Assert.AreEqual(singleTableResults.FisherExact2P, fisher2s[j], 0.0001);

    498. 

  498.             }

    499. 

  499.         }

    500. 

  500.     }

    501. 

  501. }
    502.