Description Usage Arguments Details Value Examples
Find the rejection region
1 2 3 4 5 | rejectregion.npnorm(result, alpha = 0.05)
rejectregion.npt(result, alpha = 0.05)
rejectregion(result, alpha = 0.05)
|
result |
an object class nspmix |
alpha |
the FDR controlling rate. |
Find the rejection region based on the family in the result The rejection region is calculated using the density estimate rather than data points hence robust. The rejection is based on the hypothesis is located at 0. The optimisation is done via NLopt library (The package nloptr)
a list with par is the boundary for rejection and area is the propotion of rejection
1 2 3 4 5 6 7 | data = rnorm(500, c(0, 2))
x = makeobject(data, pi0 = 0.5)
r1 = computemixdist(x)
rejectregion(r1)
x2 = makeobject(data, pi0 = 0.5, method = "nptll") # equivalent to normal
r2 = computemixdist(x2)
rejectregion(r2)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.