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 Attribution-NonCommercial-ShareAlike 4.0 License.
Σας παρακαλώ να ενημερωθείτε για κάποιους επιπλέον περιορισμούς
http://stavrakoudis.econ.uoi.gr/stavrakoudis/?iid=401.