Displaying 1 result from an estimated 1 matches for "xppaut".
Did you mean:
xppat
2005 Jun 30
0
vecortizing uniroot() for numerical solutions
...oot
}
df1
# what I'd really like to do is something like:
R12 = function(d1, d2, h1, h2) (d1/(h1^2))-(h2*(d2^2))
uniroot(R12, interval=c(-100,100), tol=0.0001, d2=c(1:5), h1=c(1:5),
h2=c(1:5))
# but R rejects the multiple values for d2.
# I could use a numerical solver such as XPPAUT, but would rather be
able to do it all in R.
# Any thoughts on vectorizing this would be greatly appreciated.
# Thanks! -Dan
--
Daniel E. Bunker
Associate Coordinator - BioMERGE
Post-Doctoral Research Scientist
Columbia University
Department of Ecology, Evolution and Environmental Bio...