文档介绍：数理统计上机报告
上机实验题目:用R软件进行回归分析
上机实验目的:
1 进一步理解回归分析的基本思想,学会使用回归进行统计推理。
2 学会利用R软件进行回归分析的方法。
一元线性回归基本理论、方法:
1 根据样本观察值对经济计量模型参数进行估计,求的回归方程。
2 对回归方程、参数估计值进行显著性检验。
3 利用回归方程进行分析、评论及预测。
P430
第十一题
上机实验步骤:
y<-c(2813,2705,11103,2590,2131,5181)
x<-c(,,,,,)
xbar<-mean(x)
L11<-sum((x-xbar)^2)
ybar<-mean(y)
Lyy<-sum((y-ybar)^2)
L1y<-sum((x-xbar)*(y-ybar))
n<-length(x)
beta_1<-L1y/L11
beta_0<-ybar-xbar*beta_1
sigma2_hat<-(Lyy-beta_1*L1y)/(n-2)
sigma_hat<-sqrt(sigma2_hat)
实例计算结果及分析:
> L11
[1]
> Lyy
[1] 59353704
> L1y
[1]
> beta_0
[1] -
> beta_1
[1]
> sigma_hat
[1]
>
P=-+4165L
P432 第十八题
上机实验步骤:
y<-c(64,60,71,61,54,77,81,93,93,51,76,96,77,93,95,54,168,99)
X<-matrix(0, nrow = 18, ncol = 4)
X[,1]<-rep(1,18)
X[,2]<-c(,,,,,,,,,,,,,,,,,)
X[,3]<-c(53,23,19,34,24,65,44,31,29,58,37,46,50,44,56,36,58,51)
X[,4]<-c(158,163,37,157,59,123,46,117,173,112,111,114,134,73,168,143,202,124)
beta<-solve(t(X)%*%X)%*%t(X)%*%y
yhat<-X%*%beta
ytidle<-y-yhat
实例计算结果及分析:
> yhat
[,1]
[1,]
[2,]
[3,]
[4,]
[5,]
[6,]
[7,]
[8,]
[9,]
[10,]