Constrained optimization with R

Question

I would like to solve a minimization problem with R \'Rsolnp\' package. My objective function is:

sum((I[i, j]*a[j] - E[i, j]*b[i])^2/(I[i, j] + E[i, j]))