- 89 Views
- Uploaded on
- Presentation posted in: General

Different Distributions

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

Different Distributions

David Purdie

Application of GEE to:

- Binary outcomes:
- logistic regression

- Events over time (rate):
- Poisson regression

- Survival data:
- Cox regression

Outcome for subject i at time j = Yij

E(Yij)=ij

Generalized linear model

g(ij)=Xi

where Xi=(xi1,…,xij) is the matrix of covariates for subject i

Outcome:

Pr(Yij = 1) = ij (probability of an event)

Pr(Yij = 0) = 1- ij.

Logit link function:

Logistic model:

where ij = E(Yij|Xi)

- Outcome: Yi= number of events in time period ti
- E(Yi): iti Var(Yiti)= iti (were i is the event rate)
- Log link function: log (i)
- Poisson model:

Parameter: tij (time to event yij)

Based on a hazard function: ht

Outcome: Tij= time till event yij

Log link function: log (ht)

Cox model:

where t is the baseline hazard rate.

If the responses are binary, it may make more sense to use a matrix of odds ratios rather than correlations.

Replace corr(Yij, Yik) with:

The ALR algorithm models ijk = log{OR(Yij,Yik)} as:

ijk=zijk

where are regression parameters and

z is fixed and needs to be specified

- E(y|u)=, var(y|u)=V(), g()=X+Zu
- Random coefficients u have dist f(u)
- y|u has the usual glm distribution
- Binary outcome:
- binomial for y|u and beta for u

- Count outcome:
- Poisson for y|u and gamma for u

- Study of bladder cancer
- All patients had superficial bladder tumours on entry which were removed
- Two randomly allocated treatments (group):
- Placebo (n=47), Thiotepa (n=38)

- Many multiple recurrences of tumours
- Month is month since treatment (1 to 53)
- Baseline covariates of number of initial tumours (number) & size of largest tumour (size)
- Lots of missing data: 3585 out of 4505 potential observations (80%) are missing
- Model missing data (yes/no) using a binomial GEE to assess if data is missing at random (logit link function)

Name in data set

procgenmod data=tumour_miss descending;

class group subject month;

model missing=group month size number

/ dist=binomial type3;

repeated subject=subject

/ type=ind corrw within=month;

estimate 'effect of thiotepa' group -11/ exp;

run;

- Response: number of new tumours (count)
- Month is month since treatment (1 to 53)
- Baseline covariates of number of initial tumours (number) & size of largest tumour (size)
- Timesince is the number of months since the last visit
- Missing data are dependent upon treatment group and time
- Model new tumour counts using a Poisson GEE to assess treatment effect (log link function)

Name in data set

Count of tumours by treatment group

procgenmod data=tumour_count;

class group subject month;

model count=group size number timesince

/ dist=poisson scale=deviance;

repeated subject=subject

/ type=exch withinsubject=month corrw;

estimate 'effect of thiotepa' group -11/ exp;

run;

*WARNING: The number of response pairs for estimating correlation is less than or equal to the number of regression parameters. A simpler correlation model might be more appropriate.

data tumour_count;

set tumour_count;

off=log(timesince+1);

run;

proc genmod data=tumour_count;

class group subject month;

model count=group size number

/ dist=poisson scale=deviance offset=off type3;

repeated subject=subject

/ type=unstr withinsubject=month;

estimate 'effect of thiotepa' group -1 1/ exp;

run;

*WARNING: The number of response pairs for estimating correlation is less than or equal to the number of regression parameters. A simpler correlation model might be more appropriate.

- Descriptive: plots of means or tables of means (percentages, etc.)
- Tables of parameter estimates and confidence intervals (odds ratios or relative risks)
- P-values for effects or interactions (possibly just in the text)
- Emphasize results from descriptive analysis and effect estimates.

- What is the distribution of the outcome?
- How were the data summarized?
- Due to the repeated nature of the data, a generalized estimated equations (GEE) approach was used to estimate parameters and test for differences between groups.
- What was the form of the correlation structure?
- What hypotheses were being tested?
- How were missing data handled?
- How were variances calculated?
- What statistical package was used?

Mean numbers of new tumours were used to summarise the data. Poisson regression was used to model tumour counts using the time between successive observations as an offset. Due to the repeated nature of the data, a generalized estimated equations (GEE) approach was used to estimate parameters and test for differences between groups. The main hypothesis being tested was whether Thiotepa affected the numbers of new tumours. The correlation between successive observations was examined and an appropriate correlation structure was specified. Drop outs and non-attendance was examined to assess for differences between the treatment groups. Robust variance estimate techniques were used to calculate standard errors and confidence intervals. All analysis were performed using SAS version 8.2.