{
tol <- 1e-6
T
}
{
#fac.design
all(dim(fac.design(rep(2,3)))==c(8,3))
}
{
#factor
class(factor(rep(1:3,4),1:3,c("dog","rat","cat")))=="factor"
}
{
#factor.names
any(names(factor.names(pigment))=="Batch")
}
{
#factor.names<-
xdf <- pigment
factor.names(xdf)<-c("batch","sample","test","moisture")
any(names(factor.names(xdf))=="test")
}
{
#factor.scope              # Support for step.glm
exists("factor.scope")
}
{
#family.glm
k <<- kyphosis[1:20,]
zkyph <- glm(Kyphosis=="present" ~ Age + Number,family=binomial,data=k)
family.glm(zkyph)$family[[1]]=="Binomial"
}
{
#family.lm
zfree <- lm(freeny.y~freeny.x)
family.lm(zfree)$family[[1]]=="Gaussian"
}
{
#family
family(zfree)$family[[1]]=="Gaussian"
}
{
#family.default
family.default(NULL)$family[[1]]=="Gaussian"
}
{
#fft
is.complex(fft(ship))
}
{
#filter
flynx <- filter(lynx,c(.3, .6, .1))
T
#var(flynx[!is.na(flynx)]) < var(lynx)
}
{
#find
is.character(find("sin"))
}
{
#find.doc                # No HELP file
T
}
{
#find.objects            
        platf <- platform()
	PAT <- if (platf == "DOS386" || platf == "WIN386") "*print*" else "print"
	a <- find.objects(pattern=PAT)
	mode(a) == "character" && length(a) > 0
}
{
#fisher.test
x <- factor(c(1, 2, 1, 1), 1:2, c("a", "b"))
y <- factor(c(2, 1, 1, 2), 1:2, c("c", "d"))
xy.tab <- table(x,y)
abs(fisher.test(xy.tab)$p.value - 1) < tol
}
{
#fitted
is.ts(fitted(zfree))
}
{
#fitted.default
is.ts(fitted.default(zfree))
}
{
#fitted.values         # Same as fitted()
T
}
{
#fix.NA                 
platform() == "DOS386" || platform() == "WIN386" || as.logical(fix.NA(1))
}
{
#floor
abs(floor(-1.9)+2) < tol 
}
{
#format
is.character(format(2^.5))
}
{
#formula
class(formula(a~b+c))=="formula"
}
{
#formula.data.frame
length(formula.data.frame(kyphosis))==3
}
{
#formula.default
length(formula.default(zkyph))==3
}
{
#formula.formula
class(formula.formula(a~b+c))=="formula"
}
{
#formula.lm
length(formula.lm(zfree))==3
}
{
#formula.terms
T
}
{
#frame.attr                 # What are these??
#frame.attr<-
#frame.attributes
T
}
{
#frequency
abs(frequency(ship)-12) < tol
}
{
#friedman.test
trmt <- factor(rep(1:3,c(4,4,4)),1:3,c("Trt 1","Trt 2","Trt 3"))
ppl <- factor(rep(1:4,3),1:4,c("Subject 1","Subject 2","Subject 3","Subject 4"))
y <- c(0.73,0.76,0.46,0.85,0.48,0.78,0.87,0.22,0.51,0.03,0.39,0.44)
abs(friedman.test(y, trmt, ppl)$p.value -  0.3678794) < tol
}
{
#gam
k <<- kyphosis[1:20,]
kyph.gam <- gam(Kyphosis ~ s(Age) + Number, family = binomial, data=k)
abs(kyph.gam$deviance - 6.678695) < 1e-2
}
{
#anova.gam
all(anova.gam(kyph.gam)$Df==1)
}
{
#summary.gam
summ <- summary.gam(kyph.gam)
abs(summ$deviance - 6.678695) < 1e-2
}
{
#predict.gam
length(predict.gam(kyph.gam))==20
}
{
#gam.control
gam.control()$maxit == 10
}
{
#gam.lo               # Support for Function lo()
#gam.match            # Support for Functions s(), lo()
#gam.nlchisq          # Support for Function gam()
#gam.random           # Support for Function random()
#gam.s                # Support for Function s()
T
}
{
#gam.scope
kyph.scope <- gam.scope(k)
class(kyph.scope$Start)=="formula"
}
{
#gam.slist
length(gam.slist)==3
}
{
#gam.sp               # Support for Function predict.gam()
T
}
{
#gam.wlist
length(gam.wlist)==2
}
{
#gamlist              # Support for Functions anova(), anova.gam()
kyph2.gam <- gam(Kyphosis ~ s(Age),family=binomial,data=k)
kyph.list <- gamlist(kyph.gam,kyph2.gam)
abs(anova.gamlist(kyph.list)$Deviance[2] + 3.675451) < 1e-2
}
{
#gamma
abs(gamma(5)-24) < tol
}
{
#gaussian
any(names(gaussian())=="family")
}
{
#get
#get.default
some.data <- c(2,4,4,5)
same.data <- get("some.data")
all(some.data==same.data)
}
{
#getenv
is.character(getenv("SHOME"))
}
{
#glim
liver.x <- cbind(injection=liver.cells,
         A.exper=ifelse(levels(liver.exper)[liver.exper]=="A", 1, 0),
         B.exper=ifelse(levels(liver.exper)[liver.exper]=="B", 1, 0),
         C.exper=ifelse(levels(liver.exper)[liver.exper]=="C", 1, 0))
liver.glim <- glim(liver.x, liver.gt[,1], error="poisson",
         link="log", intercept=F, seq=c(0,1,4))
abs(liver.glim$deviance[3] - 245.6395269) < tol
}
{
#glim.print
T
}
{
#glm
index <- 1:90 * 10
soldersmall <<- solder[index,]
solder.glm <- glm(skips ~ Opening + Solder, data = soldersmall, family = poisson)
abs(solder.glm$deviance - 377.66511221 ) < tol
}
{
#deviance.glm
abs(deviance.glm(solder.glm) - 377.665112217) < tol
}
{
#family.glm
family(solder.glm)$family[[1]]=="Poisson"
}
{
#predict.glm
length(predict.glm(solder.glm))==90
}
{
#residuals.glm
length(residuals.glm(solder.glm))==90
}
{
#glm.control
glm.control()$maxit == 10
}
{
#glm.fit            # Support for glm()
exists("glm.fit")
}
{
#glm.links          # Family Support Objects
all(dim(glm.links)==c(5,8))
}
{
#glm.null           # Support for Function gam()
exists("glm.weight")
}
{
#glm.variances      # Family Support Objects
all(dim(glm.variances)==c(3,5))
}
{
#glm.weight         # Support for Function make.family()
exists("glm.weight")
}
{
#glm.weights        # Family Support Objects
all(dim(glm.weights)==c(8,5))
}
{
#anova.glm
solder2.glm <- update(solder.glm, ~ . + Mask, data=soldersmall)
abs(anova.glm(solder.glm,solder2.glm)$Deviance[2] -  189.292775339) < tol
}
{
#glmlist
solder.list <- glmlist(solder.glm,solder2.glm)
abs(anova.glmlist(solder.list)$Deviance[2] - 189.292775339) < tol
}
{
#step.glm
solder.step <- step.glm(solder.glm,trace=F,data=soldersmall)
abs(solder.step$null.deviance - 741.7059514125) < tol
}
{
#summary.glm
summ <- summary.glm(solder.glm)
sum(summ$df)==90
}
{
#grep
platf <- platform()
if (platf != "DOS386" && platf != "WIN386") state.name [grep("Neb", state.name)]=="Nebraska" else
state.name [grep("Neb*", state.name)]=="Nebraska"
}
{
#gui.map              # Support for help.start()
platf <- platform()
if (platf != "DOS386" && platf != "WIN386") exists("gui.map") else T
}
{
#hat
h <- hat(freeny.x)
abs(sum(h) - 5) < tol
}
{
#hclust
h <- hclust(dist(freeny.x))
max(h$order)==39
}
{
#help
#help.findsum
#help.keyword
#help.name
#help.off
#help.running
#help.search
#help.send
#help.start
#help.topic
T
}
{
#history                    # Tricky!
T
}
{
#ifelse
ifelse(1>0,T,F)
}
{
#info, unix only
platf <- platform()
if (platf == "DOS386" || platf == "WIN386") T else any(names(info(file="")$Sinfo)=="SHOME")
}
{
#inherits
inherits(guayule,"design")
}
{
#integer
xi <- integer(1)
xi < tol
}
{
#interaction
length(interaction(catalyst,drop=T))==8
}
{
#interactive; we may be running LOOP with stderr redirected, or not.
#so don't rely on interactive() returning TRUE.
mode(interactive()) == "logical"
}
{
#interp
i <- interp(freeny.x[,1],freeny.x[,2],freeny.x[,3])
length(i)==3
}
{
#inverse.gaussian
inverse.gaussian()$family[[1]] == "Inverse Gaussian"
}
{
#invisible
xi <- invisible(sin(pi))
xi < tol
}
{
#is.array
is.array(iris)
}
{
#is.atomic
is.atomic(1)
}
{
#is.call
my.call <- call("sin", as.name("x"))
is.call(my.call)
}
{
#is.category
is.category(pigment$Test)
}
{
#is.character
is.character(letters)
}
{
#is.complex
is.complex(0+1i)
}
{
#is.data.frame
is.data.frame(kyphosis)
}
{
#is.dir
is.dir(search()[1])
}
{
#is.double
is.double(pi)
}
{
#is.double.factor           #Always FALSE
T
}
{
#is.expression
meow <- expression(cat)
is.expression(meow)
}
{
#is.factor
is.factor(pigment$Sample)
}
{
#is.finite
is.finite(1230)
}
{
#is.function
is.function(sin)
}
{
#is.inf
is.inf(1/0)
}
{
#is.infinite
is.infinite(1/0)
}
{
#is.infinite.warning
is.null(is.infinite.warning(3))
}
{
#is.integer
is.integer(length(lynx))
}
{
#is.integer.factor           #It is always FALSE
T
}
{
#is.language
meow <- expression(cat)
is.language(meow)
}
{
#is.list
is.list(guayule)
}
{
#is.loaded
is.loaded(symbol.For("glmfit"))
}
{
#is.logical
is.logical(2>1)
}
{
#is.matrix
is.matrix(freeny.x)
}
{
#is.na.data.frame
any(!is.na.data.frame(gun))
}
{
#is.na.factor
any(!is.na.factor(pigment$Batch))
}
{
#is.na.warning
is.null(is.na.warning(2))
}
{
#is.name
myls <- as.name("ls")
is.name(myls)
}
{
#is.nan
is.nan(0/0)
}
{
#is.nan.warning
is.null(is.nan.warning(2))
}
{
#is.not.finite.warning
is.null(is.not.finite.warning(2))
}
{
#is.not.number.warning
is.null(is.not.number.warning(2))
}
{
#is.not.orderable.warning
is.null(is.not.orderable.warning(2))
}
{
#is.null
is.null(guayule$madeup)
}
{
#is.number
is.number(3)
}
{
#is.numeric
is.numeric(3)
}
{
#is.one
is.one(1)
}
{
#is.orderable
is.orderable(4)
}
{
#is.qr
free.qr <- qr(freeny.x)
is.qr(free.qr)
}
{
#is.recursive
is.recursive(guayule)
}
{
#is.single
is.single(as.single(3))
}
{
#is.ts
is.ts(lynx)
}
{
#is.uminus
meow<-call("-","x")
is.uminus(meow)
}
{
#is.vector
is.vector(letters)
}
{
#is.zero
is.zero(0)
}