多元统计分析实验指导书——实验一-均值向量和协方差阵检验word版本

实验一SPSS软件的基本操作与均值向量和协方差阵的检验

【实验目的】

通过本次实验，了解SPSS的基本特征、结构、运行模式、主要窗口等，了解如何录入数据和建立数据文件，掌握基本的数据文件编辑与修改方法,对SPSS有一个浅层次的综合认识。同时能够掌握对均值向量和协方差阵进行检验。

【实验性质】

必修，基础层次

【实验仪器及软件】

计算机及SPSS软件

【实验内容】

1.操作SPSS的基本方法(打开、保存、编辑数据文件)

2.问卷编码

3.录入数据并练习数据相关操作

4.对均值向量和协方差阵进行检验，并给出分析结论。

【实验学时】

4学时

【实验方法与步骤】

1.开机

2.找到SPSS的快捷按纽或在程序中找到SPSS，打开SPSS

3.认识SPSS数据编辑窗、结果输出窗、帮助窗口、图表编辑窗、语句编辑窗

4.对一份给出的问卷进行编码和变量定义

5.按要求录入数据

6.练习基本的数据修改编辑方法

7.检验多元总体的均值向量和协方差阵

8.保存数据文件

9.关闭SPSS，关机。

【实验注意事项】

1.实验中不轻易改动SPSS的参数设置，以免引起系统运行问题。

2.遇到各种难以处理的问题，请询问指导教师。

3.为保证计算机的安全，上机过程中非经指导教师和实验室管理人员同意，禁

止使用移动存储器。

4.每次上机，个人应按规定要求使用同一计算机，如因故障需更换，应报指导

教师或实验室管理人员同意。

5.上机时间，禁止使用计算机从事与课程无关的工作。

【上机作业】

1.定义变量：试录入以下数据文件，并按要求进行变量定义。

表1

要求：

1）变量名同表格名，以“（）”内的内容作为变量标签。对性别（Sex）设值标签“男=0；女=1”。

2）正确设定变量类型。其中学号设为数值型；日期型统一用“mm/dd/yyyy“型号；生活费

用货币型。

3）变量值宽统一为10，身高与体重、生活费的小数位2，其余为0。

2.（1）试录入以下数据文件，保存为“数据1.sav”。

表2

（2）试录入以下数据文件，保存为“数据2.sav”。

（3）试将数据2合并到数据1，合并后的数据文件另存为“数据3.sav”。

（4）将工资进行重编码，2000以下（含2000）为1，2000-3000为2，3000-4000为 3，4000以上为4，重编码的结果保存为“工资等级”。新数据文件保存为“数据4.sav”。

（5）求出各职工刚进入公司时的年龄，保存为“初入年龄”。新数据文件保存为“数据5.sav”。

（6）试按各职员的工资数进行排秩，排秩要求工资最高的排为第一，相同数额取平均等级。排秩后的数据文件保存为“数据6.sav”。

（7）试按各职员的工资数分性别进行排序，要求先排男性，后排女性。同一性别按工资从高到低排列。排序后的数据文件保存为“数据7.sav”。

3.某航空公司38名职员性别和工资情况的调查数据

表4 某航空公司38名职员基本情况调查数据表

表5 某航空公司38名职员收入情况调查数据表

（1）将表1数据输入到SPSS的数据编辑窗口中，将gender定义为字符型变量，保存数据文件，命名为“gender.sav”。

（2）将表2数据输入到SPSS的数据编辑窗口中，将salary定义为数值型变量，保存数据文件，命名为“salary.sav”。

（3）将两个数据文件合成一个数据文件，命名为“Employee Data.sav”。

（4）要求将数据文件“Employee Data.sav”按照变量salary（收入）进行升序排序，并建立一个新数据文件“Employee Data(sorted).sav”放置排序以后的结果。

（5）要求将数据文件“Employee Data.sav”按照变量gender（性别）进行分组，对每一组的变量salary计算其算术平均数，并计算其最大观察值，并建立一个新数据文件“Employee Data(aggregate).sav”放置分类汇总以后的结果。

（6）要求以gender（性别）对数据文件“Employee Data.sav”进行拆分，并要求在以后的统计分析中可以将各拆分文件的统计分析结果放在同一表格中显示。（看不出操作结果，熟悉该操作过程即可）。

（7）要求在“Employee Data.sav”文件中，标识工资在30000元以上的员工。标识变量名设为s_ed变量标签为工资学历标识。并把文件保存为“Employee Data(counted).sav”

4.大学生的素质高低要受各方面因素的影响，其中包括家庭环境与家庭教育（x1）、学校生活环境（x2）、学校周围环境（x3）和个人向上发展的心理动机（x4）等。从某大学在校学生中抽取了 20 人对以上因素在自己成长和发展过程中的影响程度给予评分（以9分制），数据如表6所示：

表6

假定 x=(x1，x2，x3，x4)’服从四元正态分布。

试检验 H0:μ=μ0=(7，5， 4，8)，H1:μ≠μ0(α= 0.05)。

5.测量30名出生到3周岁婴幼儿的身高和体重数据（见表7），其中男女各15名，

（1）假定这两组都服从正态分布且协方差阵相等，试在显著性水平0.05下检验男女婴幼儿的这两项指标是否有显著差异？请将下表的数据输入到SPSS文件中，并进行检验。

（2）检验男性婴幼儿与女性婴幼儿的协差阵是否相等（α=0.05）。

表7

6.1992年美国总统选举的三位候选人为布什、佩罗特和克林顿。从支持三位候选人的选民中分别抽取了20 人，登记他们的年龄段（x1）、受教育程度（x2）和性别（x3）资料如下表所示：

（1）假定三组都服从多元正态分布，检验这三组的总体均值是否有显著性差异（α=0.05）。（2）检验三位候选人的协差阵是否相等（α=0.05）。

盐雾试验作业指导书

盐雾试验箱作业指导书 1．目的 使作业人员掌握该盐雾机的使用方法及维护，确保在测试过程中不因调置不当影响测试之品质，使本厂的产品满足规定要求。 2．适用范围 适用于本厂的盐雾腐蚀试验箱。 3．定义 无 4．职责 4.1．测试员负责基本的测试调节与维护，并作保养记录。 4.2．品质主管负责盐雾腐蚀试验箱的检查。 5．作业指导 5.1试验前的准备 5.1.1.工作室底部应加入自来水，以工作室的水位探测器上端面平为准。 5.1.2．箱体上部四周的密封槽试验前加入自来水，不宜过满，以关闭箱 盖后不外溢为佳。 5.1.3．压力桶内加入蒸馏水水或去离子水，水位以不少于最少水位量， 量为准。 5.1.4．把配制好的盐溶液放入试验机控制箱顶的贮水箱内，盐水由大气 压的作用靠平衡度自动流入喷雾塔内，盐水不能低于贮水箱的下限标 记。

