similar to: Symbolic derivation

The issue is that without an extensible derivative table or the proposed extensions, it is not possible to automatically produce (without manual modification of the deriv3 output) a function that avoids catastrophic cancellation regardless of the working range. Manual modification is not onerous as a one-time exercise, but can be time consuming when it must be done numerous times, for example

Hi. Unless I'm misremembering, log, exp, sin, cos, and tan are all handled in deriv3. The functions listed are specially coded slightly more accurate versions but can be substituted with native ones for which deriv/deriv3 will work automatically. I believe that if you write your functions using log(a + 1) instead of log1p(a) or log(x) / log(2) instead of log2(x) deriv3 will work fine.

Hi, 2 questions: Question 1: example of what I currently do: for(i in 1:6){sink("temp.txt",append=TRUE) dput(i+0) sink()} x=scan(file="temp.txt") print(prod(x)) file.remove("C:/R-2.5.0/temp.txt") But how to convert the output of the loop to a vector that I can manipulate (by prod or sum etc), without having to write and append to a file? Question 2: >

Dear R-developeRs, Attached follows a patch against svn 34959 that adds the printing of p-values to the TukeyHSD.aov function in stats package. I also updated the corresponding documentation file and added a 'see also' reference to the simint function of the multcomp package. As it was already brought up in a previous thread [1] in R-help, one can obtain the adjusted

The derivative table resides in the function D. In S+ that table is extensible because it is written in the S language. R is faster but less flexible, since that table is programmed in C. It would be useful if R provided a mechanism for extending the derivative table, or barring that, provided a broader table. Currently unsupported mathematical functions of one argument include expm1, log1p,

Hello, I was wondering if the functions deriv3(), deriv() etc. could be extended to handle psigamma() and its special cases (digamma(), trigamma() etc.). From the error message it seems that 'psigamma' needs to be added to the derivatives table. This might be easy since psigamma() has a deriv argument. Additionally, this error message is also obtained when requesting for the Hessian of

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

Greetings, I am running a factor analysis model in rjags and have received the following error Error parsing model file: syntax error on line 5 near "" Line 5 is model { Here is my syntax. Any insight would be greatly appreciated. require(rjags) modelstring = " model { for (i in 1 : nData) { for (j in 1 : nIndicators) { y[i,j] ~ dnorm(mu[i,j],psi[j]) ephat[i,j] <-

Dear participants to the list, this is my problem: I want to obtain an expression that represents the second derivative of one function. With "deriv3" (package "stats") it is possible to evaluate the second derivative, but I do not know how I can get the (analytical) expression of this derivative. For example: Suppose that I have a function of this form:

Hi all, I'm adjusting a nonlinear regression model for data that has a categorigal variable present. So, I can use nls() to do this considering the categorical variable, like this #------------------------------------------------------------ da <- expand.grid(tr=gl(2,1,la=c("tr")), x=1:12) da$y <- 10*da$x/(3+da$x)+rnorm(da$x,0,0.1) plot(y~x, da) n0 <-

Hello: Why is deriv3() functioning differently in R from that in Splus using library(MASS) ? For example deriv3(~(t1*log(t2)+lgamma(t1)+(1-t1)*log(y)+y/t2),c("t1","t2"),function(y,t1,t2)NULL) complains of lgamma. Mervyn

Dear Dr. Bates I replied to your email before but apparently it didn't make it so I am replying again. I would really appreciate if you could send me an example on how you suggest to build a cell means model for fixed effects using the NLME library. I am not sure whether you suggest to create a separate factor for each unique combination of levels of a factor or whether you suggest to

What does it require to run the example on page 263 of Venables/Ripley book using R 1.3.0? I got the following error and I have no clue how to correct it. Thank you very much. (If you reply to the list, please send a copy to me.) > lmix2 <- deriv3( + ~ -log(p*dnorm((x-u1)/s1)/s1 + (1-p)*dnorm((x-u2)/s2)/s2), + c("p", "u1", "s1", "u2",

So this is on Windows, then? This is documented somewhere (don't remember): If the package is not located under $R_HOME/library, then the html help will not be able to find it. Don't think there's a known work-around. Cheers, Andy > From: byandell@wisc.edu > > I did that. In fact the default uses .libPaths()[1], which I > have set to >

This may have been announced before but there are now both Debian GNU/Linux and RedHat Linux packages for R-0.61 on the CRAN archive sites. The master site for CRAN is http://www.ci.tuwien.ac.at/R/ That site provide a list of available mirror sites. Within the CRAN archive sites the pre-compiled packages for Linux are located in the bin/i386-linux directory. Under RedHat you can install

Using S+ 6 for Linux and R 1.4.1 Taking verbatim some lines from S-Plus that work perfectly, and running them on R, I get an unexpected error message: R : Copyright 2002, The R Development Core Team Version 1.4.1 (2002-01-30) > library(MASS) > data <- read.table("/usr/local/data/mcmanus.dat",header=T) > mcmanus.frm <-

I did that. In fact the default uses .libPaths()[1], which I have set to "P:/stat/yandell/public/statgen/R" The package installs fine. However, the HTML search is not done properly. Again, the break occurred after installation of the package, when make.search.html() was invoked. Brian On Thu, 11 Mar 2004, Douglas Bates wrote: > Brian Yandell <byandell@wisc.edu> writes: >

Martin Maechler <maechler@stat.math.ethz.ch> writes: > This caused 3 postings (1 to R-help, 2 to R-devel) to be aborted. > For some reason, it also seems majordomo did not save the postings in > a 'dead.letter' or something. > So could you please re-submit them? > These are the sendmail entries: > Sep 17 21:14:21 .. sendmail[9364]: VAA09363:

Hi, df <- read.table("data.txt", header=T); library(nls); fm <- nls(y ~ a*(x+d)^(-b), df, start=list(a=max(df->y,na.rm=T)/2,b=1,d=0)); I was using the following routine which was giving Singular Gradient, Error in numericDeriv(form[[3]], names(ind), env) : Missing value or an Infinity produced when evaluating the model errors. I also tried the

I have recompiled the r-base package for Debian-2.0 with readline support. I still haven't determined why the configure script does not detect the readline library on one machine but does detect it on another where I compiled this. The package is available as ftp://franz.stat.wisc.edu/pub/R/bin/i386-linux/Debian-2.0/r-base_0.61.1-3_i386.deb and should appear on the other CRAN mirrors by