The R Project SVN R

Rev

Rev 88598 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 88598 Rev 88896
Line 62... Line 62...
62
rN <- sapply(eRules, function(R) runmed(yN, 7, endrule=R))
62
rN <- sapply(eRules, function(R) runmed(yN, 7, endrule=R))
63
matlines(rN, type = "b", pch = 4, lwd = 3, lty=2,
63
matlines(rN, type = "b", pch = 4, lwd = 3, lty=2,
64
         col = adjustcolor(c("red", "orange4", "orange1"), 0.5))
64
         col = adjustcolor(c("red", "orange4", "orange1"), 0.5))
65
yN[c(1, 20:21)] <- NA # additionally
65
yN[c(1, 20:21)] <- NA # additionally
66
rN. <- sapply(eRules, function(R) runmed(yN, 7, endrule=R))
66
rN. <- sapply(eRules, function(R) runmed(yN, 7, endrule=R))
67
head(rN., 4); tail(rN.) # more NA's too, still not *so* many:
67
head(rN., 4); tail(rN.) # more NAs too, still not *so* many:
68
stopifnot(exprs = {
68
stopifnot(exprs = {
69
   !anyNA(rN[,2:3])
69
   !anyNA(rN[,2:3])
70
   identical(which(is.na(rN[,"keep"])), c(2L, 40L))
70
   identical(which(is.na(rN[,"keep"])), c(2L, 40L))
71
   identical(which(is.na(rN.), arr.ind=TRUE, useNames=FALSE),
71
   identical(which(is.na(rN.), arr.ind=TRUE, useNames=FALSE),
72
             cbind(c(1:2,40L), 1L))
72
             cbind(c(1:2,40L), 1L))