zoukankan      html  css  js  c++  java
  • ggplot2做qqplot图

    转载自http://stats.stackexchange.com/questions/12392/how-to-compare-two-datasets-with-q-q-plot-using-ggplot2

    感谢csgillespie的答案

    qqplot是可以直观反应两组数字是否属于同一分布的作图。ggplot2提供了一个qqplot的函数,但这个函数并不能对两组观测的数字进行作图。与此相对的是,R中却有原生函数qqplot来提供这个作图。

    以下是如何利用qqplot函数的方法,使用ggplot来作图。

    这是R中qqplot的原始方法:

    R> qqplot
    function (x, y, plot.it = TRUE, xlab = deparse(substitute(x)), 
        ylab = deparse(substitute(y)), ...) 
    {
        sx <- sort(x)
        sy <- sort(y)
        lenx <- length(sx)
        leny <- length(sy)
        if (leny < lenx) 
            sx <- approx(1L:lenx, sx, n = leny)$y
        if (leny > lenx) 
            sy <- approx(1L:leny, sy, n = lenx)$y
        if (plot.it) 
            plot(sx, sy, xlab = xlab, ylab = ylab, ...)
        invisible(list(x = sx, y = sy))
    }
    <environment: namespace:stats>

    这是ggplot利用同样方法进行作图的代码:

    x <- rnorm(10);y <- rnorm(20)
    
    sx <- sort(x); sy <- sort(y)
    lenx <- length(sx)
    leny <- length(sy)
    if (leny < lenx)sx <- approx(1L:lenx, sx, n = leny)$y
    if (leny > lenx)sy <- approx(1L:leny, sy, n = lenx)$y
    
    require(ggplot2)
    g = ggplot() + geom_point(aes(x=sx, y=sy))
    g
  • 相关阅读:
    博客诞生感言~
    java 字符串锁
    oracle三种表连接方式
    两张超级大表join优化
    docker安装配置gitlab详细过程
    docker安装应用
    docker安装教程-centos
    JVM参数调优
    java向word中插入Excel附件
    application.properties参数详解
  • 原文地址:https://www.cnblogs.com/yumtaoist/p/4498729.html
Copyright © 2011-2022 走看看