R语言与统计分析第四章答案.docxVIP

下载本文档

59
0
约3.56千字
约 4页
2017-07-21 发布于河南
举报
版权申诉

R语言与统计分析第四章答案.docx

1、本文档共4页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
8、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

R语言与统计分析第四章答案

#第四章习题#4.1x-rbinom(1000,100,0.3)hist(x,main=c(1000个参数为0.3的伯努利分布随机数))#4.2x-rnorm(1000,10,4)hist(x,probability=T,xlim=c(min(x),max(x)),nclass=max(x)-min(x)+1, col=lightblue,main=c(1000个正态分布随机数))lines(density(x,bw=1),col=blue,lwd=3)#4.3x-sample(c(rt(10,1),rt(10,2),rt(10,10)),1000,replace=T)hist(x,xlim=c(min(x),max(x)),probability=T,nclass=max(x)-min(x)+1,col=lightblue, main=c(3个t分布混合样本直方图))lines(density(x,bw=1),col=blue,lwd=2)#方法二k-matrix(,3,100)k[1,]=rt(100,1)k[2,]=rt(100,2)k[3,]=rt(100,10)x=c(k[1,],k[2,],k[3,]) #3个t分布混合成一个样本hist(x,xlim=c(min(x),max(x)),probability=T,nclass=max(x)-min(x)+1,col=lightblue, main=c(3个t分布混合样本直方图))lines(density(x,bw=1),col=blue,lwd=2)#4.4install.packages(DAAG)library(DAAG)data(possum)par(mfrow=c(2,2))hist(possum$age,breaks=1+(0:8)*1)hist(possum$age,breaks=0+(0:9)*1)hist(possum$age,breaks=1+(0:5)*2)hist(possum$age,breaks=0+(0:5)*2)summary(possum$age)age-possum$age[!is.na(possum$age)]summary(age)sd(age)#4.5install.packages(DAAG)library(DAAG)data(tinting)ts-table(tinting$sex,tinting$tint) #列联表barplot(ts) #联合柱状图windows() #新图op-par()layout(matrix(c(2,1,0,3),2,2,byrow=T),c(1,6),c(4,1))par(mar=c(1,1,5,1))plot(tinting$age,tinting$it)lines(lowess(tinting$age,tinting$it),lwd=2) #拟合线rug(side=2,jitter(tinting$age,5)) #细小刻度rug(side=1,jitter(tinting$it,5))par(mar=c(1,2,5,1))boxplot(tinting$age,axes=F)par(mar=c(5,1,1,2))boxplot(tinting$it,horizontal=T,axes=F)windows() #因子为tintcoplot(tinting$age~tinting$it|tinting$tint)windows() #因子为tint与sexcoplot(tinting$age~tinting$it|tinting$tint*tinting$sex)windows() #等高线图library(MASS)z-kde2d(tinting$it,tinting$csoa)contour(z,col=red,drawlabels=FALSE)windows() #matplot图d-data.frame(y1=tinting$age,y2=tinting$it,y3=tinting$csoa) matplot(d,type=l,main=matplot)#4.6data(InsectSpray