В базовой графике:
r <- 3*runif(10)
degs <- 360*runif(10)
# First you want to convert the degrees to radians
theta <- 2*pi*degs/360
# Plot your points by converting to cartesian
plot(r*sin(theta),r*cos(theta),xlim=c(-max(r),max(r)),ylim=c(-max(r),max(r)))
# Add a circle around the points
polygon(max(r)*sin(seq(0,2*pi,length.out=100)),max(r)*cos(seq(0,2*pi,length.out=100)))
Обратите внимание, что хотя бы одна из точек будет на границе круга, поэтому, если вы этого не хотите, вы должны заменить max(r)
характеристики начто-то вроде 1.1*max(r)