在使用for循环生成绘图标题时包含零

Question

我试图绘制频率与不同状态之间的转换次数（0-3）。但是，在包含状态“零”的那些地方，我得到空图。将状态等级更改为1-4可以解决空图的问题，但因此，图标题不正确。我看到this线程，其中声明循环可以从零开始。是否有任何解决方法来避免这些空白的情节？

这是代码：

n_IDs <- 150
TIME <- 0:100
set.seed(1)
(markov <- data.frame(ID=rep(1:n_IDs,each=length(TIME)),
                      TIME=rep(TIME,n_IDs),
                      STATE=sample(0:3,n_IDs*length(TIME),replace=TRUE)))

transitions <-  data.frame(ID=head(markov$ID,-1),
TIME=tail(markov$TIME,-1),FROM=head(markov$STATE,-1),TO=tail(markov$STATE,-1))

transition_table <- with(transitions,table(FROM,TO,ID))
min_n_transitions <- min(unlist(transition_table))
max_n_transitions <- max(unlist(transition_table))
max_freq <- 50  

par(mfrow=rep(4,2),mai=c(.3,.3,.3,.1))
for ( from in 0:3 ) {
  for ( to in 0:3 ) {
    foo <- transition_table[from,to,]
    hist(foo,freq=TRUE,
         breaks=seq(min_n_transitions-.5,max_n_transitions+0.5),
         xlim=c(min_n_transitions,max_n_transitions),
         ylim=c(0,max_freq),xlab="",ylab="",
         main=paste("From",from,"to",to),
         las=1,col="lightgray")
    mtext('Number of transitions', side = 1, outer = TRUE, line = 2)
    mtext('Frequency', side = 2, outer = TRUE, line = 1.5)

  }
}

先感谢您！

Answer 1

问题是你起诉from和to作为索引号并且在R索引0中不存在，要使用from和to作为行/列名称你应该将它们用作字符，如果你改变你在哪里你应该工作将foo定义为以下内容

foo <- transition_table[as.character(from), as.character(to), ]