model.matrix.pim: Create a model matrix for a probabilistic index model

In pimdevel: Probabilistic Index Models

Description

This function creates a model matrix for use in a probabilistic index model. This model matrix can be passed to pim.fit.

Usage

model.matrix(object, ...)

## S4 method for signature 'pim'
model.matrix(object, data, ...)

## S4 method for signature 'pim.formula'
model.matrix(object, data, ...)

Arguments

object	a pim.formula object that contains the formula necessary for constructing the model matrix.
...	extra arguments passed to or from other methods. This is currently only implemented in concordance with the generic model.matrix function.
data	an optional argument specifying the data frame for which the model matrix should be constructed. See also model.matrix) in the stats package.

Value

a design matrix for a pim model

Examples

data("FEVData")
# Create the "model frame"
FEVenv <- new.pim.env(FEVData, compare="unique")
# This includes the poset
pos <- poset(FEVenv, as.list=TRUE)

# create the formula and bind it to the pim.environment.
FEVform <- new.pim.formula(
  Age ~ I(L(Height) - R(Height))  ,
  FEVenv
)

# Use this formula object to construct the model matrix
MM <- model.matrix(FEVform)

# Use this formula object to construct the pseudo response
Y <- response(FEVform)

# Now pim.fit can do what it does
res <- pim.fit(MM,Y, estim = "estimator.glm", penv=FEVenv)

pimdevel documentation built on May 2, 2019, 4:58 p.m.