GuangchuangYu

f <- function(x) 1/(x^2-1)
x <- seq(-3,3, by=0.001)
y <- f(x)
d <- data.frame(x=x,y=y)

p <- ggplot()
p <- p+geom_rect(fill = "white",color="black",size=3,
                 aes(NULL, NULL,xmin=-3, xmax=3,
                     ymin=-3,ymax=3, alpha=0.1))

GuangchuangYu

require(ggplot2)

t <- seq(0,2*pi, by=0.1)
x <- 16*sin(t)^3
y <- 13*cos(t) - 5*cos(2*t) - 2*cos(3*t) - cos(4*t)

d <- data.frame(x=x,y=y, f=0)

a <- x
b <- y

GuangchuangYu

fibonacci <- function(n) {
    u <- (1+sqrt(5))/2
    (u^n - (1-u)^n) / sqrt(5)
}

GuangchuangYu

p <- ggplot(mtcars, aes(factor(cyl),mpg,
                        fill=factor(cyl),
                        colour=factor(cyl)))
 
p1 <- p+geom_violin(alpha=0.3, width=0.5) +
    geom_boxplot(width=0.2, outlier.colour=NA)
 
p2 <- p+geom_violin(alpha=0.3, width=0.5) +
    geom_dotplot(binaxis='y', stackdir='center', dotsize=0.5)