Displaying 5 results from an estimated 5 matches for "nsaid".
Did you mean:
said
2007 Mar 14
0
Logistic regression for drugs Interactions
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 ***
NSA...
2006 May 11
1
time-dependent covariate survival curves
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?
------------------------------------------
Two related papers:
Simon R, M...
2007 Jun 01
0
Metropolis code help
...the hierarchical version.
Here is the original code of the teacher with both flat prior on betas and a
hierarchical version:
www.stats.uwo.ca/faculty/murdoch/458/metropolis.r
Below is My code with a flat prior on beta only (I'd like also to have the
hierarchical version!)
X<- cbind(1,DF$nsaid,DF$diuretic,DF$diuretic*DF$nsaid)
y<- DF$Var3
Metropolis <- function(logtarget, start, R = 1000, sd
= 1) {
parmcount <- length(start)
sims <- matrix(NA, nrow=R, ncol = parmcount)
colnames(sims) <- names(start)
sims[1,] <- start
oldlogalpha <- logtarget(sta...
2007 May 03
1
Bayesian logistic regression with a beta prior (MCMClogit)
...essage (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){
sum(dbeta(beta,shape1,shape2, log=T)) }
posterior <- MCMClogit(nausea~nsaid*diuretic, data=w,
verbose=2000,burnin = 1000, mcmc = 10000,
user.prior.density=logpriorfun,shape1=1,shape2=1)
user.prior.density(beta.start) == -Inf.
Error in MCMClogit(nausea ~ nsaid * diuretic, data = w, verbose = 2000, :
Respecify and call MCMClogit() again.
--
View this message in...
2000 Sep 19
1
Re: more experience with formulas
...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 return...