[R] Coinertia randtest

Petar Milin pmilin at ff.uns.ac.rs
Mon Aug 23 21:33:19 CEST 2010


Hello!
I dunno why, but I cannot make randtes.coinertia() from ade4 package
working. I have two nice distance matrices (Euclidean):
>  dist1
            1         2         3         4         5         6         7
2  2.5776799
3  1.7892825 1.0637487
4  1.0557991 2.4270728 2.0626604
5  1.6745483 4.1505805 3.4581614 1.8900295
6  2.5045058 0.8144662 0.7157900 2.6888306 4.1708413
7  4.1367323 1.6058153 2.6451634 3.7795019 5.6220092 2.1504953
8  1.5224318 1.4444836 1.3682047 1.0048085 2.8425813 1.8464349 2.7853605
9  4.5321276 4.5859822 3.8178836 5.4277885 5.7703245 3.7731505 5.6657106
10 4.0389095 1.6702235 2.7327563 3.5495557 5.4279792 2.3611163 0.5540045
            8         9
2
3
4
5
6
7
8
9  5.1082493
10 2.5924107 6.0248981

>  dist2
            1         2         3         4         5         6         7
2  0.9364828
3  1.4333876 1.9160376
4  0.2523886 1.1884864 1.3697080
5  2.3753105 2.5185909 3.7830808 2.4140423
6  1.7928748 1.6959363 1.0841587 1.8848077 4.0827074
7  2.5117723 1.6858232 2.7455418 2.7406021 4.0115957 1.8140011
8  0.6014150 0.3748333 1.8022486 0.8514693 2.2931638 1.8062392 2.0515604
9  2.8802778 3.5524921 1.6707483 2.7134664 4.9800100 2.5480973 4.3439383
10 3.0686316 2.6196374 2.4765500 3.2133783 5.1381903 1.4228493 1.5749603
            8         9
2
3
4
5
6
7
8
9  3.3777063
10 2.8790276 3.6806114

And then, I follow from help:
>  pco1 = dudi.pco(dist1, nf=2, scannf=F)
>  pco2 = dudi.pco(dist2, nf=2, scannf=F)
>  coi = coinertia(pco1, pco2, nf=2, scannf=F)
>  testco1 = randtest.coinertia(coi, nrepet=1000)

The result is:
>  Error in randtest.coinertia(coi, nrepet = 1000) : Not yet available

Could anyone help with this?

Best,
PM



More information about the R-help mailing list