附錄嶺回歸模型代碼library(MASS)library(ridge)library(openxlsx)setwd("C:/Users/華為/Desktop/R")data4.1<-read.xlsx("data4.1.xlsx")ridge4.1<-lm.ridge(Y~.,data4.1,lambda=seq(0,2,0.1))beat<-coef(ridge4.1)beatk<-ridge4.1$lambdaplot(k,k,type="n",xlab="嶺參數k",ylab="嶺回歸系數",ylim=c(-0.6,0.6))linetype<-c(1:9)char<-c(15:22)for(iin1:9)lines(k,beat[,i],type="o",lty=linetype[i],pch=char[i],cex=0.75)ridge4.1.2<-lm.ridge(Y~.-X3,data4.1,lambda=seq(0,1.5,0.1))beat1<-coef(ridge4.1.2)plot(k,k,type="n",xlab="嶺參數k",ylab="嶺回歸系數",ylim=c(-1,1))linetype<-c(1:8)char<-c(15:21)for(iin1:8)lines(k,beat[,i],type="o",lty=linetype[i],pch=char[i],cex=0.75)ridge4.1.3<-lm.ridge(Y~.-X3-X7,data4.1,lambda=seq(0,1,0.1))beat2<-coef(ridge4.1.3)plot(k,k,type="n",xlab="嶺參數k",ylab="嶺回歸系數",ylim=c(-0.6,0.6))linetype<-c(1:7)char<-c(15:20)for(iin1:7)lines(k,beat[,i],type="o",lty=linetype[i],pch=char[i],cex=0.75)ridge4.1.4<-lm.ridge(Y~.-X3-X7-X8,data4.1,lambda=seq(0,0.8,0.1))beat3<-coef(ridge4.1.4)plot(k,k,type="n",xlab="嶺參數k",ylab="嶺回歸系數",ylim=c(-0.6,0.6))linetype<-c(1:6)char<-c(15:19)for(iin1:6)lines(k,beat[,i],type="o",lty=linetype[i],pch=char[i],cex=0.75)ridge4.1.5<-lm.ridge(Y~.-X3-X7-X8-X6,data4.1,lambda=seq(0,0.8,0.1))beat4<-coef(ridge4.1.5)plot(k,k,type="n",xlab="嶺參數k",ylab="嶺回歸系數",ylim=c(-0.6,0.6))linetype<-c(1:5)char<-c(15:18)for(iin1:5)lines(k,beat[,i],type="o",lty=linetype[i],pch=char[i],cex=0.75)summary(lm(Y~.-X3-X7-X8-X6,data4.1))GM(1,1)模型運行代碼symsab;c=[ab]';X=[8.93,9.02,9.56,9.75,10.75,11.28,11.68];n=length(X);lambda=X(1:n-1)/X(2:n)range=minmax(lambda)range1=[exp(-2/(n+1)),exp(2/n+1)]flag=1;if(range(1)>range1(1)&&range(2)<range1(2))disp('級比檢驗通過');flag=1;elseflag=0;disp('級比檢驗未通過')endB=cumsum(X)fori=2:nC(i)=(B(i)+B(i-1))/2;endC(1)=[];B=[-C;ones(1,n-1)];Y=X;Y(1)=[];Y=Y'c=inv(B*B')*B*Y;c=c';a=c(1)b=c(2)F=[];F(1)=X(1);fori=2:(n+4)F(i)=(X(1)-b