2.9. Heatmap

library(ggplot2)

## Warning: пакет 'ggplot2' был собран под R версии 4.2.1

library(AER)

## Загрузка требуемого пакета: car

## Загрузка требуемого пакета: carData

## Загрузка требуемого пакета: lmtest

## Загрузка требуемого пакета: zoo

## 
## Присоединяю пакет: 'zoo'

## Следующие объекты скрыты от 'package:base':
## 
##     as.Date, as.Date.numeric

## Загрузка требуемого пакета: sandwich

## Загрузка требуемого пакета: survival

library(dplyr)

## 
## Присоединяю пакет: 'dplyr'

## Следующий объект скрыт от 'package:car':
## 
##     recode

## Следующие объекты скрыты от 'package:stats':
## 
##     filter, lag

## Следующие объекты скрыты от 'package:base':
## 
##     intersect, setdiff, setequal, union

library(viridis)

## Загрузка требуемого пакета: viridisLite

library(gridExtra)

## 
## Присоединяю пакет: 'gridExtra'

## Следующий объект скрыт от 'package:dplyr':
## 
##     combine

data("DoctorVisits")
head(DoctorVisits)

##   visits gender  age income illness reduced health private freepoor freerepat
## 1      1 female 0.19   0.55       1       4      1     yes       no        no
## 2      1 female 0.19   0.45       1       2      1     yes       no        no
## 3      1   male 0.19   0.90       3       0      0      no       no        no
## 4      1   male 0.19   0.15       1       0      0      no       no        no
## 5      1   male 0.19   0.45       2       5      1      no       no        no
## 6      1 female 0.19   0.35       5       1      9      no       no        no
##   nchronic lchronic
## 1       no       no
## 2       no       no
## 3       no       no
## 4       no       no
## 5      yes       no
## 6      yes       no

g1 <- ggplot(data = DoctorVisits, aes(x = age, y = income)) + geom_point() +
  scale_fill_continuous(type = "viridis") + ggtitle("Диаграмма рассеивания") + theme_bw()
g1

g2 <- ggplot(data = DoctorVisits, aes(x = age, y = income)) + geom_bin2d(bins = 20) +
  scale_fill_continuous(type = "viridis") + ggtitle("Тепловая карта #1 с bins = 20") + theme_bw()
g2

g3 <- ggplot(data = DoctorVisits, aes(x = age, y = income)) + geom_bin2d(bins = 5) +
  scale_fill_continuous(type = "viridis") + ggtitle("Тепловая карта #2 с bins = 5") + theme_bw()
g3

#непрерывные величины

x <- rnorm(1000, 50, 20)
y <- rnorm(1000, 20, 50)

data <- data.frame(x = x, y = y)
g4 <- ggplot(data = data, aes(x = x, y = y)) + geom_bin2d(bins = 25) + scale_fill_viridis_c(option = "magma") + theme_bw() + ggtitle("Тепловая карта #3 для непрерывных величин с bins = 25")
g4