File: pLausen94.Rd

package info (click to toggle)
r-cran-maxstat 0.7-25-2
  • links: PTS, VCS
  • area: main
  • in suites: bookworm, bullseye, forky, sid, trixie
  • size: 332 kB
  • sloc: ansic: 331; sh: 51; makefile: 2
file content (71 lines) | stat: -rw-r--r-- 1,846 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
\name{pLausen94}
\alias{pLausen94}
\alias{qLausen94}
\title{ Approximating Maximally Selected Statistics }
\description{

Approximates the probability that a maximally selected rank statistic is
greater or equal to \code{b}.

}
\usage{
pLausen94(b, N, minprop=0.1, maxprop=0.9, m=NULL)
qLausen94(p, N, minprop=0.1, maxprop=0.9, m=NULL)
}
\arguments{
  \item{b}{ quantile. }
  \item{p}{ probability. }
  \item{N}{ number of observations.}
\item{minprop}{ at least \code{minprop}*100\% of the observations in the
first group. }
  \item{maxprop}{ not more than \code{minprop}*100\% of the observations in
the first group. }
  \item{m}{a integer vector containing the sample sizes in the first groups
for each cutpoint considered. If \code{is.null(m)} a continuous predictor is
assumed. }
}
\details{

Approximation based on an improved Bonferroni inequality.

}
\value{
The probability that, under the hypothesis of independence, a maximally
selected statistic greater equal \code{b} is observed.
}
\references{ 

	Worsley, K.J. (1982),
	An Improved Bonferroni Inequality and Applications.
	\emph{Biometrika}, \bold{69}, 297--302

	Lausen, B. (1990),
	Maximal Selektierte Rangstatistiken. Dissertation. 
	Universit\"at Dortmund 

        Lausen, B., Sauerbrei, W. & Schumacher, M. (1994).
        Classification and Regression Trees (CART) used for the
        exploration of prognostic factors measured on different scales.
        in: P. Dirschedl & R. Ostermann (Eds),
        \emph{Computational Statistics}, Heidelberg, Physica-Verlag,
        483--496

}

\examples{

p <- pLausen94(2.5, 20, 0.25, 0.75)

# Lausen 94, page 489

if (round(p, 3) != 0.073) stop("error checking pLausen94")

# the same

p2 <- pLausen94(2.5, 200, 0.25, 0.75, m=seq(from=50, to=150, by=10))

stopifnot(all.equal(round(p,3), round(p2,3)))

}

\keyword{distribution}