rで一般化最小二乗GLS（）を使用する方法[クローズ]

クローズ。この質問はオフトピックです。現在、回答を受け付けていません。

作業者が使用できる再現可能な例を追加してください。
この時点では、モデル内のグループ内の不均一分散構造についてはまだ説明していません'。 ?glsを確認し、weights引数を確認してください。それ以上のものが必要な場合は、@ gungが指摘したように実行可能な例を提供してください。
@Stefanご提案ありがとうございます。 weights引数を調べたところ、varFuncなどを指定する必要があることがわかりました。しかし、私の能力が限られているため、私の場合、具体的にそれを行う方法を理解することはできません。分散をより一定にしたいだけです。また、元のデータセットを追加しましたが、今のところ、適切な再現可能な例である可能性がありますか？

回答

この時点では、モデルのグループ内の不均一分散構造についてはまだ説明していません。?glsの例に示すように、weights=varPower()を試してください。。これにより、ケースの不均一分散が解消されます。

比較：

m1 <- gls(salary ~ age*sex) plot(m1)

m2 <- gls(salary ~ age*sex, weights=varPower()) plot(m2)

また、SとSの混合効果モデルの第5.2.1章（208ページ）を参照してください。 -さらに、Pinheiro and Bates 2000によると、nlmeの分散関数に関する多くの情報があります。この回答も役立つ場合があります：分散が等しくない回帰モデリング。