从 MarkovChain 对象的图中替换 x 轴标签答案

【问题标题】：Replace x-axis label from the plot of a MarkovChain object从 MarkovChain 对象的图中替换 x 轴标签
【发布时间】：2019-04-09 14:56:07
【问题描述】：

这是一些生成 0 阶马尔可夫链图的代码。我想用显示今年前六个月的 45 度旋转标签替换绘图的 x 轴标签（c、d、h、i、o、p）。但是，在plot 调用中使用xaxt="n" 似乎不起作用。此代码只是覆盖现有标签，而不是替换它们。如何用我想要的标签替换标签？

library(clickstream)
clickstreams <- c("User1,h,c,c,p,c,h,c,p,p,c,p,p,o",
                   "User2,i,c,i,c,c,c,d",
                   "User3,h,i,c,i,c,p,c,c,p,c,c,i,d",
                   "User4,c,c,p,c,d",
                   "User5,h,c,c,p,p,c,p,p,p,i,p,o",
                   "User6,i,h,c,c,p,p,c,p,c,d")
csf <- tempfile()
writeLines(clickstreams, csf)
cls <- readClickstreams(csf, header = TRUE)
mc <- fitMarkovChain(cls, order=0)
plot(mc, xaxt="n")
text(x=1:6, y=par()$usr[3], labels = month.name[1:6], srt=45, adj = c(1.1,1.1), xpd = TRUE, cex=.9)

【问题讨论】：

标签： r plot clickstream

【解决方案1】：

str(mc) 揭示了 S4 对象的结构。将 xy 坐标而不是整个 "MarkovChain" 对象传递到 plot 会恢复 xaxt 选项功能。使用with 使这最方便。

此外，这允许我们使用实现的 S4 plot 方法来签名 'MarkovChain' 的 clickstream 包。

library(clickstream)
with(mc@transitions[[1]], plot(states, probability, xaxt="n"))
text(x=1:6, y=par()$usr[3], labels=month.name[1:6], srt=45, adj=rep(1.1, 2), xpd=TRUE, cex=.9)

结果

数据

library(clickstream)
mc <- new("MarkovChain", states = c("h", "c", "p", "o", "i", "d"), 
    order = 0, transitions = list(structure(list(states = structure(1:6, .Label = c("c", 
    "d", "h", "i", "o", "p"), class = "factor"), frequency = c(25L, 
    4L, 5L, 7L, 2L, 17L), probability = c(0.416666666666667, 
    0.0666666666666667, 0.0833333333333333, 0.116666666666667, 
    0.0333333333333333, 0.283333333333333)), class = "data.frame", row.names = c(NA, 
    -6L))), lambda = 0, logLikelihood = -88.4241188515082, observations = 60, 
    start = structure(c(c = 0.166666666666667, h = 0.5, i = 0.333333333333333
    ), class = "table", .Dim = 3L, .Dimnames = structure(list(
        c("c", "h", "i")), .Names = "")), end = structure(c(d = 0.666666666666667, 
    o = 0.333333333333333), class = "table", .Dim = 2L, .Dimnames = structure(list(
        c("d", "o")), .Names = "")), transientStates = c("c", 
    "h", "i", "p"), absorbingStates = c("d", "o"), absorbingProbabilities = structure(list(
        state = structure(1:4, .Label = c("c", "h", "i", "p"), class = "factor"), 
        d = c(0.792201957232916, 0.864224245314581, 0.903842865008182, 
        0.517925028202586), o = c(0.207798042767084, 0.13577575468542, 
        0.0961571349918185, 0.482074971797414)), class = "data.frame", row.names = c(c = 1L, 
    h = 3L, i = 4L, p = 6L)))

【讨论】：

【解决方案2】：

使用plot.default 而不是plot 会暴露图形plot 函数，这允许使用xaxt 函数参数来删除现有的x 轴标签。由于 plot.default 不接受 MarkovChain 对象，因此需要从对象中提取绘图的 x 和 y 值。

plot.default(
   x=mc@transitions[[1]]$states, 
   y=mc@transitions[[1]]$probability,
   xaxt="n", 
   ann=FALSE, 
   pch="-", 
   cex=3
)

text(
    x = 1:6, 
    y = par()$usr[3], 
    labels = month.name[1:6], 
    srt = 45, 
    adj = c(1.1,1.1), 
    xpd = TRUE, 
    cex = .9
)

【讨论】：

with(mc@transitions[[1]], plot(states, probability, xaxt="n")) 也可以。
谢谢，您的回答比我的要好，因为它使用点击流绘图函数的原始样式。如果您想提交它作为答案，我会接受。
大功告成，乐于助人！