r - 在 R 中保持变量不变的预测边际/预测

Question

我最近切换到 R，我喜欢它。但我最想念的一件事是能够生成预测模型响应，将某些变量保持在预设水平(平均值、第 90 个百分位数等)。当试图辨别交互项、转换变量等的影响时，这非常方便。

我可以使用 adjust 在 Stata 中轻松完成此操作命令。我已经尝试并试图弄清楚如何在 R 中做到这一点，但是使用名为 R 的语言(也有一个统计 R)并搜索诸如“调整”之类的术语的一大陷阱是我可以似乎只能在调整后的 R 平方上找到命中。这是令人沮丧的。

所以，冒着问一个非常简单的问题的风险，有谁知道如何做到这一点？我研究了预测边际，这似乎至少是一种相关类型的方法，但它的实现通常涉及以某种方式标准化解释变量。

Answer 1

http://www.stata.com/support/faqs/stat/adjust.html指出:

In fact, adjust is really just a front end for predict

它的工作原理是为预测创建一个新的数据集。

R的 predict有争论 newdata用于改变预测假设，但可能不像 Stata 的 adjust 那样容易. transform可能也很有趣，请参见以下示例。

# create model
mtcars.lm <- lm(hp~disp*cyl,mtcars)

# default fit predictions

predict(mtcars.lm)
          Mazda RX4       Mazda RX4 Wag          Datsun 710      Hornet 4 Drive 
          128.60896           128.60896            80.25811           128.88296 
  Hornet Sportabout             Valiant          Duster 360           Merc 240D 
          208.48842           128.79069           208.48842            75.58796 
           Merc 230            Merc 280           Merc 280C          Merc 450SE 
           76.29995           128.63021           128.63021           197.85671 
         Merc 450SL         Merc 450SLC  Cadillac Fleetwood Lincoln Continental 
          197.85671           197.85671           222.63037           221.11516 
  Chrysler Imperial            Fiat 128         Honda Civic      Toyota Corolla 
          218.58981            83.79391            84.15593            84.71104 
      Toyota Corona    Dodge Challenger         AMC Javelin          Camaro Z28 
           78.79793           203.18519           201.41745           207.22575 
   Pontiac Firebird           Fiat X1-9       Porsche 914-2        Lotus Europa 
          213.53912            83.75770            78.77380            81.81483 
     Ford Pantera L        Ferrari Dino       Maserati Bora          Volvo 142E 
          207.35202           128.56702           201.03865            78.68933 


# predict assuming all cars have 8 cylinders

predict(mtcars.lm,newdata=transform(mtcars,cyl=8))
          Mazda RX4       Mazda RX4 Wag          Datsun 710      Hornet 4 Drive 
           183.2349            183.2349            176.6690            195.6091 
  Hornet Sportabout             Valiant          Duster 360           Merc 240D 
           208.4884            191.4423            208.4884            181.5556 
           Merc 230            Merc 280           Merc 280C          Merc 450SE 
           180.8106            184.1946            184.1946            197.8567 
         Merc 450SL         Merc 450SLC  Cadillac Fleetwood Lincoln Continental 
           197.8567            197.8567            222.6304            221.1152 
  Chrysler Imperial            Fiat 128         Honda Civic      Toyota Corolla 
           218.5898            172.9694            172.5906            172.0098 
      Toyota Corona    Dodge Challenger         AMC Javelin          Camaro Z28 
           178.1969            203.1852            201.4174            207.2257 
   Pontiac Firebird           Fiat X1-9       Porsche 914-2        Lotus Europa 
           213.5391            173.0073            178.2221            175.0402 
     Ford Pantera L        Ferrari Dino       Maserati Bora          Volvo 142E 
           207.3520            181.3409            201.0386            178.3105

至于搜索 R 相关主题，您可能会发现以下链接很有用: How to search for "R" materials?