5.1.5．检查试验箱内部漏斗与集雾器之间的连接管是否通畅完好，一定 不能影响盐雾沉降量的收集.在测试区至少放置一个干净的盐雾收集器，要求收集到的不是从测试样品或其它地方滴落的溶液。收集器的位置在样品附近。确保在连续运转至少16h时，在收集器的每80cm2的水平收集面积收集量为1.0-2.0ml/h，收集溶液的pH为6.5-7.2，且在25℃时测得比重为1.0255-1.0400才满足浓度要求。 5.1.6．开启箱盖，将样品正确放置在工作室内的样品架上，样品与样品 之间应保持一定距离，使样品暴露的表面积大于98％，再关上箱盖。 注：每批次的盐雾测试数量1-5PCS但重量不可超过2KG，测试时间按产品要求进行。 5.2控制系统 试验室温度控制器：控制试验室温度，依标准温度设定。 压力桶温度控制器：控制饱和空气桶温度。 计时器：可调式1~999hr可任意设定试验所需时间，终了自动停机，附自动归零装置。 积时器：直读式0~999hr可指示试验时间之累积。 电源开关：触屏式，控制全机之总电源。 5.3、操作说明

盐雾试验箱作业指导书

文件编号 浙江力拓新电气有限公司 版本：A 修订状态：0 LTOX-PZ-YW-02-01-01 盐雾试验箱作业指导书 共2页，第1页 一、 目的： 为了确保使用者正确的使用以及保养此仪器,以便測得正确数据并維持本仪器的寿命. 二、 范围： 品质部检验員. 三、 操作程序： 1. 盐雾试验箱组合: 2. 作业前准备: 2.1.试验室温度设定：35° 2.2.饱和空气桶温度设定：47° 2. 3. 试验机内氯化钠溶液是否足够。 2. 4. 将试样放置在支架。 2. 5. 检视试验槽内必须洁净。 编制/日期: 审核/日期: 批准/日期: 饱和桶 加水管 排雾管 排水管 饱和桶水位计 调压阀 盐水箱 总电源开关 电源 开关 喷雾 开关 周期喷雾开关 周期 计时器 计时器 实验室温度设定 实验室温度显示 压力表 饱和桶温度设定 饱和桶实际温度显示 试验室盖 操作面板 隔水槽

文件编号浙江力拓新电气有限公司版本：A修订状态：0 LTOX-PZ-QC-02-01-01盐雾试验箱作业指导书共2页，第2页 2.6. 检视水槽水位。 2.7. 检视补充液。 3.作业内容: 3.1.测试报告格式依客户要求。 3.2.设定喷雾试验时间。 3.3. 按入电源开关，键灯亮起。 3.4. 按入时间开关，键灯亮起。 3.5. 按入喷雾开关，键灯亮起。 3.6. 操作功能开始执行。 3.7. 调整压缩空气压力为1㎏f/c㎡。 3.8. 操作时间结束时，按下除雾开关，键灯亮起。 3.9. 三分钟后掀开试验槽顶盖，取出试样。 3.10.全部键按出，关闭电源。 3.11.操作期间必须每日检视执行功能，并作成记录. 四、判定： 1.试验周期与判定标准： 零件材料镀覆种类试验持续时间(h) 合格要求(主要表面) 碳钢 锌24 无白色或灰黑腐蚀剂镉96 铜+镍+铬64 铜及铜合金镍+铬96 无棕锈 镍48 无浅绿色或灰色腐蚀物锡48 无灰黑色腐蚀物 银24 无铜绿 铝和铝合金阳极氧化64 无白色腐蚀物 2.试验观察时间为4H一次. 3.以上试验要求按照GB/T10125-1997（NSS）标准执行. 五、注意事项： 1.请依指示之电压安装。 2.禁止安置在高温及有震动的场所。

多元统计分析实验报告

实验一 一、实验目的及要求 对应分析是你也降维的思想以达到减化数据结构的目的，凤的研究广泛用于定义属性变量构成的列联表利用对应分析方法分析问卷中教育程度与网上购物支付方式之间的相互关系。 二、实验环境 SPSS 19.0 window 7系统 三、实验内容及实验步骤（实践内容、设计思想与实现步骤） 实验题目： 通过分析问卷数据，绘制如下的教育程度与网上购物支付方式的交叉表，运用对应分析方法研究教育程度与网上购物所选择的支付方式之间的相关性，及揭示不同人群网上购物的特征等问题。 设计思想：原假设：H1：χ2>χα2[(n?1)(p?1)] 实现步骤： 1.在变量视窗中录入3个变量，用edu表示【教育程度】，用fangshi表示【在网上购物时采用什么样的支付方式】，用pinshu表示【频数】；如图所示：

2.先对数据进行预处理。执行【数据】→【加权个案】命令，弹出【加权个案】对话框。选中【加权个案】按钮，把【频数】放入【频率变量】框中，点击【确定】按钮完成。 3.打开主窗口，选择菜单栏中的【分析】→【降维】→【对应分析】命令，弹出【对应分析】对话框。 4.将【教育程度】导入【行】，将【在网上购物时采用什么样的支付方式】导入【列】。 5. 单击【定义范围（D）】，打开【对应分析：定义行范围】对话框； 定义行变量分类全距最小值为1，最大值为4，单击【更新】；点击【继续】，返回【对应分析】对话框；同方法打开【对应分析：定义列范围】对话框； 定义列变量全距最小值为1，最大值为5，单击【更新】； 6. 单击【统计量】打开【对应分析：统计量】对话框；选择【行轮廓表】，【列轮廓表】；单击【继续】，返回【对应分析】对话框， 7.选择【绘制】→【对应分析：图】对话框，选择【散点图】中的【行点】、【列点】选择【线图】中的【已转换的行类别】、【已转换的列类别】，单击【继续】，返回【对应分析】对话框。 8.单击【确定】按钮，完成设置并执行列联表分析。 四、调试过程及实验结果（详细记录实验在调试过程中出现的问题及解决方法。记录实验的结果） SPSS实验结果及分析： 上表显示了在32155名被调查者中，大多数消费者在网上购物时选择第三方支付和网上银行支付，在网上购物的消费人群以大学本科生相对最多。

