R 语言实践一书中的一个问题

firelife

中文版第134页

<br />
vars <- c("mpg","hp","wt")<br />
dstats <- function(x)(c(mean=mean(x), sd=sd(x)))<br />
by(mtcars[vars], mtcars$am, dstats)<br />

</p>

不能运行

Error in is.data.frame(x) :

(list) object cannot be coerced to type 'double' In addition: Warning message:

In mean.default(x) : 参数不是数值也不是逻辑值：回覆NA

这是为什么？

jasmine_007

<br />
vars <- c("mpg", "hp", "wt")<br />
dstats <- function(x) {<br />
  mean=sapply(x,mean)<br />
  sd=sapply(x,sd)<br />
  c(mean,sd)<br />
}<br />
by(mtcars[vars],mtcars$am,dstats)<br />

返回mean,sd

or

<br />
by(mtcars[vars],mtcars$am,colMeans)<br />

返回mean。
</p>

firelife

有点明白了，因为传递给dstats的是一个dataframe

那书上的代码是什么回事

jasmine_007

回复第3楼的 firelife：因为mean()的参数必须是numeric或者logical，而书上传给dstats()的是以mtcars$am分组的mtcars[vars],是一个dataframe，所以报错。对dataframe求各列均值需要用sapply()函数。假如传递的只是一列数据比如，

<br />
> dstats <- function(x)(c(mean=mean(x), sd=sd(x)))<br />
> by(mtcars$mpg, mtcars$am, dstats)<br />
mtcars$am: 0<br />
     mean        sd<br />
17.147368  3.833966<br />
-------------------------------------------------------------------------<br />
mtcars$am: 1<br />
     mean        sd<br />
24.392308  6.166504<br />

就不会出错了。因为传递给dstats()的是一列数据，R会利用降维的思想，变成numeric类型。
</p>