Wald-Wolfowitz τεστ (κριτήριο ροών) σε Octave/Matlab

Θα προσθέσω αναλυτικότερη περιγραφή αργότερα. Ένα παράδειγμα για τον υπολογισμό p-value του κριτηρίου Wald-Wolfowitz σε Octave/Matlab:

clear;
T = 10;

x  = unifrnd(0, 1, 10, 1);
w  = x>median(x);
n1 = sum(w)
n2 = T-w;

Runs  = 1;
for i = 2:T
  if ( w(i) != w(i-1) )
    ++Runs;
  end
end
Runs

muR      = 2*n1*n2/(n1+n2) + 1
sigma2R  = (muR-1) * (muR-2) / (T-1);
sigmaR   = sqrt(sigma2R);

if (Runs > muR)
  Z = (Runs - muR - 0.5) / sigmaR;
else
  Z = (Runs - muR + 0.5) / sigmaR;
end

Z
pval = normcdf(abs(Z), 0, 1)

Συνδεθείτε για περισσότερες δυνατότητες αλληλεπίδρασης,
σχολιασμοί, εξωτερικοί σύνδεσμοι, βοήθεια, ψηφοφορίες, αρχεία, κτλ.

Creative Commons License
Εκπαιδευτικό υλικό από τον Αθανάσιο Σταυρακούδη σας παρέχετε κάτω από την άδεια Creative Commons Attribution-NonCommercial-ShareAlike 4.0 License.
Σας παρακαλώ να ενημερωθείτε για κάποιους επιπλέον περιορισμούς
http://stavrakoudis.econ.uoi.gr/stavrakoudis/?iid=401.