Hi all,
I have a system of 3 equations with many defined parameters and 3 variables
I need to find solutions to. I actually know the solutions I'm aiming for
(.0061,0.287,0.070) but R tends to give me (0,0,0).
I tried the "BB" package but don't really follow how to refine my solutions
from that; these are wrong so far. Here's my code from trying that:
> f<-function(x){
+ x1<-x[1]
+ x2<-x[2]
+ x3<-x[3]
+
+ F<-rep(NA,3)
+ F[1]<-(0.075*x3+0,175*x2)*0.065*exp(-0.03*25)
+ F[2]<-(x3/0.075)*((1/50)/(0.065*5.125*0.872214*(1-(x3*(56/25))))-0.075)
+ F[3]<-(x2/0.175)*((1/70)/(0.065*0.366*0.872214*(1-(x2*(66/35))))-0.175)
+ return(F)
+ }
> p0<-c(0.1,0.1,0.1)
> dfsane(par=p0,fn=f,control=list(maxit=3000))
I have substituted estimated values for all the parameters into the above
equations.
Thanks
--
View this message in context:
http://www.nabble.com/Non-negative-solutions-to-complicated-equations-tp24622867p24622867.html
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
