3

オブジェクトの標準誤差を含む最小二乗平均は、次の関数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 )

コメントやヒントは大歓迎です。ありがとう

4

2 に答える 2

5

lme と nlme は最尤法または制限付き最尤法 (後者がデフォルト) に適合するため、結果はこれらの方法のいずれかに基づいています。

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

固定効果に切片を当てはめたので、MachineA の結果ではなく、固定効果の結果に切片項が含まれます。MachineB と MachineC の結果は切片と対比されるため、MachineB と MachineC の平均を取得するには、それぞれの値を切片の平均に追加します。しかし、標準誤差はあなたが望むものではありません。

目的の情報を取得するには、固定効果に切片項が含まれないようにモデルを適合させます (固定効果-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
于 2011-12-02T07:37:20.340 に答える
3

ダグラス・ベイツを引用するには

http://markmail.org/message/dqpk6ftztpbzgekm

「ほとんどのユーザーにとって、lsmeans の定義は「lsmeans ステートメントを使用するときに SAS から取得する数値」であると強く疑っています。そのような数値を取得するための私の提案は、SAS ライセンスを購入し、SAS を使用して適合させることです。あなたのモデル。」

于 2011-12-02T07:10:27.677 に答える