いやね、…


このCpプロットを簡単書くためにPerlを書いてただけですよ。だって、そうじゃないとモデルをもれなく、重複なく、選ばないといけないという作業やらないといけないことになってしまうんだもん。こんなものコピペでやっても書きたくないでしょ?

rss.21<-sum(lsfit(cbind(mhw,mid.displacement),log.mid.price)$residuals^2)
Cp.21<-(rss.21/unbiased.sigma.square.hat)-(length(log.mid.price)-2*3)
rss.22<-sum(lsfit(cbind(mhw,mid.mpg),log.mid.price)$residuals^2)
Cp.22<-(rss.22/unbiased.sigma.square.hat)-(length(log.mid.price)-2*3)
rss.23<-sum(lsfit(cbind(mhw,is.eu),log.mid.price)$residuals^2)
Cp.23<-(rss.23/unbiased.sigma.square.hat)-(length(log.mid.price)-2*3)
rss.24<-sum(lsfit(cbind(mhw,is.hybrid),log.mid.price)$residuals^2)
Cp.24<-(rss.24/unbiased.sigma.square.hat)-(length(log.mid.price)-2*3)
rss.25<-sum(lsfit(cbind(mhw,is.suv),log.mid.price)$residuals^2)
Cp.25<-(rss.25/unbiased.sigma.square.hat)-(length(log.mid.price)-2*3)
rss.26<-sum(lsfit(cbind(mid.displacement,mid.mpg),log.mid.price)$residuals^2)
Cp.26<-(rss.26/unbiased.sigma.square.hat)-(length(log.mid.price)-2*3)
rss.27<-sum(lsfit(cbind(mid.displacement,is.eu),log.mid.price)$residuals^2)
Cp.27<-(rss.27/unbiased.sigma.square.hat)-(length(log.mid.price)-2*3)
rss.28<-sum(lsfit(cbind(mid.displacement,is.hybrid),log.mid.price)$residuals^2)
Cp.28<-(rss.28/unbiased.sigma.square.hat)-(length(log.mid.price)-2*3)
rss.29<-sum(lsfit(cbind(mid.displacement,is.suv),log.mid.price)$residuals^2)
Cp.29<-(rss.29/unbiased.sigma.square.hat)-(length(log.mid.price)-2*3)
rss.210<-sum(lsfit(cbind(mid.mpg,is.eu),log.mid.price)$residuals^2)
Cp.210<-(rss.210/unbiased.sigma.square.hat)-(length(log.mid.price)-2*3)
rss.211<-sum(lsfit(cbind(mid.mpg,is.hybrid),log.mid.price)$residuals^2)
Cp.211<-(rss.211/unbiased.sigma.square.hat)-(length(log.mid.price)-2*3)
rss.212<-sum(lsfit(cbind(mid.mpg,is.suv),log.mid.price)$residuals^2)
Cp.212<-(rss.212/unbiased.sigma.square.hat)-(length(log.mid.price)-2*3)
rss.213<-sum(lsfit(cbind(is.eu,is.hybrid),log.mid.price)$residuals^2)
Cp.213<-(rss.213/unbiased.sigma.square.hat)-(length(log.mid.price)-2*3)
rss.214<-sum(lsfit(cbind(is.eu,is.suv),log.mid.price)$residuals^2)
Cp.214<-(rss.214/unbiased.sigma.square.hat)-(length(log.mid.price)-2*3)
rss.215<-sum(lsfit(cbind(is.hybrid,is.suv),log.mid.price)$residuals^2)
Cp.215<-(rss.215/unbiased.sigma.square.hat)-(length(log.mid.price)-2*3)

