読者です 読者をやめる 読者になる 読者になる

楽チンに課題を終わらせちゃう呪文とか

R

いや、別にこれくらい普通にやって問題ないわけだけど、この前の練習とかやってみたかったわけですよ。
jpgファイルをeqsファイルに変換させる。

perl -e '@list=();foreach(1..5){push(@list,qq/convert $_.jpg eps2:$_.eps/);}map{system $_}@list;'

7個くらいある説明変数と従属変数の散布図を書くための関数とか残差との散布図とか。

draw_cor<-function(n){
   png(paste(n,".jpg",sep=""))
   eval(parse(text=paste("plot(data$x",n,",data$y,xlab=\"X",n,"\",ylab=\"Y\")",sep="")))
   eval(parse(text=paste("abline(lm(y~x",n,",data=data))",sep="")))
   title(paste("X",n,"とYの散布図と回帰直線"))
   dev.off()
}

for(i in 1:5){
    draw_cor(i)
}

draw_res<-function(n){
   png(paste("res",n,".jpg",sep=""))
   eval(parse(text=paste("plot(data$x",n,",result$residuals,xlab=\"X",n,"\",ylab=\"Y\")",sep="")))
   title(paste("X",n,"と残差の散布図"))
   dev.off()
}

for(i in 1:7){
    draw_res(i)
}