多変量の見せ方-2 [4章]

Rによる計算機統計学

Rによる計算機統計学を読んでいる。 英語の原書はこちら。

多変量の見せ方

の続き。lattice packageを用いて、3D scatterplotを作る。irisデータを用いる。
公開されているコードから少し改変してトライ。

3D scatterplot

    library(lattice)
    attach(iris)
    #basic 3 color plot with arrows along axes
    print(cloud(Petal.Width ~ Sepal.Length * Sepal.Width,
          data=iris, groups=Species))

4C3=4つをそれぞれ配置する。

    print(cloud(Sepal.Length ~ Petal.Length * Petal.Width,
        data = iris, groups = Species, main = "A", pch=1:3,
        scales = list(draw = FALSE), zlab = "SL",
        screen = list(z = 30, x = -75, y = 0)),
        split = c(1, 1, 2, 2), more = TRUE)

    print(cloud(Sepal.Width ~ Petal.Length * Petal.Width,
        data = iris, groups = Species, main = "B", pch=1:3,
        scales = list(draw = FALSE), zlab = "SW",
        screen = list(z = 30, x = -75, y = 0)),
        split = c(2, 1, 2, 2), more = TRUE)

    print(cloud(Petal.Length ~ Sepal.Length * Sepal.Width,
        data = iris, groups = Species, main = "C", pch=1:3,
        scales = list(draw = FALSE), zlab = "PL",
        screen = list(z = 30, x = -55, y = 0)),
        split = c(1, 2, 2, 2), more = TRUE)

    print(cloud(Petal.Width ~ Sepal.Length * Sepal.Width,
        data = iris, groups = Species, main = "D", pch=1:3,
        scales = list(draw = FALSE), zlab = "PW",
        screen = list(z = 30, x = -55, y = 0)),
        split = c(2, 2, 2, 2))
    detach(iris)