The Oslo data set contains chemical concentrations of 332 samples of
different plant species collected along a 120 km transect running through
the city of Oslo, Norway. It is a subset of the
OsloTransect
data provided by the rrcov
package.
Format
A data frame with 332 observations on the following 14 variables.
site
transect site ID, a factor with levels
102
103
104
105
106
107
108
109
111
112
113
114
115
116
117
118
119
121
122
123
124
125
126
127
128
129
131
132
133
134
135
136
138
139
141
142
143
144
XC
X coordinate, a numeric vector
YC
Y coordinate, a numeric vector
forest
forest type, a factor with levels
birspr
mixdec
pine
sprbir
sprpin
spruce
weather
weather type, a factor with levels
cloud
moist
nice
rain
litho
lithological type, a factor with levels
camsed
(Cambro-Silurian sedimentary),gneis_o
(Precambrian gneisses - Oslo),gneis_r
(- Randsfjord),magm
(Magmatic rocks)altitude
altitude, a numeric vector
Cu
Copper, a numeric vector
Fe
Iron, a numeric vector
K
Potassium, a numeric vector
Mg
Magnesium, a numeric vector
Mn
Manganese, a numeric vector
P
Lead, a numeric vector
Zn
Zinc, a numeric vector
Source
Reimann, C., Arnoldussen, A., Boyd, R., Finne, T.E., Koller, F., Nordgulen, Oe., And Englmaier, P. (2007) Element contents in leaves of four plant species (birch, mountain ash, fern and spruce) along anthropogenic and geogenic concentration gradients, The Science of the Total Environment, 377, 416-433.
Details
The OsloTransect
contains 360 observations, with 9
observations per site. Only 7 chemical elements were retained from the 25
contained in the OsloTransect
data, and these were all
log-transformed, following Todorov and Filzmoser (2009).
Only complete cases on these variables were retained, and two lithological types of low frequency were removed, leaving 332 observations.
References
Todorov V. and Filzmoser P. (2009) Robust statistic for the one-way MANOVA, submitted to the Journal of Environmetrics.
Examples
data(Oslo)
table(Oslo$litho)
#>
#> camsed gneis_o gneis_r magm
#> 98 89 32 113
Oslo.mod <- lm(cbind(Cu, K, Mg, Mn, P, Zn) ~ litho, data=Oslo)
car::Anova(Oslo.mod)
#>
#> Type II MANOVA Tests: Pillai test statistic
#> Df test stat approx F num Df den Df Pr(>F)
#> litho 3 0.24655 4.8503 18 975 1.503e-10 ***
#> ---
#> Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
heplot(Oslo.mod, var=c("Cu", "Mn"))
pairs(Oslo.mod)
if (FALSE) { # \dontrun{
if(require(candisc)) {
Oslo.can <- candisc(Oslo.mod)
Oslo.can
heplot(Oslo.can)
if(requireNamespace("rgl")){
heplot3d(Oslo.can, shade=TRUE, wire=FALSE, alpha=0.5, var.col="red")
}
}
} # }