starting with the resulting expression g I need g's derivative as expression, but I get: "Function `eval' is not in the derivatives table":
#function (double gamama) h <- expression((t/d1)^a1*exp(-(t-d1)/b1)-c*(t/d2)^a2*exp(-(t-d2)/b2)) # plot it t <- seq(0, 20000, by = 100) t <- t/1000 g <- expression(eval(h)/sqrt(sum(eval(h)^2))) plot(eval(g), type = 'l') g.deriv <- D(g, "t")
> Error in D(g, "t") : Function `eval' is not in the derivatives table
is there any way one can solve this problem?
thanks a lot christoph
## -- ## parameters peak1 <- 5.4 fwhm1 <- 5.2 peak2 <- 10.8 fwhm2 <- 7.35 dip <- 0.35
b1 <- 0.9 # dispersion b2 <- 0.9 #dispersion a1 <- peak1/b1 a2 <- peak2/b2 d1 <- a1*b1 d2 <- a2*b2 c <- dip
______________________________________________ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html