Commit a6a85e97 authored by rattoma's avatar rattoma
Browse files

don't be too conservative with PCA!!

git-svn-id: file:///home/sebastien/dynare/gsa_dyn@123 f1850c17-3b45-254b-b221-fcb05880fee1
parent 02c55718
......@@ -964,7 +964,8 @@ else, % main effects analysis
% npca = min(40, max(find(cumsum(latent)./length(latent)<0.99))+1);
npca = max(find(cumsum(latent)./length(latent)<0.99))+1;
siPCA = (EET.SAcc'*abs(pcc'))';
siPCA = siPCA./(max(siPCA')'*ones(1,npT)).*(latent*ones(1,npT));
% siPCA = siPCA./(max(siPCA')'*ones(1,npT)).*(latent*ones(1,npT));
siPCA = siPCA./(max(siPCA')'*ones(1,npT));
% siPCA = sum(siPCA,1);
% siPCA = siPCA./max(siPCA);
SAcc=zeros(size(ccac,2),npT);
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment