The R Project SVN R

Rev

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

Rev 88628 Rev 88896
Line 158... Line 158...
158 
%% gives 6-line warning but does NOT break down
 158 
%% gives 6-line warning but does NOT break down
159 
lines(runmed(y, 7),      lwd = 2, col = col[3])
 159 
lines(runmed(y, 7),      lwd = 2, col = col[3])
160 
legend(length(y),max(y), c("data", "lowess(y, f = 0.3)", "runmed(y, 7)"),
 160 
legend(length(y),max(y), c("data", "lowess(y, f = 0.3)", "runmed(y, 7)"),
161 
       xjust = 1, col = col, lty = c(0, 1, 1), pch = c(1,NA,NA))
 161 
       xjust = 1, col = col, lty = c(0, 1, 1), pch = c(1,NA,NA))
162 
 
 162 
 
163 
## An example with initial NA's - used to fail badly (notably for "Turlach"):
 163 
## An example with initial NAs - used to fail badly (notably for "Turlach"):
164 
x15 <- c(rep(NA, 4), c(9, 9, 4, 22, 6, 1, 7, 5, 2, 8, 3))
 164 
x15 <- c(rep(NA, 4), c(9, 9, 4, 22, 6, 1, 7, 5, 2, 8, 3))
165 
rS15 <- cbind(Sk.3 = runmed(x15, k = 3, algorithm="S"),
 165 
rS15 <- cbind(Sk.3 = runmed(x15, k = 3, algorithm="S"),
166 
              Sk.7 = runmed(x15, k = 7, algorithm="S"),
 166 
              Sk.7 = runmed(x15, k = 7, algorithm="S"),
167 
              Sk.11= runmed(x15, k =11, algorithm="S"))
 167 
              Sk.11= runmed(x15, k =11, algorithm="S"))
168 
rT15 <- cbind(Tk.3 = runmed(x15, k = 3, algorithm="T", print.level=1),
 168 
rT15 <- cbind(Tk.3 = runmed(x15, k = 3, algorithm="T", print.level=1),