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)

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

Σχετικοί εξωτερικοί σύνδεσμοι

Αναζήτηση στο google.com για παρόμοια θέματα

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