现代统计图形 - 科学网—博客

158 附录 A 程序初步
A.1.5 函数
R中也可以自定义函数,以便编程中可以以不同的参数值重复使用一
段代码。定义函数的方式为:function(arglist) expr return(value);
其中arglist是参数列表,expr是函数的主体,return()用来返回函数值。例如
我们定义一个求峰度 n
i=1 (xi − ¯x) 4 /(n · Var(x)) − 3的函数kurtosis()如下:
1 > kurtosis = function(x, na.rm = FALSE) {
2 + # 去掉缺失值
3 + if (na.rm)
4 + x = x[!is.na(x)]
5 + return(sum((x - mean(x))^4)/(length(x) * var(x)^2) -
6 + 3)
7 + }
该函数有两个参数,数据向量x和是否删除缺失值na.rm,后者有默认
值FALSE;下面我们用均匀分布的随机数测试一下:
1 > # 理论上均匀分布的峰度应该小于1,以下为一个模拟结果
2 > kurtosis(runif(100))
[1] -0.911
注意,R中有时候不必用函数return()指定返回值,默认情况下,函数
主体的最后一行即为函数返回值。
最后,因为本书的很多作图函数都是泛型函数,我们也补充一点泛型
函数(generic function)的作用原理。实际上泛型函数是根据第一个参数
的类(class)去调用了相应的“子函数”。以plot()为例,我们用methods()查
看一下它有哪些作图方法:
1 > methods("plot")
[1] plot.Date* plot.HoltWinters*
[3] plot.POSIXct* plot.POSIXlt*
[5] plot.TukeyHSD plot.aareg*
[7] plot.acf* plot.correspondence*
[9] plot.cox.zph* plot.data.frame*
[11] plot.decomposed.ts* plot.default
[13] plot.dendrogram* plot.density
[15] plot.ecdf plot.factor*
[17] plot.formula* plot.hclust*
[19] plot.histogram* plot.isoreg*