日々の記録簿

日々のできごとの記録

日記

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

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

関連記事

WordPressの投稿一覧で【この投稿タイプの投稿を編集する権限がありません】と表示される

その他、【プラグインの管理画面にアクセスできない】【ファイルをアップロードできない(この親ディレクトリのアクセス権はサーバーによる書き込みを許可していますか?)】等のエラー表示がでていました。 一方で …

Multivariate Regression Analysis に関するメモ

Multivariate Regression Analysis (多変量回帰分析)に関するメモです。 多変量回帰にもいろいろあるようですが、ここでは「独立変数は同じだが従属変数が複数あるケース」を想 …

knitrで複数頁に渡る表を作成する

knitrで複数頁に渡る表を作成する必要があったので、メモです。 と言っても、longtable をセットするだけなのですが。 仮に、longtable.Rというファイルがあったとします。データはTI …

大規模調査の分析でSPSSの「ケースの重み付け」を使う場合の留意点

前回SPSSの「ケースの重み付け」を調べて満足していたのだが、Rのコードを見ていてとんでもないことに気づいた。 weighted.se <- function(x, w, na.rm = FAL …

beamerでQRコードを表示する

タイトルのとおりだが、beamerでプレゼンテーションをするときにQRコードを埋め込みたくなった。できるだけ簡単な方法を探ったところ、このサイトがヒット。 要は{textpos}を使って、任意の位置に …