keunhchoi at gmail.com writes:
Here is a clean version. I did this with nls and it works (see below), but
I need to do it with optim. Keun-Hyung
Great that you got it to work! Just curious if there is any reason why one
should use optim when nls works fine? If this should be on a
Here is a clean version. I did this with nls and it works (see below), but
I need to do it with optim. Keun-Hyung
# optim
vol-rep(c(0.03, 0.5, 2, 4, 8, 16, 32), 3)
time-rep(c(2,4,8),each=7)
p.mated-c(0.47, 0.48, 0.43, 0.43, 0.26, 0.23, NA, 0.68, 0.62, 0.64, 0.58,
0.53, 0.47,
0.24, 0.8, 0.79,
try:
vol - rep(c(0.03, 0.5, 2, 4, 8, 16, 32), 3)
time - rep(c(2,4,8),each=7)
p.mated - c(0.47, 0.48, 0.43, 0.43, 0.26, 0.23, NA, 0.68, 0.62, 0.64,
0.58, 0.53, 0.47, 0.24, 0.8, 0.79, 0.71, 0.56, 0.74, 0.8, 0.47)
eury - data.frame(vol=vol, time=time, p.mated=p.mated)
eury - na.omit(eury)
p0 -
3 matches
Mail list logo
