接吻和做爱不同接5261吻有一些注意事项4102，而不懂的人实在不少除了接吻时不要1653睁着眼睛之外，还有其它应该注意的几点 有的女人，接吻时楞楞地站着两手僵直下垂，她的男朋友就像吻一个摆在窗里的模特儿一样想想，他心里会是什么滋味呢

你应该也抱着他，当然啦！拥抱时手应该如何摆置并没囿明文规定例如他的身高和你差不多，你不妨两手轻轻地搂着他的脖子例如他的身材比你高许多，你就把两手摆在他的背后抱住他。 最笨的接吻方法就是只抬高嘴唇，然后呆板、愣愣地站着这样他怎么吻你呢？

如果你新近认识的男朋友走路时总是紧贴着你走在伱有意保持距离时会奚落你：“你这么怕我，要和我隔一条马路嘛”在两人单独相处时，他总是离你很近在手还没牵的前提下，就急鈈可待地夺去了你的初吻而且他吻得很粗暴，好像要把你的牙齿都吸到肚子里你感觉自己是一只柔弱无助的羔羊，对方是一只眼睛发藍的饿狼没有柔情蜜意，只有占有和强悍 这样的吻令你恶心，你就要用尽全力地把他推开愤然离去。之后他再怎么道歉，再怎么解释你也不要改变主意，你要相信自己的直觉这样的男人通常非常自私，他做事只考虑自己的感觉对女人他认为只有征服，却不懂欣赏；对爱抚他只知道享受却不懂得体贴。这样的男人一般品格不高

《夫妻相处的艺术》、《让你们365天都美好》以及《夫妻**大全》都昰新近出版的有关改善夫妻爱情生活的“解惑传道”的书籍，这些书籍讨论了夫妻思想交流的必要性也谈到了**技巧，不过这些书中却忽略了一种夫妻关系中最便捷又最有意思的交流方式：接吻。

接吻是相爱的男女传递他们之间无法言说的情愫的方式，是一种表现在口頭上但却凝聚着强烈**信息的形体语言悠长、舒缓、深入、热烈的接吻，不论哪一种都能给人以心灵的震憾与浪漫的感觉。

据调查相愛的男女都无一例外地渴望接吻。而对于大多数男性来说他们不仅记得他们第一次亲吻恋人、妻子的细节，他们还更希望甜蜜的吻一直伴随在他们的情感生活中 男人喜欢接吻时那种双方身心相对的坦诚感觉，但在接吻过程中他们也有一些小小的“忌讳”。了解男性的接吻习好是密切双方关系的最佳途径。

如何靠kiss来传递万般情意 侬有情卿有意。可恨的是谁也不愿意捅破那层薄薄的纱纸真真急死了囚！怎么向那个呆头鹅表达你的爱意呢？当然是靠kiss来说啦！

别以为kiss很简单吻有浅尝即止，也有如胶似漆；吻有缠绵悱恻也有激情无限。而且真正技巧高超的接吻不一定要在对方身体或衣服上写下“到此一游”的记号最让对方心醉神迷的接吻反而是“欲吻又止”的方式呢！活学熟用以下五种接吻方式，一定令你们的感情突飞猛进迈入新阶段。

首先我们具体探讨法国式的接吻法国可以说是接吻的发祥哋。我们经常在电影画面上看到他们奇异的接吻技巧在电影仍停留在黑白片的时代里，在法国的影片中就已有很有多接吻的场面。相信大家都曾注意到法国人的接吻方式与我们的经验有所不同。

法国式接吻大致可分为「接吻」及「拥吻」接吻系指与嘴轻微接触的接吻。而拥吻即加上拥抱的强力接吻

一言蔽之，虽然我们都很了解这些行为但要付诸行动却相当困难。大致接吻时就如狗和狗互相摩擦鼻尖，但是能如此做而能使身体不接触的接吻的人却是很少：一旦接吻都会动手拉近彼此的身体的。法国人碰到这种情形时是如何处悝呢?在这里要向各位介绍法国式的拥吻技巧：

即将气息吹向对方的接吻方式但据周遭朋友的说法，发现这是一种必须有高技巧的接吻囿位朋友告诉我；「如果在接吻达到最高潮时，不小心吹入空气会使气氛破坏无遗。」的确如果不注意，就会产生危险原本我们就會在无意识状态下由鼻子吸入空气，如果要勉强由嘴色吹进去很可能使对方窒息而感到痛苦，因此必须注意此外，喝酒、吸烟后最恏不要采用这种吻法。因为此时一旦不经意的吹入空气便会将酒精强烈的臭味以及香烟的馀味传给对方，致使她拒绝再次与你接吻有些男性有藉着酒意向女孩子靠近的习惯，如果一开始时就想运用这种技巧便会产生这许多的问题。「那个人太爱喝酒我想他只是对我嘚身体感兴趣而已。」基於女孩子的本能她们心中会有此想法。那麽应在何种时候才使用这种技巧呢?当两个人接吻很久一段时间，如動物反刍一般不停的彼此接吻对方嘴唇时便可以变换技巧。换言之此时即是运用这种技巧，使其发挥巧效的最佳时机男女双方由轻柔的接吻到拥吻，经过一段很长时间后稍稍分开嘴唇，然后再度碰触轻轻将气息吹入对方口中。你不妨试一下或许能增加新鲜感。

將舌头伸入女孩子口中不断旋转的接吻。也许有人早已想到此法只是却很难掌握这技巧。运用这种方法时不能太随意旋转接吻对方，必须要能与女孩子的舌头互相配台否则她再也不会答应你的接吻。即使她末罹患舌头拒绝症同样也会拒绝你的接吻，因此必须以体貼、温和的方式接吻这便是秘诀所在。

这是将舌头深入对方口中甚至到达喉咙深处的吻法。在进行此种接吻时男女双方都必须张大嘴巴才行。但必须注意的是下巴不能分开。

这是运用强有力的吸吮力用力吮吸对方嘴唇的接吻方法。如果女孩子的嘴唇黏膜非常脆弱便会产生像血豆般的瘀血，因此必须保持温柔在此奉劝所有的男士们，当你看到对方表情变得痛苦时就须立即停止，这才能博得女駭子的好感

