U1

软件安装、数据输入、外部文件导入、变换、子集、插值、导出

> sum(1,2,3
+ )
[1] 6

> sum(1,2,3)
[1] 6

> x=1;y=2;x+y
[1] 3

> {x=1
+ y=2
+ x+y
+ }
[1] 3

> x<-3
> x
[1] 3

直接录入

#行输入
> price=c(101,82,66,35,31,7)    #调用c函数
> price=ts(price,start=c(2005,1),frequency=12)
    #变量命名为price，start选项指定的起始读入时间
    #frequency选项指定序列每年读入的数据频率
> price
     Jan Feb Mar Apr May Jun
2005 101  82  66  35  31   7

#列输入
> price=scan()    #调用scan函数
1: 101
2: 82
3: 66
4: 35
5: 31
6: 7
7: 
Read 6 items
> price=ts(price,start=c(2005,1),frequency=12)
> price
     Jan Feb Mar Apr May Jun
2005 101  82  66  35  31   7

外部数据文件转换

> x=read.table("E:/例题data/file1.csv",sep=",",header=T)
    #sep指定制表分隔符，如果是txt格式，sep=“\t”，即文件以制表符分隔，如果是csv格式，sep=“，”，即以逗号分隔
    #header指定第一行是否包含变量名，header=TURE（或简写为“T”）表示要转换的文本文件第一行包含变量名，这是系统指定的默认值，header=FALSE（或简写为“F”）表示改文件第一行不包含变量名
> x
   Year yield
1  1884  15.2
2  1885  16.9
。。。
> x$yield    #单独考察这两个变量的特征
 [1] 15.2 16.9 15.3 14.9 15.7 15.1 16.7 16.3 16.5 13.3 16.5 15.0 15.9
[14] 15.5 16.9 16.4 14.9 14.5 16.6 15.1 14.6 16.0 16.8 16.8 15.5 17.3
[27] 15.5 15.5 14.2 15.8 15.7 14.1 14.8 14.4 15.6 13.9 14.7 14.3 14.0
[40] 14.5 15.4 15.3 16.0 16.4 17.2 17.8 14.4 15.0 16.0 16.8 16.9 16.6
[53] 16.2 14.0 18.1 17.5
>

数据变换

> y=log(x$yield)
> y
 [1] 2.721295 2.827314 2.727853 2.701361 2.753661 2.714695 2.815409
 [8] 2.791165 2.803360 2.587764 2.803360 2.708050 2.766319 2.740840
[15] 2.827314 2.797281 2.701361 2.674149 2.809403 2.714695 2.681022
[22] 2.772589 2.821379 2.821379 2.740840 2.850707 2.740840 2.740840
[29] 2.653242 2.760010 2.753661 2.646175 2.694627 2.667228 2.747271
[36] 2.631889 2.687847 2.660260 2.639057 2.674149 2.734368 2.727853
[43] 2.772589 2.797281 2.844909 2.879198 2.667228 2.708050 2.772589
[50] 2.821379 2.827314 2.809403 2.785011 2.639057 2.895912 2.862201

数据导出

> ln_yield=log(x$yield)
> x_new=data.frame(x,ln_yield)
> write.table(x_new,file="E:/例题data/file1.csv",sep=",",row.names=F)