{
tol <- 1e-6
T
}
{
#page                        # ????????
exists("page")
}
{
#param<-
myframe <- kyphosis
param(myframe, "lambda") <- 1.5
abs(attr(myframe,"parameters")$lambda - 1.5  )< tol
}
{
#param
abs(param(myframe, "lambda") - 1.5  )< tol
}
{
#parameter.names             # Support for Functions ms(), nls()
exists("parameter.names")
}
{
#parameters
abs(parameters(myframe)$lambda - 1.5  )< tol
}
{
#parameters<-
pp <- list(mu=2,lambda=1.5)
parameters(myframe)<-pp
abs(attr(myframe,"parameters")$mu - 2  )< tol
}
{
#parse
eval(parse(text="x<-T"))
}
{
#paste
x <- paste("sin","(","pi","/2",")")
y <- eval(parse(text=x))
as.logical(y)
}
{
#path.tree                   # Tested after tree()
T
}
{
#pattern                     # Support for Function alias.lm()
T
}
{
#pbeta
abs(pbeta(.5,1,1) - .5  )< tol
}
{
#pbinom
abs(pbinom(.5,10,.2) - 0.1073741824 ) < tol
}
{
#pcauchy
abs(pcauchy(.5,0,10) - 0.5159022513 ) < tol
}
{
#pchisq
abs(pchisq(12,12) - 0.5543203586 ) < tol
}
{
#peaks
pp <- peaks(lynx)
abs(sum(pp) -14 ) < tol
}
{
#pexp
abs(pexp(1) - 0.6321205588 ) < tol
}
{
#pf
abs(pf(1,4,12) - 0.5550537109 ) < tol
}
{
#pframe                      #No HELP file
T
}
{
#pgamma
abs(pgamma(1,2) - 0.26424111766 ) < tol
}
{
#pgeom
abs(pgeom(2,.5) - 0.875 ) < tol
}
{
#phyper
abs(phyper(3,4,6,7) - 0.83333333333333 ) < tol
}
{
#pipegraph                  # Support for graphics devices (e.g. dev.motif)
T
}
{
#platform
pp <- platform()
is.character(pp)
}
{
#plnorm
as.logical( 2 * plnorm(1) )
}
{
#plogis
as.logical( 2 * plogis(0) )
}
{
#plotlabels
#plotlabels.default
pp <- plotlabels.default(freeny.y)
length(pp) == 39
}
{
#plotlabels.coef
pp <- plotlabels.coef(freeny.y)
length(pp) == 39
}
{
#pmatch
as.logical(pmatch("Alab",state.name))
}
{
#pmax
pp <- pmax(freeny.x[,1],freeny.x[,2])
length(pp) == 39
}
{
#pmin
pp <- pmin(freeny.x[,1],freeny.x[,2])
length(pp) == 39
}
{
#pnbinom
abs(pnbinom(20,10,.5) - 0.97861302737 ) < tol
}
{
#pnorm
abs(pnorm(1.644854) - .95 ) < tol
}
{
#pointwise
gas.m <- loess(NOx ~ E, data = gas, span = 2/3)
attach(gas)
gas.se <- predict(gas.m, newdata = seq(min(E), max(E),length = 7),se.fit = TRUE)
pointwise(gas.se)
detach("gas")
all.equal(pmax(gas.se$lower,gas.se$upper) , gas.se$upper)
}
{
#poisson
poisson()$family[1] == "Poisson"
}
{
#poly
ppll <- poly(liver.cells,2)
all(dim(ppll) == c(52,2))
}
{
#poly.matrix                 # Support for Functions lo(), poly()
#poly.raw                    # Support for Function contr.poly()
T
}
{
#polyroot
a <- c(1.2, .5, .3, 1)
root <- polyroot(c(1, -a))
root[1] < 1
}
{
#power
quafam <- quasi(link=power(1/3))
quafam$family[1] == "Quasi-likelihood"
}
{
#ppois
abs(ppois(.5,1) - 0.3678794412  )< tol
}
{
#ppreg
round(abs(ppreg(freeny.x,freeny.y,2)$fl2 - 0.0008327241521 ), 3) - .001 < 1e-6
}
{
#prcomp
prim.pr <- prcomp(prim4)
abs(prim.pr$sdev[1] - 5.1387305260 ) < tol
}
{
#pred.tree                  # Support for predict.tree
T
}
{
#predict.gam                # Tested after gam()
#predict.glm                # Tested after glm()
#predict.lm                 # Tested after lm()
#predict.loess              # Tested after loess()
#predict.mlm                # Tested after lm()
#predict.smooth.spline      # Tested after smooth.spline()
T
}
{
#predict.smooth.spline.fit  # Support for Function predict.smooth.spline()
exists("predict.smooth.spline.fit")
}
{
#predict.tree
zfuel <- tree(formula = Mileage ~ Weight, data=fuel.frame)
length(predict.tree(zfuel))==60
}
{
#preplot.gam
gamobj <- gam(Kyphosis ~ s(Age,4), family = binomial, data=kyphosis)
preobj <- preplot.gam(gamobj)
class(preobj) == "preplot.gam"
}
{
#preplot.loess
etha.lo <- loess(NOx ~ C * E, span = 1/2, degree = 2,
parametric = "C", drop.square = "C",data=ethanol)
preobj <- preplot.loess(etha.lo)
class(preobj) == "preplot.loess"
}
{
#pretty
as.logical(max(pretty(c(0,1),5)))
}
{
#proc.time, vector has different lengths unix and dos
is.numeric(proc.time())
}
{
#process.id
platf <- platform()
if (platf == "DOS386" || platf == "WIN386") T else is.integer(process.id()$pid)
}
{
#prod
prod(1:4)==24
}
{
#profile.default             # Just a warning
#profile.ms                  # Tested right after ms()
#profile.nls                 # Tested right after nls()
T
}
{
#prog.name                   # No HELP file
is.character(prog.name())
}
{
#proj.default
gunaov.qr <- aov(Rounds ~ Method + Physique/Team, gun, qr = T)
gunaov.proj <- proj.default(gunaov.qr)
all(dim(gunaov.proj)==c(36,10))
}
{
#proj.lm
lm.object <- lm(cost ~ age + type + car.age, claims, na.action=na.omit)
tmp <- proj.lm(lm.object)
all(dim(tmp)==c(123,5))
}
#{
##prompt
#blah <<- function(x) x
#cat("\nIGNORE prompt() messages. vvv\n")
#prompt(blah)
#cat("IGNORE prompt() messages. ^^^\n\n")
#ans <- is.readable("blah.d")
#if (ans) unlink("blah.d")
#ans
#}
{
#prop.test
smokers <- c( 83, 90, 129, 70 )
patients <- c( 86, 93, 136, 82 )
ppt <- prop.test( smokers, patients )
abs(ppt$statistic - 12.6004113) < tol
}
{
#prune.tree
z.pruned <- prune.tree(zfuel,k=25)
class(z.pruned) == "tree"
}
{
#psignrank                   # No HELP file.  Support for htest functions??
abs(psignrank(78,12) - 1) < 1e-5
}
{
#pstab                       # No stable distribution function
T
}
{
#pt
abs(pt(1,2)-0.7886751346)<tol
}
{
#punif
abs(punif(.5)-.5)<tol
}
{
#pweibull
abs(pweibull(.5,2) - 0.2211992169 ) < tol
}
{
#pwilcox
abs(pwilcox(24, 4, 6) - 0.6952380952 ) < tol
}
