The problem I've encountered after trying to train neural networks isn't a new one : The fitted values I'm getting are all the same. Here's some oversimplified code as an example:
a <- c( 123, 223, 234, 226, 60)
b <- c(60, 90, 53, 54, 91)
d <- c(40,100,207,290,241)
q <- cbind(a,b,d)
nn <- neuralnet(a~b+d,data=q,hidden=2,threshold=0.01,err.fc="sse")
nn$net.result`
Previous answers I have stumbled upon suggest using nnet instead. I am getting the same results though, unless I set the decay argument to a value not equal to 0. Instead of blindly using the decay option, just because it seems to "work" though, I would appreciate understanding what goes wrong with my neuralnet model to begin with.