Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_factor_analysis.wasp
Title produced by softwareFactor Analysis
Date of computationWed, 22 Dec 2010 00:56:49 +0000
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2010/Dec/22/t1292979354v7n91q8wetirss8.htm/, Retrieved Mon, 06 May 2024 06:57:43 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=114041, Retrieved Mon, 06 May 2024 06:57:43 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact179
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Factor Analysis] [bonus - tabel 1] [2010-12-22 00:56:49] [54d0a09f418287eaabab8ba43e4b06f8] [Current]
Feedback Forum

Post a new message
Dataseries X:
1	NA	NA	3,65	15,71	15,56	6,47	3,00	5,00	3,00
2	6,30	0,69	1,50	6,91	8,79	3,74	3,00	1,00	3,00
3	NA	NA	2,64	8,13	10,70	4,09	1,00	1,00	1,00
4	NA	NA	NA	-0,08	8,65	3,22	5,00	2,00	3,00
5	2,10	0,59	4,23	14,75	15,34	6,44	3,00	5,00	4,00
6	9,10	-0,36	3,30	9,26	12,10	5,19	4,00	4,00	4,00
7	15,80	1,36	2,94	-3,77	-1,20	3,56	1,00	1,00	1,00
8	5,20	0,00	3,41	11,98	12,04	5,97	4,00	5,00	4,00
9	10,90	1,28	3,33	8,10	10,15	4,14	1,00	2,00	1,00
10	8,30	0,34	3,91	10,86	12,99	5,44	1,00	1,00	1,00
11	11,00	0,41	1,95	-0,86	8,76	4,72	5,00	4,00	4,00
12	3,20	-0,36	3,40	13,05	12,96	5,64	5,00	5,00	5,00
13	7,60	0,99	NA	-0,60	7,78	NA	2,00	1,00	2,00
14	NA	NA	3,69	12,14	12,95	5,90	5,00	5,00	5,00
15	6,30	0,74	1,25	-2,59	7,09	3,74	1,00	1,00	1,00
16	8,60	NA	3,91	8,01	10,13	3,33	2,00	2,00	2,00
17	6,60	1,41	1,79	-0,24	8,16	3,74	2,00	2,00	2,00
18	9,50	0,18	2,34	-1,61	8,52	4,79	2,00	2,00	2,00
19	4,80	0,26	3,53	7,25	9,77	NA	1,00	2,00	1,00
20	12,00	1,81	1,95	11,00	11,30	NA	1,00	1,00	1,00
21	NA	-1,20	3,33	13,18	13,43	5,99	5,00	5,00	5,00
22	3,30	-0,69	3,00	10,23	11,65	5,00	5,00	5,00	5,00
23	11,00	1,22	1,36	-2,12	6,91	2,77	3,00	1,00	2,00
24	NA	NA	3,67	12,24	12,91	5,53	1,00	4,00	1,00
25	4,70	0,41	3,71	11,35	12,69	5,74	1,00	3,00	1,00
26	NA	NA	2,79	10,50	11,69	4,14	1,00	1,00	1,00
27	10,40	1,22	2,20	-2,29	8,29	3,33	5,00	1,00	3,00
28	7,40	-0,22	2,03	6,95	8,61	4,22	5,00	3,00	4,00
29	2,10	-0,22	3,83	13,16	13,39	5,82	5,00	5,00	5,00
30	NA	NA	3,11	11,51	11,96	4,61	1,00	1,00	1,00
31	NA	NA	2,79	10,46	10,93	3,50	3,00	5,00	4,00
32	7,70	0,34	0,96	-5,30	-1,97	3,07	5,00	2,00	4,00
33	17,90	0,69	3,18	-4,61	-1,39	3,91	1,00	1,00	1,00
34	6,10	0,64	4,61	11,03	14,09	5,59	1,00	1,00	1,00
35	8,20	0,88	NA	-2,10	8,01	3,40	2,00	1,00	1,00
36	8,40	1,03	NA	7,21	9,00	3,81	3,00	1,00	3,00
37	11,90	0,26	1,16	-3,77	-0,92	2,94	4,00	1,00	3,00
38	10,80	0,69	0,69	-3,04	-1,11	3,40	4,00	1,00	3,00
39	13,80	1,72	1,61	7,44	8,75	2,48	2,00	1,00	1,00
40	14,30	1,13	1,87	8,16	9,29	4,79	2,00	1,00	1,00
41	NA	0,00	3,16	12,43	13,10	6,09	5,00	5,00	5,00
42	15,20	0,59	2,48	-0,73	9,65	4,94	2,00	2,00	2,00
43	10,00	-0,11	3,01	9,21	11,65	5,14	4,00	4,00	4,00
44	11,90	0,59	2,56	7,39	9,34	2,83	2,00	1,00	2,00
45	6,50	0,64	3,30	12,17	12,10	4,74	4,00	4,00	4,00
46	7,50	-0,11	2,89	7,82	9,40	3,43	5,00	5,00	5,00
47	NA	NA	2,62	8,36	10,58	4,14	2,00	2,00	2,00
48	10,60	0,96	1,55	-1,27	7,55	3,04	3,00	1,00	3,00
49	7,40	0,88	2,28	8,35	10,83	3,95	1,00	1,00	1,00
50	8,40	0,18	3,37	8,82	12,10	5,10	2,00	3,00	2,00
51	5,70	-0,11	1,95	-0,29	9,42	5,42	2,00	2,00	2,00
52	4,90	-0,69	1,79	8,19	9,95	5,42	3,00	2,00	3,00
53	NA	NA	2,83	9,60	11,49	5,01	5,00	5,00	5,00
54	3,20	-0,51	3,00	10,92	12,07	5,02	5,00	5,00	5,00
55	NA	NA	2,54	7,24	9,43	4,50	2,00	2,00	2,00
56	8,10	0,79	1,25	-2,81	6,91	NA	3,00	1,00	2,00
57	11,00	0,83	1,50	-0,11	7,86	4,09	2,00	1,00	2,00
58	4,90	-0,69	2,01	7,60	9,42	5,30	3,00	1,00	3,00
59	13,20	0,96	0,83	-2,26	7,82	3,83	3,00	2,00	2,00
60	9,70	-0,51	3,18	8,34	10,97	5,35	4,00	3,00	4,00
61	12,80	1,89	1,10	8,16	8,27	2,64	2,00	1,00	1,00
62	NA	NA	2,56	8,31	9,74	3,64	3,00	1,00	1,00




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time21 seconds
R Server'Sir Ronald Aylmer Fisher' @ 193.190.124.24
R Framework error message
The field 'Names of X columns' contains a hard return which cannot be interpreted.
Please, resubmit your request without hard returns in the 'Names of X columns'.

