similar to: Partial derivatives of a matrix with respect to a single element

Hi, Is there a function for raising a matrix to a power? For example if you like to compute A%*%A%*%A, is there an abbreviation similar to A^3? Atte Tenkanen > A=rbind(c(1,1),c(-1,-2)) > A [,1] [,2] [1,] 1 1 [2,] -1 -2 > A^3 [,1] [,2] [1,] 1 1 [2,] -1 -8 But: > A%*%A%*%A [,1] [,2] [1,] 1 2 [2,] -2 -5

Hello, I am trying to construct two matrices, F and V, composed of partial derivatives and then find the eigenvalues of F*Inverse(V). I have the following equations in ryacas notation: > library(Ryacas) > FIh <- Expr("betah*Sh*Iv") > FIv <- Expr("betav*Sv*Ih") > VIh <- Expr("(muh + gamma)*Ih") > VIv <- Expr("muv*Iv") I

Ryacas is an R interface to the free yacas computer algebra system. Ryacas allows one to send R expressions, unprocessed yacas strings and certain other R objects to a separate yacas process from R and get back the result. It also has facilities for manipulating yacas strings and R expressions destined for yacas processing. It can be used for exact arithmetic, symbolic math, ASCII pretty

Ryacas is an R interface to the free yacas computer algebra system. Ryacas allows one to send R expressions, unprocessed yacas strings and certain other R objects to a separate yacas process from R and get back the result. It also has facilities for manipulating yacas strings and R expressions destined for yacas processing. It can be used for exact arithmetic, symbolic math, ASCII pretty

Ryacas is an R interface to the free yacas computer algebra system. Ryacas allows one to send R expressions, unprocessed yacas strings and certain other R objects to a separate yacas process from R and get back the result. It also has facilities for manipulating yacas strings and R expressions destined for yacas processing. It can be used for exact arithmetic, symbolic math, ASCII pretty

Dear All I am doing the following: > x <- yacas("3/2") > for (i in 2:400) + x <- yacas(paste(x,"*",x)) > x expression(Inf^1.260864167e+117/Inf^6.304320836e+116) > Eval(x) [1] NaN No luck this way. However, I am successful with > y <- yacas("(3/2)^400") > y expression(7.05507910865533e+190/2.58224987808691e+120) > Eval(y) [1]

Hi all, I am trying to implement the following matlab code with Ryacas : syms U x x0 C d1=diff(U/(1+exp(-(x-x0)/C)),x); pretty(d1) d2=diff(U/(1+exp(-(x-x0)/C)),x,2); pretty(d2) solx2 = solve(d2 == 0, x, 'Real', true) pretty(solx2) slope2=subs(d1,solx2) I have tried the following : library(Ryacas) x <- Sym("x");U <- Sym("U");x0 <-

Thanks for the response. Yes, I did study the vignette but did not understand it fully. Anyway, I have tried once again now. I am happy to say that I have got what I wanted. library(Ryacas) x <- Sym("x");U <- Sym("U");x0 <- Sym("x0");C <- Sym("C") my_func <- function(x,U,x0,C) { return (U/(1+exp(-(x-x0)/C)))} FirstDeriv <-

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

Hi all, I wrote a symbollic differentiation function in R, which can be downloaded here: http://www.econ.upenn.edu/~clausen/computing/Deriv.R http://www.econ.upenn.edu/~clausen/computing/Simplify.R It is just a prototype. Of course, R already contains two differentiation functions: D and deriv. However, these functions have several limitations. They can probably be fixed, but

Have you studied the "Introduction to Ryacas" vignette that come with the package? Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Tue, Sep 19, 2017 at 2:37 AM, Vivek Sutradhara <viveksutra at gmail.com> wrote:

Hello yet again, R People: I was working with Ryacas and yacas last night and all was well. Now this morning, I keep getting the following: > a <- Sym("a") > a Error in summary.connection(x) : invalid connection > When I go to yacas from the command line, it works fine. Any suggestions, please? I'm thinking that a port might be open, but here I have: >

Hi I am currently (for pedagogical purposes) writing a simple numerical analysis library in R. I have come unstuck when writing a simple Newton-Raphson implementation, that looks like this: f <- function(x) { 2*cos(x)^2 + 3*sin(x) + 0.5 } root <- newton(f, tol=0.0001, N=20, a=1) My issue is calculating the symbolic derivative of f() inside the newton() function. I cant seem to get R to

I am having issues with Ryacas errors. I searched and found this error was reported in 2006, but nothing since. > library(Ryacas) Loading required package: XML > yacas(expression(Factor(x^2-1))) [1] "Starting Yacas!" CommandLine(1) : Expecting ) closing bracket for sub-expression, but got ^ instead > yacas("D(x)Sin(x)") CommandLine(1) : Expecting ) closing bracket

Hello again. Now I have the following: > xx [1] "Solve(1 - R/100==(2*y)/10000,R)" > yacas(xx) expression(list(R == 100 * (1 - 2 * y/10000))) > I would like to put in a value for y and obtain R. I've tried more stuff with eval and Eval, but no luck yet. Any suggestions would be much appreciated. Thanks, Erin -- Erin Hodgess Associate Professor Department of Computer

There is either something wrong with either me or is Yacas/Ryacas doing odd things. See below I ask to simpify an expression and the there's output! If this is a simplification.. well.. Do you think there is something set wrong somewhere. Thanks. > library(Ryacas) Loading required package: XML > a=Sym("a") > N=Sym("N") > b=Sym("b") >

Hi I recently gave a presentation about Automatic Differentiation (AD) and R at the Eighth Euro AD Workshop in Oxford (17/07/09). The presentation was intended as a general introduction to R and the desire for a generic AD interface for R. During the presentation I emphasised the need and the high level of interest that the R community has in developing such an interface and that input from the

Hey, what is the R function for the mathematical limit ? e.g. to calculate and return the amount that the expression X^2 +X +2 approach as X approach 2 (X-> 2) thanks hassan [[alternative HTML version deleted]]

Let's say I had a matrix like this: library(Ryacas) x<-Sym("x") m<-matrix(c(cos (x), sin(x), -sin(x), cos(x)), ncol=2) How can I use R to obtain the eigenvalues and eigenvectors? Thanks, John [[alternative HTML version deleted]]

Dear all, I am using the function 'D' in the 'stats' package to perform symbolic derivation. This works very well and it is much faster than e.g. Mathematica (at least for my purposes). First, I would like to thank the development team for this excellent function. However, I run into trouble in some cases, particularly when I am to do some operations on long expressions