最新文章专题视频专题问答1问答10问答100问答1000问答2000关键字专题1关键字专题50关键字专题500关键字专题1500TAG最新视频文章推荐1 推荐3 推荐5 推荐7 推荐9 推荐11 推荐13 推荐15 推荐17 推荐19 推荐21 推荐23 推荐25 推荐27 推荐29 推荐31 推荐33 推荐35 推荐37视频文章20视频文章30视频文章40视频文章50视频文章60 视频文章70视频文章80视频文章90视频文章100视频文章120视频文章140 视频2关键字专题关键字专题tag2tag3文章专题文章专题2文章索引1文章索引2文章索引3文章索引4文章索引5123456789101112131415文章专题3
当前位置: 首页 - 正文

《应用回归分析》课后习题部分答案何晓群版

来源:动视网 责编:小OO 时间:2025-09-23 18:42:09
文档

《应用回归分析》课后习题部分答案何晓群版

第二章一元线性回归2.14解答:(1)散点图为:(2)x与y之间大致呈线性关系。(3)设回归方程为=(4)=(5)由于服从自由度为n-2的t分布。因而也即:=可得即为:(2.49,11.5)服从自由度为n-2的t分布。因而即可得(6)x与y的决定系数(7)ANOVAx平方和df均方F显著性组间(组合)9.00024.5009.000.100线性项加权的8.16718.16716.333.056偏差.8331.8331.667.326组内1.0002.500总数10.0004由于,拒绝,说明回归
推荐度:
导读第二章一元线性回归2.14解答:(1)散点图为:(2)x与y之间大致呈线性关系。(3)设回归方程为=(4)=(5)由于服从自由度为n-2的t分布。因而也即:=可得即为:(2.49,11.5)服从自由度为n-2的t分布。因而即可得(6)x与y的决定系数(7)ANOVAx平方和df均方F显著性组间(组合)9.00024.5009.000.100线性项加权的8.16718.16716.333.056偏差.8331.8331.667.326组内1.0002.500总数10.0004由于,拒绝,说明回归
第二章  一元线性回归

 2.14 解答:(1)散点图为:

            

 (2)x与y之间大致呈线性关系。

 (3)设回归方程为

      =

(4)

        

        =

       

   

(5)由于

服从自由度为n-2的t分布。因而

也即:=

可得

即为:(2.49,11.5)

  服从自由度为n-2的t分布。因而

可得

(6)x与y的决定系数

(7)

ANOVA
x
平方和df均方F显著性
组间(组合)9.00024.5009.000.100
线性项加权的8.16718.16716.333.056
偏差.8331.8331.667.326
组内1.0002.500
总数10.0004
由于,拒绝,说明回归方程显著,x与y有显著的线性关系。

(8)       其中

      

接受原假设认为显著不为0,因变量y对自变量x的一元线性回归成立。

(9)相关系数    

                   =

小于表中的相应值同时大于表中的相应值,x与y有显著的线性关系.

(10)           

序号
111064
221013-3
3320200
442027-7
5540346
残差图为:

从图上看,残差是围绕e=0随机波动,从而模型的基本假定是满足的。

(11)当广告费=4.2万元时,销售收入

,即(17.1,39.7)

2.15 解答:

(1) 散点图为:

(2)x与y之间大致呈线性关系。

(3)设回归方程为

      =

(4) 

     

=0.2305

0.4801

(5) 由于

服从自由度为n-2的t分布。因而

也即:=

可得

即为:(0.0028,0.0044)

  服从自由度为n-2的t分布。因而

可得

(6)x与y的决定系数 =0.908

(7) 

ANOVA
x
平方和df均方F显著性
组间(组合)1231497.5007175928.2145.302.168
线性项加权的1168713.03611168713.03635.222.027
偏差62784.46104.077.315.885
组内66362.500233181.250
总数1297860.0009
由于,拒绝,说明回归方程显著,x与y有显著的线性关系。

(8)        其中

      

接受原假设认为显著不为0,因变量y对自变量x的一元线性回归成立。

(9) 相关系数    

                   =

小于表中的相应值同时大于表中的相应值,x与y有显著的线性关系.

(10)

序号
18253.5

3.07680.4232
221510.88080.1192
3107043.95880.0412
455022.0868-0.0868
548011.8348-0.8348
692033.4188-0.4188
713504.54.9688-0.4668
83251.51.27680.2232
967032.51880.4812
10121554.48080.5192