\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input & view raw input (R code)  \tabularnewline
Raw Output & view raw output of R engine  \tabularnewline
Computing time & 21 seconds \tabularnewline
R Server & 'Sir Ronald Aylmer Fisher' @ 193.190.124.24 \tabularnewline
R Framework error message & 
The field 'Names of X columns' contains a hard return which cannot be interpreted.
Please, resubmit your request without hard returns in the 'Names of X columns'.
\tabularnewline \hline \end{tabular} %Source: https://freestatistics.org/blog/index.php?pk=114041&T=0

[TABLE]
[ROW][C]Summary of computational transaction[/C][/ROW]
[ROW][C]Raw Input[/C][C]view raw input (R code) [/C][/ROW]
[ROW][C]Raw Output[/C][C]view raw output of R engine [/C][/ROW]
[ROW][C]Computing time[/C][C]21 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Sir Ronald Aylmer Fisher' @ 193.190.124.24[/C][/ROW]
[ROW][C]R Framework error message[/C][C]
The field 'Names of X columns' contains a hard return which cannot be interpreted.
Please, resubmit your request without hard returns in the 'Names of X columns'.
[/C][/ROW] [/TABLE] Source: https://freestatistics.org/blog/index.php?pk=114041&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=114041&T=0

As an alternative you can also use a QR Code:  

