小编Lin*_*ang的帖子

首先,这是mydata(121315*4)的一部分:

     LONGITUDE LATITUDE NUM_PICKUPS TOTAL_REVENUE
 1   121.6177  38.9124          21         337.0
 2   121.8069  39.0210          16         454.7
 3   121.5723  38.9645          38         696.9
 4   121.6423  38.9258         622       13609.7
 5   121.5647  38.9129         116        2016.7
 6   121.6429  38.8846         120        2417.3
 7   121.5852  38.9279         117        1975.0
 8   121.6616  38.9189          94        1712.4
 9   121.5812  38.9828          50         981.6
 10  121.6411  38.9255         225        4696.2

看到这一点,第一列和第二列是经度和纬度.

mydata[1,3]=21意思是在地方(121.6177, 38.9124),有21个拾音器.

然后,我用NUM_PICKUPSdesc 求助mydata :

LONGITUDE   LATITUDE    NUM_PICKUPS TOTAL_REVENUE
121.6019    39.0181        14243    514716
121.5382    38.9609        13244    443754.7
121.5381    38.9609 …

由于数据的隐私性,我使用mtcar数据集ggplot2来解释我的问题.

有四个图:

g1 <- ggplot(mtcars,aes(mpg,wt)) + geom_point()
g2 <- ggplot(mtcars,aes(mpg,disp)) + geom_point()
g3 <- ggplot(mtcars,aes(mpg,drat)) + geom_point()
g4 <- ggplot(mtcars,aes(mpg,qsec)) + geom_point()

我想把这四个图放在一个图中,所以我grid.arrange()在包中使用函数grid.Extra:

grid.arrange(g1,g2,g3,g4,ncol=2)

现在,我想在这些图中的每个图下面添加标题,如下图所示(我在Word中修改它,所以它不漂亮) 在询问之前,我已经搜索过SO,并且我知道如何在一个绘图下添加标题,例如,使用grid.text()或这三种方法 在ggplot2生成的绘图下方显示文本,或者element_text(vjust=-10)我不能将其应用于一个图中的四个绘图.同时,我在基础图中得到了一些结果如何为R中的每个图添加标题？或者使用par(mfrow)编译的数字面板的公共主标题,qustion是我想要进行ggplot2并且标题在每个情节下面,我该如何实现它？谢谢!

我使用rplcon()包中的函数生成一些随机变量poweRlaw

data <- rplcon(1000,10,2)

现在，我想知道哪些已知分布最适合数据。对数范数？经验？伽玛？幂律？指数截止的幂律？

所以我fitdist()在包中使用函数fitdistrplus：

fit.lnormdl <- fitdist(data,"lnorm")
fit.gammadl <- fitdist(data, "gamma", lower = c(0, 0))
fit.expdl <- fitdist(data,"exp")

由于幂律分布和具有指数截止的幂律不是根据CRAN Task View: Probability Distributions的基本概率函数，所以我根据示例 4 编写了幂律的 d,p,q 函数?fitdist

dplcon <- function (x, xmin, alpha, log = FALSE) 
{
    if (log) {
        pdf = log(alpha - 1) - log(xmin) - alpha * (log(x/xmin))
        pdf[x < xmin] = -Inf
    }
    else {
        pdf = (alpha - 1)/xmin * (x/xmin)^(-alpha)
        pdf[x < …