ggplot2做qqplot图

摘要:
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

免责声明:文章转载自《ggplot2做qqplot图》仅用于学习参考。如对内容有疑问,请及时联系本站处理。

上篇FLINK --- 写HDFSHAProxy安装文档下篇

宿迁高防,2C2G15M,22元/月;香港BGP,2C5G5M,25元/月 雨云优惠码:MjYwNzM=

随便看看

CAD转DXF怎么转换?教你三种转换方法

2.进入到CAD版本转换的界面中后,在选择“点击选择文件”,在跳转出的“打开”界面中打开需要转换的CAD图纸。...

layui 学习笔记(四) 复杂表头前台Excel导出

merges':mergeConf,'!cols':colConf,'!rows‘:rowConf}});}@...

flutter Radio单选框

单选框,允许用户从一组中选择一个选项。...

C#Win32API编程之PostMessage

本文以C#调用Win32API函数PostMessage完成指定表单的后台鼠标和键盘模拟为例,大致解释了C#调用非托管代码和Window的消息处理机制。我们可以将PostMessage用于函数。成功与否在很大程度上取决于我们传达的信息是否真实。消息表明消息是什么。请原谅我先讲故事。我希望先解释一下PostMessage函数。这是一个异步操作,如下图所示:调用...

Java 读取ANSI文件中文乱码问题解决方式[转]

Filefile=newFile(路径);InputStreamin=newjava.io.FileInputStream(文件);BufferedReader读取器=新的BufferedReader(读取);FileInputStreamin=newFileInputStream(文件);byte[]b=新字节[3];内容如下(b);...

Java 安全之:csrf攻击总结

最近,我在维护一些旧项目。在调试期间,我发现请求被反复拒绝。我仔细查看了项目的源代码,发现存在csrftoken验证。我借此机会了解了csrf攻击,并将其总结成一篇论文。受攻击的网站无法阻止攻击。在整个过程中,攻击者无法获取受害者的登录凭据,只能“冒充”。CSRF攻击成功,因为服务器将攻击者发送的请求误认为是用户的请求。服务器通过验证请求是否携带正确的令牌来...