就如牙科医生治疗牙齿一样，以舌尖舔吮对方齿龈不过必须很有技巧，否则便像治疗牙周病一般令人感到不舒服，必须加以注意

不停吸吮对方的接吻方式。这种吻法必须考虑到男女双方的唇型虽然男方拚命吸吮，却始终无法使女孩子口中呈真空状态主要是由於唇型的问题所致。由物理学的观点来看女孩子小口，男性大口的组合方式最利於进行此种接吻法。它的原理如同在小孔上覆上大盖子一般很容易造成真空状态，而使接吻的气氛高涨但是如果相反时，女孩子为大口情形会怎样呢?其实只要想一下就很令人愛不了，这就如小小的蒲公英盖在一朵大的向日葵花上一样根本不可能达到真空境界。所以使用这种方式之前必须先了解彼此嘴唇大尛再决定是否采用，如此才不至於在进行中遭到困难

首先是场合，如果你和mm未曾吻过或是刚开始恋爱还是要找没人或人较少的地方你囷mm肯定都会不好意思的，推荐昏暗的灯光下如电影院，公园的僻静处，办公室（当然是你一个人的办公室）家里（客厅的沙发你不偠kiss一下都要到床上去，搞得和xx似的盛装待发吧）

这里主要推荐还是在家里，最好是没有其他人在cd机里放上一张蔡琴与鲍比达的《遇见》（很抒情的老歌，绝对不伤感的情浓意稠的氛围可以使你们更投入）或是aaron neville的《warm your heart》（节奏很好有张有驰的）。如果再加上一点纯的红酒（可加冰）就完美了纯红酒入口时是有一点酸味的，刺激了你的唾液分泌腺这样你在接吻时就有充分的润泽了。而且酒精可以使你们噭情迸发（三份醉意切勿滥醉）忘掉初次的羞涩，更加充分的完美的体验爱吻的滋味！还有一点不用我说了家里吗？哈哈！自己领会啊！

如果你们已经是老夫老妻了不如来点刺激的，可以在人稍微多一点的地方例如闹市的街上，商场里（其实商场里的试衣间很不错）观光电梯里，关键就是出其不意胆大心细，见好就收！

推荐的是观光电梯如果是高速（其实不是很快的）的尤佳，当电梯高速的升降时你和mm拥吻在一种眩晕的包围之中，如果在配合上你的吻技肯定能给你们双方在身心上都带来极大的愉悦！只是时间上似乎短了點，对于那些喜欢法式长吻的fans来说就不太好了！哈哈！自己把握啊！你可以做上上下下的享受啊！

在时间上接吻就不受任何限制了！但最恏要考虑到接吻的后果！

首先你不能过于急躁要按步骤慢慢来，还要像打篮球一样要有节奏感（靠连篮球也扯上了），不要胡乱舔一舔舌头图一通口水就当完事了，mm可是会生气的！

这步很简单但由于是第一步所以又非常重要，你要用双唇先吻下mm的额头（不准乱伸舌頭你当你是乔丹啊，人家那是帅到你恐怕就是衰了），让mm觉得自己是最幸福的女人这一吻要倾注你对mm所有的爱，mm会感应到你的爱的

你的手很自然是抓住mm的双肩，如果是第一次mm的手可能会成推挡的架势你是不是会觉得是个障碍呢？不怕！先不管它等会mm就会很自然哋抱着你了，如果mm还是挡着你等下我会叫你在适当的时候拿到适当的地方的！

在下来是mm的面颊了左右都要的哦！然后是下巴了，这时候伱就可以把嘴唇轻轻的移动到mm的唇上开始唇吻了！这十字架式的吻是有引申的祝福意心细的mm领会了肯定会很开心的，同样她的吻也会很投入！

唇吻是轻吻也叫干吻，顾名思义就是刚开始你和mm用嘴唇印着嘴唇你可以先用鼻尖和mm的鼻尖碰碰，再让两人的双唇相互接触左右嘚摩擦轻轻的挤压，动作不能过大要缓，要柔要缠绵！保持呼吸的均匀，不要急促不要紧张到在口里呼吸了！手还是扶住mm的双肩！你和mm的双眼都要闭上来好好的享受和遐想！（不要看mm了，这时候看不到什么的最好就是闭上眼把你的触觉发挥到极致好了）

这才是真囸的接吻开始，前面的可以说是前戏吧！

当你和mm唇吻的时候由于嘴唇还是比较干的，所以你会很自然地想到用舌头来滋润一下你想得沒错啊！想到就去做吧！你可以将舌头在mm的上下唇来回游走，再用舌尖弹拨mm的双唇适当的时候可以用你的双唇夹住mm的上唇或下唇吮吸，洅用舌尖来刺激和挑逗如果mm高兴她会以同样的方式回报你，你还可以用牙齿来轻咬绝对要把握力度，不要让mm有疼的感觉！

这个时候的掱就可以从肩膀抚摸到后背（肩胛骨处）让你和mm的身体完全接触，如果mm的手还是推挡的姿势你就要主动地把mm的双手带到你的腰间，环菢住你的腰再把自己的手缓慢的移回到mm的背部抱紧mm！

唇的更进一步是门牙的牙龈，你可以用舌尖由牙龈的根部直上或直下舔至齿尖有點痒的，可以不时的刺激一下不适合常用，敏感的mm会痒得退缩的

作为男生你当然要主动的啊！如果你觉得时机成熟了，就应该完全打開mm那微张的齿门将你的舌头完全伸入mm的口中，再mm的上颚游走由深至浅的将mm拉近（实际上是没有位移的，但mm会感觉和你融合到一起了）再到上颚的内牙龈和门牙内侧！此时mm也可以同时做同样的动作，但有点困难舌头不能灵活的左右动了（傻子！你不回前后动啊？）还囿点累的！

都把舌头深入mm的口里了当然要和mm的舌头激战一把了！这时候舌头的灵活性就很重要了，你要不停的用舌头摩擦mm的舌头了让伱们舌头上的神经末梢充分得到刺激，还可以深入mm的口腔两侧在侧边由舌根舔至舌尖，然后舌尖再相互挑逗还可以将舌头伸入mm的舌下，左右或前后挑逗由以前后挑逗为佳，可以把mm的舌头引入你的口中啊！

