https://github.com/cran/gss
Tip revision: 9924457bfed29635cbc74e54959beb6a433c8123 authored by Chong Gu on 23 September 2004, 00:00:00 UTC
version 0.9-3
version 0.9-3
Tip revision: 9924457
mkran.Rd
\name{mkran}
\alias{mkran}
\title{
Generating Random Effects in Mixed-Effect Models
}
\description{
Generate entries representing random effects in mixed-effect models.
}
\usage{
mkran(formula, data)
}
\arguments{
\item{formula}{Symbolic description of the random effects.}
\item{data}{Data frame containing the variables in the model.}
}
\details{
This function generates random effects terms from simple grouping
variables, for use in nonparametric mixed-effect models as described
in Gu and Ma (2003a, b). The syntax of the formula resembles that
of similar utilities for linear and nonlinear mixed-effect models,
as described in Pinheiro and Bates (2000).
Currently, \code{mkran} takes only two kinds of formulas,
\code{~1|grp2} or \code{~grp1|grp2}. Both \code{grp1} and
\code{grp2} should be factors, and for the second formula, the
levels of \code{grp2} should be nested under those of \code{grp1}.
The Z matrix is determined by \code{grp2}. When observations are
ordered according to the levels of \code{grp2}, the Z matrix is
block diagonal of 1 vectors.
The Sigma matrix is diagonal. For \code{~1|grp2}, it has one tuning
parameter. For \code{~grp1|grp2}, the number of parameters equals
the number of levels of \code{grp1}, with each parameter shared by
the \code{grp2} levels nested under the same \code{grp1} level.
}
\value{
A list of three components.
\item{z}{Z matrix.}
\item{sigma}{Sigma matrix to be evaluated through
\code{sigma$fun(para,sigma$env)}.}
\item{init}{Initial parameter values.}
}
\note{
One may pass a formula or a list to the argument \code{random} in
calls to \code{\link{ssanova1}} or\code{\link{gssanova1}} to fit
nonparametric mixed-effect models. A formula will be converted to a
list using \code{mkran}. A list should be of the same form as the
value of \code{mkran}.
}
\author{Chong Gu, \email{chong@stat.purdue.edu}}
\references{
Gu and Ma (2003a), \emph{Optimal Smoothing in Nonparametric
Mixed-Effect Models}. Available at
\url{http://stat.purdue.edu/~chong/manu.html}.
Gu and Ma (2003b), \emph{Generalized Nonparametric Mixed-Effect
Models: Computation and Smoothing Parameter Selection}. Available
at \url{http://stat.purdue.edu/~chong/manu.html}.
Pinheiro and Bates (2000), \emph{Mixed-Effects Models in S and
S-PLUS}. New York: Springer-Verlag.
}
\examples{
## Toy data
test <- data.frame(grp=as.factor(rep(1:2,c(2,3))))
## First formula
ran.test <- mkran(~1|grp,test)
ran.test$z
ran.test$sigma$fun(2,ran.test$sigma$env) # diag(10^(-2),2)
## Second formula
ran.test <- mkran(~grp|grp,test)
ran.test$z
ran.test$sigma$fun(c(1,2),ran.test$sigma$env) # diag(10^(-1),10^(-2))
## Clean up
\dontrun{rm(test,ran.test)}
}
\keyword{models}
\keyword{regression}