stages<-c("seedling", "vegetative", "flowering")
A<-matrix(c(
0, 0, 6,
0.4, 0.6, 0.02,
0.01,0.2, 0.25
), nrow=3, byrow=TRUE,
dimnames=list(stages,stages)
)
n<-c(25,25,25)
p<-pop.projection(A,n, 15)
op <- par(ask = TRUE )
stage.vector.plot(p$stage.vectors)
stage.vector.plot(p$stage.vectors, prop=FALSE, legend=c(1,300))
## log scale
stage.vector.plot(p$stage.vectors, prop=FALSE, log='y',
legend="bottomright", ylab="Number in stage class (log scaled)")
par(op)
Run the code above in your browser using DataLab