【发布时间】:2020-12-30 15:56:41
【问题描述】:
正如问题中所说,我想连接 ggplot2 图中的点。这些点应该按照它们绘制到图形的顺序一个接一个地连接起来,如下所示:
basepop <- matrix(c(10,30),nrow = 2)
matrixA <- matrix(c(0.4,-0.325,0.3,1.2),nrow = 2)
matrixA
z <- eigen(matrixA)
z$val
eigenmatrix <- z$vec
eigenmatrix
eigeninverse <- solve(eigenmatrix)
constmatrix <- eigeninverse %*% basepop
constmatrix
popdata = NULL
for(k in 1:100){
predator <- c(constmatrix[1,1]*z$val[1]^k*z$vec[1,1] + constmatrix[2,1]*z$val[2]^k*z$vec[1,2])
prey <- c(constmatrix[1,1]*z$val[1]^k*z$vec[2,1] + constmatrix[2,1]*z$val[2]^k*z$vec[2,2])
popdata = rbind(popdata, data.frame(k,predator,prey))
}
library(ggplot2)
ggplot(popdata, aes(predator, prey, color = k)) +
geom_point()
ggplot(popdata, aes(k, predator, color = k)) +
geom_point()
ggplot(popdata, aes(k, prey, color = k)) +
geom_point()
print(popdata)
这是我计算猎物和捕食者系统流行度的代码。
【问题讨论】:
-
最后在
geom_line()!