R/calcAniSignal.R

In TIMP: Fitting Separable Nonlinear Models in Spectroscopy and Microscopy

"calcAniSignal" <- function(m, plotoptions) {
  # need to implement analysis for more than 3 datasets simult.
  perpind <- which(m[[1]]@anispec$angle == "PERP")[1]
  parind <- which(m[[1]]@anispec$angle == "PAR")[1]

  anisig <- (m[[parind]]@psi.df - m[[perpind]]@psi.df) /
    (m[[parind]]@psi.df + (2 * m[[perpind]]@psi.df))

  write.table(anisig,
    file = paste(plotoptions@makeps,
      "calculated_ani_signal.txt",
      sep = ""
    ),
    row.names = m[[1]]@x,
    col.names = m[[1]]@x2, quote = FALSE
  )
}