0

我有一个包含三个分组变量的数据集:条件、子和延迟。这是我的数据的简化版本(真实数据要长得多)

sub condition delay later_value choiceRT later_choice primeRT        cue
 10      SIZE    10          27     1832            1     888      CHILD
 10      PAST     5          11      298            0    1635      PANTS
 10      SIZE    21          13      456            0     949      CANDY
 11      SIZE   120          22      526            1    7963        BOY
 11    FUTURE   120          27      561            1    4389   CHILDREN
 11      PAST     5          13      561            1    2586     SPRING

我有一套复杂的程序来应用这些数据(细节并不重要)我编写了以下函数,当被三个分组变量分割时,它完成了我想要的。它返回 3 个我感兴趣的变量(indiff、p_intercept 和 p_lv)

 getIndiffs <- function(currdelay){
      if (mean(currdelay$later_choice) == 1) {
        indiff = 10.5
        p_intercept = "laters"
        p_lv = "laters"
      }

      else if (mean(currdelay$later_choice) == 0) {
        indiff = 30.5

        # no p-val here, code that this was not calculated
        p_intercept = "nows"
        p_lv = "nows"
      }

      else {
        F <- factor(currdelay$later_choice)

        fit <- glm(F~later_value,data=currdelay,family=binomial())
        indiff <- -coef(fit)[1]/coef(fit)[2]

        if (indiff < 10) indiff = 10.5
        else if (indiff > 30) indiff = 30.5

        p_intercept = round(summary(fit)$coef[, "Pr(>|z|)"][1],3)
        p_lv = round(summary(fit)$coef[, "Pr(>|z|)"][2], 3)
        c(indiff,p_intercept,p_lv)
      }

我正在尝试使用 ddply 将其应用于每个 3 个分组变量的数据子集:

ddply(数据,。(子,条件,延迟),getIndiffs)

但是,当我运行它时,我得到了错误

list_to_dataframe(res, attr(.data, "split_labels")) 中的错误:结果的长度不相等

奇怪的是,当我只使用 1 个分组变量但用 2+ 抛出错误时,这很好用

此外,当我自己“模拟”将数据集拆分为仅包含由 3 个分组变量拆分的子集的数据帧时,我的函数工作得很好。(注意:我尝试了返回 3 个变量甚至只返回 1 个变量的不同方法,但它也不起作用)

基本上,我想知道的是如何使用 plyr 来使用函数返回多个变量。

对我的问题的任何其他根本不同的解决方案也是受欢迎的。

4

1 回答 1

7

当我的函数应用于我的一个片段返回一个空数据框时,通常会发生该错误。在任何情况下,调试情况的一种简单方法是使用dlply而不是 ddply,并检查输出;例如

x <- dlply(data,.(sub,condition,delay),getIndiffs)
sapply(x,ncol)

检查它们是否都具有相同的列数。如果没有,请更加标准化您的功能。

看起来您的函数getIndiffs设计为在单行上运行,而不是在整个数据帧上运行。 处理由与该组匹配的观察子集组成的d*ply(x,vars,fn)整个fn()数据框。嗯,同样,该函数可以在三个不同的位置返回——在每个条件子句的末尾。我认为你的意思是放在c(indiff,p_intercept,p_lv)最后一个}(并用另一个结束你的函数})。

于 2011-09-08T05:30:24.593 回答