anesrake 錯誤：“沒有變量偏離超過 ____”

Question

我需要根據來自更廣泛人群的四個人口特征的邊際分布對樣本中的觀察結果進行加權。 我目前正在使用 package anesrake這樣做。

人口信息存儲在targets中。 這是一個包含 4 個元素的列表 - 我要根據每個受訪者屬性對樣本進行加權的一個數字向量。 每個元素的行名代表不同的類別。 我在這里創建targets ：

quota_age    <- c(0.30, 0.33, 0.37)
quota_race   <- c(0.62, 0.12, 0.17, 0.5, 0.3)
quota_gender <- c(0.52, 0.48)
quota_ed     <- c(0.41, 0.29, 0.19, 0.11)

names(quota_age)    <- c("18 to 34", "35 to 54", "55+")
names(quota_race)   <- c("White non-Hispanic", "Black non-Hispanic", "Hispanic", "Asian", "Other")
names(quota_gender) <- c("Female", "Male")
names(quota_ed)     <- c("HS or less", "Some college", "Bachelors", "Advanced")

targets <- list(quota_age, quota_race, quota_gender, quota_ed)

調查文件 ( m1b ) 是一個數據框，其中包含人口統計信息和每個受訪者的唯一 ID（ 此處鏈接到谷歌表）。 這是前幾個obs：

> head(m1b)
         ResponseId     quota_ed quota_age quota_gender         quota_race
1 R_3McITJbfcFuwc9x Some college  18 to 34       Female White non-Hispanic
2 R_2q3oeAbZgCZ5YcZ    Bachelors       55+       Female White non-Hispanic
3 R_YSVccSQ1xJ6zuDv     Advanced  35 to 54       Female White non-Hispanic
4 R_DubbKu7uJicbpQd Some college  35 to 54         Male White non-Hispanic
5 R_5zj5CNu598lCwRX    Bachelors       55+         Male              Other
6 R_21mPGFS7kHX2ELm     Advanced       55+       Female White non-Hispanic

使用anesrake package，我想構建一個名為weight的新變量，我可以在以后的分析中使用它來解釋總體和樣本邊際分布之間的差異。

但是當我像這樣調用anesrake function 時（ pctlim參數非常小，無法誇大我的觀點）：

library(anesrake)

raking <- anesrake(inputter     = targets,
                   dataframe    = m1b,
                   caseid       = m1b$ResponseId,
                   choosemethod = "total",
                   type         = "pctlim",
                   pctlim       = 0.0000001)

我收到以下錯誤：

    Error in selecthighestpcts(discrep1, inputter, pctlim) : 
      No variables are off by more than 0.00001 percent using the method you have chosen, either weighting is 
unnecessary or a smaller pre-raking limit should be chosen.

盡管這在客觀上是不正確的。 例如，考慮 quota_ed 目標：

> targets[[4]]
  HS or less Some college    Bachelors     Advanced 
        0.41         0.29         0.19         0.11 
> wpct(m1b$quota_ed)
    Advanced    Bachelors   HS or less Some college 
   0.1614583    0.3645833    0.1666667    0.3072917

任何關於我做錯了什么的想法將不勝感激。 請參閱此鏈接到 RBloggers 帖子，了解我要模擬的例程。

Answer 1

要使 anesrake function 工作，可能需要執行以下步驟：

1. 將您的權重變量轉換為因子。 確保它們不包含空級別。
2. 也從您的目標中排除空級別。 例如，假設您的數據中沒有 55 歲以上的人。 然后，您應該從 a) quota_age 變量以及 b) 從您的 m1b 數據中刪除該級別。
3. 您的列表的第一級還需要使用應該加權的特定列名命名，即在您的命令添加后： names(targets) <- c("quota_age", "quota_race", "quota_gender", "quota_ed") 。