BP神经网络
>
>
> data<-read.csv('concrete.csv')
>
>
>
> normalize<-function(x){
+ return((x-min(x))/(max(x)-min(x)))
+ }
> data_normalize<-as.data.frame(lapply(data,normalize))
>
> dat<-sample(2,nrow(data_normalize),replace = T,prob = c(0.75,0.25))
> train_data<-data_normalize[dat==1,]
> test_data<-data_normalize[dat==2,]
>
>
> library(neuralnet)
> model<-neuralnet(strength~cement+slag+ash+water+superplastic+coarseagg+fineagg+age,data=train_data,hidden = 1)
>
>
> pre_strength<-compute(model,test_data[1:8])
>
> cor(pre_strength$net.result,test_data$strength)
[,1]
[1,] 0.8574563
> mean(abs(test_data$strength-pre_strength$net.result))
[1] 0.08472698
>
>
> model2<-neuralnet(strength~cement+slag+ash+water+superplastic+coarseagg+fineagg+age,data=train_data,hidden = 5)
> plot(model2)
> pre_strength2<-compute(model2,test_data[1:8])
> cor(pre_strength2$net.result,test_data$strength)
[,1]
[1,] 0.9467139
> mean(abs(test_data$strength-pre_strength2$net.result))
[1] 0.05568677