从图上看,残差是围绕e=0随机波动,从而模型的基本假定是满足的。

(11)

(12),

即为(2.7,4.7)

近似置信区间为:,即(2.74,4.66)

(13)可得置信水平为为,即为(3.33,4.07).

2.16 (1)散点图为:

可以用直线回归描述y与x之间的关系.

(2)回归方程为:

(3) 

从图上可看出,检验误差项服从正态分布。

第三章 多元线性回归

3.11 解:(1)用SPSS算出y,x1,x2,x3相关系数矩阵:

相关性
yx1x2x3
Pearson 相关性

y1.000.556.731.724
x1.5561.000.113.398
x2.731.1131.000.547
x3.724.398.5471.000
 y..048.008.009
x1.048..378.127
x2.008.378..051
x3.009.127.051.
Ny10101010
x110101010
x210101010
x310101010
所以=

系数a

模型非标准化系数标准系数tSig.B 的 95.0% 置信区间

相关性共线性统计量
B标准 误差

试用版下限上限零阶部分容差VIF
1(常量)

-348.280176.459-1.974.096-780.06083.500
x13.7541.933.3851.942.100-.9778.485.556.621.350.8251.211
x27.1012.880.5352.465.049.05314.149.731.709.444.6871.455
x312.44710.569.2771.178.284-13.41538.310.724.433.212.5861.708
a. 因变量: y

 

(2) 

所以三元线性回归方程为    

模型汇总
模型RR 方

调整 R 方

标准 估计的误差

更改统计量
R 方更改

F 更改

df1df2Sig. F 更改

1.8a

.806.70823.44188.8068.28336.015
a. 预测变量: (常量), x3, x1, x2。

(3)

由于决定系数R方=0.708      R=0.8较大所以认为拟合度较高

(4)

Anovab

模型平方和df均方FSig.
1回归13655.37034551.7908.283.015a

残差3297.1306549.522
总计16952.5009
a. 预测变量: (常量), x3, x1, x2。

b. 因变量: y

因为F=8.283       P=0.015<0.05所以认为回归方程在整体上拟合的好

(5)

系数a

模型非标准化系数标准系数tSig.B 的 95.0% 置信区间

相关性共线性统计量
B标准 误差

试用版下限上限零阶部分容差VIF
1(常量)

-348.280176.459-1.974.096-780.06083.500
x13.7541.933.3851.942.100-.9778.485.556.621.350.8251.211
x27.1012.880.5352.465.049.05314.149.731.709.444.6871.455
x312.44710.569.2771.178.284-13.41538.310.724.433.212.5861.708
a. 因变量: y

(6)可以看到P值最大的是x3为0.284,所以x3的回归系数没有通过显著检验,应去除。

去除x3后作F检验,得:

Anovab

模型平方和df均方FSig.
1回归123.199246.60011.117.007a

残差4059.3017579.900
总计16952.5009
a. 预测变量: (常量), x2, x1。

b. 因变量: y

由表知通过F检验

继续做回归系数检验

系数a

模型非标准化系数标准系数tSig.B 的 95.0% 置信区间

相关性共线性统计量
B标准 误差

试用版下限上限零阶部分容差VIF
1(常量)

-459.624153.058-3.003.020-821.547-97.700
x14.6761.816.4792.575.037.3818.970.556.697.476.9871.013
x28.9712.468.6763.634.0083.13414.808.731.808.672.9871.013
a. 因变量: y

此时,我们发现x1,x2的显著性大大提高。

(7)x1:(-0.997,8.485)    x2:(0.053,14.149)   x3:(-13.415,38.310)

(8)

(9) 

残差统计量a

极小值极大值均值标准 偏差

N
预测值175.4748292.5545231.500038.9520610
标准 预测值

-1.4381.567.0001.00010
预测值的标准误差10.46620.19114.5263.12710
调整的预测值188.3515318.1067240.183549.8391410
残差-25.1975933.22549.0000019.1402210
标准 残差

-1.0751.417.000.81610
Student 化 残差

-2.1161.754-.1231.18810
已删除的残差-97.6152350.88274-8.6834843.4322010
Student 化 已删除的残差

-3.8322.294-.2551.65810
Mahal。 距离

