第四讲-程序设计.ppt

例如，如果变量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 若a，b和x不是等长的，则较短的循环使用 switch语句 switch是多分支语句，其用法为： switch (statement, list); statement是一个表达式，list是列表，也可以用有名定义。根据表达式与list的关系返回一个值。如果表达式返回值属于1:length(list)中的一个，则返回list中相应位置的值，否则返回NULL。 a - switch (1, 2+2, mean(1:10), rnorm(4)); a [1] 4 a - switch (2, 2+2, mean(1:10), rnorm(4)); a [1] 5.5 x-3 a - switch (x, 2+2, mean(1:10), rnorm(4)); a [1] 0.3202921 0.3057707 -0.4922125 0.1280950 a - switch (4, 2+2, mean(1:10), rnorm(4)); a [1] NULL 当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语句 循环结构中常用的是for循环，是对一个向量或列表的逐次处理，格式为 for( name in values) {expression} 如：求1到100的和 sum-0; for(i in 1:100){ sum-sum+i; i-i+1 } sum; [1] 5050 当然，如果只是要求各元素的和，只要调用sum(x)即可。所以，编写程序要精简，应尽量避免使用显式循环。 我们再举一个例子。比如，我们要计算同生日的概率。假设一共有365个生日（只考虑月、日），而