1 / 15
文档名称:

R语言常用上机命令分功能整理——时间序列分析为.doc

格式:doc   页数:15页
下载后只包含 1 个 DOC 格式的文档,没有任何的图纸或源代码,查看文件列表

如果您已付费下载过本站文档,您可以点这里二次下载

分享

预览

R语言常用上机命令分功能整理——时间序列分析为.doc

上传人:825790901 2016/3/5 文件大小:0 KB

下载得到文件列表

R语言常用上机命令分功能整理——时间序列分析为.doc

文档介绍

文档介绍:第一讲应用实例?R 的基本界面是一个交互式命令窗口,命令提示符是一个大于号,命令的结果马上显示在命令下面。?S 命令主要有两种形式: 表达式或赋值运算(用’<-’或者’=’表示)。在命令提示符后键入一个表达式表示计算此表达式并显示结果。赋值运算把赋值号右边的值计算出来赋给左边的变量。?可以用向上光标键来找回以前运行的命令再次运行或修改后再运行。?S 是区分大小写的,所以 x和X 是不同的名字。我们用一些例子来看 R 软件的特点。假设我们已经进入了 R 的交互式窗口。如果没有打开的图形窗口,在 R 中,用: > x11() 可以打开一个作图窗口。然后,输入以下语句: x1= 0:100 x2= x1*2*pi/100 y= sin(x2) plot(x2,y,type="l") 这些语句可以绘制正弦曲线图。其中, “=”是赋值运算符。 0:100 表示一个从 0到 100 的等差数列向量。第二个语句可以看出, 我们可以对向量直接进行四则运算, 计算得到的 x2是向量x1 的所有元素乘以常数 2*pi/10 0 的结果。从第三个语句可看到函数可以以向量为输入, 并可以输出一个向量,结果向量 y 的每一个分量是自变量 x2 的每一个分量的正弦函数值。 plot(x2,y, type="l",main=" 画图练****quot;,sub=" 好好练", xlab="x 轴",ylab='y 轴') 有关作图命令 plot 的详细介绍可以在 R 中输入 help(plot) 数学函数 abs , sqrt : 绝对值, 平方根 log, log10, log2 , exp : 对数与指数函数 sin , cos , tan , asin , acos , atan , atan2 :三角函数 sinh , cosh , tanh , asinh , acosh , atanh :双曲函数简单统计量 sum, mean, var, sd, min, max, range, median, IQR ( 四分位间距) 等为统计量, sort , order , ran k 与排序有关,其它还有 ave , fivenum , mad , quantile , stem 等。下面我们看一看 S 的统计功能:> marks <- c(10, 6, 4, 7, 8) > mean(marks) > sd(marks) > min(marks) > max(marks) 第一个语句输入若干数据到一个向量,函 c() 用来把数据组合为一个向量。后面用了几个函数来计算数据的均值、标准差、最小值、最大值。可以把若干行命令保存在一个文本文件中,然后用 source 函数来运行整个文件: > source("C:/") 注意字符串中的反斜杠。例: 计算 6, 4, 7,8, 10 的均值和标准差,把若干行命令保存在一个文本文件(比如 C:\ ) 中,然后用 source 函数来运行整个文件。 a<- c(10, 6, 4, 7, 8) b<-mean(a) c<-sd(a) source("C:/") 时间序列数据的输入使用函数 ts ts(1:10, frequency = 4, start = c(1959, 2)) print( ts(1:10, frequency = 7, start = c(12, 2)), calendar = TRUE) a<-ts(1:10, frequency = 4, start = c(1959, 2)) plot(a) 将外部数据读入 R 默认 header = TRUE ,也就是第一行是标签,不是数据。 默认 header = FALSE 将R 中的数据输出 write 第二讲 1. 绘制时序图、自相关图例题 d=scan("") sha=ts(d,start=1964,freq=1) (sha) # 绘制时序图 acf(sha,22) # 绘制自相关图,滞后期数 22 pacf(sha,22) # 绘制偏自相关图,滞后期数 22 corr=acf(sha,22) # 保存相关系数 cov=acf(sha,22,type = "covariance") # 保存协方差图的保存,单击选中图,在菜单栏选中“文件”,再选“另存为”。同时显示多个图:用 x11() 命令生成一个空白图,再输入作图命令。 2. 同时绘制两组数据的时序图 d=("",header=F) double=ts(d,start=1964,freq=1) plot(double,