Получение обобщенных средних наименьших квадратов для фиксированных эффектов в nlme или lme4

Question

Средние наименьших квадратов с их стандартными ошибками для объекта aov могут быть получены с помощью функции model.tables:

родовое слово

Интересно, как получить обобщенные средние по методу наименьших квадратов с их стандартными ошибками из объектов nlme или lme4:

родовое слово

Мы будем очень признательны за любые комментарии и подсказки.Спасибо

Answer 1

lme и nlme соответствуют максимальному правдоподобию или ограниченному максимальному правдоподобию (последнее используется по умолчанию), поэтому ваши результаты будут основаны на любом из этих методов

summary(fm1Machine) предоставит вам выходные данные, которые включают средние значения и стандартные ошибки:

родовое слово

Поскольку вы установили фиксированные эффекты с помощью точки пересечения, вы получаете член точки пересечения в результате фиксированных эффектов вместо результата для MachineA.Результаты для MachineB и MachineC контрастируют с перехватом, поэтому, чтобы получить средства для MachineB и MachineC, добавьте значение каждого к среднему перехвату.Но стандартные ошибки - это не те, которые вам бы хотелось.

Чтобы получить нужную информацию, подгоните модель так, чтобы в ней не было перехватывающего члена в фиксированных эффектах (см. код -1 в конце фиксированных эффектов:

родовое слово

Это даст вам необходимые средства и стандартный вывод ошибок:

родовое слово

Answer 2

Цитата Дугласа Бейтса из

http://markmail.org/message/dqpk6ftztpbzgekm

«У меня есть сильное подозрение, что для большинства пользователей lsmeans определяется как« числа, которые я получаю от SAS, когда использую оператор lsmeans ». Мое предложение для получения таких чисел - купить лицензию SAS и использоватьSAS под ваши модели ".