【Rで多変量解析#8】数量化１類

投稿者 Dken2021年9月16日2021年9月16日

記事の目的

数量化１類をRを使用して実装していきます。データの作成から実装するので、コピペで再現することが可能です。

目次

ライブラリ
データの作成
データの可視化
(データの変換)
数量化１類

1 ライブラリ

library(dplyr)
library(ggplot2)

2 データの作成

set.seed(100)
N <- 30
講義 <- c(rep("講義A", N), rep("講義B", N), rep("講義C", N))
e1 <- rnorm(N, 0, 10)
e2 <- rnorm(N, 0, 10)
e3 <- rnorm(N, 0, 10)
data <- data.frame(講義) %>% 
mutate(テスト=ifelse(講義=="講義A", 50+e1 , ifelse(講義=="講義B", 30+e2, 80+e3))) %>%
mutate(テスト=round(テスト))
data %>% sample_n(5)

3 データの可視化

data %>%
  ggplot() +
  geom_point(aes(講義, テスト)) +
  theme_classic(base_family = "HiraKakuPro-W3") +
  theme(text=element_text(size=30)) +
  labs(title="データプロット", x=NULL)

4 (データの変換)

data %>%
  mutate(講義B = ifelse(講義=="講義B", 1, 0)) %>%
  mutate(講義C = ifelse(講義=="講義C", 1, 0)) %>%
  sample_n(5)

5 数量化１類

model1 <- lm(data=data, テスト~講義)
summary(model1)