多元统计分析作业.docx
《多元统计分析作业.docx》由会员分享,可在线阅读,更多相关《多元统计分析作业.docx(15页珍藏版)》请在冰豆网上搜索。
多元统计分析作业
海洋地球化学多元统计分析作业
一、预备工作:
数据的输出管理
首先设置Fileoutputmanageroutputmanager中,选中individualwind。
AlsosendtoReportwind中,选中singlereport。
二、数据的导入
数据表(data.xls)为一个深海沉积物柱中30个样品分析结果。
第1列为样品编号,第2列为样品的采样深度(单位m),第三列起为分析的各元素含量。
将data.xls数据导入Statisticaworksheet中(操作步骤为菜单Fileopen…data.xls)
三、数据(图表)的输出
统计分析过程中生成的结果都可以输出到Word文档中(菜单FileSaveas…或PrtSc,粘贴到word中)。
对生成的图表,还可先菜单FileAddtoreport,再粘贴到word中。
本项上机实习需完成以下统计分析
一、相关及回归分析(Correlationmatrices)
1、分析两组分Co-Ni,CaO-Sr,Fe2O3-MnO,的相关关系,做出相关关系图,拟合出回归方程。
图1Co-Ni相关关系图
图2CaO-Sr相关关系图
图3Fe2O3-MnO相关关系图
2、做出三组分Cu-Pb-Zn;Sr-Cu-CaO之间的散点图(scatterplot)。
图4Cu-Co-Ni散点图
图5Sr-Cu-CaO散点图
3、计算CaO、Co、Cu、Fe2O3、MnO、Ni、Sr之间的相关关系矩阵。
表1沉积物中元素相关关系矩阵(n=30,p<0.05)
CaO
Fe2O3
MnO
Co
Cu
Ni
Sr
CaO
1.00
Fe2O3
-0.23
1.00
MnO
0.18
0.18
1.00
Co
-0.21
0.85
0.41
1.00
Cu
-0.02
-0.01
0.36
0.26
1.00
Ni
-0.10
0.96
0.24
0.88
-0.03
1.00
Sr
0.97
-0.25
0.23
-0.20
0.09
-0.13
1.00
二、聚类分析(Clusteranalysis)
1、首先将数据进行标准化(分别进行和列的标准化),得到标准化的数据集。
2、采用Treeclustering方式,Singlelinkage法,对CaO、Co、Cu、Fe2O3、MnO、Ni、Sr进行R型聚类分析,进行统计分析。
图6R型聚类分析图
注:
采用Treeclustering方式,Singlelinkage,1-Pearsonr法。
3、采用Treeclustering方式,Singlelinkage法,对各样品进行Q型聚类分析。
图7Q型聚类分析图
注:
采用Treeclustering方式,Singlelinkage,Euclideandistances法。
三、因子分析(Clusteranalysis)
对元素含量进行R型因子分析
1、步骤同聚类分析,若已对数据的列进行标准化,此步则跳过。
2、最小特征值(Eigenvalue)设置为1。
3、提取方法选用主成分提取法(Principalcomponent)
4、方差最大化正交旋转(Varimaxraw)
5、计算因子得分(Factorscore)
表2因子得分表
样品
因子1
因子2
因子3
1
-0.388923
0.28564
-1.32310
2
-0.606544
1.14243
-0.84065
3
-0.356931
0.08620
-0.80998
4
-0.049097
1.25338
-0.42622
5
-0.297932
-0.64463
-0.41052
6
0.178868
-0.66344
-2.59907
7
-0.314067
1.39519
-0.71048
8
-0.002022
1.17808
-0.47965
9
-0.353060
1.56557
0.27071
10
-0.619335
0.18952
-0.23563
11
-0.242359
0.43259
0.10879
12
-0.305155
-2.28624
-0.37271
13
-0.432760
0.47655
-0.23289
14
-0.018512
0.44646
-0.04254
15
-0.325063
-3.10883
-0.33579
16
-0.371207
0.76179
-0.34949
17
-0.510880
-0.52992
-0.60566
18
-0.567731
-0.83834
0.09357
19
0.030113
0.06677
0.19415
20
-0.000682
0.22434
1.04253
21
-0.569282
-0.66443
0.33776
22
-0.243250
-0.16570
-0.30501
23
-0.026224
-0.21573
0.22697
24
0.458336
0.53100
1.64123
25
0.069827
-0.10426
0.94841
26
-0.159663
0.22337
0.65828
27
4.997148
-0.07215
-0.98686
28
0.250117
0.46571
1.73580
29
-0.110119
-0.97015
1.42804
30
0.886389
-0.46076
2.38003
6、计算特征值(Eigenvalue)
表3特征值表
因子
特征值
变量(%Total)
累计特征值
累积变量
1
3.309234
41.36543
3.309234
41.36543
2
2.319511
28.99389
5.628746
70.35932
3
1.292839
16.16048
6.921584
86.51980
注:
提取方法选用主成分提取法。
7、绘出特征值图(Screeplot)
图8特征值图
8、绘出各因子之间的因子载荷2D及3D图(Factorloading)
图9F1–F2因子载荷2D图
图10F1–F3因子载荷2D图
图11F2–F3因子载荷2D图
图12F1–F2–F3因子载荷3D图
9、列出因子分析结果表。
表4因子分析结果表
元素
因子1
因子2
因子3
depth(m)
0.393400
-0.243300
0.730897
CaO%
-0.094961
-0.982883
-0.010813
Fe2O3%
0.966560
0.146257
0.006360
MnO
0.278851
-0.264602
0.601193
Co
0.894960
0.137751
0.331305
Cu
-0.064306
0.056627
0.917321
Ni
0.987898
0.017264
-0.009756
Sr
-0.118201
-0.975624
0.112618
Expl.Var
2.970781
2.090982
1.859820
Prp.Totl
0.371348
0.261373
0.232478
注:
提取方法选用主成分提取法,方差最大化正交旋转。
四、数据分析
结合你在海洋地球化学中学习的有关知识,将聚类分析得到的树形图谱和因子分析得到的因子分析结果进行综合分析,探讨其代表的地球化学意义,想一下为什么会形成特定的元素或样品组合,该组合反映的是什么地球化学作用过程?
1、R型聚类分析
对CaO、Co、Cu、Fe2O3、MnO、Ni、Sr进行R型聚类分析,得到了如图6所示的谱系图。
从图中可以看出:
在联接距离为6的相似性水平上,谱系图可分为两个分支:
其中,子群1为CaO和Sr的组合,子群2为Cu、Fe2O3、MnO等元素的组合。
在联接距离为4的相似性水平上,谱系图可分为三个分支:
其中,子群1为MnO,子群2为Cu和depth(m)的组合,子群3为Co、Fe2O3、Ni的组合。
2、Q型聚类分析
在Q型分析中,以欧几里德距离为相似性变量标准,得到了如图7所示的谱系图。
从图中可以看出:
在联接距离为0.3—1.2尺度的相似性水平上,谱系图可以分为两支:
其中,子群1为27号样品,子群2为其余所有样品。
27号样品的独立分支,与其他样品相区别,这反映了27号样品在来源或成因上存在显著差异。
3、因子分析
在R型因子分析中,选用Varimax正交旋转,特征值>1的因子被提取,计算结果如表2所示。
从表中可以看出:
3个因子所解释的变化信息达到了86.52%,可以反映出大部分的变量信息。
●从表4可以看出:
①因子1是由正相关变量Fe2O3、Ni、Co组成。
由于铁、镍、钴的相关性较好,体现了铁的氧化物或氢氧化物对镍、钼、钴的吸附作用。
铁的氧化物或氢氧化物是有效的吸附剂,它可以吸附或以离子置换形式容纳多种金属元素,尤以吸附Co、Ti、Ni、Sn等为特征。
一般为自生沉积,多富集在粘土中。
所以因子1可能主要代表自生沉积作用。
②因子2是由负相关变量CaO、Sr组成。
由于CaO、Sr是生物作用的标志性元素,两者在生物作用影响下,多在生源沉积物中存在。
所以因子2可以代表生物作用。
③因子3由正相关变量Cu、depth(m)组成。
由于锰的氧化物或氢氧化物也是有效的吸附剂,以吸附Cu、Zn等元素为主。
在不同的环境下,铁锰的吸附作用存在差异,在还原环境下和快速沉积的背景下,以铁的吸附作用为主。
而在氧化环境下,锰的吸附作用要较强。
此样品为南海北部陆坡沉积,水深为3000米左右,因而铜的富集量受锰的影响相对较小。
在因子3作用下,Cu、MnO表现出相近正值,可见在因子3影响下Cu、MnO表现出较好的相关性,推测因子3可能代表火山作用或与火山热液活动有关。
●从表2可以看出:
样品27上的因子得分为4.997184,而其他因子在该样品上的因子得分数值较小,所以认为样品27主要受控于因子1。
而其他样品各因子的得分差异不显著,可能同时受控于3个因子。
样品27与其他样品的因子得分有显著差异,这与图7所示的聚类分析结果是相对应的。
综合以上分析显示研究区以生物成因组合,火山源成因组合,自生源沉积组合主。
其中火山源成因组合主要为Cu、MnO等;自生源成因组合主要为Fe2O3、Ni、Co组成等;生物成因组合主要为CaO、Sr等等。
铁锰氢氧化物对于许多金属元素是最有效的吸附剂,它们可以吸附或以离子置换形式容纳多种金属元素。
Mn矿物以吸附Ni、Cu、Zn为主,Fe矿物以吸附Co、V、Ti等为主。
钙锶主要受生物作用,富集在碳酸盐沉积物中。
研究区沉积富Ca、Sr、Cu、Co,是因为有大量生物沉积和自生沉积及火山热液活动。
Ni、Cu、Co与自生沉积和火山热液活动有关,构成锰结核的重要组成。
综合以上分析内容,可粗略推测出研究区有较多的自生沉积,组成自生的铁锰结核。
表明了研究区元素具有半深海沉积物的地球化学特征。
(注:
素材和资料部分来自网络,供参考。
请预览后才下载,期待你的好评与关注!
)