Displaying 5 results from an estimated 5 matches for "lifedxm".
Did you mean:
lifed
2010 Aug 04
5
Question regarding significance of a covariate in a coxme survival model
...my dataset, I have
decided to account for this by including a random effect for family
(famid). My covariate of interest is Mother's diagnosis where a 0 is
bipolar, 1 is control, and 2 is major depression. I am trying to fit
the following model.
thorp1 <- coxme(Surv(age_sym1, sym1) ~ lifedxm + (1 | famid), data =
bip.surv)
Which provides the following output:
-------------------------------------------------
> thorp1
Cox mixed-effects model fit by maximum likelihood
Data: bip.surv
events, n = 99, 189 (3 observations deleted due to missingness)
Iterations= 10 54...
2010 Jul 23
1
Survival analysis MLE gives NA or enormous standard errors
Hi,
I am trying to fit the following model:
sr.reg.s4.nore <- survreg(Surv(age_sym4,sym4), as.factor(lifedxm),
data=bip.surv)
Where age_sym4 is the age that a subject develops clinical thought
problems; sym4 is whether they develop clinical thoughts problems (0 or
1); and lifedxm is mother's diagnosis: BIPOLAR, MAJOR DEPRESSION, or
CONTROL.
I am interested in whether or not survival differs by this...
2010 Jul 27
0
AIC from coxme
Hi,
I am running the following model:
fit1.full <- coxme(Surv(age_sym1, sym1) ~ sex + lifedxm*sex + (1|famid),
data=bip.surv)
I would like to extract the AIC from that object to calculate the AICC.
However, when I look at str(fit1.full) and summary(fit1.full) (pasted
below) I don't see anything that would allow me to get pull the AIC out
from that object.
Is there a way to retrieve t...
2011 May 18
1
Changing order of facet grid in ggplot2
Hi I am running the following code:
sym <- c(sym1,sym2,sym4)
lifedxm <- c("O-BD","O-WELL","O-UNI")
life <- c(lifedxm,lifedxm,lifedxm)
tp <- c("TP-ANY","TP-ANY", "TP-ANY", "TP-SUB", "TP-SUB", "TP-SUB", "TP-CLIN"
, "TP-CLIN", "TP-CLIN")
d...
2010 Jul 22
2
Multilevel survival model
* Please cc me if you reply as I am a digest subscriber *
Hi,
I am wondering how I can run a multilevel survival model in R? Below is
some of my data.
> head(bi0.test)
childid famid lifedxm sex age delta
1 22.02 22 CONTROL MALES 21.36893 0
2 13.02 13 MAJOR MALES 21.18001 0
3 64.02 64 CONTROL MALES 20.09377 0
4 5.02 5 CONTROL FEMALES 19.94524 0
5 183.02 183 BIPOLAR FEMALES 19.87406 0
6 37.02 37 BIPOLAR MALES 19.83847...