Boa tarde prezados, uma vez mais, solicitando a ajuda dos amigos!!

Estou necessitando gerar valores de uma normal, condicionada a parâmetros de interesse, como no exemplo reproduzível abaixo:

x<-c("X1","X2","X3","X4")
média<-c(100,120,150,180)
dp<-c(20,15,10,25)
exemplo<-data.frame(cbind(x,média,dp))

O meu objetivo é gerar de forma automática, para cada linha, uma distribuição normal com os parâmetros nessas linhas. No exemplo, para a linha que contêm X1, gerar 1000 valores de uma normal com média=100 e dp=20, para a linha que contêm X2, gerar 1000 valores de uma normal de média=120 e dp=15, e assim sucessivamente, até o último vetor, gerado pelos parâmetros da última linha (média=180 e dp=25).  Algo similar a isso:

for (i in 1:nrow(exemplo)){
  
  simu[i]<-rnorm(100,exemplo[i,2],exemplo[i,3])
  
}

Alguém teria uma solução para isso? É que o número de linhas dessa tabela "exemplo" poderá variar.

Muito obrigado!!!

Diogo Jerônimo
Bacharel em Ciências Estatísticas - ENCE/IBGE
Mestre em Metrologia - PUC-Rio/PósMQI
http://lattes.cnpq.br/8996149312896520