如何在 R 中使用樣本權重

Question

我計划對關於道德問題的觀點進行多組確認因素分析。 我將比較來自比利時瓦隆區和佛蘭德斯地區的人。 我的兩個樣本需要加權，以便在年齡、性別、教育和政黨選擇方面代表他們的人口。

我的數據集中已經提供的采樣權重。 然后我創建了一個變量wreg ，將來自瓦隆和佛蘭德斯的受訪者的權重結合起來。

我是 R 新手，閱讀了有關lavaan.survey和svydesign文檔以了解代碼。 但是，我還沒有成功地寫出正確的東西。 我總是收到有關重量部分的錯誤消息。 顯然程序無法正確讀取采樣權重變量。

這是我使用的代碼：

library(lavaan.survey)

f <- "C:/.../bges07_small.csv"
s <- read.csv(f,sep=";")

r <- s[is.na(s$flawal),]

rDesign <- svydesign(ids=~1, data=r, weights=~wreg)

model.1 <- 'ethic =~ q96_1+ q96_2 +q96_3'

fit <- cfa(model.1, data=r,ordered=c("q96_1","q96_2","q96_3"))

summary(fit, fit.measures=TRUE, modindices=FALSE,standardized=FALSE)

這是我收到的錯誤消息：

Erreur dans 1/as.matrix(weights) : 
  argument non numérique pour un opérateur binaire

關於我應該如何用 R 編寫模型的任何建議？ 非常感謝！

Answer 1

從summary(r$wreg)的結果看來，您的權重列是一個因素，而不是數字向量。 確保您已正確讀取數據，並且該列不包含任何類似字符的值。 您可以手動轉換它

r$wreg <- as.numeric(r$wreg)

在運行模型之前。 此外，那些看起來像非常大的權重值。 你確定他們是正確的嗎？