2019-03-31 17:46:41
f1 = repmat(sqrt(val),size(vec1,1),1) #val 为pcacov的特征值
1.主成分分析(PCA Principal Components Analysis)
1.1 代码段
[vec1,lamda,rate] = pcacov(V)#V = corrcoef(data)
#vel1 为V 的特征向量,
#lamda 为V的特征值,
#rate为各个主成分的贡献率
f = repmat(sign(vec1),size(vec1,1),1)
#修改符号,全为正 vec2 = vec1.*f
1.2 框架
load data.txt
data = zscore(data)
V = corrcoef(data)
[vec1,lamda,rate] = pcacov(V)
f = repmat(sign(vec1),size(vec1,1),1)
vec2 = vec1.*f
contribution = cumsum(rate)
num = input('输入的主成分数目,基于Contribute >0.85')
df = data*vec2(:,1:mun)
Comprehensive_f = df*rete(1:num)/100
[stf,ind] = sort(Comprehensive_f,'descend')
stf = stf'
ind = ind'
2. 因子分析
2.1 代码段
[lamda2,t] = rotatefactor(A,'method','varrmax') #lamda2 为旋转后的载荷矩阵 #t为变换的正交矩阵 #A 一般为lamda(:,1:num)
[B,BINT,R,RINT,STATS] = regress(X,Y)
#BINT 是95%的置信区间 of B 是回归系数
#R是残差Residues
[d,dt,e,et,stats] = regress(Tscore,[ones(n,1),y])