File: spDists.Rout.save

package info (click to toggle)
r-cran-sp 1%3A0.9-66-1
  • links: PTS, VCS
  • area: main
  • in suites: squeeze
  • size: 1,696 kB
  • ctags: 123
  • sloc: ansic: 1,475; sh: 6; makefile: 5
file content (41 lines) | stat: -rw-r--r-- 1,275 bytes parent folder | download | duplicates (6)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41

R version 2.10.0 (2009-10-26)
Copyright (C) 2009 The R Foundation for Statistical Computing
ISBN 3-900051-07-0

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(sp)
> data(meuse)
> data(meuse.grid)
> coordinates(meuse) = ~x+y
> coordinates(meuse.grid) = ~x+y
> x = coordinates(meuse)
> y = coordinates(meuse.grid)
> out = spDists(meuse,meuse.grid)
> out2 = as.matrix(dist(rbind(coordinates(meuse),coordinates(meuse.grid))))
> out2 = out2[1:155,155+1:nrow(y)]
> # should be equal:
> sum(out2 - out)
[1] 0
> summary(as.vector(out2 - out))
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
      0       0       0       0       0       0 
> 
> out = spDists(meuse.grid,meuse)
> sum(out2 - t(out))
[1] 0
> summary(as.vector(out2 - t(out)))
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
      0       0       0       0       0       0 
>