統(tǒng)計建模與r軟件課后答案(文件)
【正文】 危險因子X245有無計劃X367是否感染Y10對數(shù)據(jù)的處理，如下X1X2X3Y頻數(shù)246112460172561025602247111247087257102570034612834603034712334703356183560323571035709然后用R處理并求解模型hospital(X1=rep(c(2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3),c(1,17,0,2,11,87,+ 0,0,28,30,23,3,8,32,0,9)),X2=rep(c(4,4,5,5,4,4,5,5,4,4,4,4,5,5,5,5),+ c(1,17,0,2,11,87,+ 0,0,28,30,23,3,8,32,0,9)),X3=rep(c(6,6,6,6,7,7,7,7,6,6,7,7,6,6,7,7),+ c(1,17,0,2,11,87,0,0,28,30,23,3,8,32,0,9)),+ Y=rep(c(1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0),c(1,17,0,2,11,87,0,0,28,30,23,3,8,32,0,9))+ ) glm(Y~X1+X2+X3,family=binomial,data=hospital) summary()，則方程的系數(shù)和方程本省全部通過檢驗。39。prepredict(,(X1=pe$X1,X2=pe$X2,X3=pe$X3,X4=pe$X4,X5=pe$X5)) pexp(pre)/(1+exp(pre)) p（2） step()結(jié)果是只保留了變量X1和X4。prepredict(,(X1=pe$X1,X2=pe$X2,X3=pe$X3,X4=pe$X4,X5=pe$X5)) exp(pre)/(1+exp(pre)) 顯然經(jīng)過逐步回歸后的模型更合理。summary()Coefficients: Estimate Std. Error t value Pr(|t|) (Intercept) ***X *** summary() exp($coefficients[1,1])[1] 所以theta0=,theta1=（2） nls(Y~b0*exp(b1*X),data=peo,start=list(b0=50,b1=0)) summary()Parameters: Estimate Std. Error t value Pr(|t|) b0 ***b1 ***發(fā)現(xiàn)所求的基本上與內(nèi)在線性相同。甲39。丙39。（3） (pro$Y,pro$X) 1 2 2 3 可以看到顯著性主要有乙工廠和丙工廠造成（1） old(Y=c(20,18,19,17,15,16,13,18,22,17,26,19,26,28,+ 23,25,24,25,18,22,27,24,12,14),X=factor(rep(1:4,c(10,6,6,2)))) aov(Y~X,data=old) summary()可以發(fā)現(xiàn)影響是非常顯著的。summary()結(jié)果是A和B及其交互作用都是十分顯著的（2）首先我們要選出最優(yōu)條件組合，由（1）知影響力為ABAB。zrep(0,n)+ for( i in 1:n)+ if(x[i]==y[i]){z[i]1}else{z[i]2}+ factor(z)} pro$ABab(pro$A,pro$B)然后我們開始計算各個水平的均值，如下 Kmatrix(0,nrow=3,ncol=3,dimnames=list(1:3,c(39。,39。下面給出A3B1的點估計和區(qū)間估計。summary()可以看到影響均不顯著，那么我們干脆直接按照各因素水平的均值大小來取。密度39。此題暫且不做首先把正交試驗表的結(jié)果那一列給計算出來。A39。C39。AC39。39。plete39。centroid39。39。average39。ward39。 summary()，這牽涉到經(jīng)濟(jì)學(xué)知識，我不懂。冶金39。煤炭39。機械39。森工39。紡織39。皮革39。文教用品39。前兩個主成分的累積貢獻(xiàn)率已經(jīng)達(dá)到99%。身高39。上肢長39。體重39。胸圍39。因此第一個因子可以叫做體長因子，第二個可以叫做身寬因子。) plot(fa$scores[,1:2],type=39。 kid(39。 U(kidc[,1:3])%*%$xcoef V(kidc[,4:6])%*%$ycoef oparpar(mfrow=c(3,1)) plot(U[,1],V[,1],xlab=39。) plot(U[,2],V[,2],xlab=39。) plot(U[,3],V[,3],xlab=39。)。,ylab=39。,ylab=39。,ylab=39。,header=T) kidcscale(kid) cancor(kidc[,1:3],kidc[,4:6]) source(E:\\Documents\\R\\自編函數(shù)\\) ($cor,40,p=3,q=3)[1] 2由此可見，我們應(yīng)該從前兩對典型相關(guān)變量里來對數(shù)據(jù)進(jìn)行分析。) text(fa$scores[,1],fa$scores[,2])從圖上直觀看出，第7位同學(xué)文科強而理科弱，第8,10位同學(xué)理科強而文科弱。（2） fafactanal(X,factors=2,scores=39。胸寬39。頸圍39。下肢長39。手臂長39。 prpredict() exp$z1pr[,1] exp$z2pr[,2] lm(y~z1+z2,data=exp) summary()最后的回歸方程不再列出。 exp(x1=c(,88,131,+ ),x2=c(92,93,96,94,100,101,105,112,112,112),x3=c(,29,+ 34,40,44,49,51,53),x4=c(94,96,97,97,100,101,104,109,111,111),y=c(,+ ,)) prinp(~x1+x2+x3+x4,data=exp,cor=T) summary()。造紙39?？p紉39。食品39。建材39。化學(xué)39。電力39。下面我們只保留前四個主成分，來對這13個行業(yè)進(jìn)行分類，即聚類分析。39。centroid39。plete39。ward39。average39。39。再依據(jù)AB,應(yīng)選擇AB1，也就是說A和B應(yīng)該是同一水平。AB39。B39。summary()從分析結(jié)果可以看出，顯著性大小為BABAC，其余均不顯著下面再計算出均值，從而就可以依據(jù)顯著性來選擇最優(yōu)參數(shù)了 abfunction(x,y){+ nlength(x)。施肥量39。品種39。pro$B==1]) (pro$Y[pro$A==3amp。))) for( i in 2:4)+ for(j in 1:3)+ K[j,i1]mean(pro$Y[pro[i]==j]) K A B AB1 2 3 NaN按照影響力越大（即P值越?。?，我們首先確定AB應(yīng)選擇水平2，即A和B 不等的是最好的。,39。首先要交互作用給找出來。而國內(nèi)只有甲乙，甲丙之間存在著顯著差異。 (pro$Y[pro$X==2])。乙39。mean39。然后將得到的方程用lm回歸。更加合理一些。只有X1的系數(shù)通過了顯著性檢驗，但是也不是很理想。 prepredict(,(X1=2,X2=4,X3=6)) pexp(pre)/(1+exp(pre))。我們把這些影響變量進(jìn)行編碼。說明step函數(shù)是合理的。我們可以檢驗一下step函數(shù)去掉變量后的共線性。因此可考慮剔除5,8,9,24點再做擬合。因此應(yīng)
公司管理相關(guān)推薦
鄂ICP備17016276號-1