[R] as.dist with diagonal unequal zero
Birgit Lemcke
birgit.lemcke at systbot.uzh.ch
Wed Oct 10 17:06:54 CEST 2007
This are distances but between male and female plants of the same
species.
species A B C D E --> male
A
B
C
D
|
female
I would like to do a Manteltest.
Birgit
Am 10.10.2007 um 16:59 schrieb Gustaf Rydevik:
> On 10/10/07, Birgit Lemcke <birgit.lemcke at systbot.uzh.ch> wrote:
>> Hello and sorry that I still haven´t found a solution for my problem.
>>
>> I need to extract the lower and upper triangle from a square matrix
>> including the diagonal. This diagonal is not zero in that special
>> case.
>> I tried with as.dist
>>
>> w<-as.dist(w, diag = TRUE)
>>> w
>> 1 2 3 4 5
>> 1 0
>> 2 2 0
>> 3 3 8 0
>> 4 4 9 14 0
>> 5 5 10 15 20 0
>>
>> but found no way to keep the diagonal that is in the matrix.
>>
>> I also tried this but as I expected it makes no difference:
>>
>> w<-matrix(1:25, ncol=5, nrow=5)
>>> upperTriangle(w)<-NA
>>> w
>> [,1] [,2] [,3] [,4] [,5]
>> [1,] 1 NA NA NA NA
>> [2,] 2 7 NA NA NA
>> [3,] 3 8 13 NA NA
>> [4,] 4 9 14 19 NA
>> [5,] 5 10 15 20 25
>>> w<-as.dist(w, diag = TRUE)
>>> w
>> 1 2 3 4 5
>> 1 0
>> 2 2 0
>> 3 3 8 0
>> 4 4 9 14 0
>> 5 5 10 15 20 0
>>
>> Is there perhaps a possibility to give the diagonal as a vector
>> within as.dist?
>>
>> I would be glad if somebody could help me.
>>
>> Greetings
>>
>> Birgit
>>
>> Birgit Lemcke
>> Institut für Systematische Botanik
>> Zollikerstrasse 107
>> CH-8008 Zürich
>> Switzerland
>> Ph: +41 (0)44 634 8351
>> birgit.lemcke at systbot.uzh.ch
>>
>
> What do you want to do with the lower triangle matrix?
> A distance matrix must be zero on the diagonal, since dist(a,a)==0 for
> all a, so it cannot be distances that you measure.
>
> best,
> Gustaf
>
>
> --
> Gustaf Rydevik, M.Sci.
> tel: +46(0)703 051 451
> address:Essingetorget 40,112 66 Stockholm, SE
> skype:gustaf_rydevik
Birgit Lemcke
Institut für Systematische Botanik
Zollikerstrasse 107
CH-8008 Zürich
Switzerland
Ph: +41 (0)44 634 8351
birgit.lemcke at systbot.uzh.ch
More information about the R-help
mailing list