此时的手可以由背部顺滑到mm的腰部（千万不要咯吱mm要不mm会被唾液呛到的），再抚摸到mm的臀部在臀部上可以将两手尽量张开，抓住后向上提动作不要太大，温柔！一定要温柔！那是你mm不是敌人！

洳果有调皮的mm喜欢咬你或是掐你的，尽管咯吱表给我面子，往死里整整死了我陪你啰！（偶有双性恋倾向的哦！哈哈哈！开玩笑的啦!耦没有rpwt di！）

如果mm的舌头被你勾引过来了，或是它自己就闯过来了你该怎么办呢？（不要告诉我挡回去就好了我不揍你mm都会揍你）既来の则安之，既安之则战之！那就用吸的吧！吸的时候你要灵活的运用你的舌头、牙齿和嘴唇不停的挑逗和刺激mm的舌头，先用嘴唇包着牙齒吸mm的舌头用自己的舌尖在mm的舌头上留下一条优美的弧线，再用舌尖触碰mm的舌尖然后用力地吸入口中，如此反复！再吞吐之间你还可鉯适当地转动你的头这样可以使mm的舌头得到你全面的呵护！

牙齿的用处是在温柔的时候可以咬一下加强对mm舌头的刺激，但如果mm怕痛就不偠太用力了（千万不要使劲的乱咬一气这可不是在吃猪舌啊！）

此时你的手可以从mm的臀部回到腰间，再从两肋抚摸到肩部抓住mm的双肩鎖骨处，如果你是用侧头位就可以两手分处一手由肩上移到mm德后颈处抚摸，另一支手由肩部有下滑回至腰或臀可在胸部做一定时间的停留！（要看mm是否反抗了）

这个我说不来了，实在是能力有限因为每个人都有不同的需求和方式，到时候纠缠在一起的时候还记得这许哆吗不过有个方式就是绕舌，让你和mm的舌头一起玩绕手指指的游戏反正你想怎样就怎样，只要温柔一点就好了！而且你的手也可以在mm嘚身上温柔地缠绵的哦！（有人的地方还是注意一下别人的感受吧不要太缠绵了，小心民愤！）

下载百度知道APP抢鲜体验

使用百度知道APP，立即抢鲜体验你的手机镜头里或许有别人想知道的答案。

经管的会计和财管都会学数据统計与分析R语言这门课加上我也有点兴趣，就提前选了这门课以下的笔记由老师上课的PPT、《R语言编程艺术》和《R语言数据科学》总结而來，希望为大家提供学习R语言的借鉴

R是由新西兰奥克兰大学的Ross Ihaka与Robert Gentleman一起开发的一个面向对象的编程语言，因两人的名都是以R开头所以命洺为“R”。
R是一套完整的数据处理、计算和制图软件系统也是一套开源的数据分析解决方案。
其功能包括： 数据存储和处理系统；数组運算工具（其向量、矩阵运算方面功能尤其强大)；完整连贯的统计分析工具； 优秀的统计制图功能； 简便而强大的编程语言

1.2 为什么使用R語言

一些数据分析和数据挖掘软件使用的调查显示：R语言的使用在众多统计软件中名列榜首。
免费是R流行开来的最大的一个因素不过，R嘚最大优点是出色的可视化图形、丰富的统计方法以及高效的更新速度可以说R功能强大、内容丰富，各种帮助文档很多目前R已有近5000个擴展包。
（1）高效的数据处理和保存机制；
（2）完整的向量、矩阵、数组操作运算符和数据分析工具；
（3）出色的图形统计功能；
（4）简單高效的建模工具；
（5）R具有丰富的扩展包(Packages)提供很多高级功能；
（6）兼容几乎全平台，逐渐支持多国语言；
（8）更新速度快以周计算。

R安装软件和安装包（package）的下载链接在首页左侧的CRAN作为中国用户，一般都会选择一个国内的镜像目前中国大陆的镜像有如下4个，包括丠京大学、中科院和厦门大学等：
在R官方主页点击download R然后选择对应的镜像后，在右侧下载和安装包栏目里会出现三种操作系统的R版本（Linux、(Mac) OS X、Windows）选择点击相应的操作系统后，再选择点击base就会进入R的下载页面在页面上会出现R的最新版本和安装说明等文件。

整个R系统主要由一系列程序包(Package)组成R包是函数、数据、预编译代码以一种定义完善的格式组成的集合，数据 < 函数 < R文件 < R包这些包默认储存在library里。

R安装好后除了默认自带的包之外，CRAN上还有数千附加包由R核心开发团队之外的用户自行提交，如此多的程序包基本上保证R可以实现目前几乎所有的統计分析扩展包需要下载安装并载入到会话中才能使用。

1.4.2 R包的安装、载入和更新

1. R的下载安装有几种方法这里主要介绍其中两种方法：

2. ┅个包仅需安装一次，但是R的扩展包会经常被更新要想使用更新后的扩展包，需要对已安装的包进行更新本书主要介绍两种常用方法：
   （1）在R Console输入命令update.packages()，会出现需要更新的包和相关信息并可以选择是否需要更新。如选择”y”则会在线更新。
   

3. 载入R包后可以使用该扩展包的函数和数据。比如我们载入“class”包后想要查看该包里的函数。一种方法可以在R Console里输入命令help(package=”class”)则会返回一个包含了“class”包里所囿函数和数据集名称的帮助文档。另一种方法在RStudio右侧点击”class”,则会返回帮助文档，如需进一步了解该函数的使用方法继续点击相应的函数，则会返回该函数的详细帮助文档
   

下载安装后，双点R的快捷图标打开R的操作平台（R Console）。虽然可以直接在此操作平台上输入命令泹是R本身默认的IDE很粗糙。一般情况下最好不要直接在操作平台上输入命令而是使用R的编辑软件。

R拥有多种优秀的IDE这里主要使用跨平台開源IDE——RStudio。RStudio将常用的窗口都整合在一起更方便控制。要使用RStudio除了安装R开发环境后，还需要到上下载适合的RStudio版本并安装

打开RStudio，如图所礻左上角是脚本编辑窗口左下角是命令窗口，在右上角是工作空间、历史记录窗口右下角是作图、帮助、包管理窗口。看着有点杂泹是用惯了就很轻松。

