统计建模与R软件试验报告

本文格式为Word版，下载可任意编辑——统计建模与R软件试验报告

--话说直接copy是不道德的哈

开课学院、试验室：数学与统计学院试验时间：2023年3月日

--话说直接copy是不道德的哈

n-5;x-array(0,dim=c(n,n))for(iin1:n){+for(jin1:n){+x[i,j]-1/(i+j-1)+}+};x[,1][,2][,3][,4][,5][1,]1.00000000.50000000.33333330.25000000.2000000[2,]0.50000000.33333330.25000000.20000000.1666667[3,]0.33333330.25000000.20000000.16666670.1428571[4,]0.25000000.20000000.16666670.14285710.1250000[5,]0.20000000.16666670.14285710.12500000.1111111det(x)[1]3.749295e-12solve(x)[,1][,2][,3][,4][,5][1,]25-3001050-1400630[2,]-3004800-1890026880-12600[3,]1050-1890079380-11760056700[4,]-140026880-117600179200-88200[5,]630-1260056700-8820044100eigen(x)$values[1]1.567051e+002.085342e-011.140749e-023.058980e-043.287929e-06$vectors[,1][,2][,3][,4][,5][1,]0.76785470.6018715-0.21421360.047161810.006173863[2,]0.4457911-0.27591340.7241021-0.43266733-0.116692747[3,]0.3215783-0.42487660.12045330.667350440.506163658[4,]0.2534389-0.4439030-0.30957400.23302452-0.767191193[5,]0.2098226-0.4290134-0.5651934-0.557599950.376245545

分析：从试验结果来看。R软件在处理数据上相当确凿，便利。

教师签名年月日

--话说直接copy是不道德的哈

开课学院、试验室：数学与统计学院试验时间：2023年3月日

--话说直接copy是不道德的哈

10111213141516171819

AlfredDukeGuidoJames

M14M14M15M12

69.0112.563.5102.567.0133.057.383.062.584.059.099.572.0150.064.8128.057.585.066.5112.0

JeffreyM13JohnPhilipRobertThomasM12M16M12M11

WilliamM15

五、试验结果及实例分析student-read.table(3.7数据.txt,header=T)attach(student)cor.test(身高,体重)#Pearson相关性检验

Pearson'sproduct-momentcorrelation

data:身高and体重t=7.5549,df=17,p-value=7.887e-07alternativehypothesis:truecorrelationisnotequalto095percentconfidenceinterval:0.70443140.9523101sampleestimates:cor0.8777852其p值7.887e-070.05,拒绝原假设，所以身高与体重相关

教师签名年月日

--话说直接copy是不道德的哈

开课学院、试验室：数学与统计学院试验时间：2023年月日

--话说直接copy是不道德的哈

SSta(n1),Xta(n1)nn四、试验环境(所用软件、硬件等)及试验数据文件

X

数据见试验内容，所用软件：R2.15.1五、试验结果及实例分析在R软件中运行代码：x-c(54,67,68,78,70,66,67,70,65,69)t.test(x)#做单样本正态分布区间估计OneSamplet-testdata:xt=35.947,df=9,p-value=4.938e-11alternativehypothesis:truemeanisnotequalto095percentconfidenceinterval:63.158571.6415sampleestimates:meanofx67.4####平均脉搏点估计为67.4,95%置信度的区间估计为[63.1585,71.6415]。t.test(x,alternative=less,mu=72)#做单样本正态分布单侧区间估计OneSamplet-testdata:xt=-2.4534,df=9,p-value=0.01828alternativehypothesis:truemeanislessthan7295percentconfidenceinterval:-Inf70.83705sampleestimates:meanofx67.4p-value=0.018280.05,拒绝原假设，平均脉搏低于常人。

教师签名年月日

--话说直接copy是不道德的哈

开课学院、试验室：试验时间：2023年月日

--话说直接copy是不道德的哈

22方差122

未知，S12和S2分别是X和Y的样本方差。由统计知识可知，当H0为真时，XYT~t(n1n22)1(n11)S2(n1)S21122Sn其中Swwn12nn212因此，当T满足(成为拒绝域):双边检验：Tt(n1n22)2

单边检验I单边检验II

Tt(n1n22)Tt(n1n22)

则认为H0不成立，此方法也称为t检验法。四、试验环境(所用软件、硬件等)及试验数据文件见试验内容软件：R2.15.3R软件。五、试验结果及实例分析a-c(126,125,136,128,123,138,142,116,110,108,115,140)b-c(162,172,177,170,175,152,157,159,160,162)###正态性检验：ks.test(a,pnorm,mean(a),sd(a))

One-sampleKolmogorov-Smirnovtest

data:aD=0.1464,p-value=0.9266alternativehypothesis:two-sided

ks.test(b,pnorm,mean(b),sd(b))

One-sampleKolmogorov-Smirnovtest

data:bD=0.2222,p-value=0.707alternativehypothesis:two-sided####方差齐性检验：var.test(a,b)

--话说直接copy是不道德的哈

Ftesttocomparetwovariances

data:aandbF=1.9646,numdf=11,denomdf=9,p-value=0.32alternativehypothesis:trueratioofvariancesisnotequalto195percentconfidenceinterval:0.50219437.0488630sampleestimates:ratioofvariances1.964622####可认为a和b的方差一致。####选用方差一致模型t检验：t.test(a,b,var.equal=TRUE)

TwoSamplet-test

data:aandbt=-8.8148,df=20,p-value=2.524e-08alternativehypothesis:tr