search for: nsaids

I have the model below, for which I run a logistic regression including the interaction term (NSAID*Diuretic) ------------------------ fit1=glm(resp ~ nsaid+diuretic+I(nsaid*diuretic), family= binomial,data=w) NSAID Diuretic Present Absent 0 0 185 6527 0 1 53 1444 1 0 42 1293 1 1 25 253 Coefficients Std. Error z value Pr(>|z|) (Intercept) -3.56335 0.07456 -47.794 < 2e-16 ***

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?

Dears, I have the below code for metropolis of the GLM logit (logistic regression) using a flat prior. Can someone help me modify the prior so that the model becomes hierarchical by using a flat prior for mu and sigma, the derived density for beta ~ N(mu, sigma^2)? Actually I took my code from a teacher that posted on the internet and modified it to the GLM logit but I can't adapt it to the

Dear all, I am trying to use the logistic regression with MCMClogit (package: MCMCpack/Coda) and I want to put a beta prior on the parameters, but it's giving me error message (please see output below) no matter what shape 1 or 2 I use. It works perfect with the cauchy or normal priors. Do you know if there is a catch there somewhere? Thanks logpriorfun <- function(beta,shape1,shape2){

...tely] WSt> I do not want to call this a bug, but it bugged me. Please try the WSt> following: WSt> tform <- sqrt(RADAI) ~ sqrt(RADAI.e) + TAGE.ej + SPITAL + ARZT + DAS28 + WSt> SJC + TJC + DGA + HAQ + PGA + PAIN + YEAR.SYM + YEAR.DIA + WSt> ALTER + IndV + m.NSAIDs + m.Stereoids + DM.Cyclosp + DM.Penic WSt> WSt> as.character(tform) --> length-3 vector "~" "<LHS>" "<RHS>". However, the end of <RHS> is TRUNCATED! and that truncation of the third part is a bug in my opinion. The fact that it returns...