The R Project SVN – R – Diff – /branches/R-OBJSXP/tests/Examples/stats-Ex.Rout.save

-R Under development (unstable) (2023-09-10 r85120) -- "Unsuffered Consequences"
+R Under development (unstable) (2023-10-17 r85341) -- "Unsuffered Consequences"
 Copyright (C) 2023 The R Foundation for Statistical Computing
 Platform: x86_64-pc-linux-gnu
 R is free software and comes with ABSOLUTELY NO WARRANTY.
 You are welcome to redistribute it under certain conditions.
 +       zlab = "Pr(T <= t)")
+>
 > plot(function(x) dt(x, df = 3, ncp = 2), -3, 11, ylim = c(0, 0.32),
 +      main = "Non-central t - Density", yaxs = "i")
+>
+> ## Relation between F_t(.) = pt(x, n) and pbeta():
+> ptBet <- function(t, n) {
++     x <- n/(n + t^2)
++     r <- pb <- pbeta(x, n/2, 1/2) / 2
++     pos <- t > 0
++     r[pos] <- 1 - pb[pos]
++     r
++ }
+> x <- seq(-5, 5, by = 1/8)
+> nu <- 3:10
+> pt. <- outer(x, nu, pt)
+> ptB <- outer(x, nu, ptBet)
+> ## matplot(x, pt., type = "l")
+> stopifnot(all.equal(pt., ptB, tolerance = 1e-15))
+>
+>
+>
 > cleanEx()
 > nameEx("Tukey")
 > ### * Tukey
 > ### * ks.test
+>
 > flush(stderr()); flush(stdout())
+>
 > ### Encoding: UTF-8
->
 > ### Name: ks.test
 > ### Title: Kolmogorov-Smirnov Tests
 > ### Aliases: ks.test ks.test.default ks.test.formula
 > ### Keywords: htest
+>

The R Project SVN R