Gauss formule om P-waarde te berekenen

Het forum voor overige vragen betreffende wiskunde uit het hoger onderwijs.
Plaats reactie
Mary
Nieuw lid
Nieuw lid
Berichten: 2
Lid geworden op: 05 jun 2006, 11:48

Gauss formule om P-waarde te berekenen

Bericht door Mary » 05 jun 2006, 11:52

Hallo,

Ik ben bezig een programma te schrijven (informatica, programmeren) dat automatisch een T-waarde en een P-waarde uitrekent van een Standaard verdeling. Ik heb hiervoor een aantal berekeningen gebruikt van een ander programma op internet, wat ongeveer hetzelfde doet. Het werkt prima, alleen wil ik graag weten HOE het werkt. Ze gebruiken nogal wat wiskundige berekeningen en dat is waar ik afhaak. Kunnen jullie mij vertellen wat hier gebeurt?

Het programma voert de volgende berekeningen uit:
- T-waarde delen door wortel van n
- hier nemen ze de tangens van
- hier nemen ze vervolgens de cosinus van en de sinus van
- vervolgens berekenen ze de Z-waarde door ((1 * cosinus kwadraat * 1) : 2) + 1 te nemen. De P-waarde is dan 1- (sinus * Z-waarde)
- of (hangt af van de rest na deling bij n:2) ze berekenen de Z-waarde door ((1 * cosinus kwadraat * 2) : 3) + 1 te nemen. De P-waarde is dan 1- ((tangens * sinus * cosinus * Z-waarde) : (pi:2))

Ik heb geen flauw idee waar ik dit hele sinus en cosinus verhaal moet plaatsen in de statistiek. Berekenen ze een hoek? Of is dit een afgeleide van de Gauss formule?

Alvast bedankt voor de hulp.

Voor degenen die code kunnen lezen, hier de functie:

Code: Selecteer alles

function getPforT($t,$n) {
   $Pi= pi();                                     // Pi value
   $PiD2=$Pi / 2;
   $t=abs($t);
   $w=$t / sqrt($n);
   $th=atan($w);                                 // returns the arctangent
   
   if($n==1) {
      return 1- $th / $PiD2;
   }
   
   $sth=sin($th);                                        // returns the sinus
   $cth=cos($th);                                 // returns the cosinus
      
   if (fmod($n, 2)==1) {                           // modulo ($n%2)
      $z = getZ($cth * $cth,2,$n-3,-1);
      return 1- ($th + $sth * $cth * $z) / $PiD2;
   }
   else {
      $z = getZ($cth * $cth,1,$n-3,-1);
      return 1 - $sth * $z;
   }
}

function getZ($q,$i,$j,$b) {
   $zz=1;
   $z=$zz;
   $k=$i;
   
   while($k <= $j) {
      $zz=$zz * $q * $k / ($k-$b);
      $z=$z+$zz;
      $k=$k+2;
   }
   return $z;
}

Mary
Nieuw lid
Nieuw lid
Berichten: 2
Lid geworden op: 05 jun 2006, 11:48

Bericht door Mary » 05 jun 2006, 12:00

ohja en n is hier de degrees of freedom

Plaats reactie