二元运算符的非数字参数是我得到的错误

Question

list =ls()
set.seed(1234)
S<- numeric()
Y_bar <- numeric()
Xsq <- numeric()
T_st <- numeric()

for (i in c(1:10000)) {
  sample = rnorm(n= 20, mean=2, sd = 5) 
  s = var(sample)
  S <-append(S,s)
  
  y_bar = mean(sample)
  Y_bar<-append(Y_bar,y_bar)
  
  xsq = n-1*s**2/5**2
  Xsq <-append(Xsq,xsq)
  
  t_st = y_bar - mean/( as.numeric(s)/sqrt(n))
  T_st <-append(T_st,t_st)
  
}

mean/(as.numeric(s)/sqrt(n)) 错误：二元运算符的非数字参数

sd^2 错误：二元运算符的非数字参数

我需要所有这些计算值，如s和y_bar用于我以后的计算。

Answer 1

我的感觉是有些混乱。 您在rnorm function 中将均值、 mean和样本大小n指定为参数，但未将它们分配为全局环境中的对象。 他们没有被保存（因为他们会mean <- 5 ）。

尽可能坚持你的原始代码（虽然这可以优化），运行这个循环的更合适的方法是在全局环境中定义你的样本大小、均值和标准差，并使用它来指定参数（然后在整个过程中进行相应的替换）。

set.seed(1234)
S <- numeric()
Y_bar <- numeric()
Xsq <- numeric()
T_st <- numeric()

# specify sample size, mean, and sd
sample_n <- 20
sample_mn <- 2
sample_sd <- 5

for (i in c(1:10000)) {
  # specify parameters from those defined in the global environment
  sample = rnorm(n = sample_n, mean = sample_mn, sd = sample_sd)
  s = var(sample)
  S <- append(S, s)
  
  y_bar = mean(sample)
  Y_bar <- append(Y_bar, y_bar)
  
  # Change object names accordingly
  xsq = sample_n - 1 * s ** 2 / 5 ** 2
  Xsq <- append(Xsq, xsq)
  
  t_st = y_bar - sample_mn / (as.numeric(s) / sqrt(sample_n))
  T_st <- append(T_st, t_st)
}

请注意，我使用了像sample_mn这样的术语，因为mean <- 5不是好的做法，因为mean()已经是 function。