Learn R Programming

sse (version 0.7-17)

refine: Refining the estimation

Description

Increasing the number of iterations for estimating the sample size for the 'theta' and 'xi' as specified for the example.

Usage

refine(object, factor = 10)

Arguments

object

An object of class '>power.

factor

An integer larger than one that is multiplied with the available number of iterations to from the target number of iterations.

Value

An object of class '>power.

Examples

Run this code
# NOT RUN {
 ## takes quite some time
## defining the range of n and theta to be evaluated
psi <- powPar(theta = seq(from = 0.5, to = 1.5, by = 0.1),
              n = seq(from = 20, to = 60, by = 2))

## defining a power-function               
powFun <- function(psi){
  x <- rnorm(n(psi)/2)
  y <- rnorm(n(psi)/2) + theta(psi)
  return(wilcox.test(x = x, y = y)$p.value < 0.05)
}


## evaluating the power-function for all combinations of n and theta
calc <- powCalc(psi, powFun, n.iter = 10)

## adding example at theta of 1 and power of 0.9
pow <- powEx(calc, theta = 1, power = 0.9)

## another 900 (= 1000 - 100) iterations
refine(pow)
# }

Run the code above in your browser using DataLab