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
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
|
# $Id: space.R 625 2005-06-09 14:20:30Z nj7w $
# When there are two or more points with the same (x,y) value (or
# within x+-s[1] and x+-s[2]), space these out in the y direction so
# that the points are separated by at least distance s.
space <- function(x,y,s=1/50, na.rm=TRUE, direction="x")
{
if(direction!='x')
{
tmp <- y
y <- x
x <- tmp
}
if(na.rm)
{
ind <- is.na(x) | is.na(y)
x <- x[!ind]
y <- y[!ind]
}
if (length(s)==1) s <- c(s,s)
spacing.x <- (max(x) - min(x))*s[1]
spacing.y <- (max(y) - min(y))*s[2]
within <- function(x,y,delta) { abs(x-y) < delta }
# sort x,y so we can do the work
ord <- order(x,y)
undo <- order(ord)
x <- x[ord]
y <- y[ord]
startsame <- 1
same.x <- x[1]
same.y <- y[1]
for( i in 1:length(x) )
{
if(i>1 &&
within(x[i],same.x,spacing.x) &&
within(y[i],same.y,spacing.y) )
{
if(x[startsame] == same.x )
x[startsame] <- x[startsame]
cumrun <- i - startsame
x[i] <- x[i] + (-1)^(cumrun+1) * floor((cumrun+1) /2) * spacing.x
} else {
startsame <- i
same.x <- x[i]
same.y <- y[i]
}
}
if(direction!='x')
{
tmp <- y
y <- x
x <- tmp
}
return( list(x=x[undo], y=y[undo]) )
}
|
