For my data.frame `full` below, I'm wondering how to create a two-panel `geom_point` such that on the first panel, we have `ols.(Intercept)` (x-axis) plotted against `hlm.(Intercept)`, AND on the second panel, we have `ols.ses` (x-axis) plotted against `hlm.ses`?

```
library(lme4) 
library(tidyverse)

hsb <- read.csv('https://raw.githubusercontent.com/rnorouzian/e/master/hsb.csv')

fit <- lmer(math~ses+(ses|sch.id), data= hsb)
ch <- unique(hsb$sch.id)

ols <- map_dfr(ch,~coef(lm(math~ses, data=hsb,subset=sch.id==.)))
mlm <- coef(fit)$sch
full <- cbind(ols=ols, hlm=mlm, sch.id=ch)

head(full, n = 1)
     ols.(Intercept)  ols.ses hlm.(Intercept)  hlm.ses sch.id
1224        10.80513 2.508582        11.06002 2.504083   1224

```