Learn R Programming

scan (version 0.20)

randSC: Randomization Tests for single-case data

Description

The randSC function computes a randomization test for single or multiple baseline single-case data. The function is based on an algorithm from the SCRT package (Bulte & Onghena, 2009, 2012), but rewritten and extended for the use in AB designs.

Usage

randSC(data, statistic = "Mean B-A", number = 500, complete = FALSE, limit = 5, 
       startpoints = NA, exclude.equal = FALSE, graph = FALSE, output = "c")

Arguments

data

A single-case data frame or a list of single-case data frames. See makeSCDF to learn about this format.

statistic

Defines the statistic on which the comparison of phases A and B is based on. Default setting is statistic = "Mean B-A"). The following comparisons are possible:

  • "Mean A-B": Uses the difference between the mean of phase A and the mean of phase B. This is appropriate if a decrease of scores was expected for phase B.

  • "Mean B-A": Uses the difference between the mean of phase B and the mean of phase A. This is appropriate if an increase of scores was expected for phase B.

  • "Mean |A-B|": Uses the absolute value of the difference between the means of phases A and B.

  • "Median A-B": The same as "Mean A-B", but based on the median.

  • "Median B-A": The same as "Mean B-A", but based on the median.

number

Sample size of the randomization distribution. The exactness of the p-value can not exceed \(1/number\) (i.e., number = 100 results in p-values with an exactness of one percent). Default is number = 500. For faster processing use number = 100. For more precise p-values set number = 1000.

complete

If TRUE, the distribution is based on a complete permutation of all possible starting combinations. This setting overwrites the number Argument. The default setting is FALSE.

limit

Minimal number of data points per phase in the sample. Default is limit = 5.

startpoints

Alternative to the limit-parameter startpoints exactly defines the possible start points of phase B (e.g., startpoints = 4:9 restricts the phase B start points to measurements 4 to 9. startpoints overruns the limit-parameter.

exclude.equal

If set to exclude.equal = FALSE, which is the default, random distribution values equal to the observed distribution are counted as null-hypothesis conform. That is, they decrease the probability of rejecting the null-hypothesis (increase the p-value). exclude.equal should be set to TRUE if you analyse one single-case design (not a multiple baseline data set) to reach a sufficient power. But be aware, that it increases the chance of an alpha-error.

graph

If graph = TRUE, a histogram of the resulting distribution is plotted. It's FALSE by default.

output

If set to the default output = "C", detailed information is provided. Set output = "p", to only return the resulting p value.

Value

statistic

Character string from function call (see Arguments above).

N

Number of single-cases.

n1

Number of data points in phase A.

n2

Number of data points in phase B.

limit

Numeric from function call (see Arguments above).

startpoints

A vector defining the start points passed from the function call (see Arguments above).

p.value

P-value of the randomization test for the given data.

number

Sample size of randomization distribution from function call (see Arguments above).

complete

Logical argument from function call (see Arguments above).

observed.statistic

Test statistic observed for the given single-case data. (see statistic in the Arguments above.)

Z

Z-value of observed test statistic.

p.z.single

Probability of z-value.

distribution

Test statistic distribution from randomized data sets.

possible.combinations

Number of possible combinations under the given restrictions.

auto.corrected.number

TRUE indicates that a corrected number of combinations was used. This happens, if the number of possible combinations (under the given restrictions) undercuts the requested number of combinations.

References

Bulte, I., & Onghena, P. (2009). Randomization tests for multiple-baseline designs: An extension of the SCRT-R package. Behavior Research Methods, 41, 477-485.

Bulte, I., & Onghena, P. (2012). SCRT: Single-Case Randomization Tests. Available from: https://CRAN.R-project.org/package=SCRT

Examples

Run this code
# NOT RUN {
## Compute a randomization test on the first case of the byHeart2011 data and include a graph
randSC(byHeart2011[1], statistic = "Median B-A", graph = TRUE)

## Compute a randomization test on the Grosche2011 data using complete permutation
randSC(Grosche2011, statistic = "Median B-A", complete = TRUE, limit = 4)
# }

Run the code above in your browser using DataLab