统计软件与R语言培训课程.pptVIP

画图(时间序列) x=runif(100);x=ts(x,start=1960+(3/12),frequency=12);options(digits=2);x plot(x,main=Time series) Jan Feb Mar Apr May Jun Jul Aug Sep 1960 0.28960 0.85129 0.05823 0.46995 0.57905 0.73799 1961 0.28672 0.20442 0.41045 0.96187 0.06345 0.18318 0.73739 0.08260 0.74940 1962 0.43173 0.10568 0.44284 0.09334 0.44954 0.10872 0.05065 0.63924 0.31070 1963 0.67849 0.62165 0.21193 0.98964 0.96315 0.29877 0.95799 0.84184 0.69913 1964 0.99902 0.86650 0.83739 0.26195 0.70877 0.94352 0.56905 0.76416 0.95544 1965 0.15065 0.00066 0.87050 0.78968 0.58603 0.35572 0.04752 0.25399 0.11059 1966 0.04642 0.74758 0.42476 0.50551 0.76536 0.73733 0.48250 0.07827 0.91807 1967 0.70602 0.94552 0.32180 0.20975 0.59760 0.04847 0.46288 0.40607 0.61738 1968 0.26741 0.01976 0.44851 0.90064 0.07154 0.59455 0.51154 Oct Nov Dec 1960 0.66690 0.06929 0.10811 1961 0.16218 0.59930 0.87235 1962 0.62558 0.57628 0.68994 1963 0.23256 0.93256 0.78502 1964 0.66011 0.79481 0.32449 1965 0.91167 0.44189 0.80639 1966 0.35219 0.89851 0.32350 1967 0.40915 0.51130 0.77735 1968 第一百二十七页，共一百九十五页。 第一百二十八页，共一百九十五页。 画图(matplot) sines=outer(1:20,1:4,function(x, y) sin(x / 20 * pi * y)) matplot(sines, pch = 1:4, type = o, col = rainbow(ncol(sines))) options(digits=4) sines [,1] [,2] [,3] [,4] [1,] 1.564e-01 3.090e-01 4.540e-01 5.878e-01 [2,] 3.090e-01 5.878e-01 8.090e-01 9.511e-01 [3,] 4.540e-01 8.090e-01 9.877e-01 9.511e-01 [4,] 5.878e-01 9.511e-01 9.511e-01 5.878e-01 [5,] 7.071e-01 1.000e+00 7.071e-01 1.225e-16 [6,] 8.090e-01 9.511e-01 3.090e-01 -5.878e-01 [7,] 8.910e-01 8.090e-01 -1.564e-01 -9.511e-01 [8,] 9.511e-01 5.878e-01 -5.878e-01 -9.511e-01 [9,] 9.877e-01 3.090e-01 -8.910e-01 -5.878e-01 [10,] 1.000e+00 1.225e-16 -1.000e+00 -2.449e-16 [11,] 9.877e-01 -3.090e-01 -8.910e-01 5.878e-01 [12,] 9.511e-01 -5.878e-01 -5.878e-01 9.511e-01 [13,] 8.910e-01 -8.090e-01 -1.564