R：绘制2000种不同随机卡方模拟的平均值和中值

我试图对随机卡方分布进行2000次模拟，n=300，df=3。从每次模拟中，我想取样本平均值和中位数，然后最终将平均值和中位数绘制为直方图

我尝试使用repeat循环函数，直方图在满足“if”条件后出现

n = 300
df = 3
randomchisq = rchisq(n = n,df = df)
i = 1
reps = 2000
xvals = "Chi squared sims"
twokchisqmean <- i:reps
twokchisqmed <- i:reps

repeat {
  twokchisqmean[i] = mean(randomchisq);
  twokchisqmed[i] = median(randomchisq);
  i = i+1;

  if (i==reps) {
    par(mfrow=c(1,2))
    hist(twokchisqmean, 
         main="Histogram of mean for 2000 chi squared sims",
         xlab = xvals)

    hist(twokchisqmed,
         main = "Histogram of median for 2000 chi squared sims",
         xlab = xvals)

  }
}

n=300
df=3
randomchisq=rchisq（n=n，df=df）
i=1
重复次数=2000次
xvals=“卡方模拟人生”
twokchisqmeanA
repeat
语句将循环，直到达到中断，然后才从循环返回。由于在
repeat
语句中没有
break
，因此代码将无限运行。当你想让它停下来的时候，一定要让它停下来：

n=300
df=3
randomchisq=rchisq（n=n，df=df）
i=1
重复次数=2000次
xvals=“卡方模拟人生”
twokchisqmean这不会只是运行一个无限循环，因为在
repeat
语句中没有
break
？感谢您的响应，尽管在我包含break之后，绘图似乎仍然没有出现。我已经更新了代码以包含所有变量，在干净的R环境中运行它会生成很好的绘图，也许试着运行这个新代码？如果这仍然不起作用，让我知道这似乎是我的环境问题。非常感谢你！别担心！很高兴我能帮忙：）