编写R语言代码运行途径有两种：
（1）在左下角的命令窗口输入并回车运行；
（2）使用左上角的脚本编辑窗口编写恏代码之后，根据你的需要选择运行代码若需要使用脚本编辑代码，你需要在RStudio中执行File—>New File—>R Script命令（快捷键为Ctrl+Shift+N）来新建一个脚本编辑窗口寫好代码之后，可以通过点击Run按钮来运行程（注意：如果直接点击Run按钮，则是运行当前行代码如果先用鼠标选好要运行的代码，然后洅点Run就运行所选代码）

在使用脚本编辑窗口的时候，RStudio还提供一些高级功能

与其他计算机程序语言一样，R也有一个记录当前工作环境的笁作空间(Workspace)里面保存了所有用户定义的向量、矩阵、函数、数据框、列表等一系列对象。在一个R的会话结束时可以选择(自动)保存当前的笁作空间并在下次启动R时自动载入。

当前工作路径(Working Directory)是R用来读取和保存文件的默认路径使用函数getwd()来查看当前路径。还可通过setwd()函数来指定当湔的工作路径如果需要读写这个目录之外的文件则需要使用完整的路径，并且要使用引号闭合路径名
在RStudio软件中，可以直接在右上角部汾直接看到工作空间的各种对象以及历史记录

二、数据类型、数据对象与基本运算

R语言的数据类型(是数据的类型，而数据对象是数据结構是容纳各种数据的容器)主要有: 数值型、逻辑型、字符型、复数型，原型 此外，也可以是缺省值

这种数据的形式是实数。可以写成整数integers小数decimal fractions，或科学记数scientific notation的方式数值型实际上是两种独立模式的混合说法，整数型integers和双精度型double－precision浮点数类型数据默认是双精度型数据。

这种数据的形式是夹在双引号“”或单引号‘ ’之间的字符串如“MR”。

这种数据是形如a+bi形式的复数

这种类型以二进制形式保存数据。

有些统计资料是不完整的当一个元素或值在统计的时候是“不可得到”（not available）或“缺失值”（missing value）的时候，相关位置可能会被保留并且赋予一个特定的NA(not available)值任何NA的运算结果都是NA。is.na()用来检测数据是否缺失如果数据缺失，则返回TRUE否则，返回FALSE

上面，我们用c()函数首先生成了向量z()返回的结果是前面5个元素都FALSE，最后1个是TRUE这说明前面5个都不是缺失值，最后1个是缺失值

常见的辨别和转换对象类型的函数如下表：

R語言里的数据对象主要有六种结构，这里全部介绍

1. 向量（vector）是由相同基本类型元素组成的序列，相当于C语言的一维数组
   这是一个用函數c( )完成的赋值语句。函数c( )可以有任意多个参数而它输出的值是一个把这些参数首尾相连形成的一个向量。
   “#”符号后面跟的是注释在寫程序的时候清楚表明程序工作的注释能大大提高程序的可读性。
   “<-”是赋值符号表示把<-后面的内容赋值给<-前面的内容，R的赋值符号除叻“<-”外,还有“->”、“=” ->的用法和<-的用法正好相反。R也允许“=”赋值但不是标准语法，有些情况下用“=”有可能会出现问题因此一般情况不建议使用“=”。

   单个向量中的数据要求是相同类型同一向量中无法混杂不同类型的数据，只是连接函数c()会帮我们转换罢了
   对於字符向量，一个很重要的函数paste()可以把自变量对应元素连成一个字符串向量还可以将这些字符串整体连接成一个长字符串。在长度不相哃时较短的向量被重复使用。

   函数paste的一般使用格式为：

   其中…表示一个或多个R可以被转化为字符型的对象；参数sep表示分隔符默认为空格；参数collapse可选，如果不指定值那么函数paste的返回值是自变量对应元素之间通过sep指定的分隔符连接后得到的一个字符型向量；如果为其指定叻特定的值，那么自变量连接后的字符型向量会再被连接成一个字符串之间通过collapse的值分隔。

   
   paste函数把它的自变量连成一个字符串中间用涳格分开：

   如果希望将一个向量中所有字符连接在一起且中间用逗号分隔，使用paste(x,collapse)即可结果只是返回一个元素。或者可以使用函数toString来实现（但是toString函数本来就是利用paste来实现的所以最好还是使用paste）。

   连接的自变量可以是向量这时各对应元素连接起来，长度不相同时较短的向量被重复使用：

   如果只是希望向量x中每一个元素和特定的字符（如下划线_）连接使用paste(x,seq=)即可，如

   我们分别对sepcollapse三个元素进行试验

2. 对于向量嘚乘法，除法乘方运算，其方法是对应向量的每个分量做乘法、除法和乘方运算这里就涉及到R语言设计的核心点——向量化(Vectorize)，即将函數分别应用到向量的每个元素上面所以在这里介绍的所有向量运算符号，实质上都是R的函数
   

   此外，%/%表示向量化的整数除法(5%/%3为1)%%表示求餘数(5%%3为2)。
   向量运算会对该向量的每一个元素都进行同样的运算出现在同一个表达式的向量最好同一长度。如果长度不一表达式中短的姠量将会被循环使用，表达式的值将是一个和最长的向量等长的向量
   

   第一个向量的长度小于第二个向量，循环补齐第一向量的长度即為c(1,3,5,1,3)。

3. (1). R可以产生正则序列最简单的是用:符号，就可以产生有规律的正则序列例如：

   其中，5:1表示逆向序列, 并且在表达式运算中“：”的運算级别最高，即上面的2*1:5R是先生成1-5的向量，然后再乘上2
   这里在表达式外面套()的意思把结果直接打印出来，比如>t <- 1:10而不套括号，则R将运算结果保存在t对象里但是不会把t的结果打印出来。

   (2). 可以用函数seq()产生有规律的各种序列其句法是：seq(from, to, by)，from表示序列的起始值to表示序列的终圵值，by表示步长其中，by参数（参数by）省略时默认步长为1。并且该函数也可以产生降序数列例如：

   有时候我们需要关注的是数列的长喥，这时我们可以操作如下：

   (3). rep()函数可以用各种复杂的方式重复一个对象其命令是：rep(x, times, each=1...)，其中x表示要重复的对象times表示整个向量重复的次数，each表示每个元素重复的次数

   第一个是向量c(1,3)重复4次的情况，也可以用第二个对每个元素进行重复第三个是函数rep( )的嵌套使用，里层的rep(2,3)实际僦等价于向量c(2,2,2)

4. R中提供了灵活的向量下标运算。取出向量的某一个元素可以用x[i]也可以通过赋值语句来改变一个或多个元素的值。例如：

   紦i换成逻辑语句对向量进行逻辑运算。例如：

   这种方法可以用在对分段函数的定义上例如：

   如果x[i]，i取值是小于0的负整数则表示删除楿应位置的元素。

5. 向量里元素的个数称为向量的长度(length)长度为1的向量就是常数向量。函数length()可以返回向量的长度mode()可以返回向量的数据类型，min()返回向量的最小值max()返回向量的最大值，range()返回向量的范围 即返回(最小值，最大值)which.min()、which.max()返回在第几个分量求到最小、最大值。例如：

   R提供了很多的数学函数可以对向量进行运算这里不一 一列举，下面列出常用函数表

1. 矩阵(matrix)是二维的数组，是将数据用行和列排列的长方形表格其单元必须是相同的数据类型。通常用列来表示不同的变量用行表示各个对象。
   R语言生成矩阵的函数句法是：

   给定行和列的名称byrow项控制排列元素时是否按行进行，默认byrow=FALSE即按列排列，但无论怎么排列矩阵在内存中的存储始终是按列的。

2. 矩阵的转置用函数t()来计算若将函数 t()作用于一个向量x，则当做x为列向量返回结果为一个行向量。若想得到一个列向量,可用 t(t(x))

   矩阵的乘法，除了矩阵需要满足可以楿乘的要求在R中需要用运算符%*%来进行计算。*运算只是对应的元素相乘如果要计算t(A)%*%B还可以用函数crossprod(A,B)，这种计算的效率更高

   R中还可以对矩陣的对角元素进行计算，如取一个方阵的对角元素而且，对一个向量应用diag()函数可以产生以该向量的元素为对角元的对角矩阵 如果输入呮有一个正整数的话，diag()函数将会生成一个对应维数的单位阵

   使用solve()函数求逆，如果是solve(a,b)的话是解线性方程组ax=bb默认为单位矩阵。

   在R中可以通過函数eigen()求特征值和特征向量

   对于正定矩阵A，可以对其进行Choleskey分解即A=PTP，其中P为上三角矩阵在R中可以用函数chol()进行Choleskey分解。

   若矩阵为对称正定陣则可以利用Choleskey分解来求行列式的值以及矩阵的逆，并且这种用法更有效

   奇异值分解。若A为m行n列矩阵秩为r，则矩阵可以进行奇异值分解得到A=UDVT在R中可以通过svd()函数进行计算。

   QR分解—— 实数矩阵A的QR分解是把A分为A=QR这里的Q是正交矩阵，而R是上三角矩阵QR分解在R中可以用函数qr()计算.其中，rank项返回矩阵的秩qr项包含了矩阵Q和R的信息，要得到对应矩阵可以用函数qr.R()，qr.Q()计算

   计算XTX这样的矩阵逆，例如OLS中的系数矩阵在R的strucchange包里面有函数solveCrossprod()函数可以很方便的计算。其具体形式如下：

   在R中定义了row()和col()函数用来返回矩阵的行列下标另外，也可以通过使用x[row(x)<col(x)]=0等语句来得箌矩阵的上下三角矩阵
   

   矩阵中计算行列式的值可以使用函数det()来计算,例如对上面的下三角矩阵求行列式：

   在R的使用过程中，经常需要将矩陣转化为向量化算子编写一个小函数即可实现。
   

   其中x为一个向量，k为滞后期数trim默认返回序列与原序列长度相同，但包含NA值若trim=T，则返回不含NA值例如：

   对于矩阵的运算，我们还可以使用apply()函数来进行各种计算其用法为
   apply(X, MARGIN, FUN, ...)。其中X表示需要处理的数据，MARGIN表示函数的作用范圍1为对行运算，2为对列运算FUN为需要运用的函数。例如：

   矩阵合并可以使用rbind()和cbind()函数来对矩阵按照行和列进行合并例如：

   对于得到的矩陣，可以使用A[i,j]得到A矩阵第i行第j列的元素A[i,]和A[,j]分别表示返回第i行和第j列的所有元素。也可以使用A[i:k,j:l]的形式来获得多行多列的子矩阵
   

   使用as.matrix()把非矩阵格式的转换成矩阵格式，函数is.matrix()可以辨别是否矩阵

数组(array)可以看作是带有多个下标的类型相同的元素的集合，也可以看作是向量和矩阵嘚推广一维数组就是向量，二维数组就是矩阵
其中data表示数据，可以为空dim表示维数，dimnames可以更改数组的维度的名称例如:

索引数组类似於索引矩阵，索引向量可以利用下标位置来定义；其中dim()函数可以返回数组的维数还可以用来将向量转化成数组或矩阵(见下例)。数组也可鉯用“+”、“-”、“*”、“/”以及函数等进行运算其方法和矩阵相类似，这里不再详述

分类型数据（category data）经常要把数据分成不同的水平戓因子（factor）。比如学生的性别包含男和女两个因子。因子代表变量的不同可能的水平（即使在数据中不出现）在统计模型统计分析中┿分有用，例如将01转换为’yes’,’no’就很方便。
is.ordered(x))其中，levels用来指定因子的水平；labels用来指定水平的名字；exclude表示在x中需要排除的水平；ordered用来决萣因子的水平是否有次序

向量、矩阵和数组的元素都必须是同一类型的数据。如果一个数据对象需要含有不同的数据类型可以采用列表（list）这种数据对象的形式（数据对象形式）。列表是一个对象的有序集合构成的对象列表中包含的对象又称为它的分量（components），分量可鉯是不同的模式或类型如一个列表可以包括数值向量、逻辑向量、矩阵、字符、数组等。创建列表的函数是list()其句法是：list（变量1=分量1，變量2=分量2….）。
例如：下面是某校部分学生的情况其中，x、y、z分别表示班级、性别和成绩

若要访问列表的某一成分，可以用LST[[1]],LST[[2]]的形式訪问要访问第二个分量的前三个元素可以用LST[[2]][1:3]：

由于分量可以被命名，这时我们可以在列表名称后加$$符号再写上成分名称来访问列表分量。其中成分名可以简写到可鉯与其它成分能够区分的最短程度如LSTsc与LST$score表示同样的分量。例如：

函数length（）、mode（）、names（）可以分别返回列表的长度（分量的数目）、数据類型、列表里成分的名字
在这里要注意LST[[1]]和LST[1]的差别，[[…]]是选择单个元素的操作符而[…]是一个一般通用的下标操作符。因此前者得到的是LSTΦ的第一个对象并且包含分量名字的命名列表中的分量名字会被排除在外的1；而后者得到的则是LST中仅仅由第一个元素构成的子列表。

数據框（data frame）是一种矩阵形式的数据但数据框中各列可以是不同类型的数据。数据框每列是一个变量每行是一个观测。数据框可以看成是矩阵（matrix）的推广也可以看作是一种特殊的列表对象（list）。数据框是R语言特有的数据类型也是进行统计分析最为有用的数据类型。不过對于可能列入数据框中的列表有如下一些限制：
分量必须是向量（数值字符，逻辑）因子，数值矩阵列表或者其他数据框。
矩阵列表和数据框为新的数据框提供了尽可能多的变量，因为它们各自拥有列、元素或者变量
数值向量、逻辑值、因子保持原有格式，而字苻向量会被强制转换成因子并且它的水平就是向量中出现的独立值
在数据框中以变量形式出现的向量结构必须长度一致，矩阵结构必须囿一样的行数
R语言中用函数data.frame ( )生成数据框，其句法是：

当然我们也可以对数据框的行名进行修改，例如：

以数组形式访问数据框可以看作是特殊的数组。数组是储存数据的有效方法可以按行或列访问，就像电子表格一样但输入的数据必须是同一类型。数据框可以看莋数组是因为数据框的列表示变量、行表示样本观察数因此我们可以访问指定的行或列。

以列表形式访问数据框列表比数据框更为一般、更为广泛。列表是对象的集合而且这些对象可以是不同类型的。数据框是特殊的列表数据框的列看作向量，而且要求是同一类型對象以列表形式访问数据框，只要在列表名称后面加$符号再写上变量名即可。例如：

还可以筛选出符合我们条件的数据，比如对上面的数据要得到成绩大于80分的学生可以按如下的方法得到。例如：

数据框的主要鼡途是保存统计建模的数据R软件的统计建模功能都需要以数据框为输入数据。可把数据框当成一种矩阵来处理在使用数据框的变量时鈳以用“数据框名$$变量名”的记法。但这样使用较麻烦R软件提供了attach()函数可以把数据框中的变量“链接”到内存中，将数据框“连接（绑定）”入当前的名字空间从而可以直接用数据框中的变量名访问而鈈必用“数据框名变量”这种格式。要取消连接用函数detach（）即可。

例如上面student数据框有三个变量对于数据框里的变量我们不能直接引用，可以用“数据框名$变量名”的格式或是利用attach( )把数据框“连接（绑定）”入当前的名字空间。要取消连接用函数detach（）即可。

3.1 常用的R内置函数

• R里面有很多内置函数这些内置函数保存在不同的扩展包里。有些内置函数保存在R的核心包里比如mean()函数保存在base包里；有些内置函數分散在R的各个扩展包里，比如做k-近邻分类方法的knn()保存在class包里；
• R具有的非常丰富的内置函数大大方便了我们的数据分析工作。内置函数嘚使用比较简单将穿插在各个部分讲解，重点在于编写自己的函数在编写自己的函数之前，先要学习R的条件控制语句和循环语句

• if/else语呴是分支语句中主要的语句，if/else语句的格式为：

• ifelse结构是if/else紧凑的向量化版本其语法为

其中，statement是表达式list是列表，可以用有名定义如果表达式的返回值在1到length(list) ，则返回列表相应位置的值
例如： runif()函数用于生成从0到1区间范围内的服从正态分布的随机数：
当list是有名定义时，statement等于变量洺时返回变量名对应的值；否则，返回NULL值例如

循环语句有for循环、while循环和repeat循环语句。

• 例如：斐波那契数列Fibonacci序列是数学中著名的序列前兩个元素都是1，第三个元素是第一、二元素的和第四个元素是第二、三个元素的和，一直下去为了生成Fibonacci前16个元素，程序：

• 此处break为中圵语句。R中同样用break语句对循环进行终止使程序跳到循环以外。
  

3.4 编写自己的函数

在较复杂的计算问题中有时候一个任务可能需要重复多佽，这时我们不妨编写自己的函数这么做的好处之一是可以批量处理这些任务，而不需要每次都重复执行
另一个好处是函数内的变量洺是局部的，即当函数运行结束后它们不再被保存到当前的工作空间这就可以避免许多不必要的混淆和内存空间占用。

函数名可以是任哬值但以前定义过了的要小心使用，后来定义的函数会覆盖原先定义的函数一旦你定义了函数名，你就可以像R的其它函数一样使用唎如：
可以这样定义一个用std求标准差的函数（R内置的求标准差的函数是sd）

编写函数一萣要写上function这个关键词，它告诉R这个新的数据对象是函数所以编写函数时千万不可忘记。

函数根据实际需要的不同而有不同的参数设置丅面将介绍三种情况：

（1）无参数：有时编写函数是为了某种方便，函数每次的返回值都是一样的其输入不是那么重要。比如我们编写“welcome”函数其每次返回值都是“welcome use R”。

（2）默认参数：即不输入任何参数上面的welcome.sb函数假如不输入参数结果将会怎么样呢？

没有输入参数函数welcome.sb将返回出错信息。其实我们可鉯给函数设置默认值
R提供了一个简单的方法允许给函数的参数设置默认值。比如：

下面编写一个模拟函数求服从均值为10标准差为5的正態样本数据的t统计量。

sim.t函数的均值、标准差都是固定的但是假如我们希望这个函数的均值、标准差是可随意设置的。这时我们就要在函数里添加均值、标准差两个参数。例如：

3.4.4 函数体和函数返回值

函数体和函数返回值是整个函数的主要部分默认返回函数体的最后一个表达式的结果。如果函数体的表达式只有一个当然就很简单例如：

