1

在下面的 R 代码中，我尝试allowed.clean按因子zip_cpt（有 30 个级别）为变量创建 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)), ...)})
4

1 回答 1

0

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

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

于 2013-07-10T17:27:25.670 回答