Nothing
tt1st <- function(n,alpha,theta1,theta2,tol,itmax) {
ny <- n-1
del1 <- theta1*sqrt(n)
del2 <- theta2*sqrt(n)
err1 <- -alpha
c1 <- (del1+del2)/2
err2 <- NA
S8 <- 0
while (err1 < 0)
{ c1 <- c1-0.5
areac1_2 <- pt(c1,ny,del2)
if (is.na(areac1_2))
{ S8 <- 1
break }
h <- alpha + areac1_2
c2 <- qt(h,ny,del2)
err1 <- pt(c2,ny,del1) - pt(c1,ny,del1) - alpha
}
if (S8 == 0)
{ c1l <- c1
c1r <- c1 + 0.5
it <- 0
while (abs(err1) >= tol && it <= itmax)
{ it <- it + 1
c1 <- (c1l+c1r) / 2
h <- alpha + pt(c1,ny,del2)
c2 <- qt(h,ny,del2)
err1 <- pt(c2,ny,del1) - pt(c1,ny,del1) - alpha
if (err1 <= 0)
c1r <- c1 else
c1l <- c1
}
}
if (S8 == 1)
{ c2 <- qt(alpha,ny,del2)
c1 <- qt(1-alpha,ny,del1)
areac1_1 <- pt(c1,ny,del1)
areac2_1 <- pt(c2,ny,del1)
if (is.na(areac2_1)) areac2_1 <- 1
err1 <- areac2_1 - areac1_1 - alpha
areac1_2 <- pt(c1,ny,del2)
areac2_2 <- pt(c2,ny,del2)
if (is.na(areac1_2)) areac1_2 <- 0
err2 <- areac2_2 - areac1_2 - alpha
}
pow0 <- pt(c2,ny) - pt(c1,ny)
cat(" n =",n," alpha =",alpha," theta1 =",theta1," theta2 =",theta2,
" it =",it," c1 =",c1," c2 =",c2," ERR1 =",err1," ERR2 =",err2,
" POW0 =",pow0)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.