<?php

     if($ZZ1A[$j] != 0 AND $ZZ1A[$j] != $blink){
           $answer1[$j] = $vanZ1[$j]/$ZZ1A[$j];
           }

           if($secondZZ1A[$j] != 0 AND $secondZZ1A[$j] != $blink){
            $Sanswer1[$j] = $secondvanZA[$j]/$secondZZ1A[$j];
}

               if($thirdZZ1A[$j] != 0 AND $thirdZZ1A[$j] != $blink){
             $Tanswer1[$j] = $thirdvanZA[$j]/$thirdZZ1A[$j];
             }
                     if($fourthZZ1A[$j] != 0 AND $fourthZZ1A[$j] != $blink){
              $Fanswer1[$j] = $fourthvanZA[$j]/$fourthZZ1A[$j];
                          }

           ////// count
           $answer1a[$j] = $ZZ1A[$j];
             $Sanswer1a[$j] = $secondZZ1A[$j];
              $Tanswer1a[$j] = $thirdZZ1A[$j];
               $Fanswer1a[$j] =$fourthZZ1A[$j];

        ///////median

          $home_values2 = explode(",", $bobv1[$j]);

                 $count = $ZZ1A[$j]; //total numbers in array
    $middleval = floor(($count-1)/2); // find the middle value, or the lowest middle value
    if($count % 2) { // odd number, middle is the median
        $median1 = $home_values2[$middleval];
          $answer1b[$j]=$median1;
    } else { // even number, calculate avg of 2 medians
    if ($count>1){
        $low = $home_values2[$middleval];
        $high = $home_values2[$middleval+1];
    //    $median1 = number_format((($low+$high)/2),2);
        $median1 = number_format((($low+$high)/2),0);
          $answer1b[$j]=$median1;

     }   }
     $bobmed=trim($answer1b[$j]);
           $we1="$findQN[$Q1]$bobmed";
     //  echo "we-$we1<br/>";
  $ansAnn1b[$j] =$QansState[$we1];

//   echo "we-$we1-  $ansAnn1b[$j]<br/>";
      //   include("medianmix.php");

  /////more math

     require_once("util_calce.php");

$vals = array(3, 4, 5, 7, 8);
//echo $bobv1[$j];
$vals1[$j]=explode(",", $bobv1[$j]);
if (count($vals1[$j])>1){
$std1[$j]=stddev($vals1[$j]);
         }
$vals2[$j]=explode(",", $bobv2[$j]);
if (count($vals2[$j])>1){
$std2[$j]=stddev($vals2[$j]);
}

 $vals3[$j]=explode(",", $bobv3[$j]);
 $toop[$j]=count($vals3[$j]);

  if ($toop[$j]>1){
$std3[$j]=stddev($vals3[$j]);
                   }

$vals4[$j]=explode(",", $bobv4[$j]);
if (count($vals4[$j])>1){
$std4[$j]=stddev($vals4[$j]);
}

$stderr1[$j] = $std1[$j] / sqrt(count($vals1[$j]));
$stderr2[$j] = $std2[$j] / sqrt(count($vals2[$j]));
$stderr3[$j] = $std3[$j] / sqrt(count($vals3[$j]));
$stderr4[$j] = $std4[$j] / sqrt(count($vals4[$j]));

$topzz=count($vals1[$j])-1;
   include("confidence.php");

//now get the CI:
$CI95 = $stderr*ttable((count($vals)-1),0.025);
$CI90a[$j] = $stderr1[$j]*ttable($topzz,0.05);
$CI95a[$j] = $stderr1[$j]*ttable($topzz,0.025);

  $topzz=count($vals2[$j])-1;
   include("confidence.php");
$CI90b[$j] = $stderr2[$j]*ttable($topzz,0.05);
$CI95b[$j] = $stderr2[$j]*ttable($topzz,0.025);

  $topzz=count($vals3[$j])-1;
   include("confidence.php");
$CI90c[$j] = $stderr3[$j]*ttable($topzz,0.05);
$CI95c[$j] = $stderr3[$j]*ttable($topzz,0.025);

  $topzz=count($vals4[$j])-1;
   include("confidence.php");
$CI90d[$j] = $stderr4[$j]*ttable($topzz,0.05);
$CI95d[$j] = $stderr4[$j]*ttable($topzz,0.025);


$c1[$j] = explode(",", $bobv1[$j]);
$tB1[$j] = explode(",", $bobv2[$j]);
$tB2[$j] = explode(",", $bobv3[$j]);
$tB3[$j] = explode(",", $bobv4[$j]);

