这是我测试过的一个例子，可以参考下
> library(nnet)
> norm.data <- scale(newSet)

> nn <- nnet(y~., norm.data,size = 10,decay = 0.01, maxit = 1000, linout = T, trace = F)
> norm.preds <- predict(nn, norm.data)
> plot(norm.preds~ scale(newSet$y))
计算相对误差
> (nmse2 <- mean((norm.preds-scale(newSet$y))^2)/
+      mean((mean( scale(newSet$y))- scale(newSet$y))^2))