困难的弗里德曼事后测试

Question

我正在尝试在R中完成Friedman事后测试；到目前为止，我使用的每个软件包（PCMR，PMCMRplus，NSM3）都遇到了麻烦。我的回复数据是顺序的，我不确定是否是问题所在。有什么包装或代码建议可能有用吗？

我已经完成了弗里德曼测试，分析了24周内动物的每周行为。行为类别包括“防御”，“停止”和“逃脱”行为。

数据是0.014ppmgeneralFriedman

head(`0.014ppmgeneralFriedman`)
  crayfishid   group time behavior

Nbehavior <- data_frame(behavior = c("def","stop","escape"))

mapping <- c("def"=1,"stop"=2,"escape"=3)

orderedbehavior <- mapping[`0.014ppmgeneralFriedman`$behavior]

friedman.test(orderedbehavior~time|crayfishid, data=`0.014ppmgeneralFriedman`)

    Friedman rank sum test

data:  orderedbehavior and time and crayfishid
Friedman chi-squared = 44.152, df = 24, p-value = 0.007325

但是我无法成功完成事后！