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作为行/行，而您要起诉from和to作为索引号，并且在R中不存在索引0。列名，您应该将其用作字符，如果您将定义foo的行更改为以下行，则应该可以正常工作

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

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

1 个答案: