日々の記録簿

日々のできごとの記録

日記

surveyで相関係数を計算する

投稿日:

survey packageに相関係数を算出するコマンドがなかったのでメモ。こちらのサイトが参考になった。

library(survey)

data(api)
dstrat <- svydesign(
  id = ~1, strata = ~stype, weights = ~pw, data = apistrat, fpc = ~fpc
)

var <- svyvar(~ api00 + api99, design = dstrat)
cov2cor(as.matrix(var))

これで相関係数が算出される。標準誤差がほしい場合はややこしくなる。先のサイトでは、weights関数を使って計算している模様。

svycor <- function(model, design) {
  mf <- model.frame(model, model.frame(dstrat), na.action = na.pass)
  wts <- weights(dstrat, "sampling")

  wcors <- weights::wtd.cor(mf,
    weight = wts, bootse = TRUE, mean1 = TRUE,
    bootn = 1000, bootp = TRUE
  )
  return(wcors)
}

cors <- svycor(~api00 + api99, dstrat)

cors$correlation
cors$p.value
cors$std.err

スポンサーリンク

スポンサーリンク

-日記

執筆者:


comment

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

関連記事

X1 carbon gen 10でxf86-video-intelをインストールするといろいろおかしくなる 

タイトルで終わっているが、Thinkpad X1 carbon (gen 10)でxf86-video-intelをインストールしてはいけない。Google Chromeが固まる、VLCで動画が再生で …

Thinkpad X1 Carbon 2016 (4th gen)に、Arch Linuxをインストールする方法(その2)

Thinkpad X1 Carbon 2016 (4th gen)でLinuxが起動しないという問題の続報です。 この問題は、X1 Carbonだけでなく、Yoga260やT460sでも発生しているよ …

Thinkpad X1 Carbon 2016 (4th gen)に、Arch Linuxをインストールする方法

またまた物欲に負けて、Thinkpad X1 Carbon 2016 (4th gen)を購入してしまいました。 さっそく、Arch Linuxをインストールしようとしたのですが、起動してこない・・・ …

デスクトップ環境をXfceに変更してみた!

なんとなくGnomeやめようかな・・・と思って、デスクトップ環境をXfceに変更してみました。 ついでに、NetworkManagerとか、Ibusとか、Evinceとか、代替する方法あるのかなあ?・ …

Thinkpad X1 Carbon 2016 (4th gen)に、Arch Linuxをインストールする方法(その3)

Thinkpad X1 Carbon 2016 (4th gen)でLinuxが起動しないという問題の現状です。 最新のkernel 4.6rc3を試したところ、OneLink+ ドックを繋いだままで …