juliohm

srand(2000) # make sure we get the same data

#############
# LOGNORMAL #
#############

using SpecialFunctions

# CDF and inverse CDF of standard normal
ϕ(x) = erfc(-x/√2) / 2

juliohm

function [B,L,var,fac,E,r,issuccess] = J_CoDaPFA(X,s)
% FA Factor analysis (principal factoring technique)
%
% The reason to favour this method opposite to maximum likelihood 
% (Statistics toolbox) is no distributional assumptions. In comparison to
% SPSS program, this provides the same results. We expect identical
% numerical solution, especially advantage of singular value decomposition.
% Initial estimates of communalities are set as squared multiple
% correlations. In iteration process (convergence criterium 0.001, limit
% for iteration 75), communalities are set to one if Heywood case occurs. 

juliohm

x,y,Au,Ag,As,Ba,Be,Bi,Cd,Co,Cr,Cu,Mn,Nb,Ni,P,Pb,Rb,Sb,Sn,Ti,V,W,Zn,Zr,Si,Al,Fe,Mg,Ca,Na,K
28.09,28.306,0.8397,0.16358,12.557,640,2.566,0.41059,0.15289,10.268,54.98,19.511,586.6,14.133,24.322,578.3,39.246,115.4,1.1638,2.8783,3892.5,74.208,1.7011,68.816,321.2,6.393e+05,1.429e+05,41400,11700,15800,20200,29210
21.736,34.54,0.5309,0.097287,8.7926,610.4,2.4512,0.28667,0.11157,8.9944,48.4,17.046,543.9,11.779,19.467,598.1,23.337,110,0.78,2.287,3200.6,61.677,1.51,52.25,311.7,6.588e+05,1.224e+05,29800,12100,32800,25600,31230
30.464,26.078,2.6924,0.11951,14.807,535.3,2.4506,0.35381,0.10856,11.728,62.8,22.931,562.6,13.049,33.857,556,24.873,100.2,1.3248,2.32,3874.8,89.5,1.8328,65.283,254.6,5.702e+05,1.317e+05,49100,18100,48200,11100,23100
30.617,27.809,1.428,0.14994,15.708,565,2.4934,0.27281,0.15581,9.3256,50.9,16.907,466.8,10.346,21.71,786.3,28.206,104.8,0.91655,2.01,3111.5,63.4,1.4129,66.145,269.9,5.846e+05,1.354e+05,45900,17500,34700,13400,25100
0.76984,23.904,0.6721,0.06459,5.456,521.5,2.2487,0.2226,0.08863,5.5854,34.

juliohm

x,y,Au,Ag,As,Ba,Be,Bi,Cd,Co,Cr,Cu,Mn,Nb,Ni,P,Pb,Rb,Sb,Sn,Ti,V,W,Zn,Zr,Si,Al,Fe,Mg,Ca,Na,K
28.09,28.306,0.8397,0.16358,12.557,640,2.566,0.41059,0.15289,10.268,54.98,19.511,586.6,14.133,24.322,578.3,39.246,115.4,1.1638,2.8783,3892.5,74.208,1.7011,68.816,321.2,6.393e+05,1.429e+05,41400,11700,15800,20200,29210
21.736,34.54,0.5309,0.097287,8.7926,610.4,2.4512,0.28667,0.11157,8.9944,48.4,17.046,543.9,11.779,19.467,598.1,23.337,110,0.78,2.287,3200.6,61.677,1.51,52.25,311.7,6.588e+05,1.224e+05,29800,12100,32800,25600,31230
30.464,26.078,2.6924,0.11951,14.807,535.3,2.4506,0.35381,0.10856,11.728,62.8,22.931,562.6,13.049,33.857,556,24.873,100.2,1.3248,2.32,3874.8,89.5,1.8328,65.283,254.6,5.702e+05,1.317e+05,49100,18100,48200,11100,23100
30.617,27.809,1.428,0.14994,15.708,565,2.4934,0.27281,0.15581,9.3256,50.9,16.907,466.8,10.346,21.71,786.3,28.206,104.8,0.91655,2.01,3111.5,63.4,1.4129,66.145,269.9,5.846e+05,1.354e+05,45900,17500,34700,13400,25100
0.76984,23.904,0.6721,0.06459,5.456,521.5,2.2487,0.2226,0.08863,5.5854,34.