.45.7772.7001.55510
Cook 的距离

.0003.216.486.97610
居中杠杆值.099.2.300.17310
a. 因变量: y

所以置信区间为(175.4748,292.5545)

(10)由于x3的回归系数显著性检验未通过,所以居民非商品支出对货运总量影响不大,但是回归方程整体对数据拟合较好

3.12 解:在固定第二产业增加值,考虑第三产业增加值影响的情况下,第一产业每增加一个单位,GDP就增加0.607个单位。

     在固定第一产业增加值,考虑第三产业增加值影响的情况下,第二产业每增加一个单位,GDP就增加1.709个单位。

第四章 违背基本假设的情况

4.9 解:

系数a

模型非标准化系数标准系数tSig.
B标准 误差

试用版
1(常量)

-.831.442-1.882.065
x.004.000.83911.030.000
a. 因变量: y

由SPSS计算得:=-0.831+0.004x

残差散点图为:

(2)由残差散点图可知存在异方差性

再用等级相关系数分析:

相关系数
xt
Spearman 的 rho

x相关系数1.000.318*

Sig.(双侧)

..021
N5353
t相关系数.318*

1.000
Sig.(双侧)

.021.
N5353
*. 在置信度(双测)为 0.05 时,相关性是显著的。

P=0.021 所以方差与自变量的相关性是显著的。

(3)

模型描述
因变量y
自变量1x
权重x
幂值1.500
模型: MOD_1.

M=1.5时可以建立最优权函数,此时得到:

ANOVA
平方和df均方FSig.
回归.0061.00698.604.000
残差.00351.000
总计.00952
系数
未标准化系数标准化系数tSig.
B标准误试用版标准误
(常数)-.683.298-2.296.026
x.004.000.812.0829.930.000
所以:-0.683+0.004x

(4)

系数a

模型非标准化系数标准系数tSig.
B标准 误差

试用版
1(常量)

.582.1304.481.000
x.001.000.8059.699.000
a. 因变量: yy

4.13 解:

(1)

系数a

模型非标准化系数标准系数tSig.
B标准 误差

试用版
1(常量)

-1.435.242-5.930.000
x.176.002.999107.928.000
a. 因变量: y

=-1.435+0.176x

(2) 

模型汇总b

模型RR 方

调整 R 方

标准 估计的误差

Durbin-Watson
1.999a

.998.998.09744.663
a. 预测变量: (常量), x。

b. 因变量: y

DW=0.663  查DW分布表知:=0.95

所以DW<,故误差项存在正相关。

残差图为:

随t的变化逐次变化并不频繁的改变符号,说明误差项存在正相关。

(3)=1-0.5*DW=0.6685 计算得:

Y’        x’

7.39    44.90

7.65    45.80

6.84    40.69

8.00    48.50

7.79    46.85

8.26    49.45

7.96    48.47

8.28    50.04

7.90    48.03

8.49    51.17

7.88    47.26

8.77    52.33

8.93    52.69

9.32    54.95

9.29    55.54

9.48    56.77

9.38    55.83

9.67    58.00

9.90    59.22

模型汇总b

模型RR 方

调整 R 方

标准 估计的误差

Durbin-Watson
1.996a

.993.993.073951.344
a. 预测变量: (常量), xx。

b. 因变量: yy

系数a

模型非标准化系数标准系数tSig.
B标准 误差

试用版
1(常量)

-.303.180-1.684.110
xx.173.004.99649.011.000
a. 因变量: yy

得回归方程 =-0.303+0.173x’

即:=-0.303+0.6685+0.173(—0.6685)

(4)

模型汇总b

模型RR 方

调整 R 方

标准 估计的误差

Durbin-Watson
1.978a

.957.955.074491.480
a. 预测变量: (常量), x3。

b. 因变量: y3

系数a

模型非标准化系数标准系数tSig.
B标准 误差

试用版
1(常量)

.033.0261.273.220
x3.161.008.97819.528.000
a. 因变量: y3

△=0.033+0.161△

即:=0.033++0.161(-)

(5)差分法的DW值最大为1.48消除相关性最彻底,但是迭代法的值最小为0.07395,拟合的较好。

4.14解:(1)

模型汇总b

模型RR 方

调整 R 方

标准 估计的误差

Durbin-Watson
1.541a

.293.2329.69302.745
a. 预测变量: (常量), x2, x1。

b. 因变量: y

系数a

模型非标准化系数标准系数tSig.
B标准 误差

试用版
1(常量)

-574.062349.271-1.4.107
x1191.09873.309.3452.607.012
x22.045.911.2972.246.029
a. 因变量: y

回归方程为:=-574.062+191.098x1+2.045x2

DW=0.745残差图为:

(2)=1-0.5*DW=0.6275

模型汇总b

模型RR 方

调整 R 方

标准 估计的误差

Durbin-Watson
1.688a

.474.452257.6701.716
a. 预测变量: (常量), x22, x12。

b. 因变量: y2

系数a

模型非标准化系数标准系数tSig.
B标准 误差

试用版
1(常量)

-179.66890.337-1.9.052
x12211.77047.778.5224.432.000
x221.434.628.2692.283.027
a. 因变量: y2

此时得方程:’=-179.668+211.77x1’+1.434x2’

所以回归方程为:

(3)

模型汇总b

模型RR 方

调整 R 方

标准 估计的误差

Durbin-Watson
1.715a

.511.490283.791022.042
a. 预测变量: (常量), x23, x13。

b. 因变量: y3

系数a

模型非标准化系数标准系数tSig.
B标准 误差

试用版
1(常量)

7.69839.754.194.847
x13209.144.143.5444.755.000
x231.399.583.2742.400.020
a. 因变量: y3

此时得方程:△

所以回归方程为:

第五章 自变量选择与逐步回归

5.9 后退法:输出结果

系数a

模型非标准化系数标准系数tSig.
B标准 误差

试用版
1(常量)

1438.1202252.472.638.533
农业x1

-.626.168-1.098-3.720.002
工业x2

-.328.207-1.352-1.587.135
建筑业x3

-.383.555-.251-.691.501
人口x4

-.004.025-.014-.161.875
最终消费x5

.672.1303.7105.178.000
受灾面积x6

-.006.008-.015-.695.499
2(常量)

1079.754299.7593.602.003
农业x1

-.2.130-1.126-4.925.000
工业x2

-.303.131-1.249-2.314.035
建筑业x3

-.402.525-.263-.765.456
最终消费x5

.658.0953.6366.905.000
受灾面积x6

-.006.007-.017-.849.409
3(常量)

1083.150295.8163.662.002
农业x1

-.624.127-1.095-4.931.000
工业x2

-.373.093-1.535-3.998.001
最终消费x5

.657.0943.6276.981.000
受灾面积x6

-.005.007-.015-.758.460
4(常量)

874.604106.8698.184.000
农业x1

-.611.124-1.073-4.936.000
工业x2

-.353.088-1.454-3.994.001
最终消费x5

.637.03.5167.142.000
a. 因变量: 财政收入y

Anovae

模型平方和df均方FSig.
1回归1.365E862.274E7602.127.000a

残差528793.3191437770.951
总计1.370E820
2回归1.365E852.729E7772.734.000b

残差529767.8521535317.857
总计1.370E820
3回归1.3E843.411E7991.468.000c

残差550440.1031634402.506
总计1.370E820
4回归1.3E834.547E71355.753.000d

残差570180.9311733540.055
总计1.370E820
a. 预测变量: (常量), 受灾面积x6, 建筑业x3, 人口x4, 农业x1, 最终消费x5, 工业x2。

b. 预测变量: (常量), 受灾面积x6, 建筑业x3, 农业x1, 最终消费x5, 工业x2。

c. 预测变量: (常量), 受灾面积x6, 农业x1, 最终消费x5, 工业x2。

d. 预测变量: (常量), 农业x1, 最终消费x5, 工业x2。

e. 因变量: 财政收入y

模型汇总
模型RR 方

调整 R 方

标准 估计的误差

更改统计量
R 方更改

F 更改

df1df2Sig. F 更改

1.998a

.996.994194.34750.996602.127614.000
2.998b

.996.995187.93046.000.026114.875
3.998c

.996.995185.47913.000.585115.456
4.998d

.996.995183.13944.000.574116.460
a. 预测变量: (常量), 受灾面积x6, 建筑业x3, 人口x4, 农业x1, 最终消费x5, 工业x2。

b. 预测变量: (常量), 受灾面积x6, 建筑业x3, 农业x1, 最终消费x5, 工业x2。

