华中科技大学 陈菲 企业统计 Project.doc

企业统计项目 小组成员： 指导老师：陈非 个人报告 姓名： 学号： 指导老师： 陈非 1. 打开hitters.post1975.txt，选择数据：一个球员-赛季要求至少有100次击球和一次本垒打，即AB=100HR=1，然后对平均本垒打（HR/AB）画图，并从中选出平均本垒打（HR/AB）最高和最低的球员-赛季。 程序如下： Data-read.table(hitters.post1975.txt,header=TRUE,sep=,) data-Data[Data$AB=100Data$HR=1,] dim(Data) dim(data) 输出： AvgHR=data$HR/data$AB AvgHR hist(AvgHR,main=Histogram of AvgHR,xlab=AvgHR,ylab=Freq) Max.AvgHR=max(AvgHR) Min.AvgHR=min(AvgHR) Max.AvgHR Min.AvgHR data[which(AvgHR==Max.AvgHR),c(2,3),] data[which(AvgHR==Min.AvgHR),c(2,3),] 输出： 结果，平均本垒打最高的是bondsba01 2001赛季。平均本垒打最低的是millafe01 1975赛季。 2. 给出参数为a、b的似然函数和对数似然函数，然后用contour（）函数给出似然函数的等值线图。 程序如下： ## graphing likelihood over range of alpha and beta x-seq(0,0.15,0.000015) density.1 - dbeta(x,2,3) density.2 - dbeta(x,4,30) density.3 - dbeta(x,7,50) density.4 - dbeta(x,3,60) mindensity - min(density.1,density.2,density.3,density.4) maxdensity - max(density.1,density.2,density.3,density.4) hist(AvgHR,freq=F) lines(x,density.1,type=l,col=black,lwd=2, ,ylim=c(mindensity,maxdensity)) lines(x,density.2,type=l,col=green,lwd=2) lines(x,density.3,type=l,col=blue,lwd=2) lines(x,density.4,type=l,col=red,lwd=2) legend(0.4,25,c(Beta(1,1),Beta(1/2,1/2),Beta(1/3,1/3),Beta(0,0)),col=c(black,blue,red,green),lwd=2) #loglikBeta function n - length(AvgHR) numgrid - 100 loglikBeta - function(a,b){ loglik - -Inf if (a 0 b 0){ loglik - 0 for (i in 1:n){ loglik - loglik + a*log(AvgHR[i])+b*log(1-AvgHR[i]) loglik - loglik -log(beta(a, b)) } } loglik } alpharange - ppoints(numgrid)*0.2+1.9 # a between 1.9 and 2.1 betarange - ppoints(numgrid)*7+67 # b between 67 and 75 full - matrix(NA,nrow=numgrid,ncol=numgrid) for (i in 1:numgrid){ for (j in 1:numgrid){ full[i,j] - loglikBeta(alpharange[i],betarange[j]) } } full - exp(full - max(full)) full - full/sum(full) contour(alpharange,betarange,full,xlab=a,ylab=b,drawlabels=F,main=Contour Plot) 运行的结果