juliohm

juliohm

     32335:	binding file /usr/lib/libc.so.6 [0] to /usr/lib/libc.so.6 [0]: normal symbol `_res' [GLIBC_2.2.5]
     32335:	binding file /usr/lib/libc.so.6 [0] to /usr/lib/libc.so.6 [0]: normal symbol `stderr' [GLIBC_2.2.5]
     32335:	binding file /usr/lib/libc.so.6 [0] to /usr/lib/libc.so.6 [0]: normal symbol `error_one_per_line' [GLIBC_2.2.5]
     32335:	binding file /usr/lib/libc.so.6 [0] to /usr/lib/libc.so.6 [0]: normal symbol `realloc' [GLIBC_2.2.5]
     32335:	binding file /usr/lib/libc.so.6 [0] to /usr/lib/libc.so.6 [0]: normal symbol `__morecore' [GLIBC_2.2.5]
     32335:	binding file /usr/lib/libc.so.6 [0] to /usr/lib/libc.so.6 [0]: normal symbol `__key_encryptsession_pk_LOCAL' [GLIBC_2.2.5]
     32335:	binding file /usr/lib/libc.so.6 [0] to /usr/lib/libpthread.so.0 [0]: normal symbol `__libpthread_freeres'
     32335:	binding file /usr/lib/libc.so.6 [0] to /usr/lib/libc.so.6 [0]: normal symbol `__progname_full' [GLIBC_2.2.5]
     32335:	binding file /usr/lib/libc.so.6 [0] to /usr/lib/libc.so.6 [0]: 

juliohm

	linux-vdso.so.1 (0x00007fff7f291000)
	libclangAST.so.7 => /usr/lib/../lib/../lib/libclangAST.so.7 (0x00007f7dc2d43000)
	libclangBasic.so.7 => /usr/lib/../lib/../lib/libclangBasic.so.7 (0x00007f7dc2a0e000)
	libclangLex.so.7 => /usr/lib/../lib/../lib/libclangLex.so.7 (0x00007f7dc293f000)
	libclangSema.so.7 => /usr/lib/../lib/../lib/libclangSema.so.7 (0x00007f7dc21a4000)
	libLLVM-7.so => /usr/lib/../lib/../lib/libLLVM-7.so (0x00007f7dbe1db000)
	libstdc++.so.6 => /usr/lib/../lib/../lib/libstdc++.so.6 (0x00007f7dbe04c000)
	libc.so.6 => /usr/lib/../lib/../lib/libc.so.6 (0x00007f7dbde86000)
	libgcc_s.so.1 => /usr/lib/../lib/../lib/../lib/libgcc_s.so.1 (0x00007f7dbde6c000)
	libclangAnalysis.so.7 => /usr/lib/../lib/../lib/../lib/libclangAnalysis.so.7 (0x00007f7dbddbe000)

juliohm

9     /buildworker/worker/package_linux64/build/src/gf.c:2184; jl_apply_generic
1     /buildworker/worker/package_linux64/build/src/jltypes.c:635; typekey_compare
30140 /buildworker/worker/package_linux64/build/src/task.c:268; start_task
 30140 /buildworker/worker/package_linux64/build/src/julia.h:1536; jl_apply
  30140 /buildworker/worker/package_linux64/build/src/gf.c:2182; jl_apply_generic
   30140 /buildworker/worker/package_linux64/build/src/gf.c:1829; jl_fptr_trampoline
    30140 ./task.jl:259; (::getfield(IJulia, Symbol("##12#15")))()
     30140 /home/juliohm/.julia/packages/IJulia/DL02A/src/eventloop.jl:8; eventloop(::ZMQ.Socket)
      30140 ./essentials.jl:685; invokelatest
       30140 ./essentials.jl:686; #invokelatest#1

juliohm

C%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
C                                                                      %
C Copyright (C) 2003, Statios Software and Services Incorporated.  All %
C rights reserved.                                                     %
C                                                                      %
C This program has been modified from the one distributed in 1996 (see %
C below).  This version is also distributed in the hope that it will   %
C be useful, but WITHOUT ANY WARRANTY. Compiled programs based on this %
C code may be redistributed without restriction; however, this code is %
C for one developer only. Each developer or user of this source code   %

juliohm