多元统计分析及R语言建模（第五版）——第9章

9.4因子分析法在股价预报上的探索：在本例中为了验证因子分析方法的有效性，特意不区分行业，以上交所和深交所进行分层，然后把层内全部股票选入抽样框，以进行随机抽取。从http://www.jrj. com.cn 得到了 23 个企业在2004 年3月31日的数据，见表9-4，所考虑的指标如下：流动比率x,(<2偏低）、速动比率x2(<1偏低）、现金流动负债比x,(%)、每股收益x（元）、每股未分配利润x3（元）、每股净资产x6（元）、每股资本公积金x,（元）、每股盈余公积金x（元）、每股净资产增长率x3(%)、经营净利率x10(%)、经营毛利率x1(%)、资产利润率x(%)、资产净利率（%)、主营收入增长率x14(%)、净利润增长率x13(%)、总资产增长率x。(%)、主营利润增长率初（%)、主营成本比例x18、营业费用比例x19(%)、管理费用比例x20(%)、财务费用比率x21(%).

(1）求样本相关系数阵 R 及特征根与特征向量。
(2）确定因子的个数，并解释这些因子的含义。
(3）计算各因子得分，画出前两个因子的得分图并解释之。
(4）对因子进行旋转，比较旋转前后因子分析的结果。
(5）对这23 个上市企业财务状况进行综合评价。
 

名称	X1	X2	X3	X4	X5	X6	X7	X8	X9	X10	X11	X12	X13	X14	X15	X16	X17	X18	X19	X20	X21
深深房Ａ	1.33	0.33	-1.44	-0.02	-0.92	1.11	0.95	0.12	11.1	-12.85	8.37	-0.67	-0.67	-16.66	23.64	-11.39	-40.31	89.76	3.43	15.8	6.61
同人华塑	0.76	0.73	-4.96	0.02	-1.42	1.57	1.86	0.12	-10.38	4.91	31.44	0.6	0.37	217.62	218.61	31.3	282.07	68.22	2.86	11.85	9.59
南开戈德	1.21	1.14	-4.72	-0.02	-0.72	1.03	0.53	0.23	-46.59	-28.03	0.47	-0.62	-0.64	133.02	56.28	-32.22	-46.59	98.7	1.87	10.94	13.27
st昌源	0.79	0.6	-0.19	-0.02	-1.96	0.03	0.86	0.11	-96.94	-993.95	40.57	-0.7	-0.7	-96.89	-1841.7	-33.45	-96.94	57.86	3.8	237.64	797.55
山东巨力	0.75	0.58	7.55	-0.35	-0.36	1.78	1.02	0.11	-28.02	0.24	7.59	0.11	0.07	-13.35	28.88	27.25	24.36	92.4	3.23	3.46	1.1
一汽夏利	1.11	1.01	7.27	0.01	-0.26	1.67	0.82	0.1	9.34	0.87	7.1	0.18	0.18	10.32	-76.29	-5.76	-48.34	88.59	3.94	4.05	2.03
闽东电力	1.19	1.01	-12.73	-0.03	-0.33	4.5	3.75	0.09	-7.01	-23.16	9.73	-0.32	-0.36	-33.57	-557.39	2.44	-88.61	89.73	0.57	44.33	19.72
深本实b	1.04	0.91	-17.41	-0.3	-0.3	1.78	0.49	0.52	-14.37	1.41	34.91	0.06	0.04	-21.27	4.25	6.11	-25.12	64.79	7.75	18.72	6.44
st啤酒花	0.22	0.16	0.15	-0.03	-3.28	-2.45	0.1	0.49	-253.8	-27.05	32.74	-1.23	-0.9	-71.06	-11.4	-47.84	-44.79	55.56	26.81	33.3	11.31
云大科技	1.57	1.18	-3.77	-0.08	-0.55	1.84	1.25	0.14	-30.63	-77.83	15.55	-1.58	-1.45	-63.56	-220.84	-9.97	-81.7	83.44	23.4	40.51	30.56
中天科技	1.98	1.28	-28.33	0.01	0.16	2.71	1.4	0.15	-1.39	2.27	19.2	0.23	0.16	18.01	-4.51	10.12	-2.53	80.8	9.44	6.43	1.6
爱建股份	1.36	0.82	7.47	0.05	-0.59	3.37	2.06	0.9	-13.76	6.38	17.26	1.65	0.75	491.52	311.79	-12.49	901.91	78.15	1.37	2.65	1.3
st轻骑	0.92	0.82	-4.03	0	-1.91	0.37	1.27	0.06	122.27	0.83	8.91	0.01	0.11	38.85	229.36	86.9	-3.42	90.07	4.32	5.46	-0.3
张裕A	4.037	3.49	49.27	0.24	1.31	5.27	2.48	0.48	-7.25	17.86	53.64	5.77	3.82	22.71	25.77	5.17	30.99	38.06	20.3	7.15	-0.74
阿继电器	1.95	1.43	-72.04	0.01	0.09	1.9	0.62	0.14	-38.3	6.5	37.67	0.36	0.33	-7.19	-35.42	5.78	-17.14	62.33	7.76	20.4	2.66
广州浪奇	1.91	1.38	-16.2	0.01	-0.47	1.71	0.72	0.45	0.6	0.29	14.5	0.11	0.1	61.5	-11.09	0.15	-5.12	85.04	8.99	5.61	-0.2
浙江震元	1.47	0.93	-3.68	0.02	0.45	3.56	1.84	0.26	3.44	1.04	13.19	0.4	0.25	19.88	-32.67	1.8	-25.33	86.48	4.37	6.87	0.29
四环生物	5.47	4.25	61.7	0.04	0.02	1.1	0	0	-47.48	32.85	46.64	4.54	3.28	10.59	99.3	18.05	38.88	52.27	0.7	10.12	0.97
深宝安A	1.49	0.35	0.05	0.01	-0.59	1.19	0.7	0.35	-0.1	1.45	37.3	0.24	0.06	38.92	-32.2	4.71	28.6	59.93	13.26	15.48	10.2
深发展A	0.78	0.78	1.58	0.11	0.2	2.15	0.81	0.08	6.57	11.2	36.78	0.11	0.11	41.47	39.82	30.79	80.33	0	26.38	0	0
数码网络	85	0.63	2.52	0.04	0.1	1.6	0.26	0.23	6.63	0.12	8.86	0.08	0.03	-10.08	111.57	1.41	8.52	91.06	2.84	4.65	1.5
中色建设	2.74	2.39	-27.9	0.02	0.21	2.61	1.31	0.14	9.92	11.38	13.7	0.38	0.39	5.08	260.93	22.84	9.92	85.97	2.95	26.44	2.29
东北药	1.24	1	3	0.01	-0.02	2.92	1.92	0.02	1.25	0.67	21.45	0.16	0.11	-4.64	56.87	-11.21	-7.17	78.32	7.46	9.05	4.18

library(openxlsx)
d9.4=read.xlsx("E:/多元统计分析及R语言建模/mvexer5.xlsx",sheet ="E9.4",rowNames = T)
d9.4
X=data.frame(d9.4)
R=cor(X);R#相关系数矩阵
R.e=eigen(R,symmetric = T);R.e#特征根
R.e$vectors%*%diag(R.e$values)%*%t(R.e$vectors)#特征向量


factanal(d9.4,10,rotation = "none")#因子不旋转
#前7个的累计方差达到83.6%，选前7个

source("E:/msaR.R",encoding = "utf-8")#要添加解码方式，否则报错
pa=msa.fa(d9.4,m=7,rotation = "none")
pa$scores
pa2=msa.fa(d9.4,m=2,rotation = "none")
plot(pa2$scores,asp=1);abline(h=0,v=0,lty=2)
text(pa2$scores,labels = rownames(d9.4))
biplot(pa2$scores,pa2$loadings);abline(h=0,v=0,lty=2)


fa=msa.fa(d9.4,m=7,rotation = "varimax")#旋转
fa
#从旋转后的因子载荷矩阵可以看出，公共因子Factor1在经营净利率x10、净利润增长率x15、管理费用比例x20、财务费用比率x21上的载荷值比较大，
#公共因子Factor2在速动比率x2、现金流动负债比x3、资产利润率x12和资产净利率x13上的载荷值比较大……
#由回归法估计出因子得分以各因子的方差贡献率占四个因子总方差贡献率的比重 作为权重进行加权，汇总得出个省、市、自治区的综合得分及排名。其中排名最高的几个企业分别为张裕A、爱建股份、深发展、四环生物。
#由此可见，上市企业的财务状况是由总营业利润决定的