Mercurial > hg > octave-nkf

comparison scripts/statistics/tests/manova.m @ 20441:83792dd9bcc1

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Use in-place operators in m-files where possible. * scripts/audio/@audioplayer/set.m, scripts/audio/@audiorecorder/set.m, scripts/audio/mu2lin.m, scripts/elfun/cosd.m, scripts/general/del2.m, scripts/general/profexplore.m, scripts/general/quadl.m, scripts/general/rat.m, scripts/general/rotdim.m, scripts/help/get_first_help_sentence.m, scripts/help/private/__strip_html_tags__.m, scripts/image/cubehelix.m, scripts/io/textread.m, scripts/linear-algebra/duplication_matrix.m, scripts/linear-algebra/housh.m, scripts/linear-algebra/krylov.m, scripts/linear-algebra/logm.m, scripts/linear-algebra/normest.m, scripts/linear-algebra/onenormest.m, scripts/optimization/fminsearch.m, scripts/optimization/lsqnonneg.m, scripts/optimization/qp.m, scripts/plot/appearance/annotation.m, scripts/plot/appearance/axis.m, scripts/plot/appearance/legend.m, scripts/plot/appearance/specular.m, scripts/plot/draw/colorbar.m, scripts/plot/draw/hist.m, scripts/plot/draw/plotmatrix.m, scripts/plot/draw/private/__stem__.m, scripts/plot/util/__actual_axis_position__.m, scripts/plot/util/__gnuplot_drawnow__.m, scripts/plot/util/findobj.m, scripts/plot/util/print.m, scripts/plot/util/private/__go_draw_axes__.m, scripts/plot/util/private/__print_parse_opts__.m, scripts/plot/util/rotate.m, scripts/polynomial/pchip.m, scripts/polynomial/polyaffine.m, scripts/polynomial/polyder.m, scripts/polynomial/private/__splinefit__.m, scripts/polynomial/residue.m, scripts/signal/arch_fit.m, scripts/signal/arch_rnd.m, scripts/signal/bartlett.m, scripts/signal/blackman.m, scripts/signal/freqz.m, scripts/signal/hamming.m, scripts/signal/hanning.m, scripts/signal/spectral_adf.m, scripts/signal/spectral_xdf.m, scripts/signal/stft.m, scripts/sparse/bicgstab.m, scripts/sparse/cgs.m, scripts/sparse/private/__sprand_impl__.m, scripts/sparse/qmr.m, scripts/sparse/sprandsym.m, scripts/sparse/svds.m, scripts/specfun/legendre.m, scripts/special-matrix/gallery.m, scripts/statistics/base/gls.m, scripts/statistics/models/logistic_regression.m, scripts/statistics/tests/kruskal_wallis_test.m, scripts/statistics/tests/manova.m, scripts/statistics/tests/wilcoxon_test.m, scripts/time/datevec.m: Use in-place operators in m-files where possible.
author Rik <rik@octave.org>
date Tue, 26 May 2015 21:07:42 -0700
parents d9341b422488
children
comparison
equal deleted inserted replaced
20440:e914b5399c67 20441:83792dd9bcc1
67 error ("manova: there should be at least 2 groups"); 67 error ("manova: there should be at least 2 groups");
68 else 68 else
69 group_label = s ([1, (reshape (i, 1, k - 1) + 1)]); 69 group_label = s ([1, (reshape (i, 1, k - 1) + 1)]);
70 endif 70 endif
71 71
72 x = x - ones (n, 1) * mean (x); 72 x -= ones (n, 1) * mean (x);
73 SST = x' * x; 73 SST = x' * x;
74 74
75 s = zeros (1, p); 75 s = zeros (1, p);
76 SSB = zeros (p, p); 76 SSB = zeros (p, p);
77 for i = 1 : k; 77 for i = 1 : k;
78 v = x (find (g == group_label (i)), :); 78 v = x (find (g == group_label (i)), :);
79 s = sum (v); 79 s = sum (v);
80 SSB = SSB + s' * s / rows (v); 80 SSB += s' * s / rows (v);
81 endfor 81 endfor
82 n_b = k - 1; 82 n_b = k - 1;
83 83
84 SSW = SST - SSB; 84 SSW = SST - SSB;
85 n_w = n - k; 85 n_w = n - k;