NIRStat (version 1.0)

Slopetest: Slope statistics based Analysis for NIRS data.

Description

Estimate the slope statistics and conduct a nonparametric based test on the slope difference before transfuion and after trasfusion. If detection limit occurs at 15

Usage

Slopetest(Yvec,timevec,transfusionvec,SD_est=F,num.permu=1000)

Arguments

Yvec

The outcome of NIRS time series \(Y(t_{i})\) of length N ranging from 15 to 100.

timevec

The time index of NIRS time series \(t_{i}\) of length N.

transfusionvec

The 0/1 indicator of the transfusion status \(X(t_{i})\). \(X(t_{i})=0\) means the current time point is before transfusion and \(X(t_{i})=1\) means the current time point is after transfusion.

SD_est

Whether to estimate the SD of the SLOPE statistic for pre-transfusion and post-transfuion. Default value is FALSE.

num.permu

Number of permutation for permutation test. Default value is 1000.

Value

An R vector from Slopetest containing Slope statistics and Pvalue in the following order:

Slope.before

The estimated Slope statistic before transfusion.

Slope.after

The estimated Slope statistic after transfusion.

Slope.diff

The estimated Slope statistic difference between before transfusion and after transfusion.

Pvalue

The pvalue of testing the Slope difference to be zero or not.

SD_pre

SD of the Slope statistic for pre-transfusion. Optional, only when SD_est = TRUE.

SD_post

SD of the Slope statistic for post-transfusion. Optional, only when SD_est = TRUE.

Details

This function estimates the slope statistics before transfusion and after transfusion based on penalized regression spline method and tests the difference based on a within-band permutation approach. If there is detection limit occurs (15), it will impute the missed data based on a uniform distribution and estimate the slope statistics through a standard imputation approach. The statistical testing is conducted through a nested within-band permutation approach across all imputated datasets.

References

Guo, Y., Wang, Y., Marin, T., Kirk, E., Patel, R., Josephson, C. Statistical methods for characterizing transfusion-related changes in regional oxygenation using Near-infrared spectroscopy in preterm infants.

Examples

Run this code
# NOT RUN {
# Data Simulation
dat = data.frame(Y= rep(0,100),t=1:100,trans = c(rep(0,50),rep(1,50)))
dat$Y = apply(dat,1,function(x){rnorm(1,5*rnorm(1),6*exp(rnorm(1)))})
dat$Y = dat$Y + 15 - quantile(dat$Y,0.3) 
dat$Y[dat$Y<=15] = 15


# Estimate the Slope statistics of the NIRS data and test on the difference. 
Slopetest(dat$Y,dat$t,dat$trans,FALSE,100)
# }

Run the code above in your browser using DataCamp Workspace