similar to: Delta method using numerical derivatives

Dear R community, I've been writing simple functions for the past year and half and have come across a similar problem several times. The execution of a function within my own function produces NaN's or fails to execute as intended. My conundrum is that I can execute the function outside of my function without error, so it's difficult for me, as a novice functioneer, to figure out

Hi, I'm trying to write a simple function like case1 <- function (m, cov, Q, R) { theta <- (acos(R/sqrt(Q^3))) beta <- (-2)*sqrt(Q)*cos(theta/3)+m[1]/3 rho1 <- (-2)*sqrt(Q)*cos((theta+2*pi)/3)+m[1]/3 rho2 <- (-2)*sqrt(Q)*cos((theta-2*pi)/3)+m[1]/3 stderrb <- deltamethod( ~(-2)*sqrt(Q)*cos(theta/3)+x1/3,m,cov) stderrr1 <- deltamethod(

Stephen, thanks for your effort. The more appropriate list for this discussion is probably R-devel (as far as I understand it) so I've moved the discussion there. Related topics have already been discussed in the past. Specifically, I remember contributions by Paul Johnson ("rockchalk" package) and John Fox ("effects" and "car" package) as their packages

Hi R Team I've got the following problem I'd like to run a time series regression of the following form Regression1: At = ? + ?1 * Bt + ?2 * Bt-1 + ?3 [(Bt-2 + Bt-3 + Bt-4)/3] + ?t The B's are the input values and the A's are the output values, the subscript stands for the lag. The real Beta of this regression is ?real = ?1 + ?2 + ?3 First: How can I run the regression

I have code and documentation that don't match, but R CMD check didn't flag it. in mspath.R mspath <- function(formula, # formula with observed Markov states ~ observation times (required) qmatrix, # matrix of 1s and 0s with indices of allowed transitions (diagonal is ignored) (required) misc = FALSE, ematrix = NULL, # matrix

It's the same as PR#2993, fixed yesterday in R-patched & R-devel On Fri, 30 May 2003 rossini@blindglobe.net wrote: > > Something like the following should be done for power.t.test to make > sure that it doesn't try to evaluate abs(NULL), which results in an > error. That's a backwards patch, I presume, or you revert the fix. > --- rossini.power.t.test.R

Something like the following should be done for power.t.test to make sure that it doesn't try to evaluate abs(NULL), which results in an error. --- rossini.power.t.test.R 2003-05-30 07:24:49.000000000 -0700 +++ rossini.power.t.test.R.~1~ 2003-05-30 08:47:09.000000000 -0700 @@ -10,7 +10,7 @@ alternative <- match.arg(alternative) tsample <- switch(type, one.sample = 1,

Dear friends: I'm interested to make a stacked plot of cumulative incidence. that's, the cuminc model is fitted [fit=cuminc(time, relapse)] and cumulative incidence is in place. I'd like to stack the cuminc plots (relapse of luekemia and death free from leukemia, for example) , then the constituent ratio of leukemia relapse and treatment related mortality is very clear. Can

Dear R-help list members, I am trying to create my own axes when plotting a cumulative incidence curve using the plot.cuminc function in the CMPRSK library. The default x-axis places tick marks and labels at 0, 20, 40, 60, and 80 (my data has an upper limit of 96), whereas I want them at my own specified locations. Here is my example code: library(cmprsk) attach(MYDATA) MYCUMINC <-

Hello everyone, (This is my second question posted today on the R list). I am carrying out a competing risks analysis using the cuminc function...this takes the form: cuminc(ftime,fstatus,group) In my study, fstatus has 3 different causes of failure (1,2,3) there are also censored cases (0). "group" has two levels (0 and 1). I therefore have 6 different cumulative incidence curves:

Michael Rock napsal(a): > Hi Dave, > > I am using applydeltaiso with your delta and > CentOS-5.0-x86_64-bin-DVD.iso 4,186,786 bytes from one > of the Centos mirrors. But the block size is wrong and > does not match the image length. > > So I checked the md5 on your delta and it does not > match your md5.txt file. I downloaded it again to be > sure. > > Your

> David wrote: > Now there are right md5sums and even > x86_64 iso.delta. > David I tried using applydeltaiso and it creates a 2gb ISO but I notice it fails on applying. kdnssd-avahi-devel.i386: copying unchanged payload kernel-xen.x86_64 (???): applying delta new payload open failed __________________________________________________ Do You Yahoo!? Tired of spam? Yahoo! Mail has

btw - I installed your x86_64 rpms that has the ich9 patch directly on that box but it still will not boot without all the options. --- Michael Rock <mikerocks65 at yahoo.com> wrote: > > David wrote: > > Now there are right md5sums and even > > x86_64 iso.delta. > > David > > I tried using applydeltaiso and it creates a 2gb ISO > but I notice it fails on

David wrote: >I'll test DVD tomorrow. Today new kernel RPMS are > being built. What ptions are you talking about? > David all-generic-ide irqpoll Without irqpoll it just hangs and without all-generic-ide it get a kernel panic. So I still need both to recognize the disk and continue booting. __________________________________________________ Do You Yahoo!? Tired of spam? Yahoo!

Hi all, I've the same problem, I've download the CentOS-5.0-x86_64-bin-DVD.iso from one centos mirror, i've applied deltaiso and the created image is 2 GB, the md5sum check of iso.delta (CentOS-5.0-x86_64-bin-DVD-ich9.iso.delta) is ok, what's wrong ? deltaiso program exit on the same as Michael , maybe I've downloaded the worng DVD 64 bit iso version? do you know any good

> WELCOME TO MY WORLD :D > (Hey, are the new foobar packages done, how about the xwyz packages) > Thanks, > JOhnny Hughes When you snip the rest of my post man that last request sounds really bad on its own. Sorry guys. I saw the "uploading the new DVD's" a few days ago and thought I would ask. I really do not want to put anyone out and please do not go through any

Oops sorry, ugh another poor selection of words as I meant to say "When I saw the rest of my post snipped, I then realized that last request sounds really bad, which was not my intention. > When you snip the rest of my post man that last > request sounds really bad on its own. Sorry guys. I > saw the "uploading the new DVD's" a few days ago and > thought I would

Skipped content of type multipart/alternative-------------- next part -----= --------- A non-text attachment was scrubbed... Name: not available Type: image/gif Size: 10923 bytes Desc: not available Url : http://lists.linux-foundation.org/pipermail/virtualization/attachment= s/20070401/0fe6b7ca/attachment-0001.gif

Dear list members, I am using a cross validation of a generalised linear model (glm). The cv.glm function (from boot package) returns an error as so-called ?delta? value. I would like to get to a (cross-validated) squared q, because I want to directly compare it to the squared correlation coefficient r. I tried to find an an equation for the raw and/or adjusted cross-validation estimate of

Hi! When I use the glm.cv function I get a value called "delta" which is explained as the "raw cross-validation estimate of prediction error". I recently found a formula for that term in literature where it is defined as: alpha = 1 / N * sum over( yi - yi,pred,CV) Well it is somehow similar to the RSS for R2 and the PRESS for Q2. But this delta value increases with