Função para criar DFs separados com base no valor da coluna

Eu tenho uma tabela demográfica que consiste no registro demográfico de diferentes pacientes. Na demografia, cada paciente tem o registro daquele hospital em que o paciente está registrado (Hospital_ID) .

Estou filtrando a tabela de demografia por Hospital_ID == A434 e, em seguida, mesclando o Demography df com a prescrição de paciente por Pat_ID para cada hospital. No final, terei mesclado o registro do paciente para cada ID do hospital. Em dados reais, tenho cerca de 2.000 pacientes e 270 Hospital_ID distintos. Fazer isso manualmente para cada hospital_id não é uma boa ideia.

O que eu quero é como criar uma função para isso, para que seja executada uma vez e me dê 270 DFs como saída para cada Hospital_ID . Abaixo está o exemplo e a saída de dfs.

library(dplyr)

Demography <- read.table(text="
Pat_ID  Gender  Hospital_ID
P1  M   A434
P2  F   A232
P3  F   A331
P4  M   A434
P5  F   A21
P6  M   A232", header=TRUE)

Patient_Prescription <- read.table(text="
Pat_ID  Prescription_ID     Prescription
P1  221     Hydrolin
P2  332     ospamox
P3  331     Paxix
P4  334     Danitol
P5  441     Soapsomede
P6  3311    'Dentol hydro'", header = TRUE)

Dois DFs de saída