盐雾试验箱操作规程

盐雾试验箱操作规程 盐雾试验箱型号H/YW-90A 设备编号160111M6-3 一实验前准备 1.工作室底部加热水槽内加足蒸馏水或纯净水，以防加热时箱体开 裂、老化，准备充足的纯净水。 2.箱体上部四周水密封槽应加入适量的蒸馏水或矿泉水，不宜过满 切勿太少，以关闭箱盖时水和烟雾不能外溢为准。 3.检查储水箱和工作室内喷雾器之间的水管连接是否完好配置好。 4.检查盐雾箱后面的储物管是否连接完好以免影响盐雾排放。 二操作步骤 1.加水 （1）打开面板电源开关，箱底低水位和压力桶低水位指示灯在点亮状态，首先往试验箱加入纯净水直至试验箱低水位指示灯熄灭后，打开压力桶开关往压力桶加入纯净水直至压力桶低水位指示灯 熄灭，将开关关闭。 （2）将配置好的5%盐溶液加入盐水箱内，注意加盐水要小心不能将盐水外溅，防止对电气元器件有损坏。 2.调节压力 打开面板喷雾开关，使空气压缩机处于工作状态，拉出盐雾箱 后面的压力调节阀旋转调节压力，显示时按下旋钮。拉出面板压 力调节阀阀旋转调节压力，显示时按下旋钮。 3.设定温度

按动面板温度调节箭头设定工作湿度为35°，压力桶温度设定为47°设定好后不宜经常改变。红色显示实际温度，实际温度是不可改变的。 4.时间设定 打开计时开关计时器显示时间，上下按键可调整时间，根据实验要求设定试验时间，右边为时间累积计时器。 5.样品的摆放 打开试验箱盖将样品摆放在样品架上，样品要倾斜35°摆放，样品之间要留有一定距离，使样品暴露的表面积大于98%，盖上箱盖打开喷雾开关试验机开始工作。 6.补水 在工作的过程中，盐雾箱压力桶内的水都有一定的消耗，如发现低水位指示灯亮起时需对其进行补充纯净水直至低水位灯熄灭。压力桶补水时注意先往压力桶加满水再打开开关，防止内部热水喷出。 7.巡视检查 实验工作中巡视检查试验箱的运行情况，特别注意查看烟雾沉降量是否正常（正常每小时平均1-2毫升），若偏大或偏小随时调节挡板角度和喷雾压力直至正常为止，巡查监视水位及时补水。 8.液体排放 试验按规定的时间自然停机后，及时关闭气源和电源，使试验箱各个系统处于静止状态，实验结束时先打开除雾开关让雾气尽

盐雾测试标准

