Question

在下面的R代码中，我尝试用因子allowed.clean（有30个级别）为变量zip_cpt创建30个直方图。

对于这些直方图中的每一个，我还想添加平均值和样本大小 - 需要为因子zip_cpt的每个级别计算它们。所以我使用panel.text来做到这一点。

运行此代码之后，我在每个直方图中都有错误消息，其中显示“使用数据包21时出错......”x“缺失，有......”（我无法读取整个错误消息，因为它们不要全部出现）。我猜对象x有问题。是因为mean（x）和length（x）实际上并不适用于因子zip_cpt的每个级别的数据？

我感谢任何帮助！

histogram(~allowed.clean|zip_cpt,data=cpt.IC_CAB1,  
      type='density',
      nint=100,
      breaks=NULL,
      layout=c(10,3),  
      scales= list(y=list(relation="free"),
                   x=list(relation="free")),
      panel=function(x,...) {
                             mean.values <-mean(x)
                             sample.n <- length(x)
                             panel.text(lab=paste("Sample size = ",sample.n))
                             panel.text(lab=paste("Mean = ",mean.values))
                             panel.histogram(x,col="pink", ...) 
                             panel.mathdensity(dmath=dnorm, col="black",args=list(mean=mean(x, na.rm = TRUE),sd=sd(x, na.rm = TRUE)), ...)})

Answer 1

我在网上找到的讨论有助于在每个直方图上添加自定义文本（例如基本统计信息）：

https://stat.ethz.ch/pipermail/r-help/2007-March/126842.html

r - 使用panel.text向每个点阵直方图添加文本，但出现错误“对象x缺失”

1 个答案: