Description Slots Methods Author(s) See Also Examples
The class someMTP.object is the output of a call to
fdrOrd
. It also stores the information needed for related plots.
rej
:a logical vector indicating whenever the related hypotesis have been rejected.
p
:The vector of (raw) p-values used in the procedure.
ord
:The vector used to sort the p-values (decreasing).
idOrd
:The vector of indices used in sorting.
MTP
:The type of procedure used.
GD
:A logical value incating if the correction for General Dependence have been used or not.
q
:The level of contrelled FDR when MTP=="fdrOrd".
k
:The number of false rejection when MTP=="kfweOrd"
J
:The number of allowed Jumps when MTP=="kfweOrd"
alpha
:The significance level when MTP=="kfweOrd"
alphaprime
:The significance level of individual tests.
call
:The cal that generates the object.
someMTP.object: Prints the test results.
someMTP.object: Prints the test results (as show
).
someMTP.object: Plots results; what = c("all","ordVsP", "stepVsR")
signature(x = "someMTP.object")
: Sorts the p-values to decreasing order of ord
.
signature(x = "someMTP.object")
: The number of tests performed.
signature(x = "someMTP.object")
: Extracts the row names of the results matrix.
signature(x = "someMTP.object")
: Changes the row names of the results matrix. Duplicate names are not allowed, but see alias
.
Livio Finos: livio@stat.unipd.it
1 2 3 4 5 6 7 8 9 10 11 12 13 | # Simple examples with random data
set.seed(17)
x=matrix(rnorm(60),3,20)
x[,1:10]=x[,1:10]+2 ##variables 1:10 have tests under H1
ts=apply(x,2,function(x) t.test(x)$statistic)
ps=apply(x,2,function(x) t.test(x)$p.value)
m2=apply(x^2,2,mean)
pOrd <- fdrOrd(ps,q=.05,ord=m2)
pOrd
length(pOrd)
names(pOrd) <- paste("V",1:20,sep="")
names(pOrd)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.