$cc1[$j]= count($c1[$j]);
$ttB1[$j]= count($tB1[$j]);
$ttB2[$j]= count($tB2[$j]);
$ttB3[$j]= count($tB3[$j]);



$avgg[$j]=  avg2($c1[$j]);
//echo "<br/>0-$avg[$j] <br/>";
$avgg1[$j]=  avg2($tB1[$j]);
//echo "1-$avgg1[$j] <br/>";
$avgg2[$j]=  avg2($tB2[$j]);
//echo "2=$avgg2[$j] <br/>";
$avgg3[$j]=  avg2($tB3[$j]);
//echo "3-$avgg3[$j] <br/>";

$cCC1 =array($cc1[$j],$avgg[$j]);
$cCC2 =array($ttB1[$j],$avgg1[$j]);
$cCC3 =array($ttB2[$j],$avgg2[$j]);
$cCC4 =array($ttB3[$j],$avgg3[$j]);


//$tA = array(180, 45);
//$tB = array($ttB1, $avgg2[$j]);


$conv0=cr($cCC1);
$conv1=cr($cCC2);
$conv2=cr($cCC3);
$conv3=cr($cCC4);

//$conv=$tB;
  $qsamp=$cc1[$j];

  $bs = array(.05,0.0625, 0.0225, 0.0025);
   $bs11 = array(1.645,1.96, 2.326);
   $err11 = array(.05,.06, .07,.08,.09,.1);
$samplesize[$j]=ssize($conv0,$bs[0]);
$samplesize2[$j]=ssize2($std1[$j],$bs11[$bss11C],$err11[$errorC]); //change
 if ($qsamp<$samplesize2[$j]){
   //   $adviceSample[$j]="Increase your sample to:$samplesize2[$j]- $samplesize[$j]";
     $adviceSample[$j]="Increase your sample to: $samplesize2[$j]";
     }


$t=$tC;
$zscoreA=zscore($cCC1, $cCC3);
 // echo "s-$s , z-$z";
 if ($zscoreA >0 ){
$zscoreA="  The second sample is significant.<br/>";
 } else{
 $zscoreA=" The second sample is not significantly different. <br/>";
 }







     $avgg[$j]=number_format($avgg[$j],2);
     $avgg1[$j]=number_format($avgg1[$j],2);
     $avgg2[$j]=number_format($avgg2[$j],2);
     $avgg3[$j]=number_format($avgg3[$j],2);

     $std1[$j]=number_format($std1[$j],3);
     $std2[$j]=number_format($std2[$j],3);
    $std3[$j]=number_format($std3[$j],3);
      $std4[$j]=number_format($std4[$j],3);

        $CI90a[$j]=number_format($CI90a[$j],3);
         $CI90b[$j]=number_format($CI90b[$j],3);
          $CI90c[$j]=number_format($CI90c[$j],3);
           $CI90d[$j]=number_format($CI90d[$j],3);

               if ($operand2[$j]==29){
           $answer1[$j] = "Population1: $avgg[$j] <br/> vs Population1: $avgg1[$j] <br/>$zscoreA";
         }      //compare
//
 $not1[$j]="Population: Average: $avgg[$j], o1=$std1[$j], Sample size:$cc1[$j]<br/> CI90: $CI90a[$j]<br/> $adviceSample[$j]<br/>";

if($ttB1[$j]>1){
$not2[$j]= "Population 2 Average: $avgg1[$j], o2=$std2[$j], Sample Size:$ttB1[$j]<br/>CI90: $CI90b[$j]<br/>";
}
if($ttB2[$j]>1){
$not3[$j]= "Population 3 Average: $avgg2[$j], o3=$std3[$j],  Sample Size:$ttB2[$j]<br/>CI90: $CI90c[$j]<br/>" ;
}
if($ttB3[$j]>1){
$not4[$j]="Population 4 Average: $avgg3[$j], o4=$std4[$j],  Sample Size:$ttB3[$j]<br/>CI90: $CI90d[$j]<br/>";
         }
$tC = $ttB2[$j];
//${notes.$j} .="Control Average: $avgg[$j] ,  Sample size:$cc1[$j]<br/> Test Average: $avgg1[$j] ,  Sample Size:$ttB1[$j]<br/> ,  Test Average: $avgg2[$j] ,  Sample Size:$ttB2[$j]<br/>  Test Average: $avgg3[$j] ,  Sample Size:$ttB3[$j]<br/>  ";






?>
