R语言及在统计学上的应用.pptVIP

* 画图(image) data(volcano) x - 10*(1:nrow(volcano)) y - 10*(1:ncol(volcano)) image(x, y, volcano, col = terrain.colors(100), axes = FALSE) contour(x, y, volcano, levels = seq(90, 200, by=5), add = TRUE, col = peru) axis(1, at = seq(100, 800, by = 100)) axis(2, at = seq(100, 600, by = 100)) box() title(main = Maunga Whau Volcano, font.main = 4) * * 画图(箭头和线段) x=runif(12);y=rnorm(12);I=order(x,y);x=x[i];y=y[i];plot(x,y, main=arrows(.) and segments(.))## draw arrows from point to point : s=seq(length(x)-1)# one shorter than data arrows(x[s],y[s],x[s+1],y[s+1],col= 1:3) s=s[-length(s)];segments(x[s],y[s], x[s+2],y[s+2],col=pink) 第一部分画点,第二部分画箭头,第三部分画线段 * * 四 编写函数 * 编程基础 R可以灵活的编写程序，用户自己编写的程序可以直接调用。编程时无需声明变量的类型，这与C,C++等语言不同。 基本格式 函数名 - function(数据,参数1= 默认值,…) { 异常处理; 表达式(循环/判别); return(返回值); } 函数内部也可用#添加注释 * 函数实例 图28 函数实例 data2mat() * 程序流程控制 if if(条件) 表达式 if(条件) 表达式1 else 表达式2 举例 p = 0.03 { if(p=0.05) print(p = 0.05!) else print(p 0.05!) } * 循环 for, while for(变量 in 向量) 表达式 用法： for(i in 1:10) print(i) while(条件) 表达式 用法: i - 1 while(i10){ print(i); i - i + 1} * 返回值 返回值表示函数输出的结果。 返回值必须是一个对象。 R默认将最后一行作为返回值。 如果函数的结果需要有多个返回值，可以创建一个list()，并返回该对象。 也可以用return()函数，设定返回值。 但是一个函数的返回的对象只有一个。 * 异常处理 如数据输入不能满足要求，或者参数设定错误等等，可能造成函数给出错误的结果，则需要对函数的运行过程发出警告或终止，以提高程序的稳健性。 警告的写法 if(any(is.na(inputdata))) inputdata - na.omit(inputdata) cat(“NAs are found in the input data, and has been removed.\n) 终止的写法 if(any(is.na(xx))) stop(NAs are not allowed!\n) * 函数举例 问题:输入直角三角形的两个边长，求其斜边长。 定义函数： rcal - function(x,y){ z - x^2 + y^2 result - sqrt(z) return(result) } 调用函数： rcal(3,4) * 六 数据保存 * write.table() write.csv() save.image() sink() unlink() 若有LaTeX基础，可以用 Sweave() 函数 该函数能将脚本、程序说明和运算结果直接保存成.tex文件，用LaTeX编译成pdf文件。 数据保存 * R中的编程 编程的重要性 一个统计软件(包)应该包罗万象—所有统计方法或技术，可以完成所有的任务？？ 这是一种错误的观点！！ 太贵！！ 费时！！ SAS是一个世界上最为优秀的统计分析软件之一，但……！！当然SAS具体编程能力，而SPSS更糟！ 通过编程可以实现(开发)满足自己需要的函数或宏包 * 好的编程习惯 为了他人，更为你本人! 你的程序应该具有 可读性(