speaq (version 1.2.0)

createNullSampling: Building a null hypothesis data

Description

Create a null sampling data (N times) and write them to a file

Usage

createNullSampling(X, groupLabel, N = 100, verbose=TRUE)

Arguments

X

The spectral dataset in the matrix format in which each row contains a single sample

groupLabel

Group label of samples in the dataset.

N

The number of iteration for creating null sample distribution

verbose

A boolean value to allow print out process information.

Value

• A matrix with N rows containing the null distribution.

Details

Create a null sampling data (N times) and write them to a file

Examples

res=makeSimulatedData();
X=res$data;
groupLabel=res$label;



peakList <- detectSpecPeaks(X,
  nDivRange = c(128),                
  scales = seq(1, 16, 2),
  baselineThresh = 50000,
  SNR.Th = -1,
    verbose=FALSE
);


resFindRef<- findRef(peakList);
refInd <- resFindRef$refInd;


maxShift = 50;
Y <- dohCluster(X,
                peakList = peakList,
                refInd = refInd,
                maxShift  = maxShift,
                acceptLostPeak = TRUE, verbose=FALSE);


# find the BW-statistic
BW = BWR(Y, groupLabel);

H0 = createNullSampling(Y, groupLabel, N = 100,verbose=FALSE)