Mein Modell enthält eine Antwortvariable, fünf Prädiktoren und einen Interaktionsterm für Prädiktor_1 und Prädiktor_2. Ich möchte für jede Prädiktorvariable partielle Residuendiagramme zeichnen, die ich normalerweise mit der realisierecrPlots Funktion aus dem Paketcar. Leider beschwert sich die Funktion, dass sie mit Modellen, die Interaktionsterme enthalten, nicht funktioniert.

Gibt es eine andere Möglichkeit, das zu tun, was ich will?

EDIT: Ich habe ein kleines Beispiel erstellt, das das Problem veranschaulicht

require(car)
R <-  c(0.53,0.60,0.64,0.52,0.75,0.66,0.71,0.49,0.52,0.59)
P1 <- c(3.1,1.8,1.8,1.8,1.8,3.2,3.2,2.8,3.1,3.3)
P2 <- c(2.1,0.8,0.3,0.5,0.4,1.3,0.5,1.2,1.6,2.1)

lm.fit1 <- lm(R ~ P1 + P2)
summary(lm.fit1)
crPlots(lm.fit1) # works fine

lm.fit2 <- lm(R ~ P1*P2)
summary(lm.fit2)
crPlots(lm.fit2) # not available