表面处理用中性盐雾试验作业指导书 1.适用范围:本标准规定为各项金属底材于电镀后、有机或无机涂装后等各项表面处理用之盐水喷雾耐触性试验方法. 2.试验方法:本法是使用盐水喷雾试验机将氯化纳溶液,以雾状喷于被覆膜上之一种腐蚀试验方法.试验的主要条件如表1所示. 表1 主要的试验条件 表2 连续雾化时间表 3.试验液之配制:溶解氯化纳(化学纯、分析纯）于蒸馏水（或总溶解固体量小于200 ppm以下的水中）,调配成浓度为5±% 的试验液.此试验液在35℃喷雾后,其收集液PH值应为2.且喷雾前,此试验液不能含有浮态物. 注（1）氯化钠不能含有铜兴镍的不纯物,固体内的碘化钠含量小于.因为不纯物中可能含有腐蚀抑制剂,所以不纯物总含量须小于.

（2）在33~35℃间测量比试验液的比重应为,在25℃测量时的比重则为.此试验液的浓度亦可用硝酸银液滴定法或其它方法检定. （3）试验液须以试药级的盐酸或氢氧化钠稀溶液调整PH值,并以PH仪或其它可靠方法测量之,由于配制试验液的水中含有二氧化碳,二氧化碳在水中的溶解度随温度改变而影响溶液的PH值,故须小心控制PH值,PH值则可依下列任一方法调整: ①常温配制试验液,于35℃喷雾,因为温度的升高而使部份二氧化碳逸出溶液而升高PH值.故在常温配制试验液 时,PH值应调整在内,才可以使收集液的PH值在之间. ②PH值调整前,使试验液先煮沸再冷至35℃,或维持在35℃温度48小时.如此调整的PH值在喷雾时,将不会产生 太大的变化. ③先将水加热至35℃以上,以去除溶解的二氧化碳,而后再调制试验液并调整PH值,如此在35℃喷雾时,所调整的 PH值也不会产生太大的变化. （4）为避免喷雾嘴阻塞,此试验液须过滤或小心倾斜注入盐水桶,或于喷雾吸水管前端处装上玻璃筛过滤,盐雾测试仪每工作2000小时,需更换玻璃喷嘴. 4.设备:本试验所需之设备为喷雾嘴、盐水桶、试验片支架、喷雾液收集容器、试验室、盐水[补给桶、压力桶、压缩空气之 供给设备与排气设备等所构成,其装置如图1所示,并依照如下条件试验. 盐水喷雾试验机与其所需的管路应采用钝性材料,不能对喷雾之腐蚀试验有影响或本身被腐蚀者. 喷雾嘴不可直接将试验液喷向试样,喷雾室顶部聚集之溶液不得滴落在试验片上. 试验片滴下之试验液不可流回盐水桶,而再用于试验. 压缩空气不能含有油脂及灰尘,所以须有空气清净器.空气压力须保持在±因为压缩空气于膨胀时,有吸热现象,所以须事先有预热以增加压缩空气的温度与湿度,如附表2,以获取均一温度的喷雾. 喷雾采取器其水平采取面积为80cm2,直径约为10 cm,置于试验片附近（靠近喷嘴最近与最远之处的两个地方）. 喷雾液量以整个时间计算,在采取容器上,应每小时平均可收集至之盐水溶液.喷雾液至少应收集16小时,以其平均值表不喷雾量. 试验中盐水桶,其氯化钠溶液浓度应维持40~60g/L,试验用盐液在1星期内没用完时,需作更换,如下次试验时间超过1星期时,需将试验室内部清洗,将试验盐液,加热水槽之水排出. 压力桶之温度须保持在47±1℃,盐水桶之温度在35±1℃. 试验室的相对湿度须保持在85%以上,更高的相对湿度的要求可与客户双方协议确定. 5.试样 采取位置:试样可从制品之主要表面采取或以制品本身做为试样.但如无法以制品试验或判定时,可以用试验片代替且此试验片必须能代表被测试制品. 尺度:试验片之标准尺寸为150×70mm,或为100×65mm. 数目:试样之数目3-5PCS. 试验前处理 试样依镀层性质与清洁程度而须做适当清洗,不能使研磨剂及具有腐蚀或抑制作用之溶剂清洗,且清洗方法不得损伤表面.至于不锈钢试样使用硝酸清洗与钝化,试样经清净后(通过拔水性试验),以干净之布或吸水将水份擦干,或用无油之干燥空气吹干.不得已时,可使用氧化镁膏.此膏状物乃10g试药级氧化镁加入100mL的蒸馏水中. 若非另有规定,否则试样切口及因挂色而造成底材露出部份,或因识别记号所造成镀层缺陷处于试验时应覆以合适之保护层.如硬蜡(ceresin wax),乙烯胶带等绝缘物. 手纹污染会造成严重的不良试验结果,试样于清净后不得有任何手纹的污染. 6.试样之置放:试验中,试验室内之试样的位置符合以下列条件: 试样之主要表面与铅垂线成15度至30度之倾斜,又自试验室上方俯视时应与喷雾之主要流动方向平行.特殊零件具有很多方向的主要表面,需要同时测试时,可取多件试样置放,务必使每个主要表面能同时接受盐水的喷雾. 试样之排列,应使喷雾自由地落至全部试片上,不应防碍喷雾之自由落下. 试样不可互相接触,也不可接触到金属性导体或有毛细现象作用之物质,以及其它支持架以外之物. 盐水溶液勿从一个试样滴流至其它试样上. 试样识别记号或装配孔应覆于下方.

多元统计分析实验报告

1. 正态性检验 Kolmogorov-Smirnov a Shapir o-Wilk 统计量df Sig. 统计量df Sig. 净资产收益率.113 35 .200*.978 35 .677 总资产报酬率.121 35 .200*.964 35 .298 资产负债率.086 35 .200*.962 35 .265 总资产周转率.180 35 .006 .864 35 .000 流动资产周转率.164 35 .018 .885 35 .002 已获利息倍数.281 35 .000 .551 35 .000 销售增长率.103 35 .200*.949 35 .104 资本积累率.251 35 .000 .655 35 .000 *. 这是真实显著水平的下限。 a. Lilliefors 显著水平修正 此表给出了对每一个变量进行正态性检验的结果，因为该例中样本中n=35<2000，所以此处选用Shapiro-Wilk统计量。由Sig.值可以看到，总资产周转率、流动资产周转率、已获利息倍数及资本积累率均明显不遵从正态分布，因此，在下面的分析中，我们只对净资产收益率、总资产报酬率、资产负债率及销售增长率这四个指标进行比较，并认为这四个变量组成的向量遵从正态分布（尽管事实上并非如此）。这四个指标涉及公司的获利能力、资本结构及成长能力，我们认为这四个指标可以对公司运营能力做出近似的度量。 2. 主体间因子 N 行业电力、煤气及水的生产和供应 业 11 房地行业15 信息技术业9 多变量检验a 效应值 F 假设 df 误差 df Sig. 截距Pillai 的跟踪.967 209.405b 4.000 29.000 .000 Wilks 的 Lambda .033 209.405b 4.000 29.000 .000 Hotelling 的跟踪28.883 209.405b 4.000 29.000 .000 Roy 的最大根28.883 209.405b 4.000 29.000 .000 行业Pillai 的跟踪.481 2.373 8.000 60.000 .027 Wilks 的 Lambda .563 2.411b8.000 58.000 .025 Hotelling 的跟踪.698 2.443 8.000 56.000 .024 Roy 的最大根.559 4.193c 4.000 30.000 .008 a. 设计 : 截距 + 行业

盐雾试验作业指导书

盐雾试验作业指导书文档编制序号：[KK8UY-LL9IO69-TTO6M3-MTOL89-FTT688]

1.目的 正确使用盐雾试验箱，以检测镀层的耐腐蚀性能。 2. 适用范围 可根据客户要求对电镀零件和产品进行盐雾试验。 3. 职责 质保部是电镀产品检验和试验的归口管理部门； 电镀检验人员负责对电镀产品进行盐雾试验的整个过程。 4. 操作程序 设备: 经过鉴定符合有关标准的盐水喷雾试验机 JS-HY-SST-bo-01c 盐雾试验溶液的配制: 公升的纯净水倒入专用的塑料桶中，用PH试纸测试其PH值是否在之间； PH值若小于则加入少量的氢氧化钠，PH值若大于则加入少量的冰醋酸；加入250g氯化钠NaCl，搅拌均匀。 样品: 20个试样 操作程序: 1. 将自动加水的入水口阀门排水阀和排气阀的开关打开。

2. 将隔绝水槽加水至垫板位置。 3. 将配制好的氯化钠盐水倒入到盐水补充槽，即自动充填盐水进入试验箱内的预热槽，使盐水流至盐水预热槽。 4. 加少许水在湿球杯内，湿球温度覆盖着纱布，纱布末端置于湿球杯内。 5. 开始试验前，试样必须充分清洗，清洗方法视表面情况及污物的性质而定，不能使用任何 会侵蚀试样表面的磨料或溶剂，同时试样切口及因挂钩而造成底材露出部分，或因识别记号所造成的镀层缺陷处，试验前应用透明胶带将以覆盖。放置试样或试片于置物架上，试样在箱内放置的位置，应使受试平板试样与垂直线成15-30°角，试样的主要表面向上，并与盐雾在箱内流动的主要方向平行。特殊试样有很多的主要表面需要同时测试时，可取多件试样置放，务必使每个主要表面能同时进行盐雾试验。 6. 试验时，试样之间不得互相接触，也不与箱壁相碰，试样的间距一般不小于20mm，试样上下层必须交叉放置，试样间间隔应能使盐雾自由沉降在试样的主要表面上。一个试样上的盐水溶液不得滴在任何别的试样上。试样识别记号或装配孔应覆于下方。 7. 设定试验温度、压力和时间： 将盐水桶和试验室的温度调整至35℃，压力桶温度调整至47℃（按“+”为增加，按“—”为减少，H：时/M：分/S：秒）。喷压压力保持在±cm2，若压力

可靠性测试规范之盐雾试验作业指导书

核准： 审核： 作成：袁媛 盐雾试验作业指导书 版 本 B0 制订部门 品质部 页次 1/7 生效日期 2020.05.11 1. 0目的 指导作业，规范操作，提升试验结果的客观性及可信赖性。 2. 0范围： 所有需要盐雾测试的产品。 3. 0定义： 盐雾试验：利用盐水喷雾腐蚀来检验和鉴定电镀层封孔性之好坏，以及对镀层耐腐蚀性和对基 体保护性能的测试；或试样无表面处理时本身耐腐蚀的能力。 4. 0权责 品管部负责取样、测试、判定。 5.0设备、药品及操作条件 5.1盐水喷雾试验机 5.2氯化钠(分析纯)溶液（5%）、溶液使用纯水配制，紧急时可使用纯净水替代。 5.3操作条件 项 目 试 验 中 备 注 盐水质量百分比浓度(%) 5±0.1 盐水不得重复使用 盐水PH 值 6.5-7.2 测定收集的盐雾溶液 压缩空气压力(kgf/cm 2) 1.00±0.1 经过现场校验和认证 喷雾量(ml/80cm 2/hr) 1.5+0.5 连续不得中断，至少8H 以上 压力桶温度 47±20C 试验室温度和湿度 35±20C,90% RH 以上 样品放置角度 15°-25° 附角度参照图 试验时间(hr) 参考本文件7.0条款 6.0 试验 6.1试样准备： 在试样准备以及试验结束取样观察全过程中，不可裸手接触试样，应全程戴一次性手套或 手指套，以保护试件电镀面不被汗渍及其它外来物污染。在用手套或手指套防护下,将镍片用双面胶粘在治具上,单个产品间距不少于20mm 。 6.2试样摆放： （1）试样不应摆放在盐雾直接喷射到的位置。

核准： 审核： 作成：袁媛 盐雾试验作业指导书 版 本 B0 制订部门 品质部 页次 2/7 生效日期 2020.05.11 （2）在盐雾试验箱中被试面与垂直方向成15°~ 25°,并尽可能成20°,对于不规则的试样， 例如整个工件都是被试面，也应尽可能接近上述规定。 －带材测试：带材对折30°~50°之间，垂直放置在盐雾箱内。 20° 40° －镍片测试：借助辅助治具，确保试样被试面与垂直方向成15°~ 25° （3）试样可以摆放在试验箱不同水平面上，但不能接触箱体，也不能相互接触，单个试 样件间距不得小于20mm 。试样或其支架上的滴液不得落在其他试样上。 6.3 试验后试样的处理： （1）试验结束后取出试样,用温度不高于40℃的清洁流动水轻轻清洗以除去试样表面残留 的盐雾溶液,接着在距离试样约300mm 处用气压不超过300kPa 的空气立即吹干.或者 清洗后用无尘布轻轻吸干试件表面水份. 电镀面 错误摆放方式试件表面水渍印 试件不可纵向放 置上下间会滴液 L 型电镀面 电镀面

多元统计分析实验报告,计算协方差矩阵,相关矩阵,SAS

院系：数学与统计学学院 专业：__统计学 年级：2009 级 课程名称：统计分析 ____ 学号：____________ 姓名：_________________ 指导教师：____________ 2012年4月28日 （一）实验名称 1. 编程计算样本协方差矩阵和相关系数矩阵；

2. 多元方差分析MANOVA。 （二）实验目的 1. 学习编制sas程序计算样本协方差矩阵和相关系数矩阵; 2. 对数据进行多元方差分析。 （三）实验数据 第一题： 第二题:

（四）实验内容 1. 打开SAS软件并导入数据； 2. 编制程序计算样本协方差矩阵和相关系数矩阵; 3. 编制sas程序对数据进行多元方差分析； 4. 根据实验结果解决问题，并撰写实验报告； （五）实验体会（结论、评价与建议等） 第一题： 程序如下： proc corr data=sasuser.sha n cov; proc corr data=sasuser.sha n no simple cov; with x3 x4; partial x1 x2; run; 结果如下： （1）协方差矩阵 $AS亲坯 曲；15 Friday, Apr： I SB,沙DO COUR过程 x4 目由度=30 Xi x2x3x4x5X? -10.I9B4944-0.45E2GJ5I.3347097-G.1193E48-￡0.e75?GS

-ID. 188494669,36&Q3?9-7.22IO&OS1J5692043I5.49ee^91S.Oa97SM -8.45S2645■7,221050829.S78&S46-6.372E47I-15.3084183-21.7352376-11.5674785 1.3841097 1.G5S2M7t.3726171IJ24?17B 4.e093011 4.4C12473 2.B747CM -G. I1S3S49 1.GS92043-is.soul aa 4.B09B01I68.7978495劣』S670971S.57ai1B3 -IH.05l6l?a15.43S6569-J1.73S2376孔耶124TB27.0387097105.103225&S7.3505S7E： -2D K5752??319-11337204-1L55M7S52r9747?3i19,573118337.3S0&87E33.3SQ6452 (2) 相关系数矩阵 Pearson相关系数” N =引 当HO： Rho=0 时.Prob > |r| Xi Xi xl 1.QQ000 x2 -C.23954 0.2061 x3 -0,30459 0.0957 x4 0.18975 Q.3092 x5 '0.14157 0.4475 x6 -0.83787 0.0630 -0.49292 0.0150 x2-0.23354 1.00000-0.162750.143510.022700.181520.24438 x20.20C10.31:1?0.441?0.90350.32640.1761 x3-0.30459-0.16275 1.00000-0.06219-0.34641-0.^797-0.23674 x30.095?0.381?<.00010.0563o.oses0 JS97 x40.1S8760.14351-0.86219L000000.400540,313650.22610 x40.30920.4412<.0001 D.02EG Q.085S0.2213 x5-0J 41570.02270-0.946410.40054 1.000000.317370.26750 x50.4J750.90350.0G68Q.025&0.08130 + 1620 x6-0.33?e?0.1S162-0.397970.813650.31787LOOOOO0.82976 x60.0S300.32840.02660.08580.0813C0001辺-0.432920.24938-0.288740.22810 D.267600.92976 1.00000 x70,01500J7610.19970.22130JG20<.0001 第二题： 程序如下： proc anova data=sasuser.hua ng; class kind; model x1-x4=k ind; manova h=k ind; run; 结果如下： (1)分组水平信息 The ANNA Procedure Cla^s Level Informat ion Class Level?Values kind 3 123 Number of observatIons CO (2) x1、x2、x3、x4的方差分析

第2章 多元正态分布均值向量和协差阵的检验

第一章 多元正态分布的参数估计 一、填空题 1.设X 、Y 为两个随机向量，对一切的u 、v ，有 ，则称X 与Y 相互独立。 2.多元分析处理的数据一般都属于 数据。 3．多元正态向量()' =p X X X ,,1 的协方差阵∑是 ，则X 的各分量是相互独立的随机变量。 4．一个p 元函数() p x x x f ,,,21 能作为p R 中某个随机向量的密度函数的主要条件是 和 。 5．若p 个随机变量1X ，2X ， ，p X 的联合分布等于 ，则称1X ， 2X ， ，p X 是相互独立的。 6.多元正态分布的任何边缘分布为 。 7.若()∑,~μp N X ，A 为p s ?阶常数阵，d 为s 维常数向量，则 ~d AX + 。 8.多元正态向量X 的任何一个分量子集的分布称为X 的 。 9.多元样本中，不同样品的观测值之间一定是 。 10.多元正态总体均值向量和协差阵的极大似然估计量分别是 。 11.多元正态总体均值向量μ和协差阵∑的估计量X 、S n 1 1 -具有 、 和 。 12．设X 和S 分别是多元正态总体()∑,μp N 的样本均值向量和离差阵，则 ~X ，X 和S 。 13.若()()∑,~μαp N X ，n ,,2,1 =α且相互独立，则样本离差阵 ()()()()∑=' --=n X X X X S 1~ααα 。 14．若()∑,~i p i n W S ，k i ,,1 =，且相互独立，则~21k S S S S +++= 。 二、判断题 1.多元分布函数()x F 是单调不减函数，而且是右连续的。 2.设X 是p 维随机向量，则X 服从多元正态分布的充要条件是：它的任何组合 ()p R X ∈'αα都是一元正态分布。 3.μ是一个P 维的均值向量，当A 、B 为常数矩阵时，具有如下性质： （1）E （AX ）=AE （X ） （2）E （AXB ）=AE （X ）B 4．若P 个随机变量X 1，…X P 的联合分布等于各自边缘分布的乘积，则称X 1，… X P 是相互独立的。 5．一般情况下，对任何随机向量()'=X X X p ,,1 ，协差阵∑是对称阵，也 是正定阵。 6．多元正态向量( )' =X X X p ,,1 的任意线性变换仍然服从多元正态分布。 7．多元正态分布的任何边缘分布为正态分布，反之一样。 8．多元样本中，不同样品之间的观测值一定是相互独立的。 9．多元正态总体参数均值μ的估计量X 具有无偏性、有效性和一致性。 10． S n 1 是∑的无偏估计。 11.Wishart 分布是2 χ分布在p 维正态情况下的推广。

盐雾试验箱安全操作规程

盐雾试验箱安全操作规程 一、开机 接通烟雾机的外供电源和气源→检查压力桶和试验箱排水阀是否关闭、检查进气压力表是否达到2Kg→开启烟雾机操作面板上电源、操作开关→观察控制版面指示是否有低水位和低盐水亮红灯及报警声→出现时要进行加水和加盐水来消除亮灯及报警故障（溶解试药级氯化钠于蒸馏水（或总溶解固体量小于200ppm以下的水中），调配成浓度为5±1%的试验液）→将需测试物品放置到试验箱内（试样之主要表面与铅垂线成15度至30度之倾斜）、同时调整试验产品需测试时间→调整计时器→开启控制版面计时开关、开启喷雾开关→控制版面压力表迅速动作（喷雾压力需设定为1Kg）→盐雾机喷雾工作开始； 二、设备运行 设备运行时：保障电源、水源、气源正常； 设备运行时：压力桶之温度须保持在47±1℃，盐水桶之温度在35±1℃； 设备运行时：试验室的相对温度须保持在85%； 设备运行时：喷雾液量以整个时间计算，在采取容器上，应每小时平均可收集1.0至2.0ml之盐水溶液。喷雾液至少应收集16小时，以其平均值表求喷雾量。 三、其他 停机：先关闭控制版面喷雾开关→开启除雾开关、待实

验箱雾出除干净后才能打开实验箱盖； 保养：每月清洗一次压力桶、实验箱、盐水桶及喷嘴；1周以上时间停用机器需全面清洗干净；（喷嘴拆卸时要轻拿轻放，并用清水冲洗干净） 盐雾试验箱安全操作规程 四、开机 接通烟雾机的外供电源和气源→检查压力桶和试验箱排水阀是否关闭、检查进气压力表是否达到2Kg→开启烟雾机操作面板上电源、操作开关→观察控制版面指示是否有低水位和低盐水亮红灯及报警声→出现时要进行加水和加盐水来消除亮灯及报警故障（溶解试药级氯化钠于蒸馏水（或总溶解固体量小于200ppm以下的水中），调配成浓度为5±1%的试验液）→将需测试物品放置到试验箱内（试样之主要表面与铅垂线成15度至30度之倾斜）、同时调整试验产品需测试时间→调整计时器→开启控制版面计时开关、开启喷雾开关→控制版面压力表迅速动作（喷雾压力需设定为1Kg）→盐雾机喷雾工作开始； 五、设备运行 设备运行时：保障电源、水源、气源正常； 设备运行时：压力桶之温度须保持在47±1℃，盐水桶之温度在35±1℃；

多元统计分析实验报告

多元统计分析实验报告 Document serial number【UU89WT-UU98YT-UU8CB-UUUT-UUT108】

1. 正态性检验 Kolmogorov-Smirnov a Shapir o-Wilk 统计量df Sig.统计量df Sig. 净资产收益 .11335.200*.97835.677 率 总资产报酬 .12135.200*.96435.298 率 资产负债率.08635.200*.96235.265 总资产周转 .18035.006.86435.000 率 流动资产周 .16435.018.88535.002 转率 已获利息倍 .28135.000.55135.000 数 销售增长率.10335.200*.94935.104 资本积累率.25135.000.65535.000 *. 这是真实显着水平的下限。 a. Lilliefors 显着水平修正 此表给出了对每一个变量进行正态性检验的结果，因为该例中样本中 n=35<2000，所以此处选用Shapiro-Wilk统计量。由Sig.值可以看到，总资产周转率、流动资产周转率、已获利息倍数及资本积累率均明显不遵从正态分布，因此，在下面的分析中，我们只对净资产收益率、总资产报酬率、资产负债率及销售增长率这四个指标进行比较，并认为这四个变量组成的向量遵从正态分布（尽管事实上并非如此）。这四个指标涉及公司的获利能力、资本结构及成长能力，我们认为这四个指标可以对公司运营能力做出近似的度量。 2. 主体间因子 N

行业电力、煤气及水的 生产和供应业 11 房地行业15 信息技术业9 多变量检验a 效应值F假设 df 误差 df Sig. 截距Pillai 的跟 踪 .967.000 Wilks 的 Lambda .033.000 Hotelling 的跟踪 .000 Roy 的最大 根 .000 行业Pillai 的跟 踪 .481.027 Wilks 的 Lambda .563.025 Hotelling 的跟踪 .698.024 Roy 的最大 根 .559.008 a. 设计 : 截距 + 行业 b. 精确统计量 c. 该统计量是 F 的上限，它产生了一个关于显着性级别的下 限。 上面第一张表是样本数据分别来自三个行业的个数。第二张表是多变量检验表，该表给出了几个统计量，由Sig.值可以看到，无论从哪个统计量来看，三个行业的运营能力（从净资产收益率、总资产报酬率、资产负债率及销售增长率这四个指标的整体来看）都是有显着差别的。 3. 主体间效应的检验

盐雾试验作业指导书

1.目的 正确使用盐雾试验箱，以检测镀层的耐腐蚀性能。 2. 适用范围 可根据客户要求对电镀零件和产品进行盐雾试验。 3. 职责 3.1 质保部是电镀产品检验和试验的归口管理部门； 3.2 电镀检验人员负责对电镀产品进行盐雾试验的整个过程。 4. 操作程序 4.1设备: 经过鉴定符合有关标准的盐水喷雾试验机JS-HY-SST-bo-01c 4.2盐雾试验溶液的配制: 4.2.1. 调制方法：将4.75公升的纯净水倒入专用的塑料桶中，用PH试纸测试其PH值是否在6.5-7.2之间； 4.2.2. PH值若小于6.5则加入少量的氢氧化钠，PH值若大于7.2则加入少量的冰醋酸； 4.2.3. 加入250g氯化钠NaCl，搅拌均匀。 4.3样品: 20个试样 4.4操作程序: 1. 将自动加水的入水口阀门排水阀和排气阀的开关打开。 2. 将隔绝水槽加水至垫板位置。 3. 将配制好的氯化钠盐水倒入到盐水补充槽，即自动充填盐水进入试验箱内的预热槽，使盐水流至盐水预热槽。 4. 加少许水在湿球杯内，湿球温度覆盖着纱布，纱布末端置于湿球杯内。 5. 开始试验前，试样必须充分清洗，清洗方法视表面情况及污物的性质而定，不能使用任何 会侵蚀试样表面的磨料或溶剂，同时试样切口及因挂钩而造成底材露出部分，或因识别记号所造成的镀层缺陷处，试验前应用透明胶带将以覆盖。放置试样或试片于置物架上，试样在箱内放置的位置，应使受试平板试样与垂直线成15-30°角，试样的主要表面向上，并与盐雾在箱内流动的主要方向平行。特殊试样有很多的主要表面需要同时测试时，可取多件试样置放，务必使每个主要表面能同时进行盐雾试验。 6. 试验时，试样之间不得互相接触，也不与箱壁相碰，试样的间距一般不小于20mm，试样上下层必须交叉放置，试样间间隔应能使盐雾自由沉降在试样的主要表面上。一个试样上的盐水溶液不得滴在任何别的试样上。试样识别记号或装配孔应覆于下方。 7. 设定试验温度、压力和时间： 将盐水桶和试验室的温度调整至35℃，压力桶温度调整至47℃（按“+”为增加，按“—”为减少，H：时/M：分/S：秒）。喷压压力保持在1.00±0.01kgf/cm2，若压力不在范围内，可利用调压阀将压力调整至规定范围（顺时针为增加，逆时针为减少）。测试时间一般为24小时（按“△”为增加，按“▽“为减少），若客户有特殊要求则可另行设定，测试时间一般可设定为8、16、24、48、96、168、336、672小时，在规定的试验周期内喷雾不得中断，只有当需要短暂观察试样时才能打开盐雾箱，开箱检查的时间和次数应尽可能少。

盐雾测试标准-

？ 表面处理用中性盐雾试验作业指导书 1.适用范围:本标准规定为各项金属底材于电镀后、有机或无机涂装后等各项表面处理用之盐水喷雾耐触性试验方法. 2.试验方法:本法是使用盐水喷雾试验机将氯化纳溶液,以雾状喷于被覆膜上之一种腐蚀试验方法.试验的主要条件如表1所示. 表1 主要的试验条件 表2 连续雾化时间表

3.试验液之配制:溶解氯化纳(化学纯、分析纯）于蒸馏水（或总溶解固体量小于200 ppm以下的水中）,调配成浓度为5 ±%的试验液.此试验液在35℃喷雾后,其收集液PH值应为~.且喷雾前,此试验液不能含有浮态物. 注（1）氯化钠不能含有铜兴镍的不纯物,固体内的碘化钠含量小于%.因为不纯物中可能含有腐蚀抑制剂,所以不纯物总含量须小于% . （2）在33~35℃间测量比试验液的比重应为~,在25℃测量时的比重则为~.此试验液的浓度亦可用硝酸银液滴定法或其它方法检定. （3）试验液须以试药级的盐酸或氢氧化钠稀溶液调整PH值,并以PH仪或其它可靠方法测量之,由于配制试验液的水中含有二氧化碳,二氧化碳在水中的溶解度随温度改变而影响溶液的PH值,故须小心控制PH值,PH值则可依下列任一方法调整: ①常温配制试验液,于35℃喷雾,因为温度的升高而使部份二氧化碳逸出溶液而升高PH值.故在常温配制试验液 时,PH值应调整在内,才可以使收集液的PH值在~之间. ②PH值调整前,使试验液先煮沸再冷至35℃,或维持在35℃温度48小时.如此调整的PH值在喷雾时,将不会产生 太大的变化. \ ③先将水加热至35℃以上,以去除溶解的二氧化碳,而后再调制试验液并调整PH值,如此在35℃喷雾时,所调整的 PH值也不会产生太大的变化. （4）为避免喷雾嘴阻塞,此试验液须过滤或小心倾斜注入盐水桶,或于喷雾吸水管前端处装上玻璃筛过滤,盐雾测试仪每工作2000小时,需更换玻璃喷嘴. 4.设备:本试验所需之设备为喷雾嘴、盐水桶、试验片支架、喷雾液收集容器、试验室、盐水[补给桶、压力桶、压缩空气 之供给设备与排气设备等所构成,其装置如图1所示,并依照如下条件试验. 盐水喷雾试验机与其所需的管路应采用钝性材料,不能对喷雾之腐蚀试验有影响或本身被腐蚀者. 喷雾嘴不可直接将试验液喷向试样,喷雾室顶部聚集之溶液不得滴落在试验片上. 试验片滴下之试验液不可流回盐水桶,而再用于试验. 压缩空气不能含有油脂及灰尘,所以须有空气清净器.空气压力须保持在±cm2,因为压缩空气于膨胀时,有吸热现象,所以须事先有预热以增加压缩空气的温度与湿度,如附表2,以获取均一温度的喷雾. 喷雾采取器其水平采取面积为80cm2,直径约为10 cm,置于试验片附近（靠近喷嘴最近与最远之处的两个地方）. { 喷雾液量以整个时间计算,在采取容器上,应每小时平均可收集至之盐水溶液.喷雾液至少应收集16小时,以其平均值表不喷雾量. 试验中盐水桶,其氯化钠溶液浓度应维持40~60g/L,试验用盐液在1星期内没用完时,需作更换,如下次试验时间超过1星期时,需将试验室内部清洗,将试验盐液,加热水槽之水排出. 压力桶之温度须保持在47±1℃,盐水桶之温度在35±1℃. 试验室的相对湿度须保持在85%以上,更高的相对湿度的要求可与客户双方协议确定. 实验工作中,实验盖边口有实验雾气漏出时,需向密封水槽内加水,试验期间如超过1个月,需将热水槽内水更换. 5.试样 采取位置:试样可从制品之主要表面采取或以制品本身做为试样.但如无法以制品试验或判定时,可以用试验片代替且此试验片必须能代表被测试制品. 尺度:试验片之标准尺寸为150×70mm,或为100×65mm.

多元统计分析实验报告doc

多元统计与程序设计》课程实验报告 项目名称： 学生姓名： 学生学号： 指导教师： 完成日期：

1 实验内容 2 模型建立与求解 2.1聚类分析的形成思路 2.2.1类平均法 2.2.2谱系图的形成 2.3.快速聚类法 (以上内容见课本) 3 实验数据与实验结果 3.1实验数据 设有20个土壤样品分别对5个变量的观测数据如表5.16所示，试利用 聚类法对其进行样品聚类分析 样品号 含沙量1X 淤泥含量2X 粘土含量3X 有机物4X PH 值5X 1 77.3 13.0 9.7 1.5 6.4 2 82.5 10.0 7.5 1.5 6.5 3 66.9 20.0 12.5 2.3 7.0 4 47.2 33.3 19.0 2.8 5.8 5 65.3 20.5 14.2 1.9 6.9 6 83.3 10.0 6.7 2.2 7.0 7 81.6 12.7 5.7 2.9 6.7 8 47.8 36.5 15.7 2.3 7.2 9 ４8.6 37.1 14.3 2.1 7.2 10 61.6 25.5 12.6 1.9 7.3 11 58.6 26.5 14.9 2.4 6.7 12 69.3 22.3 8.4 4.0 7.0 13 61.8 30.8 7.4 2.7 6.4 14 67.7 25.3 7.0 4.8 7.3 15 57.2 31.2 11.6 2.4 6.3 16 67.2 22.7 10.1 33.3 6.2 17 59.2 31.2 9.6 2.4 6.0 18 80.2 13.2 6.6 2.0 5.8

19 82.2 11.1 6.7 2.2 7.2 20 69.7 20.7 9.6 3.1 5.9 3.2实验过程及结果 Case Processing Summary(a) Cases Valid Missing Total N Percent N Percent N Percent 20 100.0% 0 .0% 20 100.0% a Squared Euclidean Distance used 上表是接近度矩阵，计算距离使用的是平方欧氏距离，所以样品间距离越大，样品越相异，由表中矩阵可以看出样品8号和样品9号的距离是最小的，因此它们最先聚为一类。 Average Linkage (Between Groups) Agglomeration Schedule Stage Cluster Combined Coefficient s Stage Cluster First Appears Next Stage Cluster 1 Cluster 2 Cluster 1 Cluster 2 1 8 9 .153 16

多元统计分析实验报告

多元统计分析实验报告 1、实验内容 根据课本习题3-12做相关分析。 2、实验目的 （1）检验H0：；H1：协方差阵不全相等。 （2）检验H0： U1=U2 ; H1：U1≠U2; （3）检验H0： U1=U2 =U3 ; H1：U1,U2,U3不全等; （4）检验三种化学成分相互独立。 3、实验方案分析 （1）这是关于判断三个3元正态总体的协方差阵是否相等的问题； (2)均值是否相等，在两个协方差阵相等的情况下均值是否相等的问题； （3）比较三组的3项指标是否有差异的问题，就是多总体均值向量是否相等的检验问题； （4）检验 是否独立相当于检验任意2个子向量的协方差阵是否为零矩阵； 4、实验原理及操作过程，结果如下： （1）SAS 代码实现过程如下： data d3121; input y1-y3 group @@; cards; 47.22 5.06 0.10 1 1 23 ==∑∑∑

47.45 4.35 0.15 1 47.52 6.85 0.12 1 47.86 4.19 0.17 1 47.31 7.57 0.18 1 54.33 6.22 0.12 2 56.17 3.31 0.15 2 54.40 2.43 0.22 2 52.62 5.92 0.12 2 43.12 10.33 0.05 3 42.05 9.67 0.08 3 42.50 9.62 0.02 3 40.77 9.68 0.04 3 ; proc iml; n1=5;n2=4;n3=4; n=n1+n2+n3;k=3; p=3; use d3121(obs=5); xa={y1 y2 y3 }; read all var xa into x1; print x1; use d3121(firstobs=6 obs=9); read all var xa into x2; print x2; use d3121(firstobs=10 obs=13); read all var xa into x3; print x3; xx=x1//x2//x3; ln={[5] 1} ; x10=(ln*x1)/n1; print x10; mm1=i(n1)-j(n1,n1,1)/n1; mm=i(n)-j(n,n,1)/n; a1=x1`*mm1*x1; print a1; ln={[4] 1} ; x10=(ln*x2)/n2; print x20;