[R] manipulate a matrix2
Jon Hak
Jon_Hak at natureserve.org
Mon Jul 21 16:22:24 CEST 2008
Thanks Jim, that was exactly what I was after.
On a second note, do you have any insight into pulling out the
duplicates in this type of matrix?
I thought that was what the upper=FALSE is in:
csv.dis <- vegdist(csv.m, method='jaccard', binary=FALSE, diag=FALSE,
upper=FALSE). I just need either the lower or upper portion, with the
zeros (,3 & ,3) being the dividing line.
[,3] [,5] [,6] [,9] [,11]
[3,] 0 2 3 4 5
[5,] 2 0 8 9 10
[6,] 3 8 0 14 15
[9,] 4 9 14 0 20
[11,] 5 10 15 20 0
Thanks again,
Jon
-----Original Message-----
From: jim holtman [mailto:jholtman at gmail.com]
Sent: Friday, July 18, 2008 9:56 AM
To: Jon Hak
Cc: r-help at r-project.org
Subject: Re: [R] manipulate a matrix2
Is this what you want:
> x
[,3] [,5] [,6] [,9] [,11]
[,3] 1 6 11 16 21
[,5] 2 7 12 17 22
[,6] 3 8 13 18 23
[,9] 4 9 14 19 24
[,11] 5 10 15 20 25
> library(reshape)
> melt(x)
X1 X2 value
1 [,3] [,3] 1
2 [,5] [,3] 2
3 [,6] [,3] 3
4 [,9] [,3] 4
5 [,11] [,3] 5
6 [,3] [,5] 6
7 [,5] [,5] 7
8 [,6] [,5] 8
9 [,9] [,5] 9
10 [,11] [,5] 10
11 [,3] [,6] 11
12 [,5] [,6] 12
13 [,6] [,6] 13
14 [,9] [,6] 14
15 [,11] [,6] 15
16 [,3] [,9] 16
17 [,5] [,9] 17
18 [,6] [,9] 18
19 [,9] [,9] 19
20 [,11] [,9] 20
21 [,3] [,11] 21
22 [,5] [,11] 22
23 [,6] [,11] 23
24 [,9] [,11] 24
25 [,11] [,11] 25
>
On Fri, Jul 18, 2008 at 11:10 AM, Jon Hak <Jon_Hak at natureserve.org>
wrote:
> Building upon Jim's answer below (Thanks Jim, that helped a lot), I
need
> to pickup where this thread left off. I'm using Vegan to calculate
the
> Jaccard's Index and the Row.Names and column names are represented in
my
> matrix as seen here.
>
> [,3] [,5] [,6] [,9] [,11]
>
> [3,] 0 6 11 16 21
>
> [5,] 2 0 12 17 22
>
> [6,] 3 8 0 18 23
>
> [9,] 4 9 14 0 24
>
> [11,] 5 10 15 20 0
>
>
>
> When I use the command;
> xy <- cbind(row=as.vector(row.names(x)), col=as.vector(colnames(x)),
> value=as.vector(x))
>
> I get the list (the column value is the issue);
>
> row col value
> [1,] 3 1 0
> [2,] 5 1 2
> [3,] 6 1 3
> [4,] 9 1 4
> [5,] 11 1 5
> [6,] 3 2 6
> [7,] 5 2 0
> [8,] 6 2 8
> [9,] 9 2 9
> [10,] 11 2 10
> [11,] 3 3 11
> [12,] 5 3 0
>
> I would really like the col value to equal the actual name, not the
> column number. What am I missing? The analysis is very large, 6k x6k
> matrix so automating the process is a high priority.
>
> Thanks,
> Jon
>
>
>
> From: jim holtman <jholtman_at_gmail.com
>
<mailto:jholtman_at_gmail.com?Subject=Re:%20%5BR%5D%20manipulate%20a%20m
> atrix> >
> Date: Mon, 25 Jun 2007 12:39:46 -0400
>
> Is this what you want?
>
>> x
>
> [,1] [,2] [,3] [,4] [,5]
>
> [1,] 1 6 11 16 21
> [2,] 2 7 12 17 22
> [3,] 3 8 13 18 23
> [4,] 4 9 14 19 24
> [5,] 5 10 15 20 25
>
>
>
> cbind(row=as.vector(row(x)), col=as.vector(col(x)),
value=as.vector(x))
>
>
> row col value
> [1,] 1 1 1
> [2,] 2 1 2
> [3,] 3 1 3
> [4,] 4 1 4
> [5,] 5 1 5
> [6,] 1 2 6
> [7,] 2 2 7
> [8,] 3 2 8
> [9,] 4 2 9
> [10,] 5 2 10
> [11,] 1 3 11
> [12,] 2 3 12
> [13,] 3 3 13
> [14,] 4 3 14
> [15,] 5 3 15
> [16,] 1 4 16
> [17,] 2 4 17
> [18,] 3 4 18
> [19,] 4 4 19
> [20,] 5 4 20
> [21,] 1 5 21
> [22,] 2 5 22
> [23,] 3 5 23
>
>
> [24,] 4 5 24
> [25,] 5 5 25
>>
>
> On 6/25/07, Jon Hak <Jon_Hak_at_natureserve.org> wrote:
>>
>> I have read everything I can find on how to manipulate a results
> matrix in
>
<http://tolstoy.newcastle.edu.au/R/e2/help/07/06/19875.html#19887qlink1>
>
>> R and I have to admit I'm stumped. I have set up a process to extract
> a
>> dataset from ArcGIS to compute a similarity index (Jaccards) in
Vegan.
> The
>> dataset is fairly simple, but large, and consists of rows = sample
> area, and
>> columns = elements. I've been able to view the results in R, but I
> want to
>> get the results out to a database and a matrix that is 6000-rows x
>> 6000-columns can be very difficult to manipulate in Windows XP. I
> would to
>> rotate the matrix so that the output would look like the old
condensed
>
>> format in programs like Conoco. Ideally, I would like format to look
>> something like this;
>>
>>
>> Site-row Site-col Jaccard
>> 1 1 1
>> 1 2 .9
>> 1 3 .6
>> 2 1 .9
>> 2 2 1
>> 2 3 .75
>>
>> Thanks for any help,
>>
>>
>>
>>
>> ***********************************************************
>> John Hak
>> Senior GIS Analyst/Sr. Ecologist
>> NatureServe
>> 4001 Discovery Drive
>> Boulder, CO 80303
>> (703) 797-4809
>>
>> There is perhaps no better demonstration of the folly of human
> conceits
>> than this distant image of our tiny world. To me, it underscores our
>> responsibility to deal more kindly with one another, and to preserve
> and
>> cherish the pale blue dot, the only home we've ever known. --Carl
> Sagan
>>
>>
>> ______________________________________________
>> R-help_at_stat.math.ethz.ch mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide
>> http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.
>>
>
> --
> Jim Holtman
> Cincinnati, OH
> +1 513 646 9390
>
> What is the problem you are trying to solve?
>
> [[alternative HTML version deleted]]
>
>
>
>
>
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>
--
Jim Holtman
Cincinnati, OH
+1 513 646 9390
What is the problem you are trying to solve?
More information about the R-help
mailing list