第四讲R语言语句组、循环和条件控制

第四讲：语句组，循环，条件控制
事件处理的逻辑模式
顺序执行
按条件作出是否执行的判断 根据某种条件变化反复执行 根据不同的情况作出相应的执行
一、组合表达式
为实现某一功能而使各种命令表达式组合在
一起形成一个复合表达式，形式一般为： { expr_1 expr_2 … expr_n }
其他循环形式
repeat循环 repeat { expr
#组合表达式
break Байду номын сангаас达到结束循环的条件,跳处循环 }
while循环 while (condition) { expr }
#while循环中首先执行condition语句， 当condition语句的执行结果为TRUE 时，执行expr，否则结束循环
例子：将容量为256M的内存的价格普调增长10%
> > + + + + + + + > count<-1 repeat { if(MemInfo$Mem[count]=="256") MemInfo$Price[count]<-MemInfo$Price[count]*1.1 count=count+1 if (count>length(MemInfo$Price)) break } MemInfo Vender ProduceName Type Mem Frequency Price 1 Vender1 256MB DDR 400 DDR 256 400 346.5 2 Vender2 256MB DDR 333 DDR 256 333 330.0 3 Vender3 512MB DDR 400 DDR 512 400 650.0 4 Vender4 256MB DDR 400 DDR 256 400 308.0 5 Vender5 256MB DDR 333 DDR 256 333 291.5 6 Vender6 256MB SDRAM 266 SDRAM 256 266 291.5 7 Vender7 256MB DDR 400 DDR 256 400 319.0 8 Vender8 256MB DDR 400 DDR 256 400 379.5 9 Vender9 512MB DDR 400 DDR 512 400 575.0
[,1] [,2] [1,] 2 8 [2,] 4 10 [3,] 6 12
可以改变计算函数FUN:+,-,*,/,^
> outer(Y,X,FUN="+") , , 1, 1 >X [,1] [,2] [1,] 1 3 [2,] 2 4 >Y [,1] [,2] [1,] 1 4 [2,] 2 5 [3,] 3 6 [,1] [,2] [1,] 2 5 [2,] 3 6 [3,] 4 7 , , 2, 1 [,1] [,2] [1,] 3 6 [2,] 4 7 [3,] 5 8 , , 1, 2 [,1] [,2] [1,] 4 7 [2,] 5 8 [3,] 6 9 , , 2, 2 [,1] [,2] [1,] 5 8 [2,] 6 9 [3,] 7 10
二、条件控制语句
if (condition1) statement1 else if (condition2) statement2 else statement3
condition语句返回一逻辑值，为TRUE时执行，为FALSE时跳过
逻辑运算符： && #逻辑与 || #逻辑或 优先对向量的第一个元素进行运算 & #逻辑与 | #逻辑或 ！ #逻辑非 按照逐个元素的方式进行计算 == #恒等 > #大于 >= #大于等于 < #小于 <= #小于等于
for (j in 1:100) { min.n<-0 norm<-rnorm(100,1,1) min.norm[j]<-100 for (i in 1:100) { if (norm[i]<=min.norm[j]) { min.norm[j]=norm[i] min.count[j]=i } } } plot(min.count,min.norm)
例子：
X<-1:9 #x<-paste(c("X","Y"),1:10,sep="") If (!is.numeric(x) && !plex(x) && !is.logical(x)) { warning("argument is not numeric or logical: returning NA") return(NA_real_) }
例子：分析outer函数的计算过程 两个数组的外积
> X<-matrix(1:4,nrow=2) > Y<-array(1:6,dim=c(3,2)) >X [,1] [,2] [1,] 1 3 [2,] 2 4 >Y [,1] [,2] [1,] 1 4 [2,] 2 5 [3,] 3 6
> outer(Y,X) #Y%o%X , , 1, 2 , , 1, 1 [,1] [,2] [1,] 1 4 [2,] 2 5 [3,] 3 6 , , 2, 1 [,1] [,2] [1,] 3 12 [2,] 6 15 [3,] 9 18 , , 2, 2 [,1] [,2] [1,] 4 16 [2,] 8 20 [3,] 12 24
outer(X,Y,FUN=“*”)：提取两个矩阵的属性
outer(X,Y,FUN=“*”)：两个矩阵的外积计算过程
定义一个数据框
> Vender<-paste(c("Vender"),1:9,sep="") > ProduceName<-c("256MB DDR 400","256MB DDR 333","512MB DDR 400", + "256MB DDR 400","256MB DDR 333","256MB SDRAM 266", + "256MB DDR 400","256MB DDR 400","512MB DDR 400") > Type<-rep(c("DDR"),9) > Type[6]<-"SDRAM" > Mem<-c(256,256,512,256,256,256,256,256,512) > Frequency<-c(400,333,400,400,333,266,400,400,400) > Price<-c(315,300,650,280,265,265,290,345,575) > MemInfo<-data.frame(Vender,ProduceName,Type,Mem,Frequency,Price) > MemInfo Vender ProduceName Type Mem Frequency Price 1 Vender1 256MB DDR 400 DDR 256 400 315 2 Vender2 256MB DDR 333 DDR 256 333 300 3 Vender3 512MB DDR 400 DDR 512 400 650 4 Vender4 256MB DDR 400 DDR 256 400 280 5 Vender5 256MB DDR 333 DDR 256 333 265 6 Vender6 256MB SDRAM 266 SDRAM 256 266 265 7 Vender7 256MB DDR 400 DDR 256 400 290 8 Vender8 256MB DDR 400 DDR 256 400 345 9 Vender9 512MB DDR 400 DDR 512 400 575
练习4.2
找出以上向量的最大值及其出现的位置。 找处100次产生的随机向量的最小值于出现位
置的关系。(用plot(min.count,min.norm)表示， 其中min.count为位置向量，min.norm为最小 值向量)
norm<-rnorm(100,1,1) max.norm<--100 for (i in 1:100) { if (norm[i]>=max.norm) { max.norm=norm[i] max.count=i } }
for (count in 1:length(MemInfo$Price)) { if (MemInfo$Frequency[count]=="266") case<-1 else if(MemInfo$Frequency[count]=="333") case<-2 else case<-3 switch(case, { MemInfo$Price[count]<-MemInfo$Price[count]*1.1 }, { MemInfo$Price[count]<-MemInfo$Price[count]*1.15 }, { MemInfo$Price[count]<-MemInfo$Price[count]*1.12 }) }
例子：定义一个长度为100、mean=1，sd=1 的正态分布随机向量，并找出该向量的最小值 及其出现的位置
norm<-rnorm(100,1,1) #随机向量 min.norm<-100 #设定初始最大值 for (i in 1:100) { if (norm[i]<=min.norm) #第i个值与现有最大值做比较， { 如果比最大值大则更新最大值， 并记录位置，否则不做处理 min.norm=norm[i] min.count=i } }