I'm tinkering around in Ryacas trying to find an easy way to get the first and second partial derivatives of mu and sigma from a normal distribution (actually a bitterly ugly likelihood but this example works for now).
I've done all of the work in Mathematica, but I then need to manually write R code whereas I think Ryacas will give me the expression such that I can export for Tex and produce code I can use for a function. Is there a way I can use Ryacas with built in functions, like dnorm? For example, I can define the normal and then use the D function as: ndist<-expression( (1/(s*sqrt(2*pi))) * exp(- ((x-u)^2/(2*s^2)))) D(ndist, "u") But, I don't see in the vignette a way to use dnorm? ______________________________________________ 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.
