similar to: Question on passing in parameter to Cox hazard

... and the mailing list is picky about attachments... whatever you attached did not conform to the stringent requirements mentioned in the Posting Guide. Pasting the code right into the email is usually safest, though you DO have to post using plain text (as the Posting Guide indicates) or your code may get mangled by the automatic html format removal. On May 15, 2018 7:04:31 AM PDT, Bert Gunter

Unless there is good reason not to, always cc the list -- there are lots of smarter folks than I on it who can help. I may or may not have time to look at this. Hopefully someone else will. -- 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

I am trying to use t.test on the following data: date type INTERVAL nCASES MTF SDF MTO SDO nFST MF nOBS MO MB BIASCV BIASEV ME MAE RMSE CRCF 2001-06-15 avn GE1.00 4385 0.246 0.300 1.502 0.556 1367 1.373 4385 1.502 1.471 0.285 0.164 -1.256 1.266 1.399 0.056 2001-06-15 avn

OK, Let's try this again! Here is the reproducible script; it is long because I had to copy the panel dataset here. My question is related to systemfit; I don't know how to get the result for the entire panel. #Reproducible script Empdata<- read.csv("/Users/ngwinuiazenui/Documents/UPLOADemp.csv") View(Empdata) install.packages("systemfit")

Hi all, simple Q: how do I extract the upper and lower CI for predicted probabilities directly for a glm - I'm sure there's a one line to do it but I can't find it. the predicted values I get with the predict (.. "response") Thanks Gerard ********************************************************************************** The information transmitted is intended only for

Hi all, I have a multiple queing situation I'd like to simulate to get some idea of the distributions - waiting times and allocations etc. Does R has a package available for this - many years ago there used to be a language called "simscript" for discrete event simulation and I was wondering if R has an equivalent (or hopefully with graphics, something better!). Apologies if there

Dear all, I've replicated the cheese tasting example on p175 of GLM's by McCullagh and Nelder. This is a 4 treatment (rows) by 9 ordinal response (cols) table. Here's my simple code: #### cheese library(MASS) options(contrasts = c("contr.treatment", "contr.poly")) y = c(0,0, 1, 7, 8,8,19, 8,1, 6,9,12,11, 7,6, 1, 0,0, 1,1, 6, 8,23,7,

Sadly you failed to set your email program to send plain text and the data is corrupted at my end. I also think you need to reduce the size of the data set... the intent here is to increase your understanding, not debug your particular analysis. I will say that I am having a very challenging time understanding what you are trying to accomplish though. What are the equations that you think need

Hi all, apologies if this is obvious - but I can't see it and would appreciate some quick help! the matrix mhouse is 26x3 and I'm computing odds ratios. The simple code below "should" compute the odds vector for every pair (325) i.e. 26C2 in cols 1 and 2. On the first i=1 outer loop the inner j loop runs from 2 to 26 ok and then I get the error (Error: subscript out of bounds)

Hi, In my analysis of impacts of insecticide-treated bednets on malaria, I look at the relationship between malaria incidence and mosquito behaviors. The condensed data set is copied here. Ordinary regression (lm) shows that Incidence was negatively related to Mortality. This makes sense because the latter reflected the strength of killing mosquitoes by insecticide-treated nets. Since the

Hello everyone, I am a new user of R. Does anybody know how hazard ratios are extracted for each factor level in a stratified Cox proportional hazard regression model? I have a cancer data set where the variable ?differentiation? is a factor with three levels: poor, intermediate and good. I would like to extract the hazard ratio for each grade level and relate it to another prognostic factor.

He, I used the "coxph" function, with four covariates. Let's say something like that > model.1 <- coxph(Surv(Time,Event)~X1+X2+X3+X4,data=DATA) So I obtain the 4 coefficients B1,B2,B3,B4 such that h(t) = h0(t) exp(B1*X1+ B2*X2 + B3*X3 + B4*X4). When I use the function on the same data > predict.coxph(model.1,type="lp") how it works in making the prediction?

I am using a syntax like coxph(Surv(start, stop, event) ~ X, data) to estimate the effect of X, which may change at each measurement (every 6 months). Is there anyone who knows a way to test the proportional hazard assumption in that case? Thank you in advance Jean-François Boudreau Sherbrooke University [[alternative HTML version deleted]]

-- begin included message --- I am hoping for some advice regarding obtaining the values for the hazard function in a cox regression that I have undertaken. I have a model in the following form, analysed with the package survival (v. 2.34-1) and a log-log plot obtained using Design (v. 2.1-2). For two variables, the lines in the survival curves crossed. The statistician I been obtaining

Dear R users, Here is the coxme output I obtain on my survival dataset having 3 strains and 2 infection status (i: infected, ni: non infected) coxme(Surv(lay) ~ infection*strain, data=datalay, random= ~1 |block) Cox mixed-effects model fit by maximum likelihood Data: datalay n= 1194 Iterations= 3 77 NULL Integrated Penalized Log-likelihood -7270.028 -7223.859 -7218.175

Dear all, Is there a way to calculate the non-cumulative hazard (instantaneous hazard), which is the product of baseline hazard and exp{beta*covariate} ? I knew in survfit, we can get the estimator of cumulative baseline hazard, but how can we get the non-cumulative one? Thank you very much! Koshihaku -- View this message in context:

Hi, Here is the summary-output of the Coxph-model I used (the output is based on the best final model i.e. all significant explanatory variables and their interactions are included): coxph(formula = Y ~ LT + Food + Temp2 + LT:Food + LT:Temp2 + Food:Temp2 + LT:Food:Temp2) n= 555 coef exp(coef) se(coef) z

What is the best way to report the standard error when publishing Kaplan-Meier plots? In my field (Vascular Surgery), practitioners loosely refer to the "10% error" cutoff as the point at which to stop drawing the KM curve. I am interpreting this as the *standard error of the cumulative hazard*, although I'm having a difficult time finding some guidelines about this (perhaps I am

Hi all, I wish to do build discrete choice model to analyse habitat selection of wolverines. This can be done with a 'tricked' stratified Cox proportional hazard model. For each individual animal each selected position, and possible alternative non-used available positions are combined into strata. This means that one stratum contains a set of 1 used position and several positions which

Dear Forum, I have a question about interaction estimate in the Cox model: why the hazard ratio of the interaction is not produced in the summary of the model? (Instead, the estimate of the coefficient is given in the print of the model.) # Example: modINT <-cph( Surv(T_BASE, T_FIN,STATUS)~ NYHA + ASINI + RFP + FE_REC + XX_PR*XX_DISF) print(modINT) coef se(coef) z