我有下面的脚本为5个粒子做布朗运动,这将不会输出 colored颜色 组.使用BASE R可以很好地完成这项工作,但我正在try 处理ggploy.
5个单独粒子执行10个布朗步骤的简单脚本(这些步骤的正确性无关紧要),并在Tibble col3中附加了一个群因子:
#library(tidyverse) (if needed)
######### Construct string num_points starting from zero and moving -1,0,+1 on each step
num_points <- 10
series <- function (num_points) { cumsum(c(0,replicate(num_points-1,(sample(c(-1,0,+1),1)))) ) }
######### Construct dataset of num_sets groups of random walks. Group ID is in col 'z'
num_particles <- 5
data_set <- function (num_points) {
q <- NULL
for (i in 1:num_particles) {q_tem <- tibble(x = series(num_points),y = series(num_points),group=as.factor(i))
q <- bind_rows(q,q_tem)
}
return(q)
}
z <- data_set(num_points)
group.colors <- c("#333BFF","#CC6600","#9633FF","#E2FF33","#E3DB71")
ggplot(z, aes(x=x,y=y, group = group, fill = group)) + # Properly adding colors
geom_point() +
scale_color_manual(values = group.colors)
这将打印出所有50个点,全部为黑色,而不是指定的 colored颜色 .我想定义我的调色板.
我试过很多例子,上面的代码基本上都是从帮助站点复制过来的.但预期的 colored颜色 并没有出现.最后,我想做更多的颗粒和更大的样本,用叠印打印,我已经想好了怎么做.