r - 如何在 R 中编写 summarySE 函数

Question

我是 R 的新用户，想用 std 制作折线图使用 ggplot .
为此，我使用 summarySE功能。

vrn1_summary<-summarySE(data,measurevar = "vrn1",groupvars = c("genotype","treatment","time"))

由于我无法创建自己的函数，因此我从 R Cookbook 中复制如下:

summarySE <- function(data=NULL, measurevar, groupvars=NULL, na.rm=FALSE,
                      conf.interval=.95, .drop=TRUE) {
    library(plyr)

    # New version of length which can handle NA's: if na.rm==T, don't count them
    length2 <- function (x, na.rm=FALSE) {
        if (na.rm) sum(!is.na(x))
        else       length(x)
    }

    # This does the summary. For each group's data frame, return a vector with
    # N, mean, and sd
    datac <- ddply(data, groupvars, .drop=.drop,
      .fun = function(xx, col) {
        c(N    = length2(xx[[col]], na.rm=na.rm),
          mean = mean   (xx[[col]], na.rm=na.rm),
          sd   = sd     (xx[[col]], na.rm=na.rm)
        )
      },
      measurevar
    )

    # Rename the "mean" column    
    datac <- rename(datac, c("mean" = measurevar))

    datac$se <- datac$sd / sqrt(datac$N)  # Calculate standard error of the mean

    # Confidence interval multiplier for standard error
    # Calculate t-statistic for confidence interval: 
    # e.g., if conf.interval is .95, use .975 (above/below), and use df=N-1
    ciMult <- qt(conf.interval/2 + .5, datac$N-1)
    datac$ci <- datac$se * ciMult

    return(datac)
}

当我调用这个函数时，它给出 sd (标准偏差)，se (标准错误)和 ci (置信区间)，但显示 NA为 mean ，在我的数据中是“vrn1”，也显示警告。
据我了解，这可能是 summarySE 中的问题。功能，但我无法找出在哪里/为什么。

谁能帮我解决这个问题。或者指导我如何写一个简单的summarySE功能？

Answer 1

请检查您的数据

str(data)

我猜你的measurevar vrn1 是 Chrome .