第四讲程序设计.pptVIP

第四讲-程序设计;程序控制结构 ;分支结构 (if/else、switch语言)if/else语句;例如，如果变量lambda为缺失值就给它赋一个缺省值，可用： if(is.na(lambda)) lambda - 0.5; 又比如要计算向量x的重对数，这只有在元素都为正且对数都为正时才能做到，因此需要先检查： if (all(x0) all(log(x))0) { y - log(log(x)); print(cbind(x,y)); } else { cat(Unable to comply\n); } 注意“”表示“与”，它是一个短路运算符，即第一个条件为假时就不计算第二个条件。如果不这样此例中计算对数就可以有无效值。 在条件中也可以用“||”（两个连续的竖线符号）表示“或”，它也是短路运算符，当第一个条件为真时就不再计算第二个条件。;在用R编程序时一定要时刻牢记R是一个向量语言，几乎所有操作都是对向量进行的。而R中的if语句却是一个少见的例外，它的判断条件是标量的真值或假值。比如，我们要定义一个分段函数f(x)，当x为正时返回1，否则返回0，马上可以想到用if语句实现如下： if(x0) 1 else 0 当x是标量时这个定义是有效的，但是当自变量x是一个向量时，比较的结果也是一个向量，这时条件无法使用。所以，这个分段函数应该这样编程： y= numeric(length(x)) y[x0] - 1 y[x=0] - 0 # 这句可以省略 ;对于条件语句是向量形式，R提供了ifelse()函数，它的使用格式为 ifelse(condition, a, b) 它表示，当condition[i]成立的时候，对应的取值为a[i],否则取值为b[i]。最终返回一个和参数向量同长的向量 x-c(55, 60, 85, 90); a-pass； b-fail ifelse(x=60, a, b) [1] fail pass pass pass a-c(a1, a2, a3, a4) b-c(b1, b2,b3,b4) ifelse(x=60, a, b) [1] b1 a2 a3 a4;有多个if语句时else与最近的一个配对。可以使用if ... else if ... else if ... else ...的多重判断结构表示多分支。多分支也可以使用switch()函数。 x=0; if (x0) { value -x } else if (x0) { value - -x } else if (x==0) { value=ZERO} value [1] ZERO ;switch语句;当list是有名定义时，statement等于变量名时，返回该变量名所对应的值，否则返回NULL x - fruit; a - switch (x, fruit=apple, drink=coffee, vegetable=broccoli); a [1] apple x - meat; a - switch (x, fruit=apple, drink=coffee, vegetable=broccoli); a [1] NULL ;循环结构(for/while/repeat)for语句 ;当然，如果只是要求各元素的和，只要调用sum(x)即可。所以，编写程序要精简，应尽量避免使用显式循环。 我们再举一个例子。比如，我们要计算同生日的概率。假设一共有365个生日（只考虑月、日），而且各生日的概率是相等的（这里忽略了闰年的情况以及可能存在的出生日期分布的不均匀）。设一个班有n个人，当n大于365时至少两个人生日相同是必然事件。当n小于等于365时，我们可以计算P{至少有两人同生日}= 1 - P{n个人生日彼此不同}，这时，n个人的生日可取值数为365n ，而n个人彼此不同的可能数为365中取n个的排列数，彼此不同的概率为P3n65 。因此，为了计算n＝1，2，...，365的情况下的同生日概率，可以用如下循环实现：; x = numeric(365) for(i in 1:365) { x[i]= 1