Встроенное квадратное преобразование, дающее результаты, отличные от преобразованных значений столбца - PullRequest
Новая
       14

Встроенное квадратное преобразование, дающее результаты, отличные от преобразованных значений столбца

1 голос
/ 06 апреля 2019

При использовании встроенного квадратного преобразования (т. Е. С помощью ~ параметра ^ 2 в формуле модели) результаты будут отличаться от результатов, когда я сначала преобразовываю параметр в другой столбец и использую этот столбец в модели (т. Е. Parameter2 = параметр ^ 2и затем используя ~ Parameter2 в формуле) 

Worst.Perimeter2 = Worst.Perimeter^2

m0x_ <- glm(formula = Diagnosis ~ Worst.Perimeter2, family = binomial, 
    data = cancerdata)
summary(m0x_)

results:
Coefficients:
                   Estimate Std. Error z value Pr(>|z|)    
(Intercept)      -1.090e+01  1.132e+00  -9.634   <2e-16 ***
Worst.Perimeter2  8.974e-04  9.792e-05   9.165   <2e-16 ***



m0x1_ <- glm(formula = Diagnosis ~ Worst.Perimeter^2, family = binomial, 
    data = cancerdata)
summary(m0x1_)

Results:
Coefficients:
                 Estimate Std. Error z value Pr(>|z|)    
(Intercept)     -20.81003    2.22116  -9.369   <2e-16 ***
Worst.Perimeter   0.18931    0.02068   9.156   <2e-16 ***

Я бы ожидал, что они дадут те же результаты для коэффициентов, но, очевидно, это не так.Как мне интерпретировать эти разные результаты?

1 Ответ

2 голосов
/ 06 апреля 2019

^ - оператор формулы.Вы можете использовать его, чтобы указать формулу, которая включает взаимодействия, как в этом примере: 

# Formula `(a + b)^2` means `1 + a + b + a*b`
head(model.matrix(~ (mpg + wt)^2, data = mtcars))
#>                   (Intercept)  mpg    wt mpg:wt
#> Mazda RX4                   1 21.0 2.620 55.020
#> Mazda RX4 Wag               1 21.0 2.875 60.375
#> Datsun 710                  1 22.8 2.320 52.896
#> Hornet 4 Drive              1 21.4 3.215 68.801
#> Hornet Sportabout           1 18.7 3.440 64.328
#> Valiant                     1 18.1 3.460 62.626

Если вы хотите, чтобы ^ интерпретировался как арифметический оператор, используйте оператор I (AsIs) какв этом примере: 

head(model.matrix(~ I((mpg + wt)^2), data = mtcars))
#>                   (Intercept) I((mpg + wt)^2)
#> Mazda RX4                   1        557.9044
#> Mazda RX4 Wag               1        570.0156
#> Datsun 710                  1        631.0144
#> Hornet 4 Drive              1        605.8982
#> Hornet Sportabout           1        490.1796
#> Valiant                     1        464.8336

Создано в 2019-04-06 пакетом представлением (v0.2.1)

...