統(tǒng)計(jì)建模與r軟件課后答案-展示頁(yè)
2025-06-29 02:20本頁(yè)面
【正文】 出的結(jié)果看是相關(guān)的。) windows() persp(x,y,z,theta=30,phi=30,expand=,col=39。yseq(1,7,) ffunction(x,y)+ x^42*x^2*y+x^22*x*y+2*y^2+9*x/24*y+4 zouter(x,y,f)contour(x,y,z,levels=c(0,1,2,3,4,5,10,15,20,30,40,50,60,80,100),col=39。39。從上圖可知所有的圖中偶讀沒(méi)有這樣的趨勢(shì)，故均不相關(guān)。x339。x239。x139。x2c(5,6,8,5,10,7,12,12,6,6)。pnorm39。blue39。red39。利用data_outline函數(shù)。運(yùn)算成功39。要求輸入一個(gè)正整數(shù)39。39。39。39。5239。39。4939。4239。15939。16239。15739。16539。15639。1539。1439。1639。1539。1439。女39。男39。女39。男39。女39。丁一39。趙六39。王五39。李四39。張三39。E(4) FA[1:3,1:3](5) GB[,3] xc(rep(1,5),rep(2,3),rep(3,4),rep(4,2))。C(2) DA%*%B。z2 [,1] [,2] [,3][1,] 4 5 6[2,] 8 10 12[3,] 12 15 18(1) Amatrix(1:20,nrow=4)。z[1] 7 10 13 z1crossprod(x,y)。第二章 xc(1,2,3)。yc(4,5,6) ec(1,1,1) z2*x+y+e。z1 [,1][1,] 32 z2outer(x,y)。Bmatrix(1:20,nrow=4,byrow=T) CA+B。D(3) EA*B。x Hmatrix(nrow=5,ncol=5) for (i in 1:5)+ for(j in 1:5)+ H[i,j]1/(i+j1)（1） det(H)（2） solve(H)（3） eigen(H) studentdata(姓名=c(39。,39。,39。,39。,39。)+ ,性別=c(39。,39。,39。,39。,39。),年齡=c(39。,39。,39。,39。,39。),+ 身高=c(39。,39。,39。,39。,39。),體重=c(39。,39。,39。,39。,39。)) (studentdata,file=39。) (studentdata,file=39。)countfunction(n){if (n=0)print(39。)else{repeat{if (n%%2==0)nn/2elsen(3*n+1)if(n==1)break}print(39。)}}第三章首先將數(shù)據(jù)錄入為x。如下 data_outline(x) hist(x,freq=F) lines(density(x),col=39。) ymin(x):max(x) lines(y,dnorm(y,),col=39。) plot(ecdf(x),verticals=T,=F) lines(y,pnorm(y,)) qqnorm(x) qqline(x) stem(x) boxplot(x) fivenum(x) (x) (x,39。,) Onesample KolmogorovSmirnov testdata: xD = , pvalue = alternative hypothesis: twosidedWarning message:In (x, pnorm, , ) : ties should not be present for the KolmogorovSmirnov test這里出現(xiàn)警告信息是因?yàn)閗s檢驗(yàn)要求樣本數(shù)據(jù)是連續(xù)的，不允許出現(xiàn)重復(fù)值x1c(2,4,3,2,4,7,7,2,2,5,4)。x3c(7,11,6,6,7,9,5,5,10,6,3,10) boxplot(x1,x2,x3,names=c(39。,39。,39。),vcol=c(2,3,4))windows()plot(factor(c(rep(1,length(x1)),rep(2,length(x2)),rep(3,length(x3)))),c(x1,x2,x3)) rubber(x1=c(65,70,70,69,66,67,68,72,66,68),+x2=c(45,45,48,46,50,46,47,43,47,48),x3=c(,)) plot(rubber)具體有相關(guān)關(guān)系的兩個(gè)變量的散點(diǎn)圖要么是從左下角到右上角（正相關(guān)），要么是從左上角到右下角（負(fù)相關(guān)）。（1） student(39。) attach(student) plot(體重~身高)（2） coplot(體重~身高|性別)（3） coplot(體重~身高|年齡)（4） coplot(體重~身高|年齡+性別)只列出（4）的結(jié)果，如下圖 xseq(2,3,)。blue39。red39。具體結(jié)果不再列出 df(39。) stars(df)然后按照G的標(biāo)準(zhǔn)來(lái)畫出星圖 attach(df) df$G1(SC+LC+SMS+DRV+AMB+GSP+POT)/7 df$G2(FL+EXP+SUIT)/3 df$G3(LA+HON+KJ)/3 df$G4AA df$G5APP ascale(df[,17:21]) stars(a)這里從17開(kāi)始取，是因?yàn)樵赿f中將ID也作為了一列使用P159已經(jīng)編好的函數(shù)unison，接著上題，直接有 unison(a)第四章（1）先求矩估計(jì)。因此我們有。在R中實(shí)現(xiàn)如下 xc(,) (2*mean(x)1)/(1mean(x))[1] （2）采用極大似然估計(jì)首先求出極大似然函數(shù)為L(zhǎng)a。x=nlna+1+aln(i=1nxi最后求導(dǎo)?lnL(a。我們知道泊松分布中的參數(shù)λ，既是均值又是方差。故在R中實(shí)現(xiàn)如下 xc(54,67,68,78,70,66,67,70,65,69) mean(x)[1] ，如下 (x) (x,alternative=39。) (x,alternative=39。)此時(shí)我們只需要區(qū)間估計(jì)的結(jié)果。下面我們來(lái)做是否低于72的均值假設(shè)檢驗(yàn)。greater39。因此這10名患者的平均脈搏次數(shù)比正常人要小。但得出的結(jié)論是不影響的（他們的期望差別很大） Ac(,) Bc(,) (A,B) xc(140,137,136,140,145,148,140,135,144,141) yc(135,118,115,140,128,131,130,115,131,125) (x,y) (x,y,=F)泊松分布的參數(shù)就等于它的均值也等于方差。 xc(rep(0,7),rep(1,10),rep(2,12),rep(3,8),rep(4,3),rep(5,2)) mean(x)[1] (x)正態(tài)總體均值用樣本均值來(lái)估計(jì)。greater39。首先我們要把該正態(tài)分布的均值和方差給估計(jì)出來(lái)，這個(gè)就利用樣本即可。如下 xc(1067,919,1196,785,1126,936,918,1156,920,948) pnorm(1000,mean(x),sd(x))[1] [1] 這是檢驗(yàn)兩個(gè)總體是否存在差異的問(wèn)題。兩種方法實(shí)現(xiàn)如下 xc(113,120,138,120,100,118,138,123) yc(138,116,125,136,110,132,130,110) (sum(xy),length(x))pvalue = 1 (x,y,exact=F)pvalue = ，故接受原假設(shè)，他們無(wú)差異（1）采用w檢驗(yàn)法xc(,2,4,3,6,)yc(,5,6,2,2,2) (x) (y)采用ks檢驗(yàn)法 (x,39。,mean(x
公司管理相關(guān)推薦
鄂ICP備17016276號(hào)-1