R 多变量模型的dummy.coef（）

我想知道

dummy.coef（）

中的一些东西，它将方差分析模型中的估计参数（对比度）转换为原始参数。它只适用于单变量模型，但使其适用于多变量模型所需的更改似乎很小。在

dummy.coef.lm（）中

：

---

• 第52行

  coef这只解决了第二个问题（关于代码第20-22行的
  if（）x else y
  语句的问题）
  
  首先，尝试将这两个块剪切并粘贴到R会话中：
  
  test <- TRUE
  
  # Block 1 -- Doesn't work
  if(test) 
      cat("test is TRUE\n")
      else
      cat("test is FALSE\n")
  
  # Block 2 -- Works
  {
  if(test) 
      cat("test is TRUE\n")
      else
      cat("test is FALSE\n")
  }
  
  test我们有一个“类似”的补丁，我现在正在看它——它看起来不错，尽管一次改变的东西比我预期的要多
  
  当我在谷歌上搜索dummy.coef的其他问题时，我得到了这个旧的SO线程。
  让我回答一下（作为团队成员），如果代码（和文档）变化不大，我们对多元模型的扩展感兴趣。如果您有兴趣帮忙，请回复。
  我会接受这个答案，但我还是很想知道第一个问题的答案……谢谢。这是你的第一个问题，我想不出你不应该这么做的理由，而且因为R是开源的，谢天谢地，你很容易就能做到。您可能还需要对print.dummy_coef进行一些编辑，以获得完整的功能，但这也不太难。感谢您的关注！我的建议是，快速登记。请注意，这会将
  dummy.coef.lm（）
  返回的列表的组件更改为始终为矩阵。此外，在多变量情况下打印对象的方法可能需要更多考虑。是否可以/应该将其扩展到
  dummy.coef.aovlist（）
  ？如果它更适合您的工作流程，我可以用补丁打开正式的错误报告。在建议的
  dummy.coef.lm（）
  版本中，唯一修改的行是71、78-80，而建议的
  print.dummy\u coef（）版本
  与原始版本的差异更大。我的dummy.coef.R源文件的工作版本与在线版本相比已经有了很大的变化，因为我已经解决了上面提到的（和链接的）错误报告。因此，我将手动“diff”并应用我看到的更改。同样感谢您提供的示例。事实上，我认为添加官方bug报告（尽管没有补丁，请参见上文）或添加到[现有的一个][1]并添加注释是有意义的。[1] ：谢谢，我在这里打开了一个错误报告：
  for (i in vars) args[[i]] <- if (nxl[[i]] == 1) {
      rep.int(1, nl)
  } else { factor(rep.int(xl[[i]][1L], nl), levels = xl[[i]]) }
  
  test <- TRUE
  
  # Block 1 -- Doesn't work
  if(test) 
      cat("test is TRUE\n")
      else
      cat("test is FALSE\n")
  
  # Block 2 -- Works
  {
  if(test) 
      cat("test is TRUE\n")
      else
      cat("test is FALSE\n")
  }