R-5_编程.pptVIP

写函数:中位数和均值的比值函数用该函数计算一万个(均值为1的)指数分布的中位数和均值的比 median.mean.ratio - function(x) { return(median(x)/mean(x)) } z - rexp(10000) median.mean.ratio(z) 函数 函数之例:中位数和均值的比值函数用该函数计算一万个(均值为1的)指数分布的中位数和均值的比 median.mean.ratio - function(x){return(median(x)/mean(x))} z - rexp(10000) median.mean.ratio(z) 在函数中的{}中间,可以有多个语句;语句之间用换行或分号(;)分开. 函数的编辑 fix(median.mean.ratio) 这时会出现一个记事本编辑器 编辑并试图存了之后如果发生错误(不会存),往往会让你编辑改过但有语法错误的fix(); 如果还使用fix(median.mean.ratio)则是编辑以前的版本. 函数 函数可以有缺省值,例如 Expo=function(y,x=2){z=y^x;z} y=seq(0,1,length=100);plot(y,Expo(y),type=l) 函数可以不写return,这时最后一个值为return的值.为了输出多个值最好使用list. liti5_1.R: n=1:10; x=sin(n*pi/10); cat(x=,x) 例1 对n=1,2,…,10,求xn= 的值 source(E:/liti5_1.R) x=0.309017 0.5877853 0.809017 0.9510565 1 0.9510565 0.809017 0.5877853 0.309017 1.224606e-16 Liti5_2.R: money=10000 years=0 while (money20000){ years=years+1; money=money*(1+2.52/100) } 例2 设银行年利率为2.52%。将10000元钱存入银行，问至少多长时间会连本带利翻一番？4.14 呢？ source(E:/liti5_2.R) years=28 此时 money= 20074.31 先建立R文件liti5_3.R定义函数f(x)，再在命令窗口输入fun1(2),fun1(-1)即可。 source(E:/liti5_3.R) fun1- function(x) { y=numeric(length(x));y[x1]=x[x1]^2+1; y[x=1]=2*x[x=1];cat(fun1=,y)} fun1(c(-1,2,3)) fun1=-2 5 10 例3: 作业:对以下问题,编写R文件: (1) 编程求出矩阵最大值及其所在的位置. 要求: 编写通用函数，然后对不同的矩阵调用输出相应的结果，比如 A=matrix(floor(rnorm(100,0,4)),4,25) 提示：考虑用which函数 (2)有一函数 ,写一程序,输入自变量的值,输出函数值. (3) 有缺失数据矩阵的插补。对缺失矩阵每一列的缺失值用该列所有观测到数据的平均值代替，得到新的无缺失值矩阵。 要求同(1), 例如 A=matrix(floor(rnorm(100,0,4)),4,25) A[2,5]=NA; A[4,8]=NA; A[3,15]=NA (4) 用二分法求函数 在(0,1)上的根，要求写出程序并给出结果. 课上留的作业 if (all(x0) all(log(x))0) { y - log(log(x)); print(cbind(x,y)); } else { cat(Unable to comply\n); } 试编写一个函数，对向量求重对数(log(log(x))) 当对X的分量重对数计算有意义时给出值，没 有意义时输出对应的下标及分量的值。 以下内容自学 R的对象 R是一种面向对象的语言。一般说来，R的对象包含了若干个元素作为其数据，这些元素的个数叫做此对象的长度（length），这些元素的共同的类型叫做此对象的模式（mode）。另外，对象还可以包含一些特殊数据，称为属性（attribute），如列表的每一个成员（元素）都可以有变量名，这些变量名组成的字符型向量为此列表的names属性。 R的面向对象能力依