r - 在为分类预测器绘制 gam 预测时，有没有办法使用基本 R 地毯函数

Question

我在 r 中运行了一些游戏，并且正在绘制 y 轴上的预测响应和 x 上的预测变量。我使用 rug() 来指示所有带有连续预测变量的图上的数据点数量。我现在正在做相同的预测图，但使用分类（因子）预测器。rug() 函数现在不起作用。我有以下数据和代码：

BFRT_FRC_L  <-  c('0',  '1',    '2',    '3',     
'4',    '5',    '6')
fit1    <-  c(4.794081, 4.024341,   3.747784,    
3.690677,   1.991833,   5.28066,    6.82015)
se.fit1 <-  c(0.5903663,    0.5822973, 
0.5809198,  0.5825511,  0.5999183,  0.5852767,   
0.5937017)
b   <-  c(22,   136,    145,    76, 36, 18, 3)
plot(y=exp(fit1), x=(as.factor(BFRT_FRC_L)), type="l", lwd=2, 
xlab="BFRT_FRC",
ylab="Predicted Attendance", ylim = c(0,500))
points(y=exp(fit1+se.fit1), x=BFRT_FRC_p, lwd=1, lty=2)
points(y=exp(fit1-se.fit1), x=BFRT_FRC_p, lwd=1, lty=2)
rug(x = as.integer(b))

从 gam 部分图中它应该看起来像这样：

任何让地毯工作的帮助将不胜感激。

标签： rplot

Answer 1

rug需要比预期大小的（小）向量多一点。也许这就足够了？

set.seed(2)
rug(jitter(rep.int(seq_along(b), b)))