《R语言与应用统计分析实验指导》实验2（R语言程序设计）.pptx

第一章R语言入门实验2R语言程序设计

2.1实验目的2.2实验过程目录

1. 掌握建立R函数的方法；2. 掌握R函数的顺序结构；3. 掌握R函数的分支结构；4. 掌握R函数的循环结构。2.1实验目的

2.2实验过程1.建立R函数文件R语言虽然包含了丰富的程序包，能够满足基本的数据分析应用，但是我们在使用的过程中，往往有一些个性化的需求，此时就须要编写自己的函数，以实现这些个性化的功能。R语言的函数编写，可以在R编辑器中进行。在主界面单击“文件→新建程序脚本”，打开“R编辑器”，在该窗口下编写程序。当程序编写完成后，单击“文件→保存”(或是对应的快捷键)，打开“保存程序脚本为”窗口，输入文件名后，点击“保存”按钮，可以把编写好的程序保存在当前目录下

2.函数的定义定义函数的格式为：funname-function(arg1,arg2,...)expression注：R语言定义的函数与其它语言的不一样，在定义的时候没有输出的变量，放在expression中最后的对象(数值、向量、矩阵、数组、列表或是数据框)为函数的返回值。funname：表示函数名；function：表示定义函数的关键字；arg1,arg2等：表示输入参数；expression：表示函数的主体部分，包含函数要执行的语句。

3.调用函数?

?t.stat-function(x,y){n1-length(x);n2-length(y);xb-mean(x);yb-mean(y);sx2-var(x);sy2-var(y);s-((n1-1)*sx2+(n2-1)*sy2)/(n1+n2-2);(xb-yb)/sqrt(s*(1/n1+1/n2))}

例2.1：现有两个样本如下：A：79.98,80.04,80.02,80.04,80.03,80.03,80.04,79.97,80.05,80.03,80.02,80.00,80.02B：80.02,80.03,80.02,80.00,80.02计算两个样本的t统计量。X-c(79.98,80.04,80.02,80.04,80.03,80.03,80.04,79.97,80.05,80.03,80.02,80.00,80.02);Y-c(80.02,80.03,80.02,80.00,80.02);source(t.stat.R)t.stat(X,Y)

4.分支结构If/else函数函数格式为：if(condition)expr1elseexpr2其中条件表达式condition必须返回一个逻辑值，操作符和||经常被用于条件部分。和|与和||的区别在于，和|按照逐个元素的方式进行计算，和||对向量的第一个元素进行运算，只有在必需的时候才对第二个参数求值。if/else结构的向量版本是函数ifelse，其形式为ifelse(condition,a,b)，产生函数结果的规则是：如果condition[i]为真，对应a[i]元素；反之对应的是b[i]元素。根据这个原则函数返回一个由a，b中相应元素组成的向量，向量长度与其最长的参数等长。

例2.2：输入一个成绩，如果大于等于90，则输出“优”；如果大于等于80且小于90，则输出“良”，如果大于等于70且小80，则输出“中”；如果大于等于60且小于70，则输出“及格”；如果小于60，则输出“不及格”。

chengji-function(x){if(x=90){print(优)}elseif(x=80){print(良)}elseif(x=70){print(中)}elseif(x=60){print(及格)}elseprint(不及格)}source(chengji.R)chengji(50)测试结果注：R的分支结构if/else与其它语言有很大区别，else要放在前一个if语句的{}后，否则运行将会出错，读者可自行测试。

switch语句switch适用于多分支的情况，其使用方法为：switch(expr,list)其中，expr为表达式，list为列表。如果expr的取值在1~length(list)之间，则函数返回列表相应位置的值；如果expr超出范围，则无返回值。

例2.3：用switch语句实现例2.2中的问题。chengji2-function(x){G-switch(x%/%10+1,不及格,不及格,不及格,不及格,