tpaired1t.R

# Raw data
n     <- as.integer(round(runif(1, min=20, max=45)))
Id    <- 1:n
T0    <- round(runif(n, min=1, max=10))
EOT   <- round(T0 + runif(n, min=-4, max=1))
EOT   <-ifelse(EOT < 0, 1, EOT)
data  <- data.frame(Id, T0, EOT)

# Summary statistics
mu    <- round(runif(1, min=0.5, max=1.5), 1)
t0    <- mean(data$T0)
s_t0  <- sd(data$T0)
eot   <- mean(data$EOT)
s_eot <- sd(data$EOT)

d1    <- mean(data$T0 - data$EOT)
# d is a range of possible values
d     <- call('...', floor(d1*10)/10, ceiling(d1*10)/10)

s1_d  <- sd(data$T0 - data$EOT)
s_d   <- call('...', floor(s1_d*10)/10, ceiling(s1_d*10)/10)

tails <- "one-tailed"
alpha <- 0.05
incr <- FALSE            #if an increment is good (TRUE, e.g., in the WAIS) or bad (FALSE, e.g., in the HDRS)

# Export data for download
download <- function(fname)
{
    write.csv2(data, fname, row.names=FALSE)
}