sasefg {lazy.sasef} | R Documentation |
General Form of SAS Estimable Functions
sasefg(obj = NULL, print = 0)
obj |
A design object |
print |
= 1 to print the result |
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
G = matSwp( t(X)%*%X ) %*% t(X)%*%X
Note that all-zero rows are included
A list of G matrix, Gc and the location of non-zero rows of G,
where Gc is the G matrix in SAS format.
data <- du232[,1:3] resnf <- design_mat( data, type=-1, pattern="a", sort=0 ) temp <- sasefg( resnf, print=1 )