2 -1 2
r2-2r1
1 1 1
0 -3 0
-->
1 0 1
0 1 0
所以所求向量为 (1,0,-1)
R语言常用函数整理本篇是基础篇，即R语言自带的函数。 vector：向量 numeric：数值型向量 logical：逻辑型向量 character；字符型向量 list：列表 data.frame：数据框 c：连接为向量或列表 length：求长度 subset：求子集 seq，from:to，sequence：等差序列 rep：重复 NA：缺失值 NULL：空对象 sort，order，unique，rev：排序 unlist：展平列表 attr，attributes：对象属性 mode，class，typeof：对象存储模式与类型 names：对象的名字属性 字符型向量 nchar：字符数 substr：取子串 format，formatC：把对象用格式转换为字符串 paste()、paste0()不仅可以连接多个字符串，还可以将对象自动转换为字符串再相连，另外还能处理向量。 strsplit：连接或拆分 charmatch，pmatch：字符串匹配 grep，sub，gsub：模式匹配与替换 complex，Re，Im，Mod，Arg，Conj：复数函数 factor：因子 codes：因子的编码 levels：因子的各水平的名字 nlevels：因子的水平个数 cut：把数值型对象分区间转换为因子 table：交叉频数表 split：按因子分组 aggregate：计算各数据子集的概括统计量 tapply：对“不规则”数组应用函数 dev.new() 新建画板 plot()绘制点线图,条形图,散点图. barplot( ) 绘制条形图 dotchart( ) 绘制点图 pie( )绘制饼图. pair( )绘制散点图阵 boxplot( )绘制箱线图 hist( )绘制直方图 scatterplot3D( )绘制3D散点图. par()可以添加很多参数来修改图形 title( ) 添加标题 axis( ) 调整刻度 rug( ) 添加轴密度 grid( ) 添加网格线 abline( ) 添加直线 lines( ) 添加曲线 text( ) 添加标签 legend(） 添加图例 +, -, *, /, ^, %%, %/%：四则运算 ceiling，floor，round，signif 1、round() #四舍五入 例：x <- c(3.1416, 15.377, 269.7) round(x, 0) #保留整数位 round(x, 2) #保留两位小数 round(x, -1) #保留到十位 2、signif() #取有效数字(跟学过的有效数字不是一个意思) 例：略 3、trunc() #取整 floor() #向下取整 ceiling() #向上取整 例：xx <- c(3.60, 12.47, -3.60, -12.47) trunc(xx) floor(xx) ceiling(xx) max，min，pmax，pmin：最大最小值 range：最大值和最小值 sum，prod：向量元素和，积 cumsum，cumprod，cummax，cummin：累加、累乘 sort：排序 approx和approx fun：插值 diff：差分 sign：符号函数 abs，sqrt：绝对值，平方根 log, exp, log10, log2：对数与指数函数 sin，cos，tan，asin，acos，atan，atan2：三角函数 sinh，cosh，tanh，asinh，acosh，atanh：双曲函数 beta，lbeta，gamma，lgamma，digamma，trigamma，tetragamma，pentagamma，choose ，lchoose：与贝塔函数、伽玛函数、组合数有关的特殊函数 fft，mvfft，convolve：富利叶变换及卷积 polyroot：多项式求根 poly：正交多项式 spline，splinefun：样条差值 besselI，besselK，besselJ，besselY，gammaCody：Bessel函数 deriv：简单表达式的符号微分或算法微分 array：建立数组 matrix：生成矩阵 data.matrix：把数据框转换为数值型矩阵 lower.tri：矩阵的下三角部分 mat.or.vec：生成矩阵或向量 t：矩阵转置 cbind：把列合并为矩阵 rbind：把行合并为矩阵 diag：矩阵对角元素向量或生成对角矩阵 aperm：数组转置 nrow, ncol：计算数组的行数和列数 dim：对象的维向量 dimnames：对象的维名 rownames，colnames：行名或列名 %*%：矩阵乘法 crossprod：矩阵交叉乘积（内积） outer：数组外积 kronecker：数组的Kronecker积 apply：对数组的某些维应用函数 tapply：对“不规则”数组应用函数 sweep：计算数组的概括统计量 aggregate：计算数据子集的概括统计量 scale：矩阵标准化 matplot：对矩阵各列绘图 cor：相关阵或协差阵 Contrast：对照矩阵 row：矩阵的行下标集 col：求列下标集 solve：解线性方程组或求逆 eigen：矩阵的特征值分解 svd：矩阵的奇异值分解 backsolve：解上三角或下三角方程组 chol：Choleski分解 qr：矩阵的QR分解 chol2inv：由Choleski分解求逆 ><，>，<=，>=，==，!=：比较运算符 !，&，&&，|，||，xor()： 逻辑运算符 logical： 生成逻辑向量 all， any：逻辑向量都为真或存在真 ifelse()：二者择一 match， %in%：查找 unique：找出互不相同的元素 which：找到真值下标集合 duplicated：找到重复元素 optimize，uniroot，polyroot：一维优化与求根 if，else， ifelse， switch： 分支 for，while，repeat，break，next： 循环 apply，lapply，sapply，tapply，sweep：替代循环的函数。 function：函数定义 source：调用文件 ’ call：函数调用 . C，.Fortran：调用C或者Fortran子程序的动态链接库。 Recall：递归调用 browser，debug，trace，traceback：程序调试 options：指定系统参数 missing：判断虚参是否有对应实参 nargs：参数个数 stop：终止函数执行 on.exit：指定退出时执行 eval，expression：表达式计算 system.time：表达式计算计时 invisible：使变量不显示 menu：选择菜单（字符列表菜单） 其它与函数有关的还有： delay， delete.response， deparse， do.call， dput， environment ， formals， format.info， interactive， is.finite， is.function， is.language， is.recursive ， match.arg， match.call， match.fun， model.extract， name， parse 函数能将字符串转换为表达式expression deparse 将表达式expression转换为字符串 eval 函数能对表达式求解 substitute， sys.parent ， warning， machine cat，print：显示对象 sink：输出转向到指定文件 dump，save，dput，write：输出对象 scan，read.table，readlines, load，dget：读入 ls，objects：显示对象列表 rm, remove：删除对象 q，quit：退出系统 .First，.Last：初始运行函数与退出运行函数。 options：系统选项 ?，help，help.start，apropos：帮助功能 data：列出数据集 head()查看数据的头几行 tail()查看数据的最后几行 每一种分布有四个函数： d―density（密度函数），p―分布函数，q―分位数函数，r―随机数函数。 比如，正态分布的这四个函数为dnorm，pnorm，qnorm，rnorm。下面我们列出各分布后缀，前面加前缀d、p、q或r就构成函数名： norm：正态， t：t分布， f：F分布， chisq：卡方（包括非中心） unif：均匀， exp：指数， weibull：威布尔， gamma：伽玛， beta：贝塔 lnorm：对数正态， logis：逻辑分布， cauchy：柯西， binom：二项分布， geom：几何分布， hyper：超几何， nbinom：负二项， pois：泊松 signrank：符号秩， wilcox：秩和， tukey：学生化极差 sum, mean, var, sd, min, max, range, median, IQR（四分位间距）等为统计量， sort，order，rank与排序有关， 其它还有ave，fivenum，mad，quantile，stem等。 R中已实现的有chisq.test，prop.test，t.test。 cor，cov.wt，var：协方差阵及相关阵计算 biplot，biplot.princomp：多元数据biplot图 cancor：典则相关 princomp：主成分分析 hclust：谱系聚类 kmeans：k-均值聚类 cmdscale：经典多维标度 其它有dist，mahalanobis，cov.rob。 ts：时间序列对象 diff：计算差分 time：时间序列的采样时间 window：时间窗 lm，glm，aov：线性模型、广义线性模型、方差分析 quo()等价于quote() enquo()等价于substitute(）案例：水稻播种机穴盘育秧播种装置
1.水稻播种机穴盘育秧播种装置的试验设计 随着栽培技术的不断更新，高效、节本、高产的抛秧栽培法获得了迅速发展和推广。为了改善原有播种装置中窝眼辊轮结构，我们研制成功了穴盘育秧播种装置，它不仅解决了手工操作进行育秧培育的劳动强度大，工作效率低等问题，而且能大幅度地提高播种量的稳定性和播种的均匀性，使水稻播种机械更趋实用与完善。
（1）试验目的 考虑影响播种性能的主要因素对水稻播种机穴盘育秧播种装置播种性能的影响程度，以达到优化设计参数。
（2）试验条件 种子品种：杂交稻（协优46号）
种子状况：经过脱芒、浸种、催芽露白、去杂质
秧盘规格：600mm×340mm，561穴
种子千粒重：26.9g
试验盘数：100盘
秧盘运行速度与排种胶带线速度严格一致。
（3）试验因素
选用三个可变因素：
生产率（盘/小时）、播种量（粒/穴）、投种高度（mm）。
A.可变因素
B.可变的水平数每个因素分别取三个水平数
C.实验因素与水平
为了研究生产率、播种量及投种高度对播种性能的影响，特安排了三因素三水平的正交试验，试验因素与水平见下表所示。
2.正交试验方案与试验结果分析 （1）正交试验方案与试验结果 选用L9(34)正交表进行试验设计，试验方案与试验结果见下表所示。其数据采集方法为：在每种工况（每个试验号）条件下进行随机抽样3盘测定，测定播种合格率时，每盘随机连片抽样100穴。最后，把3次测定的各项数据的平均值记入试验结果。
（2）试验结果分析如下表所示
注：
（1）T为因素试验结果之和，如T1 = 93.0 + 91.0 + 89.0 = 273.0。
（2）t为因素试验结果之和的均值，如。
（3）R为t值中的大数-小数。
（4）播种合格率：每盘随机测定的100穴，其中种子粒数合格的穴数所占的百分比（种子粒数合格范围为：杂交稻(1-3粒/穴，常规稻3-6粒/穴）。
（5）播种变异系数V
x——每盘播种量；
——平均盘播种量（g）；n——试验盘数，s——标准差。
（6）空穴率：每盘随机测定的100穴，其中空穴数所占的百分比。
由上面两表得出影响3项指标的主次因素和较优水平为：播种合格率C1A1B3；播种变异系数C1B3A1；空穴率C1B3A2。
考虑到水稻播种的实际需要，经综合分析，选取各试验因素的较优水平组合为：A1B3C1、A2B3C1、A1B3C1。因为在上述正交试验中未出现过A1B3C1以及A2B3C1，为此专门安排了单因素（生产量）三水平试验，试验结果见下表所示。
从上表可知，最佳组合为A2B3C1，播种合格率96.0%，播种变异系数1.9%，空穴率0.5%。
3.试验结论 （1）400盘/小时是该播种装置杂交稻播种的临界生产率，高出此值，则各项性能指标受重大影响。
（2）播种量越大，各项性能指标越好。
（3）投种高度对播种质量的影响十分显著，投种高度越低，播种质量越好。
