sasef {lazy.sasef}R Documentation

Calculation of SAS Estimable Functions


Calculation of SAS Estimable Functions


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



A design object


The type of EF: I, II, or III


= 1 to print the result in SAS format
= 2 to add EF in numerical format
= 3 to add EF in fraction format


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

X is the design matrix
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


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.


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 )

[Package lazy.sasef version 0.1.2 Index]