Re: PCA with Octave

[Mahmood Naderan]

Finally I successfully installed the statistics package and now I can see princomp(). Here is what I did

octave:1> X=[1 0.7; 2 0.4; 3 0.55; 4 0.2; 5 0.63; 6 0.83; 7 0.48; 8 0.91; 9 0.73]; 
octave:5> [pc,score,l,t] = princomp(X)
 warning: XXX FIXME XXX Tsq return from princomp fails some tests
 pc =

    -0.999358  -0.035833
    -0.035833   0.999358

 score =

    3.9940e+00   2.3994e-01
    3.0054e+00  -9.5704e-02
    2.0006e+00   1.8367e-02
    1.0138e+00  -3.6724e-01
   -9.5555e-04   2.6650e-02
   -1.0075e+00   1.9069e-01
   -1.9943e+00  -1.9492e-01
   -3.0091e+00   1.9897e-01
   -4.0020e+00  -1.6747e-02

 l =

    7.509591
    0.039609

 t =

    3.577651
    1.433997
    0.541503     
    3.541846
    0.017931
    1.053196
    1.488853
    2.205235
    2.139788

As I look into the figures (from websites) I see two axis which show the direction of the data points.
How can I plot the to axis?

Regards,
Mahmood

---

From: Juan Pablo Carbajal <address@hidden>
To: Jordi Gutiérrez Hermoso <address@hidden>
Cc: "address@hidden" <address@hidden>
Sent: Monday, July 8, 2013 6:31 PM
Subject: Re: PCA with Octave

On Mon, Jul 8, 2013 at 3:17 PM, Jordi Gutiérrez Hermoso
<address@hidden> wrote:
> Octave has a princomps function that nobody ever seems to find because
> it's not named "pca".
>
> You don't need to reimplement PCA yourself.
>
> - Jordi G. H.

It is in the statistics package...time to rename it...
http://octave.sourceforge.net/statistics/function/princomp.html
_______________________________________________
Help-octave mailing list
address@hidden
https://mailman.cae.wisc.edu/listinfo/help-octave