Running 'tseqimp.R' [169s]
Running the tests in 'tests/tseqimp.R' failed.
Complete output:
  > ## VT::15.09.2013 - this will render the output independent
  > ##  from the version of the package
  > suppressPackageStartupMessages(library(rrcovNA))
  > 
  > alpha <- 0.55
  > 
  > data(phosphor); x <- y <- phosphor[,1:2]; x[10,2] <- NA; x[15,1] <- NA
  > mcdc <- CovMcd(y)                                           # mcdc - complete
  > ximp <- impSeq(x); mcds <- CovMcd(ximp)                     # mcd - sequential imputation + MCD
  > ximp <- impSeqRob(x, alpha=alpha); mcd <- CovMcd(ximp$x)    # mcd - robust sequentioal imputation + MCD
  > mcdna <- CovNAMcd(x)                                        # mcdna - norm + MCD
  > as.vector(which(mcdc@wt==0))
  [1]  1  6 10
  > as.vector(which(mcds@wt==0))
  [1] 1 6
  > as.vector(which(mcd@wt==0))
  [1]  1  4  6  7 16
  > as.vector(which(mcdna@wt==0))
  integer(0)
  > ##cbind(mcdc@wt, mcdna@wt, mcd@wt)
  > 
  > data(heart); x <- y <- heart; x[10,2] <- NA; x[2,1] <- NA
  > mcdc <- CovMcd(y)
  > ximp <- impSeq(x); mcds <- CovMcd(ximp)                     # mcd - sequential imputation + MCD
  > ximp <- impSeqRob(x, alpha=alpha); mcd <- CovMcd(ximp$x)
  > mcdna <- CovNAMcd(x)
  > as.vector(which(mcdc@wt==0))
  [1]  3  8  9 10
  > as.vector(which(mcds@wt==0))
  [1]  3  8  9 10
  > as.vector(which(mcd@wt==0))
  [1]  3  8  9 10
  > as.vector(which(mcdna@wt==0))
  integer(0)
  > ##cbind(mcdc@wt, mcdna@wt, mcd@wt)
  > 
  > data(starsCYG); x <- y <- starsCYG; x[10,2] <- NA; x[2,1] <- NA; x[33,1] <- NA; x[41,1] <- NA
  > mcdc <- CovMcd(y)
  > ximp <- impSeq(x); mcds <- CovMcd(ximp)                     # mcd - sequential imputation + MCD
  > ximp <- impSeqRob(x, alpha=alpha); mcd <- CovMcd(ximp$x)
  > mcdna <- CovNAMcd(x)
  > as.vector(which(mcdc@wt==0))
  [1]  7  9 11 14 20 30 34
  > as.vector(which(mcds@wt==0))
  [1]  2  7  9 11 14 20 30 34
  > as.vector(which(mcd@wt==0))
  [1]  7  9 11 14 20 30 34
  > as.vector(which(mcdna@wt==0))
  integer(0)
  > ##cbind(mcdc@wt, mcdna@wt, mcd@wt)
  > 
  > data(stackloss); x <- y <- stack.x; x[10,2] <- NA; x[6,1] <- NA; x[13,3] <- NA
  > mcdc <- CovMcd(y)
  > ximp <- impSeq(x); mcds <- CovMcd(ximp)                     # mcd - sequential imputation + MCD
  > ximp <- impSeqRob(x, alpha=alpha); mcd <- CovMcd(ximp$x)
  > mcdna <- CovNAMcd(x)
  > as.vector(which(mcdc@wt==0))
  [1]  1  2  3 15 16 17 18 19 21
  > as.vector(which(mcds@wt==0))
  [1]  7  8 10 14 17 18 19 20
  > as.vector(which(mcd@wt==0))
  [1]  1  2  3 15 16 17 18 19 21
  > as.vector(which(mcdna@wt==0))
  integer(0)
  > ##cbind(mcdc@wt, mcdna@wt, mcd@wt)
  > 
  > data(coleman); x <- y <- data.matrix(subset(coleman, select = -Y)); x[5,2] <- NA; x[8,4] <- NA; x[13,3] <- NA
  > mcdc <- CovMcd(y)
  > ximp <- impSeq(x); mcds <- CovMcd(ximp)                     # mcd - sequential imputation + MCD
  > ximp <- impSeqRob(x, alpha=alpha); mcd <- CovMcd(ximp$x)
  > mcdna <- CovNAMcd(x)
  > ##as.vector(which(mcdc@wt==0))
  > ##as.vector(which(mcds@wt==0))
  > ##as.vector(which(mcd@wt==0))
  > ##as.vector(which(mcdna@wt==0))
  > ##cbind(mcdc@wt, mcdna@wt, mcd@wt)
  > 
  > data(salinity); x <- y <- data.matrix(subset(salinity, select = -Y)); x[1,2] <- NA; x[8,3] <- NA; x[13,3] <- NA
  > mcdc <- CovMcd(y)
  > ximp <- impSeq(x); mcds <- CovMcd(ximp)                     # mcd - sequential imputation + MCD
  > ximp <- impSeqRob(x, alpha=alpha); mcd <- CovMcd(ximp$x)
  > mcdna <- CovNAMcd(x)
  > as.vector(which(mcdc@wt==0))
  [1]  5 16 23 24
  > as.vector(which(mcds@wt==0))
  [1]  3  5 11 16 23 24
  > as.vector(which(mcd@wt==0))
  [1]  3  5 11 16 23 24
  > as.vector(which(mcdna@wt==0))
  integer(0)
  > ##cbind(mcdc@wt, mcdna@wt, mcd@wt)
  > 
  > data(wood); x <- y <- data.matrix(subset(wood, select = -y)); x[1,2] <- NA; x[10,3] <- NA; x[13,4] <- NA
  > mcdc <- CovMcd(y)
  > ximp <- impSeq(x); mcds <- CovMcd(ximp)                     # mcd - sequential imputation + MCD
  > ximp <- impSeqRob(x, alpha=alpha); mcd <- CovMcd(ximp$x)
  > mcdna <- CovNAMcd(x)
  > as.vector(which(mcdc@wt==0))
  [1]  4  6  7  8 11 16 19
  > as.vector(which(mcds@wt==0))
  [1]  1  7  9 11 12 16 17
  > as.vector(which(mcd@wt==0))
  [1]  4  6  8 19
  > as.vector(which(mcdna@wt==0))
  integer(0)
  > ##cbind(mcdc@wt, mcdna@wt, mcd@wt)
  > 
  > data(hbk); x <- y <- data.matrix(subset(hbk, select = -Y)); x[30,2] <- NA; x[40,3] <- NA; x[17,3] <- NA
  > mcdc <- CovMcd(y)
  > ximp <- impSeq(x); mcds <- CovMcd(ximp)                     # mcd - sequential imputation + MCD
  > ximp <- impSeqRob(x, alpha=alpha); mcd <- CovMcd(ximp$x)
  > mcdna <- CovNAMcd(x)
  > as.vector(which(mcdc@wt==0))
   [1]  1  2  3  4  5  6  7  8  9 10 11 12 13 14
  > as.vector(which(mcds@wt==0))
   [1]  1  2  3  4  5  6  7  8  9 10 11 12 13 14
  > as.vector(which(mcd@wt==0))
   [1]  1  2  3  4  5  6  7  8  9 10 11 12 13 14
  > as.vector(which(mcdna@wt==0))
  integer(0)
  > ##cbind(mcdc@wt, mcdna@wt, mcd@wt)
  > 
  > proc.time()
     user  system elapsed 
     1.18    0.17    1.36