correctGeno {MPR} | R Documentation |
~~ A concise (1-5 lines) description of what the function does. ~~
correctGeno(geno.data, base.position = as.numeric(sub("[^0-9]*([0-9]*)[^0-9]*", "\1", rownames(geno.data))), correct.FUN = correctFUNHMM, minInterval = 1, verbose = TRUE, ...)
geno.data |
~~Describe geno.data here~~ |
base.position |
~~Describe base.position here~~ |
correct.FUN |
~~Describe correct.FUN here~~ |
minInterval |
~~Describe minInterval here~~ |
verbose |
~~Describe verbose here~~ |
... |
~~Describe ... here~~ |
~~ If necessary, more details than the description above ~~
~Describe the value returned If it is a LIST, use
comp1 |
Description of 'comp1' |
comp2 |
Description of 'comp2' |
...
....
~~further notes~~
~Make other sections like Warning with section{Warning }{....} ~
~~who you are~~
~put references to the literature/web site here ~
~~objects to See Also as help
, ~~~
##---- Should be DIRECTLY executable !! ---- ##-- ==> Define data, use random, ##-- or do help(data=index) for the standard data sets. ## The function is currently defined as function (geno.data, base.position = as.numeric(sub("[^0-9]*([0-9]*)[^0-9]*", "\1", rownames(geno.data))), correct.FUN = correctFUNHMM, minInterval = 1, verbose = TRUE, ...) { i <- 0 geno.data.cr <- apply(geno.data, 2, function(x, ...) { if (verbose) cat("\r", i <<- i + 1) x.nna <- which(!is.na(x)) ids <- sort(unique(x.nna[c(1, which(diff(x[x.nna]) != 0 | diff(base.position[x.nna]) >= minInterval) + 1)])) x.cr <- correct.FUN(x[ids], base.position = base.position[ids], ...) x[x.nna] <- NA x[ids] <- x.cr x }, ...) if (verbose) cat("\tDone.\n") geno.data.cr }