// Selecione o arquivo PCA , cole no notepad e salve em algum dir:
// por exemplo : "d:\bin\scilab-2.6\quimiometria"
// com o nome "pca.sci", para isto selecione, salvar como no notepad(bloco de notas)
// em tipo do arquivo selecion *.*, ou todos os arquivos

// Para executar os comandos deste  tutorial copie os comandos entre:
//---------------------------------------------------------------//

//---------------------------------------------------------------//

// na janela do SCILAB
// Obs.: Para colar no  SCILAB use o botão da direita e escolhe"paste"
////////////////////////////////////////////////////////////////////////////
// Matriz de Dados
// Descr. no artigo: M.M.C. Ferreira et al.,Food Chemistry 69 (2000) 259
//---------------------------------------------------------------//
X=[
31.8   46.1    21.78   0.37 43.59 20.31  2.24 207.56 1327.11  210.62  1.28
34.13  35.54   26.93   1.33  6.44 17.82  0.86 396.96  935.39  225.57  2.32
33.94  35.55   29.23   1.04  4.02 18.57  1.14 329.67  934.56  250.42  3.63
33.61  32.99   32.22   1.17  2.96 17.16  0.95 310.97  714.51  219.83  2.51
32.49  43.72   22.45   0.7   3.53 19.03  0.49 337.95  995.12  225.41  0.85
32.08  32.08   32.17   2.44  4.12 25.22  0.7  361.32  788.1   299.37  1.37
38.34  38.85   22.21   0.42 16.55 25.38  1.12 208.21  680.97  309.87  1.52
41.07  40.92   17.65   0.27  7.57 18.94  0.88 123.56  688.34  237.8   1.31
];
//---------------------------------------------------------------//
////////////////////////////////////////////////////////////////////////////
// Para carregar o arquivo salvo em disco
// use por exemplo:
//---------------------------------------------------------------//
caminho_arquivo=['d:\work2\data\carnes\carnes.dat']);
mtlb_load(caminho_arquivo)
X=carnes;
//---------------------------------------------------------------//
//
// caminho_arquivo - é o caminho do diretório onde esta gravado o arquivo carnes.dat
// o comando mtlb_load carrega o arquivo carnes.dat
// pode se usar  o comando:
// "X=excel2sci(caminho_arquivo,[,sep])
// para carregar arquivos do excel.
// [,sep] é o caracter usado pelo excel para separa dados de colunas
// diferentes.
//
///////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////
//Para carregar os arquivos do PCA faça:
//---------------------------------------------------------------//
caminho=['d:\bin\scilab-2.6\quimiometria']
getf([caminho+'\pca.sci']);
 

///////////////////////////////////////////////////////////////////
// Nome das amostras(labx) e variáveis(laby)
//---------------------------------------------------------------//
labx=['Frankfurter';'Roule';'Ham';'Smoked-ham';'Blanquet';'Smoked-chest';'Hamburger';'Meatball'];
laby=['Sat-fat';'Mono-fat';'Polyunsat-fat';'Non-identi(r)ed';'Calcium';'Magnesium';'Iron';'Phosphorus';'Sodium';'Potassium';'Zinc'];
//---------------------------------------------------------------//
///////////////////////////////////////////////////////////////////

///////////////////////////////////////////////////////////////////
// Para fazer gráfico das amostras
//---------------------------------------------------------------//
amostras=[1,2,3,4];
plot2d(X(amostras',:))
//---------------------------------------------------------------//
// Para fazer gráfico das Variáveis
//---------------------------------------------------------------//
variaveis=[1,2,3,4,5];
plot2d(X(:,variaveis)')
//---------------------------------------------------------------//
///////////////////////////////////////////////////////////////////

/////////////////////////////- PCA -////////////////////////////////
//---------------------------------------------------------------//
[t,p,Var]=pca(X,labx,laby);
//---------------------------------------------------------------//
///////////////////////////////////////////////////////////////////
 

///////////////////////////////////////////////////////////////////
//                          Gráficos
///////////////////////////////////////////////////////////////////

///////////////////////////////////////////////////////////////////

// Observação para Escores e Loadinsg 3D: Após terem sido escolhidos
// os Componentes Principais para os eixos x, y e z, aparecerá uma janela
// com a opção FIM e outra com o gráfico. Este gráfico pode ser girado
// usando a o botão "3D Rot" na janela do gráfico.
// Quanto opção FIM (+ok no windows) for  escolhida será pedido os ângulos
// "alfa" e "theta" mostrados na parte inferior do gráfico.
//  Após a entrada deste ângulos um gráfico com nomes e eixos é gerado.
// Caso queira mudar o ângulo lhe será dada uma nova opção.
// ATENÇÃO "USE APENAS O GRÁFICO SEM NOMES ou NÚMEROS PARA ESCOLHER
// NOVOS ÂNGULOS SE O GRÁFICO FINAL (COM NOMES OU NÚMEROS) FOR GIRADO
// OS NOMES DAS AMOSTRAS NÃO ESTARAM NA POSIÇÃO CORRETA".