cross$ export LD_LIBRARY_PATH=/APPS/ggobi/lib
cross$ R
> library("Rggobi")
>
> data(mtcars)
>
> class(mtcars)
[1] "data.frame"
> names(mtcars)
 [1] "mpg"  "cyl"  "disp" "hp"   "drat" "wt"   "qsec" "vs"   "am"   "gear"
[11] "carb"
> help(mtcars)
> ggobi(mtcars)
[1] "ggobi reference (1)"

> help(identify)
> n <- 20
> x <- rnorm(n)
> y <- x + rt(n, 4)
> plot(x, y)
> identify(x, y)
warning: no point with 0.25 inches
[1] 10 15 16 18
> i <- identify(x, y)
> i
[1]  3  7 11
> plot(x[-i], y[-i])
> points(x[i], y[i], col = "red")
> points(x[i], y[i], col = "red", pch = 1)
> points(x[i], y[i], col = "red", pch = 2)

 plot(x[-i], y[-i], xlab="x", ylab="y", main="some outliers deleted")

> plot(c(0, 1), c(0,1), type = "n")
> legend(locator(1), as.character(0:18), pch = 0:18)
> plot(c(0, 1), c(0,1), type = "n")
> legend(locator(2), as.character(0:18), pch = 0:18)

> density(x)

Call:
        density(x = x)

Data: x (300 obs.);     Bandwidth 'bw' = 0.5253

       x                y
 Min.   :-4.829   Min.   :2.932e-05
 1st Qu.:-1.624   1st Qu.:7.112e-03
 Median : 1.581   Median :6.322e-02
 Mean   : 1.581   Mean   :7.792e-02
 3rd Qu.: 4.786   3rd Qu.:1.178e-01
 Max.   : 7.991   Max.   :2.548e-01
> plot(density(x))
> plot(density(x, bw = 3))
> plot(density(x, bw = 1 / 3))
> plot(density(x, bw = 1 / 10))
> plot(density(x, bw = "ucv"))
> par(mfrow = c(2, 1))
> plot(density(x))
> plot(density(x, bw = "ucv"))

> par(mar = c(5,4,1,1)+.1)
> plot(density(x), main = "")
> plot(density(x, bw = "ucv"), main = "")

> plot(density(x, ker = "e"))