2023年のベネッセの模試(?)の結果について、5教科平均のランキングがこちらにありました。
いろいろと調べるのに、ここからテキストエディターにコピペし、エクセルの北海道から順番に並んだ表に書きこんでみました。
まあ「各都道府県でベネッセの(駿台模試の?)テストを受けた人だけのデータですが。
PISA のデータがあればいいのにな・・・
エクセルの表をダウンロードできるようにしましたが、ダウンロードすると「.crdownload」という拡張子がついてきます。その部分を削除して「.xlsx」に書き換えてやると、エクセルのデータとして使えます。
なお、ウィルスチェックとかはようしませんので、ダウンロードは自己責任で。
で、2列目に人口、3列目に得点を入れ R で
> (wd <- getwd())
> library(psych)
> x <- read.csv(paste0(wd,"/jinkou&seiseki.csv"))
> library(psych)
> x <- read.csv(paste0(wd,"/jinkou&seiseki.csv"))
> names(x)
[1] "todoufuken" "jinkou" "seiseki"
[1] "todoufuken" "jinkou" "seiseki"
> with(data = x,plot(x[ ,2],x[ ,3]))
とやると
となって相関係数は
> corr.test(x$jinkou, x$seiseki)
Call:corr.test(x = x$jinkou, y = x$seiseki)
Correlation matrix
[1] 0.71
Sample Size
[1] 47
These are the unadjusted probability values.
The probability values adjusted for multiple tests are in the p.adj object.
[1] 0
Call:corr.test(x = x$jinkou, y = x$seiseki)
Correlation matrix
[1] 0.71
Sample Size
[1] 47
These are the unadjusted probability values.
The probability values adjusted for multiple tests are in the p.adj object.
[1] 0
ってことでした。