- float varA = 4 / Atotal * (sumPcubed - a * a);
- float varB = 4 / Btotal * (sumQcubed - b * b);
- float varD = sumPQsq / Atotal + sumPsqQ / Btotal - d * d * (1/Atotal + 1/Btotal);
- float covAD = 2 / Atotal * (sumPsqQ - a * d);
- float covBD = 2 / Btotal * (sumPQsq - b* d);
+ double varA = 4 / Atotal * (sumPcubed - a * a);
+ double varB = 4 / Btotal * (sumQcubed - b * b);
+ double varD = sumPQsq / Atotal + sumPsqQ / Btotal - d * d * (1/Atotal + 1/Btotal);
+ double covAD = 2 / Atotal * (sumPsqQ - a * d);
+ double covBD = 2 / Btotal * (sumPQsq - b* d);