Displaying 1 result from an estimated 1 matches for "pp_s2".
Did you mean:
pp_o2
2008 Aug 14
3
extending the derivs table/fools rushing in
...day's
R-devel src/deriv/main.c is available at
http://www.zoology.ufl.edu/bolker/deriv_patch.txt .
I pretty much followed the framework of the other symbols;
here was my incantation
- } else if (CAR(expr) == PlogisSymbol) {
- ans = simplify(TimesSymbol,
- PP_S(TimesSymbol,
- PP_S2(ExpSymbol,
- PP_S2(MinusSymbol,CADR(expr))),
- PP_S(PowerSymbol,
- PP_S(PlusSymbol,
- Constant(1.),
- PP_S2(ExpSymbol,
- PP_S2(MinusSymbol,CADR(expr)))),
- Constant(-2.))),
- PP(D(CADR(expr),var)));
- UNPROTECT(8);
It seems to work:...