2017-02-21 4 views
1

При запуске команды data.frame в R (показано ниже - обратите внимание, что «Макро» - моя переменная, представляющая интерес для модели), я получаю вывод для моей переменной fit, se, lower и upper. Я знаю, что каждый вывод говорит мне, кроме fit.Что такое «fit» в команде data.frame в пакете эффектов в R?

> data.frame(effect(c("Macro"), model)) 
    Macro fit  se   lower  upper 
    1  C 45.30041 5.650558 34.14164 56.45918 
    2  R 33.73317 4.394917 25.05406 42.41229 

Когда я запускаю команду эффекта (который я первоначально думал, дает мне мое стандартное отклонение) я получаю те же номера, как приступе:

> effect(c("macro"), model) 
    Macro effect 
    Macro 
      C   R 
    45.30041 33.73317 

Является ли это действительно стандартное отклонение, или годен более представитель среднего? И, конечно же, всегда есть вариант, что я полностью не согласен с обеими этими потенциальными интерпретациями.

+0

Это среднее значение для параметров C и R. Это не вопрос статистики, а только на очень базовом уровне. – Carl

+0

Я голосую, чтобы закрыть этот вопрос не по теме, потому что это больше вопрос интерпретации выходных данных, чем вопрос статистики. – Carl

+0

Это неоднозначно, ИМО, @ Карл. Мы обычно разрешаем вопросы, связанные с интерпретацией, даже несмотря на то, что они несколько специфичны для программного обеспечения. OTOH, это кажется менее похожим на вопросы такого типа, которые мы оставляем открытым. – gung

ответ

2

fit представляет ваши установленные или прогнозируемые значения с учетом вашей модели регрессии. В вашем случае с категорических предсказателей, это значит:

library(effects) ## to access the effect() function 

m1 <- lm(weight ~ group, data = PlantGrowth) 

data.frame(effect(c("group"), m1)) 
    group fit  se lower upper 
1 ctrl 5.032 0.1971284 4.627526 5.436474 
2 trt1 4.661 0.1971284 4.256526 5.065474 
3 trt2 5.526 0.1971284 5.121526 5.930474 

# CALCULATE MEANS 
aggregate(weight ~ group, data = PlantGrowth, mean) 
    group weight 
1 ctrl 5.032 
2 trt1 4.661 
3 trt2 5.526 

Не знаю, почему вы думали, что effect() даст вам стандартное отклонение. Посмотрите на ?effect, чтобы узнать, что вы получите при использовании функции.

Смежные вопросы