sasef {lazy.sasef}

R Documentation

Calculation of SAS Estimable Functions

Description

Calculation of SAS Estimable Functions

Usage

sasef(obj = NULL, type = 3, print = 0)

Arguments

obj	A design object
type	The type of EF: I, II, or III
print	= 1 to print the result in SAS format = 2 to add EF in numerical format = 3 to add EF in fraction format

Details

A design object is a list of the design matrix (X), infomation list (info) and parameter list (param) where

X is the design matrix
and
info is a list whose length is the number of effects.

info[[i]] contains the following:
ename name of effect[i]
order order of effect[i]: 0, 1, ..., maxorder
elevels # of levels of effect[i]
df degrees of freedom of effect[i]
erange range of effect[i] in the (full rank) estimable functions
range=cbind(from,to) where
from is the starting independent colmn of effect[i] in X
to is the ending independent column of effect[i] in X
rangef=cbind(from,to) where
from is the starting colmn of effect[i] in X
to is the ending column of effect[i] in X
menum vector of main effect number involved in effect[i] contained vector of effect numbers which contains effect[i]
where effect number is defined according to the order of column of X.
contains vector of effects which effect[i] contais.
param is a list consisting of

type, drop, maxorder, sort, pattern

Information of estimable function may be added to this object later.

Estimable function in general form is the rows of
matSwp( t(X)%*%X ) %*% t(X)%*%X

Note that all-zero rows are included

Value

A design object with the following list attached to it:
ef=list(type, Q, ef, locnz)
where Q is the EF in numeric form
and ef is the EF in SAS format.

Examples

data <- du232[,1:3]
resnf <- design_mat( data,  type=-1, pattern="a", sort=0 )
temp <- sasef( resnf, type=1, print=1 )
temp <- sasef( resnf, type=2, print=1 )
temp <- sasef( resnf, type="III", print=1 )

---