Springen naar inhoud

Gauss formule om P-waarde te berekenen: help wanted!


  • Log in om te kunnen reageren

#1

StarMadeMary

    StarMadeMary


  • 0 - 25 berichten
  • 2 berichten
  • Gebruiker

Geplaatst op 05 juni 2006 - 12:03

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:

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;

}

Dit forum kan gratis blijven vanwege banners als deze. Door te registeren zal de onderstaande banner overigens verdwijnen.

#2

StarMadeMary

    StarMadeMary


  • 0 - 25 berichten
  • 2 berichten
  • Gebruiker

Geplaatst op 05 juni 2006 - 12:04

ohja en n is hier de degrees of freedom

#3

Albert E

    Albert E


  • 0 - 25 berichten
  • 2 berichten
  • Gebruiker

Geplaatst op 05 juni 2006 - 13:20

hmmm..... k ga even voor je kijken





0 gebruiker(s) lezen dit onderwerp

0 leden, 0 bezoekers, 0 anonieme gebruikers

Ook adverteren op onze website? Lees hier meer!

Gesponsorde vacatures

Vacatures