plot.neLht: Confidence interval plots for linear hypotheses in natural...

Description Usage Arguments Details See Also Examples

View source: R/neLht.R

Description

Confidence interval plots for linear hypotheses in natural effect models.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
## S3 method for class 'neEffdecomp'
plot(x, level = 0.95, transf = identity, ylabels,
  yticks.at, ...)

## S3 method for class 'neLht'
plot(x, level = 0.95, transf = identity, ylabels, yticks.at,
  ...)

## S3 method for class 'neLhtBoot'
plot(x, level = 0.95, ci.type = "norm",
  transf = identity, ylabels, yticks.at, ...)

Arguments

x

an object of class neLht.

level

the confidence level required.

transf

transformation function to be applied internally on the (linear hypothesis) estimates and their confidence intervals (e.g. exp for logit or Poisson regression). The default is identity (i.e. no transformation).

ylabels

character vector containing the labels for the (linear hypothesis) estimates to be plotted on the y-axis.

yticks.at

numeric vector containing the y-coordinates (from 0 to 1) to draw the tick marks for the different estimates and their corresponding confidence intervals.

...

additional arguments.

ci.type

the type of bootstrap intervals required (see type argument in neModel-methods).

Details

This function is an adapted version of plot.glht from the multcomp package and yields confidence interval plots for each of the linear hypothesis parameters.

See Also

neModel, neLht, neEffdecomp

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
data(UPBdata)

impData <- neImpute(UPB ~ att * negaff + gender + educ + age, 
                    family = binomial, data = UPBdata)
neMod <- neModel(UPB ~ att0 * att1 + gender + educ + age, 
                 family = binomial, expData = impData, se = "robust")

lht <- neLht(neMod, linfct = c("att0 = 0", "att0 + att0:att1 = 0", 
                               "att1 = 0", "att1 + att0:att1 = 0", 
                               "att0 + att1 + att0:att1 = 0"))

## all pairs return identical output
plot(confint(lht), transf = exp)
plot(lht, transf = exp)

plot(neEffdecomp(neMod), transf = exp)
plot(neMod, transf = exp)

medflex documentation built on Sept. 15, 2017, 9:03 a.m.