Description Usage Format Details Source References See Also Examples
The Chaetocnema dataset includes anatomical measures of 74 males of fleabeetles, of three different genera: Ch. Concinna (A, 21 cases), Ch. Heikertingeri (B, 31 cases) and Ch. Heptapotamica (C, 22 cases).
1 |
Data frame with 74 observations on the following 8 variables.
X10
width of the first joint of the first tarsus in microns (the sum of measurements for both tarsi)
X12
the same for the second joint
X14
the maximal width of the aedeagus in the fore-part in microns
X18
the front angle of the aedeagus (1 unit = 7.5 degrees)
X40
the maximal width of the head between the external edges of the eyes in 0.01mm
X48
the aedeagus width from the side in microns
species
species of the fleabeetles. A factor with
levels A
: Ch. Concinna, B
: Ch. Heikertingeri,
C
: Ch. Heptapotamica
area
colletion area of specimens. A factor with
levels A
: Environs of Uljianovsk, B
: Khvalynsk,
the Volga, C
: Pern, D
: Environs of Leningrad,
E
: the Ukraine, F
: Ashkhabd, Turkmenistan,
G
: France, H
: Ustj-Zilma, I
: Gabra, Abkhazia,
L
: Ussuri district, M
: Yakutsk district,
N
: Khabarovsk, O
: Germany,
P
: Environs of Lake Issyk-Kul, Kirghizia,
Q
: Alma-ata, Kazachstan, R
: Environs of Frunza, Kirghizia
The aedeagus is the median lobe of the male genital capsule that is surrounded by the phallobase and through which the sperm duct traverses.
A.A. Lubischew (1962) On the use of discriminant functions in taxonomy. Biometrics, 18:4, 455-477.
Hand, D.J., et al. (1994) A Handbook of Small Data Sets, London: Chapman & Hall, 254-255.
The dataset is also present in lubisch
.
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 37 38 39 | ## Not run:
data(chaetocnema)
lab <- chaetocnema[,7]
chaetocnema <- chaetocnema[,4:5]
set.seed(1234)
chaetosd <- apply(chaetocnema, 2, sd)
chaetocnema <- chaetocnema+cbind(rnorm(nrow(chaetocnema), 0,
chaetosd[1]/100), rnorm(nrow(chaetocnema), 0, chaetosd[2]/100))
gx <- seq(min(chaetocnema[,1]), max(chaetocnema[,1]), length=60)
gy <- seq(min(chaetocnema[,2]), max(chaetocnema[,2]), length=60)
y <- as.matrix(expand.grid(gx, gy))
tauvolume <- quantile.localdepth(chaetocnema, probs=0.20, use='volume')
system.time(depthcontours <- localdepth(x=chaetocnema, y=y,
tau=tauvolume, use='volume', method='simplicial'))
system.time(depthcontoursapprox <- localdepth(x=chaetocnema, y=y,
tau=tauvolume, use='volume', method='simplicial', type='approx'))
par(mfrow=c(2,2))
plot(chaetocnema, pch=as.character(lab), main='Exact Simplicial Depth')
contour(x=gx, y=gy, z=matrix(depthcontours$depth/depthcontours$max.depth,
nrow=length(gx), ncol=length(gy)), add=TRUE)
plot(chaetocnema, pch=as.character(lab), main='Approx Simplicial Depth')
contour(x=gx, y=gy,
z=matrix(depthcontoursapprox$depth/depthcontoursapprox$max.depth,
nrow=length(gx), ncol=length(gy)), add=TRUE)
plot(chaetocnema, pch=as.character(lab), main='Exact Simplicial Local Depth, Tau=0.2')
contour(x=gx, y=gy, z=matrix(depthcontours$localdepth/depthcontours$max.localdepth,
nrow=length(gx), ncol=length(gy)), add=TRUE)
plot(chaetocnema, pch=as.character(lab), main='Approx Simplicial Local Depth, Tau=0.2')
contour(x=gx, y=gy, z=matrix(depthcontoursapprox$localdepth/depthcontoursapprox$max.localdepth,
nrow=length(gx), ncol=length(gy)), add=TRUE)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.