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 16:59:57 +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/t1293037074rezabzoav28y2tm.htm/, Retrieved Sun, 05 May 2024 23:36:15 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=114406, Retrieved Sun, 05 May 2024 23:36:15 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact152
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Factor Analysis] [bonus - varimax] [2010-12-22 16:59:57] [e569a00cc6e8044e6afea1f18dd335a0] [Current]
Feedback Forum

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




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time3 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135
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 & 3 seconds \tabularnewline
R Server & 'Gwilym Jenkins' @ 72.249.127.135 \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=114406&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]3 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Gwilym Jenkins' @ 72.249.127.135[/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=114406&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=114406&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 time3 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135
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.686-0.452
Pslog-0.427-0.725
Llog-0.669-0.61
wblog0.8780.066
wbrlog0.9020.01
tglog0.8060.283
P-0.0820.937
S0.530.71
D 0.1540.955

\begin{tabular}{lllllllll}
\hline
Rotated Factor Loadings \tabularnewline
Variables & Factor1 & Factor2 \tabularnewline
SWS & -0.686 & -0.452 \tabularnewline
Pslog & -0.427 & -0.725 \tabularnewline
Llog & -0.669 & -0.61 \tabularnewline
wblog & 0.878 & 0.066 \tabularnewline
wbrlog & 0.902 & 0.01 \tabularnewline
tglog & 0.806 & 0.283 \tabularnewline
P & -0.082 & 0.937 \tabularnewline
S & 0.53 & 0.71 \tabularnewline
D
 & 0.154 & 0.955 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=114406&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.686[/C][C]-0.452[/C][/ROW]
[ROW][C]Pslog[/C][C]-0.427[/C][C]-0.725[/C][/ROW]
[ROW][C]Llog[/C][C]-0.669[/C][C]-0.61[/C][/ROW]
[ROW][C]wblog[/C][C]0.878[/C][C]0.066[/C][/ROW]
[ROW][C]wbrlog[/C][C]0.902[/C][C]0.01[/C][/ROW]
[ROW][C]tglog[/C][C]0.806[/C][C]0.283[/C][/ROW]
[ROW][C]P[/C][C]-0.082[/C][C]0.937[/C][/ROW]
[ROW][C]S[/C][C]0.53[/C][C]0.71[/C][/ROW]
[ROW][C]D
[/C][C]0.154[/C][C]0.955[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=114406&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=114406&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.686-0.452
Pslog-0.427-0.725
Llog-0.669-0.61
wblog0.8780.066
wbrlog0.9020.01
tglog0.8060.283
P-0.0820.937
S0.530.71
D 0.1540.955



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')