日々の記録簿

日々のできごとの記録

日記

SPSSの「ケースの重み付け」とは何なのか

投稿日:

前回の続き。だいたい納得したのだが、SPSSの「ケースの重み付け」が一体何をしているのかもう少し調べてみたくなった。「ケースの重み付け」を設定したとき、分散や標準偏差はどうやって計算されているのだろう?

weighted.var <- function(x, w, na.rm = FALSE) {
  if (na.rm) {
    w <- w[!is.na(x)]
    x <- x[!is.na(x)]
  }
  sum(w * (x - weighted.mean(x, w))^2) / (sum(w) - 1)
}
weighted.se <- function(x, w, na.rm = FALSE) {
  if (na.rm) {
    w <- w[!is.na(x)]
    x <- x[!is.na(x)]
  }
  v <- sum(w * (x - weighted.mean(x, w))^2) / (sum(w) - 1)
  sqrt(v) / sqrt(sum(w))
}

d <- data.frame(
  ans = c(1, 2, 2, 3, 3, 4, 4, 1, 2, 3),
  gender = c(rep(1, 7), rep(2, 3)),
  w = c(rep(5 / 7, 7), rep(5 / 3, 3))
)

weighted.mean(d$ans, d$w) # 平均
sqrt(weighted.var(d$ans, d$w)) # 標準偏差
weighted.se(d$ans, d$w) # 標準誤差

というわけでRのコード。例題は前回と同じくこちらのサイトのものを流用。これならウェイトが整数でないにもかかわらず、それっぽい(誤った)数値が算出される。na.rmのオプションも付けたので、欠測があっても大丈夫。もっとも使うことはないだろうが。

スポンサーリンク

スポンサーリンク

-日記

執筆者:


comment

メールアドレスが公開されることはありません。 が付いている欄は必須項目です

関連記事

Rを勉強するためにいろいろ寄り道をしていた

Rをもっと勉強したいなと思って、『R言語徹底解説』を買ってみました。・・・が、ちんぷんかんぷん。クラスとか、オブジェクト指向とか、遅延評価とか、さっぱりわかりません。 これはいかんと思い、ここ2ヶ月ほ …

Nexus 6を購入し、Cyanogenmod 12を導入した!

Google Playの在庫不足で、まったく入手できなかったのですが、ようやくNexus 6が手に入りました。 購入したのは32GB版です。3月11日にGoogle Playで注文し、「4〜5週間かか …

nftablesが停止してしまう?

今朝まで気づいていなかったのだが、起動時にftablesが停止している旨のログがでていて、心底驚いた。「systemctl is-active nftables」で確認しても停止していると表示される。 …

XfburnでDVDドライブが認識されない!

XfburnでDVDを焼こうと思ったら、DVDドライブが認識されていないというエラーに遭遇したので、その対処法です。 症状は、「DVDドライブは接続済みにもかかわらず、メッセージがでて、ドライブが認識 …

Rで日本の校区地図を描く

Rで日本の校区地図を描く方法です。 データの準備 国土数値情報 ダウンロードサービスより、校区地図のもとになるデータ(shapeファイル)をダウンロードする。 データは、リンク先ページの、JPGIS2 …