＃

Question

我尝试在r＆＃34中编写程序;使用函数＆＃34;从任何分布生成随机样本。但它在hist.default（xbars）中显示＆＃34;错误：＆＃39; x＆＃39;必须是数字＆＃34; 我的程序在这里

sim.clt <- function(n, ran.func,..., simsize,...)
{
  xbars<-vector()
  for(i in 1:simsize=simsize)
  {
    x<-function(ran.func)

      xbars[i]<-mean(x)
 }
par(mfrow=c(2,1))
hist(xbars)
qqnorm(xbars)
return(xbars)
}
sim.out<-sim.clt(n=20,ran.func="rexp",simsize=5000)
shapiro.test(sim.out)

＃

我是r编程的新手，所以无法弄明白，如何解决问题。感谢...

Answer 1

这里有很多不妥之处。

for(i in 1:simsize=simsize)

应该抛出错误：

> for(i in 1:simsize=simsize) { print(i)}
Error: unexpected '=' in "for(i in 1:simsize="

更好的是

for(i in seq_len(simsize))

然后

x <- function(ran.func)

没有做你想象的那样;它返回一个以xbars[i]<-mean(x)为主体的函数，如：

> x <- function(ran.func)
+ 
+       xbars[i]<-mean(x)
> x
function(ran.func)

      xbars[i]<-mean(x)
> is.function(x)
[1] TRUE

我想你想打电话给ran.func所以你可能需要

FUN <- match.fun(ran.func)
x <- FUN()

但是这会失败，因为您似乎没有传递ran.func的任何参数，即使只是使用n的示例中的rexp。

错误消息源于最后一点。您将xbars定义为空vector()，默认情况下会创建一个空逻辑向量：

> xbars <- vector()
> xbars
logical(0)
> is.numeric(xbars)
[1] FALSE

现在，如果您在定义x时没有犯错，那就不会有问题了（回忆xbars[i]<-mean(x)现在位于函数x的正文中，并且从未明确调用过），这意味着xbars仍然是一个空的逻辑向量。由于这不是数字，hist会抛出您看到的错误。

另一个错误是你不能在函数定义中使用...两次。您是否尝试将第一个包含参数传递给ran.func而将第二个...包含在其他内容中。你不能在R中做到这一点。

这是你想要的吗？

sim.clt <- function(n, ran.func, simsize, ...) {
  ## ... passed to ran.func for other parameters of distribution
  xbars <- numeric(simsize)
  for(i in seq_len(simsize)) {
    FUN <- match.fun(ran.func)
    x <- FUN(n = n, ...)
    xbars[i] <- mean(x)
  }
  ## plot
  op <- par(mfrow=c(2,1))
  on.exit(op)
  hist(xbars)
  qqnorm(xbars)
  xbars
}

> sim.out<-sim.clt(n=20,ran.func="rexp",simsize=5000)
> shapiro.test(sim.out)

    Shapiro-Wilk normality test

data:  sim.out
W = 0.9867, p-value < 2.2e-16

R：x中的错误必须是数字

＃

1 个答案: