R语言作业的完成过程

实验五数据结构（向量、因子、列表）

实验目的：

1.熟悉R中数据基本结构特点及使用方法 2.熟悉R中向量的编辑方法 3.熟悉R中因子的设置方法 4.熟悉R中列表的编辑方法

实验要求：完成每个实验内容的同时，在每题后面附上程序代码。要求独立完成。

实验内容：

一、编写一个函数，具有如下功能： 1.输入参数为一个向量或者数据框

2.计算向量（或数据框中的列向量）的均值、方差、标准差、中位数 3.计算向量（或数据框中的列向量）的标准化得分

4.将2与3中的计算结果存储到一个列表中作为函数返回对象 > f<-function(x){ + mean<-mean(x) + var<-var(x) + sd<-sd(x)

+ median<-median(x) + m<-(x-mean(x))/sd(x)

+ j<-list(f.mean=mean,f.var=var,f.sd=sd,f.median=median,f.m=m) + j + }

二、有三组学生使用了不同的数学辅导教材，期末考试成绩如下。问这三组学生的成绩有无显著差异？（使用stack函数将三组向量合并成一个向量以及一个平行因子，然后使用aov函数或者oneway.test函数进行方差分析。）

99 88 99 89 94 90 第一组

55 50 67 67 56 56 第二组

79 56 89 99 70 89 第三组

> x1<-c(99,88,99,89,94,90) > x2<-c(55,50,67,67,56,56) > x3<-c(79,56,89,99,70,89) > f<-stac

泊松过程的模拟和检验

对保险人而言，资产和负债是影响保险人稳定经营至关重要的因素。资产和负债的差额称为盈余，简记作：

其中A(t)A(t)表示时刻tt的资产，L(t)L(t)表示时刻tt的负债，t=0t=0时刻的盈余被称为初始盈余，简记为uu，即U(0)=uU(0)=u。对这个初步的理论模型进行简化并根据实际情况设置一些假定情况，会得出很多不同的盈余过程模型，最经典的有Sparre Andersen的古典盈余过程模型：

这是一个以uu为初值，以时间tt为指标集的随机过程。其中

称为总理赔过程，满足：

N(t)N(t)表示[0,t][0,t]内的总理赔次数，XiXi表示[0,t][0,t]内第ii次理赔的金额。

根据这个古典盈余过程模型可以引出破产模型，在这个盈余过程模型中，一方面有连续不断的保费收入并以速度c进行积累，另一方面则是不断会有理赔需要支付，因此这是一个不断跳跃变化的过程。从保险人的角度来看，当然希望ct?S(t)ct?S(t)恒大于0，否则就有可能出现U(t)<0U(t)<0的情况，这种情况可以定义为理论意义上的破产，以示与实际中的破产相区分，本文中后面出现的“破产”在没有特殊说明的情况下都是指这种理论情况。从研究保险人破产角度出发，可以把这个盈余过程模型看做是一个特殊的破产模型。

一、 泊松过程的模拟

理论基础：泊松过程构造定理 具体步骤： 1、 2、 3、 4、

即满足泊松过程

生成一定数量的满足指数分布的随机数，用()表示 ()表示第n次事件到达的

01. R语言的安装与使用

R是统计分析软件，在数据分析、数据挖掘领域有优势；与SAS比最大的好处就是免费，其次是有各种现成的程序包可用，再次是具有极好的制图功能（数据可视化）。

一、下载安装

到R语言官网或镜像站：

https://www.r-project.org/

https://mirrors.tuna.tsinghua.edu.cn/CRAN/ 下载最新版R-3.3.1安装即可，安装时可选择32位或64位版本。

建议再接着安装相应版本的R语言集成开发环境Rstudio，地址：

https://www.rstudio.com/products/rstudio/download/

升级R：

在其它目录安装R，再将旧版本保留的library目录下的文件（程序包packages）拷贝到新版本library目录下，然后update.packages()

二、操作界面

图1R-3.3.1操作界面

类似Matlab的命令窗口，输入命令，按回车得到运行结果。 【程序包】菜单可以实现载入程序包、设定CRAN镜像、安装及更新程序包等功能。

图2 RStudio操作界面

代码编辑窗口：新建或打开脚本可调出，用来编写代码，打开R脚本，单击Run运行代码

Stata如何做cubic spine图横轴：Number of nodes examined；纵轴：Log hazard ratio Number of nodes examined Survival months Sensored Age Gender 18 11 6 16 12 11 8 12 2 3 2 5 4 1 4

13 8 5 6 4 16 30 26 11 3 5 15 68 16 46 1 1 1 1 1 1 1 1 1 0 1 1 1 1 1 33 63 73 63 68 65 51 74 72 72 69 66 59 54 54 Male Female Male `Male Male Male Male Male Female Male Male Male Male Female Male

rcspline.plot(x,y,model=c(\\\xrange, event, nk=5,knots=NULL, show=c(\adj=NULL, xlab, ylab,ylim, plim=c(0,1), plotcl=TRUE, showknots=TRUE, add=FALSE,subset, lty=1, nop

基于R语言的七种多重比较方法

一花视界

百家号10-1403:18

多重比较的方法很多，根据试验设计的目的不同有不同的应用。

若试验设计之初，便明确要比较某几个组均数间是否有差异，称为事前比较。常用的事前比较方法有LSD、Bonferroni和Dunnett法。

若研究目的是方差分析有统计学差异后，想知道哪些组间的均数有差异，便是事后比较。事后比较的常用方法有SNK、Turkey、Scheffe和Bonferroni法。 本文仅介绍7种方法及R语言函数，可解决绝大部分多重比较问题。 1.LSD法

LSD法即最小显著差法；该法一般用于计划好的多重比较。它其实只是t检验的一个简单变形，并未对检验水准做出任何校正，只是为所有组的均数统一估计了一个更为稳健的标准误。

LSD法比较效果较为灵敏，在R语言中可利用agricolae包中的LSD.test函数实现，其调用格式为：

LSD.test(y, trt, DFerror, MSerror, alpha = 0.05, p.adj=c("none","holm","hommel", "hochberg", "bonferroni", "BH", "BY", "fdr"), …)

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，typeof：对象存储模式与类型 names：对象的名字属性

二、字符串处理 character：字符型向量 nchar：字符数 substr：取子串

format，format C：把对象用格式转换为字符串 paste，strsplit：连接或拆分 charmatch，pmatch：字符串匹配 grep，sub，gsub：模式匹配与替换

三、复数

complex，Re，Im，Mod，Arg，Conj：复数函数

四、因子 factor：因子 codes：因子的编码 levels：因子的各水平的名字 nlevels：因子的水

作业讲解 - 用R语言绘制犯罪率增长高于20%的美国州地图

A criminology student has collected data about the number of violent assaults in each of the US states in 2014 and 2015. For his upcoming seminar, he is interested in identifying the states in which there has been at least a 20% increase in assaults.

Please help this student by identifying the relevant states, presenting them in a list or table, and creating a map of the USA where the relevant states are colored in red. This will require finding a relevant R package and using it to plot the map

###1.R introduction #### #语言历史和特点 #1.1 R语言历史####

#R是S语言的一种实现。S语言是由AT&T贝尔实验室开发的一种用来进行数据探索、 #统计分析、作图的解释型语言。最初S语言的实现版本主要是S-PLUS.

#S-PLUS是一个商业软件，它基于S语言，并由MathSoft公司的统计科学部进一步完善。 #后来Auckland大学的Robert Gentleman和Ross Ihaka及其他志愿人员开发了一个R系统 #R的使用与S-PLUS有很多类似之处，两个软件有一定的兼容性

#1.2 R的特点####

#1.有效地数据处理和保存机制 #2.拥有一整套数组和矩阵的操作运算符 #3.一系列连贯而又完整的数据分析中间工具

#4.图形统计可以对数据直接进行分析和显示，可用于多种图形设备。

#5.一种相当完善、简洁和高效的程序设计语言。包括条件语句、循环语句、用户自定义的递归函数以及输入输出接口。

#6.R语言是彻底面向对象的统计编程语言

#7.R语言和其他编程语言、数据库之间有很好的接口

#8.R语言是自由软件，可以放心大胆地使用，但其功能却不必任何其他同类软件差 #9.R语言具有丰富的网上资源

###

实验4 R绘图（一）

一、实验目的：

1. 掌握描述性统计分析中常用的统计量；

2. 掌握R语言绘制直方图、密度估计曲线、经验分布图和QQ图的方法； 3. 掌握R语言绘制茎叶图、箱线图的方法；

4. 掌握W检验方法和K-S检验方法完成数据的正态分布检验。

二、实验内容：

练习： 要求：①完成练习并粘贴运行截图到文档相应位置(截图方法见下)，并将所有自己输入文字的字体颜色设为红色（包括后面的思考及小结），②回答思考题，③简要书写实验小结。④修改本文档名为“本人完整学号姓名1”，其中1表示第1次实验，以后更改为2,3,...。如文件名为“1305543109张立1”，表示学号为1305543109的张立同学的第1次实验，注意文件名中没有空格及任何其它字符。最后连同数据文件、源程序文件等（如果有的话，本次实验没有），一起压缩打包发给课代表，压缩包的文件名同上。 截图方法：

法1：调整需要截图的窗口至合适的大小，并使该窗口为当前激活窗口（即该窗口在屏幕最前方），按住键盘Alt键（空格键两侧各有一个）不放，再按键盘右上角的截图键（通常印有“印屏幕”或“Pr Scrn”等字符），即完成截图。再粘贴到word文档的相应位置即可。 法2：利用QQ输入法的

统计软件及应用课程设计

-----虫情危害预测

学 院： 班 级： 学 号： 姓 名： 指导老师：

目录

一．背景与意义.................................................................................................... 2 二．问题重述......................................................................................................... 2 三．方法简介......................................................................................................... 2

(一).判别分析................................................................................................ 2 (二)