similar to: Derivative of a function

I have 3 basic questions (sorry, but R is a bit new to me) First) The following code will work on C, but something happens with R: > x<-c(0,1,1,0,0,1,0,0,1) > dim(x)<-c(3,3) > types<-numeric(3) > for (i in 1:3) { + j<-1 + while (x[i,j]!=1) j<-j+1 + if (j==4) types[i]<-0 + else types[i]<-1 + } Error: subscript out of bounds Any ideas about what is wrong?

Hi, all. Thanks Timo for answering my previous post. It's starting fine now :) I must be doing something wrong again but I can't connect from any imap client. If I telnet to the server I get: $ telnet carpa 11143 Trying 143.107.209.25... Connected to carpa.ciagri.usp.br. Escape character is '^]'. * OK dovecot ready. a001 login marcelo ##### a001 OK Logged in. a002 select INBOX

Hi. Sorry for the basic question but... how do I start the server? From inetd or in standalone mode (which binary)? Thanks

Hello everyone, I have the following issue: one function generates a very big array (can be more than 1 Gb) and returns a few variables, including this big one. Memory allocation is OK while the function is running, but the final steps make some copies that can be problematic. I looked for a way to return the values without copying (even tried Rmemprof), but without success. Any ideas? The code

Hi, I'm an Apache James committer and we are curious to use imaptest in order to validate our IMAP protocol implementation. I'm using the nightlybuild : imaptest-20170506 I follow the examples given in the /T/ /est IMAP server compliancy/ examples https://imapwiki.org/ImapTest/Examples. And I wanted to analyse the errors reported, here are some: Error: user1 at james.org[11]: seq too

Hi R-users! I'd like to know how to avoid automatic coercion of numeric variables passed to array() to character, when there is one categorical variable between others. I have tried I() function, but it is just functional in data.frame. Best Regards! Eduardo Dutra de Armas

Hello list, I would like to know how can I use contrats in statistical analysis. I have a simple model, "one-way analysis of variance" "or two-way..." (balanced case). I've already used the command aov (model) to get the analysis of variance, but i know that my treatment could be organized in groups, so i would like to test these groups. Example: I have 4 treatments t1,

Hello everyone, I'm using the very good imaptest [0] tool to test my little imap server implementation. I've tried to use the dovecot-crlf [1] file, but it looks like there are some major issues : $ grep -n "In-Reply-To.*;" tests/data/dovecot-crlf 479:In-Reply-To: <20020806175441.GA7148 at linux.taugt.net>; from rueckert at informatik.uni-rostock.de on Tue, Aug 06, 2002

While writing a function that includes placing grid lines at the same position as the axis ticks, I found that the axis* functions don't return anything. Thus I have had to copy the appropriate function, removing the call to axis() and adding a line to return the tick positions. Is there a more elegant way to determine the tick positions on an axis? Thanks. Jim (normally bitwrit at

I fit my CDF to sum of exponentials and now I want to take the numerical derivative of this function to obtain probability density.I will really appreciate your help reagrding the error messages I am getting which I don't understand. * * > fitterma <- function(xtime) { a <- -0.09144115 b <- -0.01335756 c <- -2.368057 d <- -0.00600052

Hello, I am wondering if someone can help me. I have the following function that I derived using nls() SSlogis. I would like to find its derivative. I thought I had done this using deriv(), but for some reason this isn't working out for me. Here is the function: asym <- 84.951 xmid <- 66.90742 scal <- -6.3 x.seq <- seq(1, 153,, 153) nls.fn <- asym/((1+exp((xmid-x.seq)/scal)))

Dear All, I tried following for getting derivative of a polynomial in R i<- -10:10 x<-i*i*i+3*i*i+2 fun_spline<-splinefun(i,x) plot(x,type="l") lines(x,fx_spline(x, deriv=1), col='green') lines(x,fx_spline(x, deriv=2), col='green') Now when I plot 3*i*i + 6*i and 6*i + 6 the plot was not same for first deivative. where as the 2nd derivative was same Is this a

Dear R experts. 'D()' function recognizes some of the analitical functions, such as sin, cos, etc. But I'd like to take analytical derivatives from asin, atan etc. functions. Are there any R packages providing that features? Thanks. -- Timur.

Hello, I have been trying to take the derivative of a quadratic B-spline obtained by using the COBS library. What I would like to do is similar to what one can do by using fit<-smooth.spline(cdf) xx<-seq(-10,10,.1) predict(fit, xx, deriv = 1) The goal is to fit the spline to data that is approximating a cumulative distribution function (e.g. in my example, cdf is a 2-column matrix with x

Hello, I am trying to find a hessian matrix that involves log(gamma(1/p)) second derivative, p being one of the parameters of the function. I am using a function "deriv" with the hessian=TRUE option, but psigamma is not on the list of derivative functions. I know that it is possible to use 'psigamma(p,deriv)', but it doesn't work with 1/p. Does anybody can help with this?

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

Readers, A data set comprises A B C 10 6 .2 20 7 .4 30 8 .16 40 9 .0256 My requirement is to obtain the derivative for values of A with respect to B, create a function in R and plot this derivative against another variable (e.g. values in column C). Have searched the mailing list and found reference to a function 'D', but the help accessed via '?D' and '?deriv' does not

Hi all, I am trying to read multiple out variables for a sensitivity analysis. Currently using one output value as follows: Y<-(E1) However I need to run analysis against 12 values of Y. So E1-E12. My matrix will be: Inputs are Column=4, Rows = 40 i.e. 40 rows of 4 input variables in different combinations. These will be analysed against 40 rows of output variables for 12 columns.

Hello, everyone. I have the following problem with TclTk: I create some windows and want to change their position with geometry manage (sometimes they will be centered, sometimes not). If the toplevel is created and its dimensions are gathered via 'tkwinfo', I get (usually) correct values. However, if this window is created by a function (in the following example, by 'ask.format')

Is there a way of calculating the derivative of a function returned by splinefun()? Such a function is a cubic spline, whence it has a calculable derivative, but is there a (simple) way of getting at it? One workaround that I have thought of is to take a fine grid of points, evaluate the function returned by splinefun() at these points, put an interpolating spline through these points using