Computation of Scores
This function can be used to compute several scores for a data vector.
## Default S3 method:
cscores(y, type=c("Data", "Wilcoxon", "NormalQuantile",
"AnsariBradley", "Median", "Savage", "ConSal"), int=FALSE,
maxs=length(y), ... )
## S3 method for class 'factor'
cscores(y, ...)
## S3 method for class 'Surv'
cscores(y, type="LogRank", int=FALSE, maxs=nrow(y), ...)y |
a numeric, factor or logical vector or an object of class
|
type |
a character string which specifies the type of the scores to be
computed. |
int |
a logical, forcing integer valued scores. |
maxs |
an integer defining the maximal value of the scores
if |
... |
additional arguments, not passed to anything at the moment. |
This function will serve as the basis for a more general framework of rank and permutation tests in future versions of this package. Currently, it is only used in the examples.
The logrank scores are computed as given in Hothorn & Lausen (2002).
If integer valued scores are requested (int = TRUE), the
scores are mapped into integers by
round(scores*length(scores)/max(scores)). See dperm for
more details.
type is self descriptive, except for ConSal which implements
scores suggested by Conover & Salsburg (1988).
A vector of scores for y with an attribute scores indicating
the kind of scores used is returned.
Torsten Hothorn <Torsten.Hothorn@rzmail.uni-erlangen.de>
Torsten Hothorn & Berthold Lausen (2003), On the exact distribution of maximally selected rank statistics. Computational Statistics \& Data Analysis, 43(2), 121-137.
William J. Conover & David S. Salsburg (1988), Locally most powerful tests for detecting treatment effects when only a subset of patients can be expected to "respond" to treatment. Biometrics, 44, 189-196.
y <- rnorm(50) # v.d. Waerden scores nq <- cscores(y, type="Normal", int=TRUE) # quantile for m=20 observations in the first group qperm(0.1, nq, 20)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.