c. 预测变量: (常量), 受灾面积x6, 农业x1, 最终消费x5, 工业x2。

d. 预测变量: (常量), 农业x1, 最终消费x5, 工业x2。

回归方程为:

   逐步回归法:输出结果

模型汇总
模型RR 方

调整 R 方

标准 估计的误差

更改统计量
R 方更改

F 更改

df1df2Sig. F 更改

1.994a

.9.988285.68373.91659.441119.000
2.996b

.992.991247.77768.0037.258118.015
3.998c

.996.995183.13944.00415.948117.001
a. 预测变量: (常量), 最终消费x5。

b. 预测变量: (常量), 最终消费x5, 农业x1。

c. 预测变量: (常量), 最终消费x5, 农业x1, 工业x2。

Anovad

模型平方和df均方FSig.
1回归1.354E811.354E81659.441.000a

残差1550688.6541981615.192
总计1.370E820
2回归1.359E826.794E71106.637.000b

残差1105088.0031861393.778
总计1.370E820
3回归1.3E834.547E71355.753.000c

残差570180.9311733540.055
总计1.370E820
a. 预测变量: (常量), 最终消费x5。

b. 预测变量: (常量), 最终消费x5, 农业x1。

c. 预测变量: (常量), 最终消费x5, 农业x1, 工业x2。

d. 因变量: 财政收入y

系数a

模型非标准化系数标准系数tSig.相关性
B标准 误差

试用版零阶部分
1(常量)

710.37290.17.816.000
最终消费x5

.180.004.99440.736.000.994.994.994
2(常量)

1011.912136.9017.392.000
最终消费x5

.311.0491.7186.374.000.994.832.135
农业x1

-.414.154-.726-2.694.015.987-.536-.057
3(常量)

874.604106.8698.184.000
最终消费x5

.637.03.5167.142.000.994.866.112
农业x1

-.611.124-1.073-4.936.000.987-.767-.077
工业x2

-.353.088-1.454-3.994.001.992-.696-.062
a. 因变量: 财政收入y

回归方程为:

5.10

(1)

模型汇总
模型RR 方

调整 R 方

标准 估计的误差

1.908a

.824.736625.88326
2.000b

.000.0001217.15945
a. 预测变量: (常量), x6, x3, x2, x4, x5。

b. 预测变量: (常量)

Anovac

模型平方和df均方FSig.
1回归1.830E753660971.6839.346.002a

残差3917298.52210391729.852
总计2.222E715
2回归.0000.000..b

残差2.222E7151481477.129
总计2.222E715
a. 预测变量: (常量), x6, x3, x2, x4, x5。

b. 预测变量: (常量)

c. 因变量: y

系数a

模型非标准化系数标准系数tSig.
B标准 误差

试用版
1(常量)

5922.8272504.3152.365.040
x24.82.507.6771.940.081
x32.374.842.7822.818.018
x4-817.901187.279-1.156-4.367.001
x514.539147.078.050.099.923
x6-846.867291.634-.9-2.904.016
2(常量)

7542.938304.29024.7.000
a. 因变量: y

回归方程为:

(2)后退法:输出结果

模型汇总
模型RR 方

调整 R 方

标准 估计的误差

1.908a

.824.736625.88326
2.907b

.824.759597.04776
a. 预测变量: (常量), x6, x3, x2, x4, x5。

b. 预测变量: (常量), x6, x3, x2, x4。

Anovac

模型平方和df均方FSig.
1回归1.830E753660971.6839.346.002a

残差3917298.52210391729.852
总计2.222E715
2回归1.830E744575257.66912.835.000b

残差3921126.262113566.024
总计2.222E715
a. 预测变量: (常量), x6, x3, x2, x4, x5。

b. 预测变量: (常量), x6, x3, x2, x4。

c. 因变量: y

系数a

模型非标准化系数标准系数tSig.
B标准 误差

试用版
1(常量)

5922.8272504.3152.365.040
x24.82.507.6771.940.081
x32.374.842.7822.818.018
x4-817.901187.279-1.156-4.367.001
x514.539147.078.050.099.923
x6-846.867291.634-.9-2.904.016
2(常量)

