Sitka spruce data

Description

The sitka data frame contains measurements of log-size for 79 Sitka spruce trees grown in normal or ozone-enriched environments. Within each year, the data are organised in four blocks, corresponding to four controlled environment chambers. The first two chambers, containing 27 trees each, have an ozone-enriched atmosphere, the remaining two, containing 12 and 13 trees respectively, have a normal (control) atmosphere.

Usage

data(sitka)

Format

This data frame contains the following columns:

id.num

identification number of tree.

order

time order ranking within each tree.

days

time in days since 1st January, 1988.

log.size

tree size measured on a logarithmic scale.

ozone

indicator ozone treatment: 0=control,1=ozone.

Source

Diggle, P.J., Heagerty, P., Liang, K.-Y. and Zeger, S.L. (2002). Analysis of Longitudinal Data, Second Edition, Oxord: Oxford University Press.

References

Ruppert, D., Wand, M.P. and Carroll, R.J. (2003)
Semiparametric Regression Cambridge University Press.
http://stat.tamu.edu/~carroll/semiregbook/

Examples

library(SemiPar)
data(sitka)
attach(sitka)
library(lattice)
ozone.char <- rep("control",nrow(sitka))
ozone.char[ozone==1] <- "ozone"
xyplot(log.size~days|ozone.char,data=sitka,groups=id.num,type="b")

Results

R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(SemiPar)
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/SemiPar/sitka.Rd_%03d_medium.png", width=480, height=480)
> ### Name: sitka
> ### Title: Sitka spruce data
> ### Aliases: sitka
> ### Keywords: datasets
> 
> ### ** Examples
> 
> library(SemiPar)
> data(sitka)
> attach(sitka)
> library(lattice)
> ozone.char <- rep("control",nrow(sitka))
> ozone.char[ozone==1] <- "ozone"
> xyplot(log.size~days|ozone.char,data=sitka,groups=id.num,type="b")
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>