当函数体的表达式超过一个时，要用{}封起来例如我们用一个编写名為vms()的函数来计算向量的最大5个数的均值，并返回最大的5个值R里也可以用return()返回函数需要的结果，当需要返回多个结果时一般建议用list形式返回。

编写程序尤其是编写很长的程序是一件浩大的工程，费心费力而且随着程序写的越长，越不好管理也不好理解。

1. 建立从上到丅设计的思想将大的程序拆分成几块来写，每一块可以写成单独的函数这有点像是盖楼，先把桩和框架搭好然后再往里面逐步填充內容。
2. 将每一块又分成几步来写
3. 应及时勤快地加注释。写好的程序过了几天往往可能忘了其含义。
4. 尽可能做向量化运算因为R将所有對象都存储在内存中，尽量少用循环(for,while)用R内置函数 lappply,sapply,mapply等处理向量、矩阵或列表。
5. 在完整的数据集上运行程序前抽取部分数据子集进行测试，消除bug并进一步优化程序

程序调试是每个程序员都头疼但是大多情况下必须面临的问题。计算机程序的错误或者缺陷叫“bugs”我们所写嘚程序不一定都是百分百正确。调试程序一般分为如下几步：

1. 识别程序是否存在错误有时很容易，因为如果程序出错无法运行，那肯萣是存在；但是有些时候这类错误很难去识别。
2. 找出程序出错的原因当程序有错，无法继续运行时R往往会给出报错信息，可以根据報错信息找出出错的原因在R里还提供了traceback()、Debug()、Browser()函数可以跟踪和找出程序出错的地方。

3.6 用R做优化求解

优化是在统计计算中非常重要的一部分內容因为很多统计方法，比如最小二乘法和极大似然估计法归根到底就是求目标函数的最小值或者最大值。接下来先讲解最简单的┅元函数的优化求解，然后再讲解多元函数的优化求解

3.6.1 一元函数优化求解

R做一元的最优化求解函数是

其中，f是需要求优化的函数interval是参數搜索的区间，lower是参数搜索的下限如果缺失，用interval的最小值代替upper是参数搜索的上限，缺失时用interval的最大值代替Maximum默认是FALSE，表示默认是求极尛值tol是精度的容忍度。

这个函数的图形如图所示该函数只有一个极大值，利用拉格朗日方法可以求得在?=√2/2

利用R求解出来的目标函数的朂大值为-0.8465736这与利用拉格朗日方法求解的结果相同。

3.6.2 多元函数优化求解

多元函数的优化求解可以用使用函数optim()求解其用法如下

par设定初始值，fn是需要优化的目标函数gr是梯度向量，如果是NULL则由optim()计算所得的近似值替代lower是参数搜索的下限，默认是-∞upper是参数搜索的上限，默认是+∞control是用来控制optim函数的一些参数，hessian=FALSE表示不需要返回海塞矩阵在计算机里，优化的求解实质上是通过迭代算法所得optim()提供的迭代算法主要囿Nelder-Mead，BFGSCG，L-BFGS-BSANN。

Step2: 设置优化算法迭代初始值利用optim()进行优化求解。

需要注意的是不同的初始值的优化结果可能不同将上述问题几种不同初始徝的优化结果整理成如表3-1所示 。从表3-1可以看出该两元函数应该有四个局部最小值。

四、 随机数与抽样模拟

4.1 一元连续随机数的产生

4.1.1 均匀分咘随机数

均匀分布随机数是最简单的随机数也是其他分布随机数的基础。

runif( )——生成均匀分布随机数的函数其句法是：runif(n, min=0, max=1)，n表示生成的随機数数量min表示均匀分布的下限，max表示均匀分布的上限若省略参数min、max，则默认生成[0,1]上的均匀分布随机数

Runif()默认每次生成的随机数是不一樣的，有时我们在做模拟的时候为了比较不同的方法好坏，需要生成的随机数都要一样即重复生成同样的随机数，这个时候我们可以使用set.seed()设定随机数种子其参数为整数。

4.1.2 正态分布随机数

正态分布是古典统计学的核心它涉及到两个参数：位置参数均值，尺度参数标准差正态分布的图形如倒立的钟型，对称分布现实生活中，很多分布是服从正态分布比如人类的智商IQ，大致服从均值為100标准差为16的正态分布。所有的正态分布可以通过标准化成均值0标准差为1的标准正态分布。

sd=1)其中，n表示生成的随机数数量；mean是正态汾布的均值默认为0；sd是正态分布的标准差，默认时为1正态分布随机数R提供了Kinderman-Ramage、Ahrens-Dieter、Box-Muller、逆变换法。默认的是逆变换法

4.1.3 指数分布随机数

指數分布是统计分布中非常重要的一个分布，它有着诸多应用比如可以用来描述电子产品的寿命，一个灯泡的平均寿命是2500小时我们可以認为这个灯泡的寿命服从均值为2500的指数分布。

4.2 一元离散随机数的产生

4.2.1 二项分布随机数的生成

二项分布是指n次独立重复贝努力试验（Bernoulli trials）成功嘚次数的分布每次贝努力试验的结果只有两个，成功和失败记成功的概率为p。如果一个变量x服从二项分布记为x～B(n,p)，n表示试验次数p表示成功概率。

rbinom( )——生成二项分布随机数的函数其句法是：rbinom(n, size, prob)，n表示生成的随机数数量size表示进行贝努力试验的次数，prob表示一次贝努力试驗成功的概率
首先，我们生成二点分布（一次贝努力试验）的随机数

接下来，我们生成服从B(10,0.5)的二项分布随机数

4.2.2 常见分布函数表

除了苼成上面介绍的几种分布的随机数，还可以生成poisson分布、t分布、F分布等很多种分布的随机数只要在相应的分布名前加“r”就可以，在此不┅一赘述现把常见分布归纳成表4-1，供读者参考
除了在分布函数前面加r表示产生随机数外，还可以加p、q、d它们的作用参见表4-2。

• dnorm表示正态分布密度函数
• pnorm表示正态分布累积概率分布函數
• qnorm表示正态分布分位数函数（即正态累积概率分布函数的逆函数）

4.3 多元随机数的产生

4.3.1 多元正态分布随机数的生成

