lme.value: Retrieve items from an mixed-model, a glm, or from a lm...

Description Usage Arguments Value Author(s) See Also Examples

Description

Use these functions in Sweave to retrieve individual factor contrast from a summary.lme, summary.glm or summary.lm object.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
## S3 method for class 'value'
lme(x, factor, signif = 2)

## S3 method for class 'stderr'
lme(x, factor, signif = 2)

## S3 method for class 'DF'
lme(x, factor)

## S3 method for class 'p'
lme(x, factor, signif = 2)

## S3 method for class 'ps'
lme(x, factor, signif = 2)

lm.value(x, factor, signif = 2)

lm.stderr(x, factor, signif = 2)

lm.t(x, factor, signif = 2)

lm.p(x, factor, signif = 2)

lm.ps(x, factor, signif = 2)

glm.value(x, factor, signif = 2)

glm.stderr(x, factor, signif = 2)

glm.z(x, factor, signif = 2)

glm.p(x, factor, signif = 2)

glm.ps(x, factor, signif = 2)

latexSNdouble(val)

Arguments

x

an object of class summary.lme.

factor

name of factor or interaction in contrast table. Shorthand for (Intercept)

signif

number of significant digits, default = 2.

val

a floating point value

Value

Returns a scalar rounded for use with Sweave and \Sexpr. Functions lme.ps, lm.ps and glm.ps return a string in latex math mode of the form $p=0.34$.

Author(s)

Dieter Menne, dieter.menne@menne-biomed.de

See Also

latex.lme, summary.lme

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
# Pinheiro/Bates page 47
library(nlme)
fm1Oats.sum <- summary(lme(yield~ordered(nitro)*Variety, data=Oats,
  random = ~1|Block/Variety))
lme.value(fm1Oats.sum,"VarietyMarvellous", signif = 3)
lme.stderr(fm1Oats.sum,"VarietyMarvellous", signif = 4)
lme.p(fm1Oats.sum,"VarietyMarvellous")

# For lm
ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14)
trt <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69)
sex <- as.factor(rep(c("f","m"),10))
group <- gl(2,10,20, labels=c("Ctl","Trt"))
weight <- c(ctl, trt)
x <- summary(lm(weight ~ group*sex))
lm.value(x,"sexm", signif = 1)
lm.stderr(x,"sexm")
lm.p(x,"sexm", signif = 4)
lm.ps(x,"sexm", signif = 4)

latexSNdouble(1.342E-12)

dmenne/dlatex documentation built on May 15, 2019, 9:32 a.m.