5

R言語を使ってネストされたロジットモデルを推定したい. そのような問題のための標準的なパッケージとして "mlogit" を使用しています. ここで、複数の段階でモデルを推定したいと思います。問題は次のとおりです。

  1. 段階: 人々は、米国に移住するかどうかを決定します。
  2. 段階: 移住を決めたすべての人は、アメリカのどの地域に行きたいかを決める必要があります (アメリカは 6 つの地域に分けられます)。
  3. 段階: 地域内で、どのような地域に住みたいかを決定します。都会と田舎

mlogit.data() を使用して既にデータを変換しました。

data <- mlogit.data(data = usa_canada_uk, choice = "migrant"))

これが「データ」の外観です。

                sex   marst numperhh_cat age_cat   famsize                               nchlt5 labour_code sample                               nchlt10
3888.no  female married          1-5     16+         1          no chiled aged 5 or younger not working   8262          no chiled aged 10 or younger
12874.no female married          1-5     16+ 2 or more at least one child aged 5 or younger   ancillary   8262 at least one child aged 10 or younger
13084.no female married          1-5     16+ 2 or more at least one child aged 5 or younger   ancillary   8262 at least one child aged 10 or younger
9359.yes female married          1-5     16+         1          no chiled aged 5 or younger     service   8262          no chiled aged 10 or younger
7569.no  female married          1-5     16+         1          no chiled aged 5 or younger     service   8262          no chiled aged 10 or younger
5778.no  female married          1-5     16+         1          no chiled aged 5 or younger not working   8262          no chiled aged 10 or younger
         perwt        labforce age migrant country_of_birth       region citypop urban work.prob.home work.prob.abroad migration.prob stay.prob  chid
3888.no      1     in labforce  26    TRUE   United Kingdom home_country      NA urban             NA               NA             NA        NA  3888
12874.no     1     in labforce  47    TRUE   United Kingdom home_country      NA rural             NA               NA             NA        NA 12874
13084.no     1     in labforce  22    TRUE   United Kingdom home_country      NA urban             NA               NA             NA        NA 13084
9359.yes     1     in labforce  28   FALSE   United Kingdom home_country      NA urban             NA               NA             NA        NA  9359
7569.no      1     in labforce  32    TRUE   United Kingdom home_country      NA urban             NA               NA             NA        NA  7569
5778.no      1 not in labforce  38    TRUE   United Kingdom home_country      NA rural             NA               NA             NA        NA  5778
         alt
3888.no   no
12874.no  no
13084.no  no
9359.yes yes
7569.no   no
5778.no   no

これが私の考えです。どのようにコーディングしたいのですが、うまくいきません：

mlog <- mlogit(migrant ~ 1  | age + numperhh_cat + sex + famsize + work.prob.home,
             nests = list(home = c("home_country"),
                         foreign = c(region_1 = c("rural", "urban"),
                                     region_2 = c("rural", "urban"),
                                     region_3 = c("rural", "urban"),
                                     region_4 = c("rural", "urban"),
                                     region_5 = c("rural", "urban"),
                                     region_6 = c("rural", "urban"))
                                     ),
           reflevel = "yes",
           weights = perwt,
           data = data)

ご覧のとおり、最初のステージ (移行しないことを決定 = "home_country") の 1 つのネストが縮退されています。

誰かが私を助けることができれば、それは素晴らしいことです.

幸運をお祈りしています、

クリス

4

0 に答える 0