mrbin (version 1.9.4)

atnv: A function replacing negative values.

Description

This function replaces (column-wise) negative values by a small positive number. The number is calculated as an affine transformation to the range of the lowest positive number to 0,01*the lowest positive number (of this column). Ranks stay unchanged. Positive numbers are not altered. If sample-wise noise levels are available, the median noise level of samples with negative values is calculated and replaces the lowest positive number in case it is smaller. If no noise data is available, the 1 positive values in the data set is used as an estimate. It is recommended to us this function AFTER noise removal and other data clean-up methods, as it may alter (reduce) the noise level.

Usage

atnv(NMRdata, noiseLevels = NULL, verbose = TRUE, errorsAsWarnings = FALSE)

Value

An invisible matrix or mrbin object containing NMR data without negative values.

Arguments

NMRdata

A matrix or mrbin object containing NMR data. For matrix: columns=frequencies,rows=samples

noiseLevels

A vector (can be omitted if NMRdata is an mrbin object)

verbose

Should a summary be displayed if NMRdata is an mrbin object

errorsAsWarnings

If TRUE, errors will be turned into warnings. Should be used with care, as errors indicate undocumented changes to the data.

Examples

resetEnv()#clean up previous data from the package environment 
results<-mrbin(silent=TRUE,
                   parameters=list(verbose=TRUE,dimension="1D",PQNScaling="No",
                   binwidth1D=0.005,signal_to_noise1D=1,PCA="No",binRegion=c(9.5,7.5,10,156),
                   saveFiles="No",referenceScaling="No",noiseRemoval="No",
                   fixNegatives="No",logTrafo="No",noiseThreshold=.05,NMRvendor="mrbin",
                   example=TRUE,#only used for the package examples
                   NMRfolders=c(system.file("extdata/2.mr1",package="mrbin"),
                              system.file("extdata/3.mr1",package="mrbin"))
                   ))
sum(results$bins<=0)
exampleNMRpositive<-atnv(NMRdata=results$bins, noiseLevels=results$parameters$noise_level_adjusted)
sum(exampleNMRpositive<=0)