الحصول على المربعات الصغرى المعممة يعني للآثار الثابتة في نلم أو لم 4

Question

المربعات الصغرى تعني مع أخطائها القياسية ل aov يمكن الحصول على الكائن مع model.tables وظيفة:

npk.aov <- aov(yield ~ block + N*P*K, npk)
model.tables(npk.aov, "means", se = TRUE)

أتساءل عن كيفية الحصول على المربعات الصغرى المعممة تعني أخطائها القياسية من nlme أو lme4 الكائنات:

library(nlme)
data(Machines)
fm1Machine <- lme(score ~ Machine, data = Machines, random = ~ 1 | Worker )

أي تعليق وتلميح سيكون محل تقدير كبير.شكرا

Answer 1

لم و نلم تناسب من خلال أقصى احتمال أو الحد الأقصى من احتمال الحد (وهذا الأخير هو الافتراضي) ، لذلك سوف تستند النتائج الخاصة بك على أي من تلك الأساليب

summary(fm1Machine) سوف توفر لك مع الإخراج الذي يتضمن الوسائل والأخطاء القياسية:

....irrelevant output deleted
Fixed effects: score ~ Machine 
               Value Std.Error DF  t-value p-value
(Intercept) 52.35556  2.229312 46 23.48507       0
MachineB     7.96667  1.053883 46  7.55935       0
MachineC    13.91667  1.053883 46 13.20514       0
 Correlation: 
....irrelevant output deleted

نظرا لأنك قمت بتركيب التأثيرات الثابتة مع اعتراض ، تحصل على مصطلح اعتراض في نتيجة التأثيرات الثابتة بدلا من نتيجة للآلة.نتائج ماشينب و ماشينيك هي تناقضات مع اعتراض ، وذلك للحصول على وسائل ماشينب و ماشينيك ، إضافة قيمة كل إلى متوسط اعتراض.لكن الأخطاء القياسية ليست هي الأخطاء التي تريدها.

للحصول على المعلومات التي تبحث عنها ، قم بتناسب النموذج بحيث لا يحتوي على مصطلح اعتراض في التأثيرات الثابتة (انظر -1 في نهاية الآثار الثابتة:

fm1Machine <- lme(score ~ Machine-1, data = Machines, random = ~ 1 | Worker )

سيعطيك هذا بعد ذلك الوسائل وإخراج الخطأ القياسي الذي تريده:

....irrelevant output deleted
Fixed effects: score ~ Machine - 1 
            Value Std.Error DF  t-value p-value
MachineA 52.35556  2.229312 46 23.48507       0
MachineB 60.32222  2.229312 46 27.05867       0
MachineC 66.27222  2.229312 46 29.72765       0
....irrelevant output deleted

Answer 2

الاقتباس من دوجلاس بيتس من

http://markmail.org/message/dqpk6ftztpbzgekm

"لدي شك قوي في أن تعريف lsmeans لمعظم المستخدمين هو" الأرقام التي أحصل عليها من SAS عندما أستخدم عبارة lsmeans ". اقتراحي للحصول على هذه الأرقام هو شراء ترخيص SAS واستخدامهSAS لتناسب نماذجك. "