关于正交旋转,请高人指点

shiyiming

我刚接触SAS，兴趣很高,能力却不够.最近要处理一些样品数据,需使用方差最大的正交旋转方法对因子载荷矩阵进行正交旋转(方法从论文上看到的),却不知道如何操作.望高人不吝指教,感激!

shiyiming

先进行主成分分析，以初步估计因子个数
[code:3dbad]proc factor /* 调用因子分析模块 */
method=p  /* 分析方法为主成分因子分析法 */
scree     /* 输出碎石图，以分析因子的个数 */
outstat=fac1;  /* 因子分析结果输出，以供下一步分析用 */
var  要分析的指标 ;
run; [/code:3dbad]

[code:3dbad]proc factor
data=fac1  /* 从主成分分析结果开始分析 */
method=p  /* 分析方法为主成分因子分析法 */
priors=input
/* 指定各指标的公因子方差为起先主成分分析的 */
n=2 /* 指定因子个数为2 */
outstat=fac2;  /* 因子分析结果输出，以供下一步分析用 */
run; [/code:3dbad]

因子载荷矩阵旋转：各因子保持正交、使因子间方差最大
[code:3dbad]proc factor data=fac2 /* 从主因子分析结果开始分析 */
rotate=quartimax
/* 对主因子分析结果作quartimax正交四次方最大旋转 */
score; /* 输出因子得分系数矩阵 */
run;[/code:3dbad]
载荷矩阵旋转的目的：增强公因子解释能力

shiyiming

非常感谢你大哥。 <!-- s:D --><img src="{SMILIES_PATH}/icon_biggrin.gif" alt=":D" title="Very Happy" /><!-- s:D --> 我回去试一下。非常感谢！！！！！！！！！！！！

shiyiming

　　我要分析的是17组样品监测13种金属.因子个数为3.要分析因子对各变量(金属元素)的载荷来判断样品的来源及其贡献．
　　使用样本相关系数矩阵来进行因子分析，参数估计方法使用主分量法，公因子数的确定准则为特征根大于1，使用方差最大的正交旋转方法对因子载荷矩阵进行正交旋转，得到最终因子载荷矩阵．
　　感觉您的代码都是很对应的，可出来的结果有问题，与期望值差别很大．只有再麻烦您了．你看看哪些地方有问题．谢谢了．

shiyiming

两个问题：

1、因子的个数确定. 考虑到累积贡献率否？

2、试试其他的旋转方式，“因子命名是一门艺术”， <!-- s:) --><img src="{SMILIES_PATH}/icon_smile.gif" alt=":)" title="Smile" /><!-- s:) -->

shiyiming

[color=olive:06805]我所做的是17个样品监测各样品中13种金属.因子数是3.[/color:06805]
是通过分析各因子对各变量(金属元素)的载荷来判断样品的来源及其贡
献.使用样本相关系数矩阵来进行因子分析，参数估计方法使用主分量法，公因子数的确定准则为特征根大于1，使用方差最大的正交旋转方法对因
子载荷矩阵进行正交旋转.得到最终因子载荷矩阵.
    感觉您的方法没有问题,可出来的结果与期望值差别很大.是不是第3条代码中有东西需要改一下?
[color=brown:06805]我没办法.只有再麻烦你了.谢谢了![/color:06805]
     也可以依照相关矩阵解出特征值,由此得出对应的特征向量,然后使用累计方差贡献率来判别因子数.并由此算出初始因子载荷矩阵,然后采用方差最大正交旋转法得到最终因子载荷矩阵.
    [color=red:06805]两种方法好象差不多.应用范围也挺广.无奈我才上大3,知识太少,学的是地球科学.很无奈.希望您也能给我推荐一些学习方法.感觉SAS要学的东西很多.
  谢谢! <!-- s:P --><img src="{SMILIES_PATH}/icon_razz.gif" alt=":P" title="Razz" /><!-- s:P --> [/color:06805]