Announcement

Collapse
No announcement yet.
X
  • Filter
  • Time
  • Show
Clear All
new posts

  • how to get adjusted annulized mean if I have linear mixed effect regression model

    Dear all,

    May I kindly have your guidance on the following please:

    I have a small sample size data set of 5 participants and 5 time points for each of them (25 observations in total). I'm reporting the change in several continuous scores over time to show the progression in the disease.
    The time points I'm having are 0, 1, 3 , 9 and 15 months.
    because I wanted to have an annualized rate of change ( to compare my results to other studies), I have fitted the following model:

    Code:
    mixed motscore yearsfscr  || participantstudyid: yearsfscr
    ( motscore is continuous dependent outcome, yearsfrscr is a variable I generated to represent the number of years from the first time point)

    from this model I get coeff, however, I want to adjust my results for age, gender, and allele size

    I have applied the following model:
    Code:
    mixed motscore age sex allele yearsfscr  || participantstudyid: yearsfscr
    ( not sure if the code is the right one, I have also tried something like :
    Code:
    mixed motscore age sex alle##yearsfscr  || participantstudyid: yearsfscr
    but I get the following error msg:
    yearsfscr: factor variables may not contain noninteger values


    the other question is: Other studies that I want to relate my results to, have presented adjusted annualized change in the mean, Is there a way to get adjusted annualized mean from my data to present my results similary rather than presenting the Coeff?



    Code:
    * Example generated by -dataex-. To install: ssc install dataex
    clear
    input byte participantstudyid float age byte(sex allele2l visnum) float yearsfscr byte motscore
    1 54 1 43 1          0 15
    1 54 1 43 2  .07664303 14
    1 55 1 43 3  .29014865 15
    1 55 1 43 4   .7664304 15
    1 56 1 43 5  1.4397942 22
    2 54 0 40 1          0  3
    2 54 0 40 2  .07390578 11
    2 55 0 40 3   .3038349 15
    2 55 0 40 4   .7636931  9
    2 56 0 40 5  1.4343196 12
    3 59 0 40 1          0  7
    3 59 0 40 2  .07664303 13
    3 59 0 40 3   .3038349 11
    3 60 0 40 4   .7691676 12
    3 60 0 40 5  1.4890647 13
    4 44 0 44 1          0 31
    4 44 0 44 2 .071168534 26
    4 44 0 44 3   .3010976 24
    4 45 0 44 4   .7609559 28
    4 45 0 44 5   1.568445 40
    5 49 0 44 1          0 31
    5 49 0 44 2  .06295678 33
    5 50 0 44 3  .27372512 21
    5 50 0 44 4   .7527441 25
    5 51 0 44 5  1.4753784 36
    end
    N.B I'm using Stata 15.1
    Many thanks in advance
Working...
X