<>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

友情链接
KaDraw流程图
API参考文档
OK工具箱
云服务器优惠
阿里云优惠券
腾讯云优惠券
华为云优惠券
站点信息
问题反馈
邮箱:ixiaoyang8@qq.com
QQ群:637538335
关注微信