rss.31<-sum(lsfit(cbind(mhw,mid.displacement,mid.mpg),log.mid.price)$residuals^2)
Cp.31<-(rss.31/unbiased.sigma.square.hat)-(length(log.mid.price)-2*4)
rss.32<-sum(lsfit(cbind(mhw,mid.displacement,is.eu),log.mid.price)$residuals^2)
Cp.32<-(rss.32/unbiased.sigma.square.hat)-(length(log.mid.price)-2*4)
rss.33<-sum(lsfit(cbind(mhw,mid.displacement,is.hybrid),log.mid.price)$residuals^2)
Cp.33<-(rss.33/unbiased.sigma.square.hat)-(length(log.mid.price)-2*4)
rss.34<-sum(lsfit(cbind(mhw,mid.displacement,is.suv),log.mid.price)$residuals^2)
Cp.34<-(rss.34/unbiased.sigma.square.hat)-(length(log.mid.price)-2*4)
rss.35<-sum(lsfit(cbind(mhw,mid.mpg,is.eu),log.mid.price)$residuals^2)
Cp.35<-(rss.35/unbiased.sigma.square.hat)-(length(log.mid.price)-2*4)
rss.36<-sum(lsfit(cbind(mhw,mid.mpg,is.hybrid),log.mid.price)$residuals^2)
Cp.36<-(rss.36/unbiased.sigma.square.hat)-(length(log.mid.price)-2*4)
rss.37<-sum(lsfit(cbind(mhw,mid.mpg,is.suv),log.mid.price)$residuals^2)
Cp.37<-(rss.37/unbiased.sigma.square.hat)-(length(log.mid.price)-2*4)
rss.38<-sum(lsfit(cbind(mhw,is.eu,is.hybrid),log.mid.price)$residuals^2)
Cp.38<-(rss.38/unbiased.sigma.square.hat)-(length(log.mid.price)-2*4)
rss.39<-sum(lsfit(cbind(mhw,is.eu,is.suv),log.mid.price)$residuals^2)
Cp.39<-(rss.39/unbiased.sigma.square.hat)-(length(log.mid.price)-2*4)
rss.310<-sum(lsfit(cbind(mhw,is.hybrid,is.suv),log.mid.price)$residuals^2)
Cp.310<-(rss.310/unbiased.sigma.square.hat)-(length(log.mid.price)-2*4)
rss.311<-sum(lsfit(cbind(mid.displacement,mid.mpg,is.eu),log.mid.price)$residuals^2)
Cp.311<-(rss.311/unbiased.sigma.square.hat)-(length(log.mid.price)-2*4)
rss.312<-sum(lsfit(cbind(mid.displacement,mid.mpg,is.hybrid),log.mid.price)$residuals^2)
Cp.312<-(rss.312/unbiased.sigma.square.hat)-(length(log.mid.price)-2*4)
rss.313<-sum(lsfit(cbind(mid.displacement,mid.mpg,is.suv),log.mid.price)$residuals^2)
Cp.313<-(rss.313/unbiased.sigma.square.hat)-(length(log.mid.price)-2*4)
rss.314<-sum(lsfit(cbind(mid.displacement,is.eu,is.hybrid),log.mid.price)$residuals^2)
Cp.314<-(rss.314/unbiased.sigma.square.hat)-(length(log.mid.price)-2*4)
rss.315<-sum(lsfit(cbind(mid.displacement,is.eu,is.suv),log.mid.price)$residuals^2)
Cp.315<-(rss.315/unbiased.sigma.square.hat)-(length(log.mid.price)-2*4)
rss.316<-sum(lsfit(cbind(mid.displacement,is.hybrid,is.suv),log.mid.price)$residuals^2)
Cp.316<-(rss.316/unbiased.sigma.square.hat)-(length(log.mid.price)-2*4)
rss.317<-sum(lsfit(cbind(mid.mpg,is.eu,is.hybrid),log.mid.price)$residuals^2)
Cp.317<-(rss.317/unbruals^2)
Cp.61<-(rss.61/unbiased.sigma.square.hat)-(length(log.mid.price)-2*7)

もうコピペでやるか、プログラムにやってもらうか説明するまでもないですね。僕はこんなの手でやってたら蕁麻疹出てきますw。