Thanks for the tip, Nikhil. However, i need only one matrix as input for
another to compute for non-bipartite matching which minimizes pairwise
distances between observations. As such, I need the georeference (id) of the
observations for subsequent processing. Below is an illustration.
> #generate data
> x <- as.matrix(runif(5))
> Sx <- var(x)
>
> #generate id
> set.seed(1)
> id1 <- sample(1:2,5, replace=T)
> id2 <- c(1:5)
> rownames(x) <- paste(id1, id2)
>
> #generate distance
> dist <- as.matrix(
+ apply(x,1,function(i){
+ mahalanobis(x,i,Sx)
+ }
+ )
+ )
>
> #print matrices
> x
[,1]
1 1 0.2059746
1 2 0.1765568
2 3 0.6870228
2 4 0.3841037
1 5 0.7698414
> dist
1 1 1 2 2 3 2 4 1 5
1 1 0.00000000 0.01165534 3.11660015 0.4273402 4.28210082
1 2 0.01165534 0.00000000 3.50943798 0.5801450 4.74056406
2 3 3.11660015 3.50943798 0.00000000 1.2358255 0.09237602
2 4 0.42734018 0.58014499 1.23582554 0.0000000 2.00395492
1 5 4.28210082 4.74056406 0.09237602 2.0039549 0.00000000
The geo-id is composed of two references, the first digit for the region and
the next for the observation itself. What I'm thinking of is for pairwise
distance between observations of different regions, say site-11 and site-23
or site-24 to be replaced by a large number, say 999999. I need the id for
future processing, though.
Maybe I can stack the matrices generated using your tip to form a block
diagonal matrix, but then I do not have my ids? Im really sorry. Im a bit
lost.
Cheers,
Michael
On Mon, Jul 19, 2010 at 10:10 PM, Nikhil Kaza <nikhil.l...@gmail.com> wrote:
>
> replace dist with mahalanobis distance in the following example.
>
> a <- cbind(runif(10), sample(1:3, 10, replace=T))
> a.L <- split(a,a[,2])
> dist.L <- lapply(a.L, dist)
>
>
>
> Nikhil Kaza
> Asst. Professor,
> City and Regional Planning
> University of North Carolina
>
> nikhil.l...@gmail.com
>
>
> On Jul 19, 2010, at 9:24 AM, Michael Ralph M. Abrigo wrote:
>
> Hi! I am trying to implement non-bipartite matching. I have around 500
>> sites
>> which can be clustered by 10 regions. I am able to calculate pairwise
>> Mahalanobis distances between sites (thanks to another post in the forum).
>> However, I want to constrain my match to sites within the same region.
>> Thus
>> I want to replace elements of the distance matrix with a high value, say
>> 999999, for sites not of the same region so that the pair will not be
>> matched.
>> In the original data file I have information on which sites belong to what
>> region. However, when I compute for pairwise Mahalanobis distances, I only
>> use a subset of the file, which, naturally, does not include the
>> georeference of the sites. How should I do this? Any hint will be most
>> appreciated.
>> Btw, I am relatively new in using R. I may export the matrix to another
>> program and replace the elements there, but that is a very very dirty and
>> rough trick that I would rather not do given better options.
>> Many thanks in advance.
>>
>> Cheers,
>> Michael
>>
>> --
>> "I am most anxious for liberties for our country... but I place as a prior
>> condition the education of the people so that our country may have an
>> individuality of its own and make itself worthy of liberties... " Jose
>> Rizal,1896
>>
>> [[alternative HTML version deleted]]
>>
>> ______________________________________________
>> R-help@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.
>>
>
>
--
"I am most anxious for liberties for our country... but I place as a prior
condition the education of the people so that our country may have an
individuality of its own and make itself worthy of liberties... " Jose
Rizal,1896
[[alternative HTML version deleted]]
______________________________________________
R-help@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.