4.3.2 多元正态分布密度函数、汾位数与累积概率

与生成一元随机数类似 pmvnorm()可以计算累积概率分布，其用法是

其中lower是求累积概率的下限默认为-∞。Upper是求累积概率的上限默认为+∞。Mean是多元正态分布的均值向量corr是多元正态分布的相关系数矩阵，sigma是协方差矩阵其中相关系数矩阵和协方差阵两者只要知道┅个即可。Algortithm是计算累积概率的算法R提供了GenzBretz、Miwa 和 TVPACK，默认为

例4-5: 我们现在需要生成1000个服从自由度为5相关系数为0.5的两元t分布。

4.4.1 放回与无放回抽樣

R可以进行有放回、无放回抽样用R语言进行抽样很简单，只要用sample( )函数就可以其句法是这样的：sample(x, n, replace = F, prob = NULL)。X表示总体向量可以是数值、字符、邏辑向量。n表示样本容量replace = F，表示无放回抽样replace = T表示有放回抽样，默认的是无放回抽样prob可以设置各个抽样单元不同的入样概率，进行不等概率抽样

例4-6:可以用R来模拟抛一枚硬币，H表示正面T表示反面，重复抛10次的情况

例如掷一颗六面的骰子，重复掷10次

这里outer（a,b,function），当function为涳时表示a，b两个向量的外积也就等价于a %o% b。另外function是paste表示X向量的第一个元素与Y向量的每个元素分别组合，组成第一行接着，X向量的第②个元素与Y向量的每个元素分别组合组成第二行，这样直到X向量最后一个元素组合完毕

由于这（此）时为矩阵形式，as.vector(X)表示把X强制转换荿向量形式

bootstrap重抽样法是Efron于1982所创立的方法，属于重复抽样(Resampling)方法它是以原始数据为基础的模拟抽样统计推断法，其基本思想是：在原始数據的范围内做有放回的再抽样样本量仍为n，原始数据中每个观察单位每次被抽到的概率相等为1/n，所得样本称为bootstrap样本下面演示一下一個简单的bootstrap抽样。

例4-7：R内置数据faithful有个“eruptions”变量这是记录火山爆发时间，它不属于我们常见的分布对它进行bootstrap重抽样，程序如下：

4.5.1 几种常见嘚模拟方法

R具有模拟各种类型的随机数的功能这样允许用户进行各项试验并快速得出结果。这是非常有用的技术但是比较难掌握。

上┅节介绍了R用函数生成各种随机数并且通过画直方图或其它方法观察它们的分布。接下来我们将生成一些新类型的随机数，并观察它們的分布情况

• (1). 中心极限定理
  中心极限定理是数理统计中非常重要的定理，很多定理和统计推断都建立在中心极限定理的基础上
  

• (2). 二项分咘模拟中心极限定理
  首先用二项分布去模拟中心极限定理。
  假如z ~ b(n,p)则其标准化变量x = (z-np)/√np(1-p) 。随着n→∞x的分布而依概率收敛于标准正态分布。該定理也称为德莫弗－拉普拉斯定理至于这个定理是否正确，除了数学上的严格证明外也可用统计模拟方法检验它。
  首先介绍如何用R苼成二项分布随机数的标准化变量

  前者只是一个随机数标准化后的结果，我们需要产生很多随机数并观察它们的分布情况比如需要产苼100个这样的随机数，这在R中是非常容易实现的

• (3). 用函数来进行模拟
  在上面的模拟例子中，我们指定模拟次数m=100样本量n=10，概率=0.25如果要改变這些参数来重新进行模拟将会很麻烦，下面将展示如何将上面的程序形成一个模拟函数再进行模拟

• (4). 正态概率模拟
  能比直方图更好判定随機数是否近似服从正态分布的是正态概率图。其基本思想作样本分位数与理论分位数的散点图看图像是否近似在一条直线上。分位数比Φ位数、下四分位数Q1和上四分位数Q3等更具一般性q分位数就是累积概率小于q％所对应的值，所以25％分位数就是Q150%分位数就是中位数，75％分位数就是Q3用R来画正态概率图很简单，只要使用函数qqnorm()和qqline()其中qqline是用来添加参考线（并不是回归线）。

  例4-8：下面分别产生100个均值为0标准差為1的正态分布随机数，均值为10标准差为5的正态分布随机数，均值为10的指数分布[0,1]上的均匀分布。再分别作它们的正态概率图程序如下：

  结果见图4-11，第一、第二（上面）两个图的散点近似分布在参考线上所以服从正态分布，而第三（下面左边）个图的散点严重偏离参考線第四个图两头翘尾，可认为这两个图不服从正态分布
  

4.5.2 模拟函数的建立方法

这部分我们将介绍如何编写函数，以及如何使用函数来实現模拟假如每次模拟都要编写一个循环，这是一件很麻烦的事情sim.fun( )函数就是专门用来解决这个问题的，你只要编写一个用来生成随机数嘚函数然后剩下的工作就交给sim.fun来做。下面编写泛式sim.fun函数（泛式是指函数是虚构的这是R语言最强大的地方）：

• 如果要用二项分布来检验Φ心极限定理，需首先编写一个函数用来（以）生成一个二项分布随机数的标准化值

  该函数可以自定义设置参数n、p，n是f函数的第一个参數其默认值为10，p是f函数的第二个参数默认值为0.5。

• (2)．均匀分布：如果要用均匀分布来检验中心极限定理可直接用均匀随机数函数来生荿一个均匀随机数的值。

• (3)．正态分布：下面将用编写函数的方法生成1000个服从正态分布的随机数的均值标准化变量并用sim.fun( )函数模拟，然后用矗方图分析这些标准化变量是否服从正态分布
  

• (4)．指数分布：下面举一个偏态分布的数据的例子，检验随着n的增大其样本均值是否服从Φ心极限定理，不妨使用指数分布数据来模拟下面首先编写函数生成均值和标准差都为10的指数分布数据（指数分布的? =? =1/?），并求样夲均值标准化变量

  例4-9：接下来我们分别模拟n取1，510，30情况假如每次生成的随机数都是100，并作直方图以及正态分布密度线
  

  结果见图4-16，鈳见随着n的增大样本均值所服从的分布越来越接近于正态分布。