小编Arn*_*rne的帖子

我对整个R-thing仍然很陌生.

我有以下目标; 我有一个正弦函数描述了钙粒子数随时间的变化:类似于y = a*sin(b*t)+ c

实际上,在随机事件中描述了钙的产生和去除,我想在我的函数中添加一个随机噪声项(最好是平均噪声幅度可扩展).

比如z = y + random*Amplitude

你能帮我吗？

最好

我用 R 的格子包可视化了一个矩阵，以获得 10x10 的网格。

不幸的是，我确实在微调方面遇到了一些问题。

颜色条需要一个解释性标题（垂直写在它旁边）。

这是一个示例代码，可以看看它现在的样子以及一张图片。

library(lattice)

#Build the horizontal and vertical axis information

hor=c("0.0005", "0.001", "0.005", "0.01", "0.05", "0.1", "0.5", "1", "5", "10")
ver=c("1000","2000","3000","4000","5000","6000","7000","8000","9000","10000")

nrowcol=length(ver)
cor = matrix(runif(nrowcol*nrowcol, min=0.4), nrow=nrowcol, ncol=nrowcol, dimnames = list(hor, ver))
for (i in 1:nrowcol) cor[i,i] = 1

rgb.palette <- colorRampPalette(c("blue", "yellow"), space = "rgb")
levelplot(cor, col.regions=rgb.palette(120), cuts=100, at=seq(0,1,0.01),
          xlab=expression("DAG depletion rate k"[B49] *" [ s"^"-1"*" ]"),
          ylab=expression("PKC activation rate k"[D5] *" [ l / (mol*s) ]"))

你猜我该如何解决这些小问题吗？

最好，非常感谢！

我正在尝试使用 ggplot2 进行数据绘图。

纯粹出于肤浅的原因，我想在我的图例周围画一条线，以便更好地将它与情节区分开来（即图例框周围的黑色轮廓）。我在任何论坛上都找不到这个问题的答案，但也许您有提示？

library(ggplot2)

Res = matrix(ncol = 3, nrow = 500)
Res[,1] = 1:500
Res[,2] = sin((2*pi)/100*Res[,1])
Res[,3] = cos((2*pi)/100*Res[,1])
Res = as.data.frame(Res)
colnames(Res) = c("X", "Y1", "Y2")

ggplot(Res, aes(X)) + 
  geom_line(aes(y = Y1, colour = "1"), size = 2) +
  geom_line(aes(y = Y2, colour = "2"), size = 2) +
  scale_color_discrete(name = "Y's", labels = c(" sine", " cosine")) +
  theme(legend.position=c(0.9, 0.7))