小编Suk*_*ngh的帖子

目前,绘图图例的运行范围为0.5到-0.25,我想将此值增加到1到-1,每0.25个间隔后中断.我怎样才能做到这一点.

如果我使用scale_fill_continous(),它将覆盖我之前填充的z值.所以,解决方案是

+ scale_fill_continuous(场所= SEQ(-1,1,通过= 0.25),低= '白色',高= '钢青')

但是有一个错误说Error in unit(tic_pos.c, "mm") : 'x' and 'units' must have length > 0 我尝试了另一个例子,我提到范围从0到1,但它仍然从0.555开始,具体取决于最大值.和分钟.z的值

应增加图例的尺寸,使颜色也应与之对应,即低为-1,高为+1.

整个想法是,值大约-0.25(中间)的瓷砖不应该是白色的,因为它们不是最低的,最低的是-1,-0.25的颜色应该是浅蓝色并且无关紧要,如果图中显示的值低于该值,则提供自定义限制.

示例图:

用于生成它的代码

pp <- function (n,r=4) {
x <- seq(-r*pi, r*pi, len=n)
df <- expand.grid(x=x, y=x)
df$r <- sqrt(df$x^2 + df$y^2)
df$z <- cos(df$r^2)*exp(-df$r/6)
df
}
ggplot(pp(20),aes(x=x,y=y))+geom_tile(aes(fill=z))

有没有办法使用命令行下载研究论文。我无法从家里访问，但我可以登录我部门的服务器，该服务器位于 Intranet 中。所以，理想情况下我应该能够做到这一点，但我不知道如何获得 pdf 的直接链接，然后获取它，或者我们可以在远程服务器中卷曲或lynx纸质网页并进行模式匹配或手动搜索链接。

可能有人已经做到了。论文示例：http://www.nature.com/nprot/journal/v7/n3/full/nprot.2012.016.html

谢谢

考虑这两个df示例

df1=data.frame(names=c('a','b','c'),value=1:3)
df2=data.frame(names=c('c','a','b'),value=1:3)

以便

> df1
  names value
1     a     1
2     b     2
3     c     3
> df2
  names value
1     c     1
2     a     2
3     b     3

现在,我想df1按照与names列相同的顺序进行排序df2,以获得

names value
c     3
a     1
b     2

我怎样才能做到这一点？

考虑测试数据集,

dat=data.frame(name=c('A','A','B','C','C','C'),val=c(1,1,2,2,3,2))

name val
A   1
A   1
B   2
C   2
C   3
C   2

获得此输出的最有效方法是什么

name val
A   1
A-1   1
B   2
C   2
C-1   3
C-2   2

因此,只需使用自定义标识符标记重复项.我可以考虑使用公共标识符来标记它们paste(dat[which(duplicated(dat$name)),1],"-1",sep=''),但是这只会在所有重复的标记之前加上"-1".我想如果该项目第3次出现,请用"-2"标记,依此类推.

干杯

我被困在一个小东西里.我在R中有一个像这样的数据框

chrom exonCount
chr1         3
chr1         4
chr1         5
chr1         5
chr1         9
chr1        10
chr2         7
chr2        11
chr2        13
chr3         7
chr4         7

我只想输出为

chr1        36
chr2        31
chr3         7
chr4         7

我假设聚合函数可以做到这一点,但我在使用中丢失了.

谢谢

我想head <(cat file.txt)在R中使用系统运行,但转义<是一个问题.

system(paste("head <(cat file.txt)"))
sh: -c: line 0: syntax error near unexpected token `('
sh: -c: line 0: `head <(cat file.txt)'

I've tried escaping it, but its not working

system(paste("head /<(cat file.txt)"))
head: cannot open `/<(cat file.txt)' for reading: No such file or directory

有人可以建议替代方案.

干杯