Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_factor_analysis.wasp
Title produced by softwareFactor Analysis
Date of computationTue, 08 Aug 2023 15:52:52 +0200
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2023/Aug/08/t1691502873r9bf2pnherm8xfy.htm/, Retrieved Thu, 10 Sep 2026 18:43:26 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=319930, Retrieved Thu, 10 Sep 2026 18:43:26 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsheavy metals, assessment, risk,
Estimated Impact124
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Factor Analysis] [HM assrssment] [2023-08-08 13:52:52] [d41d8cd98f00b204e9800998ecf8427e] [Current]
Feedback Forum

Post a new message
Dataseries X:
"Fe"   "Zn" "Rb" "Sr" "Zr""Ti" "Mn" "Co" "Cu"  "Pb"
37822	68	27	41	218	5524	554	584	50	26
45117	39	0.02	64	288	4913	181	773	55	29
37879	45	18	12	192	4205	210	409	54	19
35803	34	13	40	263	2696	257	499	65	32
5602	30	13	18	495	2798	224	87	57	29
12379	57	19	27	396	3874	386	486	61	27
50080	69	11	76	264	4216	245	294	58	22
30203	59	17	43	254	4305	566	401	48	27
23236	30	14	57	337	4667	149	723	38	22
24517	44	9	88	138	3161	408	88	27	30
15779	42	20	65	373	4702	177	614	35	21
9167	62	23	73	303	5803	422	82	43	30
8107	0.02	14	39	190	3533	334	299	22	34
28424	45	0.02	36	192	4201	285	92	40	19
35612	26	5	48	325	2291	410	409	24	15
45190	22	24	10	218	2434	317	733	38	22
6864	23	20	24	498	3853	190	453	59	20
40127	49	26	53	221	4271	164	282	53	19
47680	49	17	35	187	4367	199	80	53	21
32799	46	18	96	323	2414	302	500	80	14
5124	39	7	70	205	3451	204	297	43	19
11424	50	14	46	308	4283	295	300	50	20
34425	0.02	17	45	284	4969	288	281	47	22
35127	21	8	37	271	3589	421	80	60	17
19617	53	23	100	391	4397	327	580	44	11
28243	53	9	36	247	3668	107	556	41	21
9112	0.02	20	63	125	5117	206	232	33	13
14032	20	10	34	413	5206	161	321	15	25
12204	21	10	38	320	5905	502	103	29	21
10619	0.02	18	75	402	4313	232	300	44	22
23650	0.02	7	16	363	2354	259	548	32	25
10049	49	13	9	253	4405	420	337	33	10
11419	35	9	43	297	3623	253	356	35	15
15813	20	10	42	508	2306	143	340	24	24
16166	27	12	57	225	3881	273	517	49	20
29586	61	18	24	429	4348	272	575	70	19
21461	46	24	23	460	2854	512	290	50	14
46132	38	16	65	333	3737	241	291	70	9
4661	44	19	61	275	4314	189	403	22	20
8465	29	8	31	225	4195	216	448	41	23
22138	0.02	16	60	337	3452	364	398	23	9
37945	19	15	11	387	4681	309	70	29	10
53799	0.02	10	13	280	3051	145	385	39	21
22337	35	28	36	462	5125	325	83	38	16
28520	30	13	52	357	2524	493	606	27	17
21647	0.02	8	66	156	3964	515	88	18	15
13632	37	10	38	328	4215	269	388	25	21
17514	30	0.02	94	341	3909	303	397	21	16
23289	33	17	79	352	1772	163	450	25	20
52695	41	15	34	236	2310	455	250	21	21
17284	17	13	60	325	4184	196	391	30	10
'




Summary of computational transaction
Raw Input view raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R ServerBig Analytics Cloud Computing Center

\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input view raw input (R code)  \tabularnewline
Raw Outputview raw output of R engine  \tabularnewline
Computing time1 seconds \tabularnewline
R ServerBig Analytics Cloud Computing Center \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=319930&T=0

[TABLE]
[ROW]
Summary of computational transaction[/C][/ROW] [ROW]Raw Input[/C] view raw input (R code) [/C][/ROW] [ROW]Raw Output[/C]view raw output of R engine [/C][/ROW] [ROW]Computing time[/C]1 seconds[/C][/ROW] [ROW]R Server[/C]Big Analytics Cloud Computing Center[/C][/ROW] [/TABLE] Source: https://freestatistics.org/blog/index.php?pk=319930&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=319930&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 Input view raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R ServerBig Analytics Cloud Computing Center



Parameters (Session):
par1 = 4 ;
Parameters (R input):
par1 = 4 ;
R code (references can be found in the software module):
par1 <- '2'
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]
print(y)
fit <- principal(y, nfactors=par1, rotate='varimax')
fit
fs <- factor.scores(y,fit)
fs
bitmap(file='test2.png')
plot(fs$scores,pch=20)
text(fs$scores,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')