6007.3202245.4812.675.022
x25.0681.360.7063.727.003
x32.308.486.7604.750.001
x4-824.261167.776-1.165-4.913.000
x6-862.699232.4-.916-3.711.003
a. 因变量: y

(3)逐步回归

模型汇总
模型RR 方

调整 R 方

标准 估计的误差

1.498a

.248.1941092.83206
2.697b

.485.406937.95038
3.811c

.657.572796.60909
a. 预测变量: (常量), x3。

b. 预测变量: (常量), x3, x5。

c. 预测变量: (常量), x3, x5, x4。

Anovad

模型平方和df均方FSig.
1回归5502210.09015502210.0904.607.050a

残差1.672E7141194281.918
总计2.222E715
2回归1.079E725392697.5546.130.013b

残差1.144E713879750.910
总计2.222E715
3回归1.461E734869041.5067.673.004c

残差7615032.41812634586.035
总计2.222E715
a. 预测变量: (常量), x3。

b. 预测变量: (常量), x3, x5。

c. 预测变量: (常量), x3, x5, x4。

d. 因变量: y

系数a

模型非标准化系数标准系数tSig.
B标准 误差

试用版
1(常量)

5161.2591142.7444.517.000
x31.511.704.4982.146.050
2(常量)

472.2982150.138.220.830
x33.188.9131.0503.492.004
x5212.32586.3.7372.451.029
3(常量)

1412.8071865.912.757.4
x33.440.7821.1334.398.001
x5348.72992.2201.2103.782.003
x4-415.136169.163-.587-2.454.030
a. 因变量: y

(4)两种方法得到的模型是不同的,回退法剔除了x5,保留了x6, x3, x2, x4作为最终模型。而逐步回归法只引入了x3。说明了方法对自变量重要性的认可不同的,这与自变量的相关性有关联。相比之下,后退法首先做全模型的回归,每一个变量都有机会展示自己的作用,所得结果更有说服力

第六章 多重共线性的情形及其处理

6.6

解:由下表我们可以看出

系数a

模型非标准化系数标准系数tSig.共线性统计量
B标准 误差

试用版容差VIF
1(常量)

-6381.5752736.958-2.332.035
x1-.593.279-1.040-2.127.052.003318.536
x2.549.1992.2602.753.016.0017.470
x3-.756.911-.495-.830.420.002472.951
x4.080.031.2812.590.021.015.706
x5.006.006.038.918.374.4342.305
x6-.010.014-.027-.750.466.5741.742
a. 因变量: y

方差扩大因子最大的为VIF2=7.470,故首先应剔除变量x2.将剩下变量继续进行回归得下表:

系数a

模型非标准化系数标准系数tSig.共线性统计量
B标准 误差

试用版容差VIF
1(常量)

-2677.4222858.846-.937.3
x1-.053.237-.092-.221.828.006160.620
x31.433.533.9372.690.017.009112.478
x4.036.032.1271.137.274.08711.509
x5.006.008.041.822.424.4342.303
x6.002.015.006.157.878.71.545
a. 因变量: y

此时,有最大的方差扩大因子VIF1=160.620,且此时x1系数为负,故x1也应被剔除,继续将剩下变量进行回归得:

系数a

模型非标准化系数标准系数tSig.共线性统计量
B标准 误差

试用版容差VIF
1(常量)

-2214.1291888.503-1.172.258
x31.318.109.86212.068.000.1995.023
x4.031.019.1071.586.132.2214.523
x5.006.007.041.841.412.4342.302
x6.003.015.008.209.837.6711.4
a. 因变量: y

此时,所有方差扩大因子都小于10,故回归方程如下:

=-2214.129+1.318x3+0.031x4+0.006x5+0.003x6

文档

《应用回归分析》课后习题部分答案何晓群版

第二章一元线性回归2.14解答:(1)散点图为:(2)x与y之间大致呈线性关系。(3)设回归方程为=(4)=(5)由于服从自由度为n-2的t分布。因而也即:=可得即为:(2.49,11.5)服从自由度为n-2的t分布。因而即可得(6)x与y的决定系数(7)ANOVAx平方和df均方F显著性组间(组合)9.00024.5009.000.100线性项加权的8.16718.16716.333.056偏差.8331.8331.667.326组内1.0002.500总数10.0004由于,拒绝,说明回归
推荐度:
  • 热门焦点

最新推荐

猜你喜欢

热门推荐

专题
Top