R/linlogtics.R

In TIMP: Fitting Separable Nonlinear Models in Spectroscopy and Microscopy

"linlogtics" <- function (x, mu, alpha) 
{
    maxorigx <- max(x)
    minorigx <- min(x)
    ticsl <- c(-alpha)
    tics <- c(-alpha)
    cntmin <- -alpha
    
    while (cntmin > minorigx) {
        cntmin <- cntmin * 10
        ticsl <- append(ticsl, cntmin)
        tics <- append(tics, -alpha - (alpha * log10(-cntmin/alpha)))
    }
    ticsl <- append(sort(ticsl), c(0, alpha))
    tics <- append(sort(tics), c(0, alpha))
    cntmax <- alpha
    while (cntmax < maxorigx) {
        cntmax <- cntmax * 10
        ticsl <- append(ticsl, cntmax)
        tics <- append(tics, alpha + (alpha * log10(cntmax/alpha)))
    }
    ## new x values as column 1
    ## new x labels as colum 2
    ret<-cbind(tics, ticsl)

    ret
}