小编Ben*_*son的帖子

我想在x轴上放一个度数符号,但结果有一个额外的空间,我似乎无法摆脱它.文本应为"温度(*C)",而不是"温度(*C)".我尝试了两种不同的解决方案,但似乎无法摆脱空间.

ggdat<-data.frame(x=rnorm(100),y=rnorm(100))

#neither of these approaches work
xlab <- expression(paste('Temperature (',~degree,'C)',sep=''))
xlab <- expression('Temperature ('*~degree*C*')')

ggplot(data=ggdat,aes(x=x,y=y)) +
  geom_point() +
  labs(x=xlab)

任何帮助表示赞赏!

本

我的数据包含多个人的几个数据点,每个人都来自特定的研究站点.我想绘制所有点,为每个人绘制95%的椭圆,但然后按研究地点为椭圆绘制颜色.不幸的是,当我指定按站点着色时,会为聚合组绘制椭圆.

数据如下所示:

dat1 <- data.frame(X=rnorm(21),Y=rnorm(21),indiv_id=rep(c(1,2,3),7),group_id=rep(1,21))
dat2 <- data.frame(X=rnorm(21,5),Y=rnorm(21,5),indiv_id=rep(c(4,5,6),7),group_id=rep(2,21))
dat3 <- data.frame(X=rnorm(21,10),Y=rnorm(21,10),indiv_id=rep(c(7,8,9),7),group_id=rep(3,21))
ggdat <- rbind(dat1,dat2,dat3)
ggdat$indiv_id <- as.factor(ggdat$indiv_id)
ggdat$group_id <- as.factor(ggdat$group_id)

如果我按个人绘制省略号,我可以分别看到所有的省略号:

ggplot(ggdat) +
  geom_point(aes(x=X, y=Y,color=indiv_id),size=1) + # 
  stat_ellipse(aes(x=X, y=Y,color=indiv_id),type = "norm")

但如果我按小组绘制,它每组只有一个椭圆:

ggplot(ggdat) +
  geom_point(aes(x=X, y=Y,color=indiv_id),size=1) + # 
  stat_ellipse(aes(x=X, y=Y,color=group_id),type = "norm") + #, linetype = 2
  theme(legend.position='none')

如何绘制所有9个椭圆,但按组颜色进行着色？谢谢您的帮助!

我正在尝试使用 sf 包和管道/tidyverse 工作流来根据另一列中定义的组生成边界框。我认为它应该像下面那样工作，但似乎 st_bbox 不尊重群体。

我希望收到三个多边形记录，它们代表 a、b 和 c 的 proints 周围的边界框，但我收到三个多边形记录，它们代表所有点的边界框。

library(dplyr)
library(sf)

a <- data.frame(group=rep('a',100), lon=rnorm(100,11,.2), lat=rnorm(100,53,.2))
b <- data.frame(group=rep('b',100), lon=rnorm(100,11.5,.2), lat=rnorm(100,53.5,.2))
c <- data.frame(group=rep('c',100), lon=rnorm(100,12,.2), lat=rnorm(100,54,.2))
dat <- rbind(a,b,c)

pts <- dat %>% st_as_sf(coords=c('lon','lat'),crs=4326) 

pts %>%
  group_by(group) %>%
  summarize(geometry = st_as_sfc(st_bbox(geometry)))

这将返回：

Simple feature collection with 3 features and 1 field
geometry type:  POLYGON
dimension:      XY
bbox:           xmin: 10.34313 ymin: 52.43993 xmax: 12.54254 ymax: 54.54012
epsg (SRID):    4326
proj4string:    +proj=longlat +datum=WGS84 +no_defs
# A tibble: 3 …

我在Ubuntu 12.04上使用pgAdmin III 1.18.1.每当我重启机器时,pgAdmin都不记得我的服务器连接.如果我进入服务器连接然后关闭pgadmin3并重新打开,则保留服务器连接,只有当我重新启动计算机时才会丢失连接.

在我创建服务器连接并关闭pgAdmin后,我看到.pgadmin3文件创建得很好:

benc@Ubuntu-T7400:~$ ls -al | grep .pgadmin3
-rw-rw-r--   1 benc benc  2088 Mar 28 11:57 .pgadmin3

然后重新启动pgadmin3将显示我的服务器连接.但是,如果我重启我的机器,我的.pgadmin3文件设置为0字节.

重新启动Ubtuntu之后:

benc@Ubuntu-T7400:~$ ls -al | grep .pgadmin3
-rw-rw-r--   1 benc benc     0 Mar 28 12:10 .pgadmin3

知道这里可能会发生什么吗？

我认为dplyr太棒了.但是我最近更新了这个软件包,似乎总结不再为我分组.类似于以下用于在升级之前按组汇总的代码:

iris_tdt <- tbl_dt(iris)
iris_tdt %.% group_by(Species) %.% summarise(avg_petal_width = mean(Petal.Width))

  avg_petal_width
1        1.199333

这用于输出具有Species和avg_petal_width的表.现在,物种列似乎被删除,avg_petal_width聚合为单个值.group_by似乎正在工作,所以我猜这是一个总结的问题.

grp <- group_by(iris_tdt,Species)
groups(grp)

[[1]]
Species

甚至Vignette的例子都没有正常工作.

hflights_df <- tbl_df(hflights)
planes <- group_by(hflights_df, TailNum)
delay <- summarise(planes,
  dist = mean(Distance, na.rm = TRUE),
  delay = mean(ArrDelay, na.rm = TRUE))

delay
      dist    delay
1 787.7832 7.094334

任何建议将不胜感激.

packageDescription("dplyr")$Version #--> 0.1.2
R.version.string #--> "R version 3.0.2 (2013-09-25)"