Desenhar círculos em R

Estou usando o pacote plotrix para desenhar círculos.

E eu não entendi o que está errado com o meu código ... :-(

Eu tenho três pontos. O primeiro ponto(1,1) deve ser o centro do círculo. Os dois pontos seguintes(1,4) e(4,1) tem a mesma distância / raio para o centro. Então o círculo na trama deve passar por esses pontos, certo?

E eu não sei porque o círculo parece errado. Existe uma explicação?

  p1 <- c(1,1)
  p2 <- c(4,1)
  p3 <- c(1,4)
  r <- sqrt(sum((p1-p2)^2))

  plot(x=c(p1[1], p2[1], p3[1]),
       y=c(p1[2], p2[2], p3[2]), 
       ylim=c(-5,5), xlim=c(-5,5))
  draw.circle(x=p1[1], y=p1[2], radius=(r))
  abline(v=-5:5, col="#0000FF66")
  abline(h=-5:5, col="#0000FF66")

Dê uma olhada na saída produzidaAqui

questionAnswers(1)

yourAnswerToTheQuestion