R语言统计入门第二章R语言环境——2.3R编程

2.3R编程

x<-runif(50,0,2)
hist.with.normal<-function(x,xlab=deparse(substitute(x)),...)
{h<-hist(x,plot=F,...)
s<-sd(x)
m<-mean(x)
ylim<-range(0,h$density,dnorm(0,sd=s))
hist(x,fre=F,ylim = ylim,xlab=xlab,...)
curve(dnorm(x,m,s),add=T)}

hist.with.normal(x)

2.3.1流程控制

y<-12345
x<-y/2
while(abs(x*x-y)>1e-10)x<-(x+y/x)/2#前面的条件是真，则对后面表达式求值，abs为绝对值函数
x

[1] 111.1081

x^2

[1] 12345

x<-y/2
repeat{
  x<-(x+y/x)/2
  if(abs(x*x-y)<1e-10)break#条件是真则break结束循环
}
x

[1] 111.1081

x<-seq(0,1,.05)
plot(x,x,ylab = "y",type = "l")
for (j in 2:8)lines(x,x^j)#for循环,j为循环变量，执行lines(x,x^j)语句循环其中j分别为2到8步长为1，绘制一组幂曲线

2.3.2类和类函数

t.test(bmi,mu=22.5)$p.value#提取假设检验的p值

[1] 0.7442183