search for: chd

...ged was signaled while in a function called from GDB. GDB remains in the frame where the signal was received. To change this behavior use "set unwindonsignal on" Evaluation of the expression containing the function (rb_p) will be abandoned. (gdb) eval "caller()" [weavere at c13-chd-app1 tools]$ Program received signal SIGINT, Interrupt. 0x00000000004132d9 in rb_thread_schedule () at eval.c:10618 10618 in eval.c The program being debugged was signaled while in a function called from GDB. GDB remains in the frame where the signal was received. To change this behavior use &quo...

Dear r-users, Does anyone know how to draw time-dependent survival curves? Example: Event outcome: CHD Time-dependent covariate: NSAID use, which changes over time for each subject I'm interested in survival curves stratified by NSAID use. I'd like to implement Simon & Makuch (1984) method. Is there a R package/function to draw this graph? ------------------------------------------...

Hi, We're using stats4 for a logistic regression. The code is chdreg.logit2 <- glm(chd ~ age + sex, family = binomial) summary(chdreg.logit2) oddsratios <- coef(chdreg.logit2) exp(oddsratios) # Calculate model predicted values pred <- predict(chdreg.logit2,type="response") The glm part runs fine, and up to now so has the predict function. H...

...and want to calculate an odds ratio between 2 sets of predictors. It is easy to obtain the difference in the predicted logodds using the predict() function, and thus get a point-estimate OR. But I can't see how to obtain the confidence interval for such an OR. For example: model <- glm(chd ~age.cat + male + lowed, family=binomial(logit)) pred1 <- predict(model, newdata=data.frame(age.cat=1,male=1,lowed=1)) pred2 <- predict(model, newdata=data.frame(age.cat=2,male=0,lowed=0)) OR <- exp(pred2-pred1) Thanks [[alternative HTML version deleted]]

Dear Listers, Apologize to pile up on the 'tz' issue in POSIX objects. I have a 'simple' thing on which I must make up my mind but cannot do it from the existing R-help threads. I am currently working on dog telemetry in China, and download time information from GPS collars. I would like to set up the corresponding POSIXxx variables in R to a given time zone. Eg Pekin