The GUIDs for individual cells are displayed in the table below:

Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time21 seconds
R Server'Sir Ronald Aylmer Fisher' @ 193.190.124.24
R Framework error message
The field 'Names of X columns' contains a hard return which cannot be interpreted.
Please, resubmit your request without hard returns in the 'Names of X columns'.







Rotated Factor Loadings
VariablesFactor1Factor2
SWS-0.573-0.475
PSln-0.399-0.738
Lln0.850.034
wbln0.8840.125
wbrln0.8850.082
tgln0.7990.322
P-0.120.944
S0.5290.734
D 0.1180.967

\begin{tabular}{lllllllll}
\hline
Rotated Factor Loadings \tabularnewline
Variables & Factor1 & Factor2 \tabularnewline
SWS & -0.573 & -0.475 \tabularnewline
PSln & -0.399 & -0.738 \tabularnewline
Lln & 0.85 & 0.034 \tabularnewline
wbln & 0.884 & 0.125 \tabularnewline
wbrln & 0.885 & 0.082 \tabularnewline
tgln & 0.799 & 0.322 \tabularnewline
P & -0.12 & 0.944 \tabularnewline
S & 0.529 & 0.734 \tabularnewline
D
 & 0.118 & 0.967 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=114041&T=1

[TABLE]
[ROW][C]Rotated Factor Loadings[/C][/ROW]
[ROW][C]Variables[/C][C]Factor1[/C][C]Factor2[/C][/ROW]
[ROW][C]SWS[/C][C]-0.573[/C][C]-0.475[/C][/ROW]
[ROW][C]PSln[/C][C]-0.399[/C][C]-0.738[/C][/ROW]
[ROW][C]Lln[/C][C]0.85[/C][C]0.034[/C][/ROW]
[ROW][C]wbln[/C][C]0.884[/C][C]0.125[/C][/ROW]
[ROW][C]wbrln[/C][C]0.885[/C][C]0.082[/C][/ROW]
[ROW][C]tgln[/C][C]0.799[/C][C]0.322[/C][/ROW]
[ROW][C]P[/C][C]-0.12[/C][C]0.944[/C][/ROW]
[ROW][C]S[/C][C]0.529[/C][C]0.734[/C][/ROW]
[ROW][C]D
[/C][C]0.118[/C][C]0.967[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=114041&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=114041&T=1

As an alternative you can also use a QR Code:  

The GUIDs for individual cells are displayed in the table below:

Rotated Factor Loadings
VariablesFactor1Factor2
SWS-0.573-0.475
PSln-0.399-0.738
Lln0.850.034
wbln0.8840.125
wbrln0.8850.082
tgln0.7990.322
P-0.120.944
S0.5290.734
D 0.1180.967



Parameters (Session):
par1 = 2 ;
Parameters (R input):
par1 = 2 ;
R code (references can be found in the software module):
library(psych)
par1 <- as.numeric(par1)
x <- t(x)
nrows <- length(x[,1])
ncols <- length(x[1,])
y <- array(as.double(x[1:nrows,2:ncols]),dim=c(nrows,ncols-1))
colnames(y) <- colnames(x)[2:ncols]
rownames(y) <- x[,1]
y
fit <- principal(y, nfactors=par1, rotate='varimax')
fit
fs <- factor.scores(y,fit)
fs
bitmap(file='test1.png')
fa.diagram(fit)
dev.off()
bitmap(file='test2.png')
plot(fs,pch=20)
text(fs,labels=rownames(y),pos=3)
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Rotated Factor Loadings',par1+1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Variables',1,TRUE)
for (i in 1:par1) {
a<-table.element(a,paste('Factor',i,sep=''),1,TRUE)
}
a<-table.row.end(a)
for (j in 1:length(fit$loadings[,1])) {
a<-table.row.start(a)
a<-table.element(a,rownames(fit$loadings)[j],header=TRUE)
for (i in 1:par1) {
a<-table.element(a,round(fit$loadings[j,i],3))
}
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable.tab')