1、读取数据

rm(list=ls())
a=read.table("D:/统计数据表.txt",header=T)
round(a[1:10,],2)

2、计算均值、最小值、中位数、最大值、标准差

a1=a[a$年份==2005,-1]
Mean=sapply(a1,mean)
Min=sapply(a1,min)
Median=sapply(a1,median)
Max=sapply(a1,max)
Sd=sapply(a1,sd)
cbind(Mean,Min,Median,Max,Sd)

3、各变量相关分析
round(cor(a),3)

5、相关系最大的两个变量的散点图
plot(a1$CPI,a1$GDP)

这几项都没有问题

为什么最小二乘估计的时候出错呢？？？？？？

lm1=lm(CPI~GDP+价格+涨幅+收入,date=a1)
summary(lm1)

> rm(list=ls())
> a=read.table("D:/统计数据表.txt",header=T)
> round(a[1:10,],2)
   年份 价格 涨幅 收入  GDP CPI
1  2005 1000  700 1000 4000 0.8
2  2005 1020  500 2000 3000 0.1
3  2008 1300  510 3000 2000 0.2
4  2005 1000  600 5000 1000 0.3
5  2005 4000  400 6000 5000 0.4
6  2005 6000  510 5000 5000 0.6
7  2009 4000  600 2000 1000 0.4
8  2005 2000  700 4000 2000 0.1
9  2015 1100  560 7000 1000 0.2
10 2005 1000  520 8000 2000 0.1
>
> a1=a[a$年份==2005,-1]
> Mean=sapply(a1,mean)
> Min=sapply(a1,min)
> Median=sapply(a1,median)
> Max=sapply(a1,max)
> Sd=sapply(a1,sd)
> cbind(Mean,Min,Median,Max,Sd)
             Mean   Min Median    Max           Sd
价格 2387.1428571 1e+03 1020.0 6000.0 1852.6962334
涨幅  561.4285714 4e+02  520.0  700.0  103.9430934
收入 4432.6530612 1e+03 5000.0 8200.0 2224.2402968
GDP  3142.8571429 1e+03 3000.0 5000.0 1471.9601444
CPI     0.3428571 1e-01    0.3    0.8    0.2581989
> round(cor(a),3)
       年份   价格   涨幅   收入    GDP    CPI
年份  1.000 -0.159  0.004  0.203 -0.522 -0.180
价格 -0.159  1.000 -0.293  0.025  0.508  0.445
涨幅  0.004 -0.293  1.000 -0.422 -0.370  0.203
收入  0.203  0.025 -0.422  1.000 -0.086 -0.381
GDP  -0.522  0.508 -0.370 -0.086  1.000  0.564
CPI  -0.180  0.445  0.203 -0.381  0.564  1.000
> plot(a1$CPI,a1$GDP)
> lm1=lm(CPI~GDP+价格+涨幅+收入,date=a1)
错误于eval(expr, envir, enclos) : 找不到这个对象"CPI"
> summary(lm1)
错误于summary(lm1) : 找不到这个对象"lm1"
>