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
|
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/multipletesting.R
\name{closed.testing}
\alias{closed.testing}
\alias{p.correct}
\title{Closed testing procedure}
\usage{
closed.testing(
object,
idx = seq_along(coef(object)),
null = rep(0, length(idx)),
...
)
}
\arguments{
\item{object}{estimate object}
\item{idx}{Index of parameters to adjust for multiple testing}
\item{null}{Null hypothesis value}
\item{...}{Additional arguments}
}
\description{
Closed testing procedure
}
\examples{
m <- lvm()
regression(m, c(y1,y2,y3,y4,y5,y6,y7)~x) <- c(0,0.25,0,0.25,0.25,0,0)
regression(m, to=endogenous(m), from="u") <- 1
variance(m,endogenous(m)) <- 1
set.seed(2)
d <- sim(m,200)
l1 <- lm(y1~x,d)
l2 <- lm(y2~x,d)
l3 <- lm(y3~x,d)
l4 <- lm(y4~x,d)
l5 <- lm(y5~x,d)
l6 <- lm(y6~x,d)
l7 <- lm(y7~x,d)
(a <- merge(l1,l2,l3,l4,l5,l6,l7,subset=2))
if (requireNamespace("mets",quietly=TRUE)) {
p.correct(a)
}
as.vector(closed.testing(a))
}
|
