R语言软件入门指导.pptVIP

* 第一种方式：从控制台输入数据 数据量较少时可以从控制台直接输入： height-c(1.75, 1.80, 1.65, 1.90, 1.74, 1.91) weight-c(60, 72, 57, 90, 95, 72) sq.height-height^2 ratio-weight/sq.height t.test(ratio, mu=22.5) * 第二种方式 从外部读取数据 数据量较大时用read.table函数从外部txt文件读取 第1步 将Excel中的数据另存为.txt格式（制表符间隔）或.csv格式。 第2步 用read.table()或read.csv()函数将数据读入R工作空间，并赋值给一个对象。 * 图14 在Excel中将数据存为txt文件 * 例：t检验(续) 一般从txt文档读取数据。每一行作为一个观测值。每一行的变量用制表符，空格或逗号间隔开。 read.table(”位置”, header=T) read.csv(”位置”,header=T) #从外部读取数据 data1-read.table(d:/t.test.data.txt,header=T) bmi- data1$weight/data1$height^2 t.test(bmi, mu=22.5) #t检验 * 例2：单因素方差分析 将三种不同菌型的伤寒病毒a,b,c分别接种于10，9，和11只小白鼠上，观察其存活天数，问三种菌型下小白鼠的平均存活天数是否有显著差异。 a菌株：2, 4, 3, 2, 4, 7, 7, 2, 5, 4 b菌株：5, 6, 8, 5, 10, 7, 12, 6, 6 c菌株：7,11,6, 6, 7, 9, 5, 10, 6, 3, 10 * 准备数据表 图15 数据表的准备 day和type 各为一列 * 例2：方差分析(续) #数据读取，将test1.txt中的内容保存到bac中， header=T表示保留标题行。 bac-read.table(d:/trees.txt,header=T) #将ba数据框中的type转换为因子(factor) bac$type-as.factor(bac$type) ba.an-aov(lm(day~type, data=bac)) summary(ba.an) boxplot(day~type,data=bac,col=red) * 图16 三种菌型对小白鼠影响的箱线图 * 例3-线性回归 转速rpm 20 22 24 26 28 30 32 34 36 38 40 42 杂质率% 8.4 9.5 11.8 10.4 13.3 14.8 13.2 14.7 16.4 16.5 18.9 18.5 对一批涂料进行研究，确定搅拌速度对杂质含量的影响，数据如下，试进行回归分析 表3 搅拌速度对涂料中杂质的影响 * rate-c(20, 22, 24, 26, 28, 30, 32, 34, 36, 38, 40, 42) impurity -c(8.4, 9.5, 11.8, 10.4, 13.3, 14.8, 13.2, 14.7, 16.4, 16.5, 18.9, 18.5) plot(impurity~rate) reg-lm(impurity~rate) abline(reg,col=red) summary(reg) * 四 简单编程 * 编程基础 R可以灵活的编写程序，用户自己编写的程序可以直接调用。 基本格式 函数名 - function(数据,参数1= 默认值,…) { 表达式(循环/判别); return(返回值); } 函数内部也可用#添加注释 * 程序流程控制 if if(条件) 表达式 if(条件) 表达式1 else 表达式2 举例 { 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()函数，设定返回值。 但是一个函数的返回的对象只有一个。 * 简单函数举例