[R] Help Interpreting Linear Mixed Model

Joshua Dixon joshuamichaeldixon at gmail.com
Tue Apr 28 22:26:31 CEST 2015


*Edit*

Where "F" position are in the same AgeGr as well.

Thanks,

Josh

On Tue, Apr 28, 2015 at 9:25 PM, Joshua Dixon <joshuamichaeldixon at gmail.com>
wrote:

> *John* -  Lot's of missing data for height unfortunately.  Which is
> needed for BMI calculation.
>
> How would I look compare very specific parts of the data, i.e. comparing
> YoYo outcomes between "F" and "M" position that are both in the
> PREMIER_LEAGUE Level?
>
> Still can't figure it out!
>
> Josh
>
> On Tue, Apr 28, 2015 at 2:39 AM, John Kane <jrkrideau at inbox.com> wrote:
>
>>
>> Looks great.  How come so many NA's in Height and BMI? Just no data
>> available?
>>
>>  str(dat1)
>> 'data.frame':   100 obs. of  8 variables:
>>  $ Id      : int  7451 148 10393 10200 1961 10428 10541 10012 9895 10626
>> ...
>>  $ Level   : Factor w/ 5 levels "CHAMPIONSHIP",..: 1 1 1 1 1 1 1 1 1 1 ...
>>  $ AgeGr   : int  14 16 10 10 13 10 10 10 10 10 ...
>>  $ Position: Factor w/ 4 levels "D","F","GK","M": 4 1 1 2 3 3 2 3 1 1 ...
>>  $ Height  : int  NA NA NA NA NA NA NA NA NA NA ...
>>  $ Weight  : num  63 64 36 46 67 40 25 30 36 33 ...
>>  $ BMI     : num  NA NA NA NA NA NA NA NA NA NA ...
>>  $ YoYo    : int  80 80 160 160 160 160 160 160 160 160 ...
>>
>> John Kane
>> Kingston ON Canada
>>
>> -----Original Message-----
>> From: joshuamichaeldixon at gmail.com
>> Sent: Mon, 27 Apr 2015 23:35:13 +0100
>> To: jrkrideau at inbox.com
>> Subject: Re: [R] Help Interpreting Linear Mixed Model
>>
>> Thanks John!
>>
>> This ok?
>>
>> > dput(head(data, 100))
>>
>> structure(list(Id = c(7451L, 148L, 10393L, 10200L, 1961L, 10428L,
>>
>> 10541L, 10012L, 9895L, 10626L, 1151L, 8775L, 10083L, 6217L, 90L,
>>
>> 10168L, 10291L, 8549L, 3451L, 10003L, 5907L, 10136L, 6182L, 6315L,
>>
>> 10015L, 9956L, 2040L, 4710L, 10747L, 6787L, 1222L, 10757L, 2892L,
>>
>> 117L, 10328L, 10503L, 768L, 2979L, 1961L, 10520L, 10498L, 3018L,
>>
>> 10335L, 2448L, 9027L, 362L, 8499L, 10603L, 9489L, 2124L, 707L,
>>
>> 8501L, 4908L, 9905L, 3000L, 2819L, 9973L, 10550L, 9921L, 10639L,
>>
>> 8771L, 10121L, 32L, 9935L, 9299L, 3246L, 682L, 10325L, 6741L,
>>
>> 3295L, 5270L, 727L, 8500L, 50L, 4705L, 3018L, 787L, 2953L, 1391L,
>>
>> 3682L, 7974L, 5023L, 652L, 727L, 679L, 10212L, 9488L, 9987L,
>>
>> 10039L, 5025L, 250L, 2539L, 787L, 3000L, 1151L, 8946L, 6177L,
>>
>> 3296L, 250L, 498L), Level = structure(c(1L, 1L, 1L, 1L, 1L, 1L,
>>
>> 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
>>
>> 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
>>
>> 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
>>
>> 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
>>
>> 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
>>
>> 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), .Label =
>> c("CHAMPIONSHIP",
>>
>> "CONFERENCE", "LEAGUE_ONE", "LEAGUE_TWO", "PREMIER_LEAGUE"), class =
>> "factor"),
>>
>>     AgeGr = c(14L, 16L, 10L, 10L, 13L, 10L, 10L, 10L, 10L, 10L,
>>
>>     14L, 10L, 10L, 10L, 12L, 10L, 10L, 12L, 10L, 10L, 10L, 10L,
>>
>>     12L, 10L, 10L, 10L, 10L, 10L, 10L, 15L, 10L, 10L, 10L, 12L,
>>
>>     10L, 10L, 13L, 10L, 13L, 11L, 11L, 13L, 12L, 11L, 12L, 14L,
>>
>>     13L, 13L, 13L, 13L, 12L, 11L, 15L, 11L, 14L, 13L, 11L, 11L,
>>
>>     11L, 12L, 14L, 12L, 13L, 11L, 13L, 15L, 11L, 13L, 13L, 13L,
>>
>>     14L, 13L, 13L, 12L, 13L, 13L, 13L, 14L, 12L, 14L, 13L, 13L,
>>
>>     13L, 13L, 13L, 12L, 13L, 14L, 13L, 14L, 13L, 14L, 13L, 14L,
>>
>>     14L, 13L, 14L, 13L, 13L, 13L), Position = structure(c(4L,
>>
>>     1L, 1L, 2L, 3L, 3L, 2L, 3L, 1L, 1L, 1L, 2L, 4L, 3L, 2L, 3L,
>>
>>     4L, 3L, 4L, 2L, 4L, 2L, 3L, 1L, 1L, 2L, 4L, 4L, 2L, 4L, 4L,
>>
>>     2L, 1L, 4L, 1L, 1L, 2L, 4L, 3L, 1L, 4L, 1L, 2L, 3L, 3L, 1L,
>>
>>     1L, 3L, 1L, 3L, 4L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 3L, 1L,
>>
>>     2L, 2L, 2L, 4L, 4L, 2L, 4L, 2L, 3L, 3L, 4L, 4L, 1L, 1L, 1L,
>>
>>     2L, 2L, 4L, 1L, 1L, 1L, 2L, 4L, 1L, 3L, 4L, 4L, 4L, 4L, 2L,
>>
>>     2L, 2L, 1L, 1L, 4L, 1L, 4L, 2L, 2L), .Label = c("D", "F",
>>
>>     "GK", "M"), class = "factor"), Height = c(NA, NA, NA, NA,
>>
>>     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
>>
>>     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
>>
>>     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 151L, NA,
>>
>>     154L, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 156L, NA,
>>
>>     147L, NA, NA, NA, NA, NA, 138L, 172L, NA, NA, 150L, NA, NA,
>>
>>     NA, NA, NA, NA, NA, 140L, 153L, NA, NA, NA, NA, NA, NA, NA,
>>
>>     158L, NA, NA, NA, NA, NA, NA, NA, NA, NA, 156L), Weight = c(63,
>>
>>     64, 36, 46, 67, 40, 25, 30, 36, 33, 61, 31, 29, 34, 47, 38,
>>
>>     32, 44, 32, 32, 30, 34, 51, 34, 28, 27, 33, 31, 28, 44, 37,
>>
>>     46, 26, 42, 32, 32, 43, 31, 72, 27, 30, 55, 53, 50, 51, 55,
>>
>>     48.6, 49, 48, 64, 35, 32, 55, 32, 50, 61, 42, 33, 37, 45,
>>
>>     45, 50, 36, 33, 49, 59, 42, 43, 35.1, 66.9, 52, 47, 40, 38,
>>
>>     45, 53, 44, 54, 39, 62, 33, 53.8, 42, 46, 39, 48, 39, 54,
>>
>>     40, 42.4, 50, 48, 46, 52, 58, 40, 46, 51, 54, 42), BMI = c(NA,
>>
>>     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
>>
>>     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
>>
>>     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
>>
>>     21.2, NA, 20.24, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
>>
>>     NA, 18.49, NA, 16.66, NA, NA, NA, NA, NA, 18.57, 22.61, NA,
>>
>>     NA, 17.77, NA, NA, NA, NA, NA, NA, NA, 16.84, 22.86, NA,
>>
>>     NA, NA, NA, NA, NA, NA, 16.9, NA, NA, NA, NA, NA, NA, NA,
>>
>>     NA, NA, 17.26), YoYo = c(80L, 80L, 160L, 160L, 160L, 160L,
>>
>>     160L, 160L, 160L, 160L, 160L, 160L, 160L, 160L, 160L, 160L,
>>
>>     160L, 160L, 160L, 160L, 160L, 160L, 160L, 160L, 160L, 160L,
>>
>>     160L, 160L, 160L, 160L, 160L, 160L, 160L, 160L, 160L, 160L,
>>
>>     160L, 160L, 200L, 200L, 200L, 200L, 200L, 200L, 200L, 200L,
>>
>>     200L, 200L, 200L, 200L, 200L, 200L, 200L, 200L, 200L, 200L,
>>
>>     200L, 200L, 200L, 200L, 200L, 200L, 200L, 200L, 200L, 200L,
>>
>>     200L, 240L, 240L, 240L, 240L, 240L, 240L, 240L, 240L, 240L,
>>
>>     240L, 240L, 240L, 240L, 240L, 240L, 240L, 240L, 240L, 240L,
>>
>>     240L, 240L, 240L, 240L, 240L, 240L, 240L, 240L, 240L, 240L,
>>
>>     240L, 240L, 240L, 240L)), .Names = c("Id", "Level", "AgeGr",
>>
>> "Position", "Height", "Weight", "BMI", "YoYo"), row.names = c(NA,
>>
>> 100L), class = "data.frame")
>>
>> On Mon, Apr 27, 2015 at 10:43 PM, John Kane <jrkrideau at inbox.com> wrote:
>>
>>  Hi Josh,
>>
>>  Just a sample  is usually  fine. As long as it cover a representative
>> (must be time for dinner---I was going to type reprehensibe) sample of the
>> data then something like dput(head(mydata, 100) ) works well.
>>
>>  Kingston ON Canada
>>
>>  -----Original Message-----
>>  From: joshuamichaeldixon at gmail.com
>>
>> Sent: Mon, 27 Apr 2015 21:30:39 +0100
>>  To: lists at dewey.myzen.co.uk
>>  Subject: Re: [R] Help Interpreting Linear Mixed Model
>>
>>  Apologies for my ignorance!
>>
>>  Thierry - thank you for the reading.  I'll look into those ASAP!
>>
>>  John - The data set I have is quite large, when using the dput() command
>> I'm unsure if it actually fits the whole output into the console.  I can't
>> scroll up far enough to see the actual command.  I can paste what is there
>> if that may help?  The bottom line:
>>
>>  Names = c("Id", "Level", "AgeGr", "Position", "Height", "Weight", "BMI",
>> "YoYo"), class = "data.frame", row.names = c(NA, -9689L))
>>
>>  Michael - Essentially, I'm looking for differences between "YoYo"
>> outcome for "Positions", "Levels" and accounting for repeated measures
>> using "Id" as a random factor.  So I was able to figure out points 2 and 3.
>>
>>  I've searched for definitions of "Scaled residuals", "Random
>> effects", "Fixed effects", "Correlation of Fixed Effects".  However, I'm
>> confused at the different interpretations I've found.  Or quite possibly,
>> I'm just confused...  What should I be looking out for in these variables?
>>
>>  I've tried to take my analysis smaller, and just look at specifics, to
>> make it simpler.  Such as, comparing YoYo (outcome score) for a
>> Premier_League (Level), 22 (AgeGr) F (Position) with a Premier_League
>> (Level), 22 (AgeGr) M (Position).  How do I convert these into a factors
>> for analysis?
>>
>>  Simple question maybe, but it's not when you can't find the answer!
>>
>>  Thank you,
>>
>>  Josh
>>
>>  On Mon, Apr 27, 2015 at 4:10 PM, Michael Dewey <lists at dewey.myzen.co.uk>
>> wrote:
>>
>>          Dear Joshua
>>
>>   It would also help if you told us what your scientific question was. At
>> the moment we know what R commands you used and have seen the head of your
>> dataset but not why you are doing it.
>>
>>   I would summarise what you have given us as
>>
>>   1 - most ID only occur once
>>   2 - goal keepers do worse than outfield players
>>   3 - older people (presumably in fact age is in years as a continuous
>> variable) do better
>>
>>   On 27/04/2015 12:42, John Kane wrote:
>>
>>   John Kane
>>   Kingston ON Canada
>>
>>           -----Original Message-----
>>   From: joshuamichaeldixon at gmail.com
>>   Sent: Mon, 27 Apr 2015 08:54:51 +0100
>>   To: thierry.onkelinx at inbo.be
>>   Subject: Re: [R] Help Interpreting Linear Mixed Model
>>
>>   Hello Thierry,
>>
>>   No, this isn't homework. Not that young unfortunately.
>>
>>   A few years ago a friend of mine and her daughter were neck-in-neck on
>> who got their Ph.D first. What's this "not that young" business?
>>
>>   BTW, a better way to supply sample data is to use the dput() command.
>>
>>   Do a dput(mydata), copy the results into the email and you have
>> supplied us with an exact copy of your data.
>>
>>   It is possible for many reasons that I will not read in your data, as
>> you supplied it, in the format you have it in.  This can lead to real
>> confusion.
>>
>>           Josh
>>
>>           On 27 Apr 2015, at 08:06, Thierry Onkelinx <
>> thierry.onkelinx at inbo.be>
>>   wrote:
>>
>>   Dear Josh,
>>
>>   Is this homework? Because the list has a no homework policy.
>>
>>   Best regards,
>>
>>   ir. Thierry Onkelinx
>>   Instituut voor natuur- en bosonderzoek / Research Institute for Nature
>>   and Forest
>>   team Biometrie & Kwaliteitszorg / team Biometrics & Quality Assurance
>>   Kliniekstraat 25
>>   1070 Anderlecht
>>   Belgium
>>
>>   To call in the statistician after the experiment is done may be no more
>>   than asking him to perform a post-mortem examination: he may be able to
>>   say what the experiment died of. ~ Sir Ronald Aylmer Fisher
>>   The plural of anecdote is not data. ~ Roger Brinner
>>   The combination of some data and an aching desire for an answer does not
>>   ensure that a reasonable answer can be extracted from a given body of
>>   data. ~ John Tukey
>>
>>   2015-04-27 2:26 GMT+02:00 Joshua Dixon <joshuamichaeldixon at gmail.com>:
>>
>>           Hello!
>>
>>   Very new to R (10 days), and I've run the linear mixed model, below.
>>   Attempting to interpret what it means...  What do I need to look for?
>>   Residuals, correlations of fixed effects?!
>>
>>   How would I look at very specific interactions, such as PREMIER_LEAGUE
>>   (Level) 18 (AgeGr) GK (Position) mean difference to CHAMPIONSHIP 18
>>   GK?
>>
>>   For reference my data set looks like this:
>>
>>   Id Level AgeGr   Position Height Weight BMI YoYo
>>   7451 CHAMPIONSHIP 14 M NA 63 NA 80
>>   148 PREMIER_LEAGUE 16 D NA 64 NA 80
>>   10393 CONFERENCE 10 D NA 36 NA 160
>>   10200 CHAMPIONSHIP 10 F NA 46 NA 160
>>   1961 LEAGUE_TWO 13 GK NA 67 NA 160
>>   10428 CHAMPIONSHIP 10 GK NA 40 NA 160
>>   10541 LEAGUE_ONE 10 F NA 25 NA 160
>>   10012 CHAMPIONSHIP 10 GK NA 30 NA 160
>>   9895 CHAMPIONSHIP 10 D NA 36 NA 160
>>
>>   Many thanks in advance for time and help.  Really appreciate it.
>>
>>   Josh
>>
>>           summary(lmer(YoYo~AgeGr+Position+(1|Id)))
>>
>>   Linear mixed model fit by REML ['lmerMod']
>>   Formula: YoYo ~ AgeGr + Position + (1 | Id)
>>
>>   REML criterion at convergence: 125712.2
>>
>>   Scaled residuals:
>>        Min      1Q  Median      3Q     Max
>>   -3.4407 -0.5288 -0.0874  0.4531  4.8242
>>
>>   Random effects:
>>     Groups   Name        Variance Std.Dev.
>>     Id       (Intercept) 15300    123.7
>>     Residual             16530    128.6
>>   Number of obs: 9609, groups:  Id, 6071
>>
>>   Fixed effects:
>>                 Estimate Std. Error t value
>>   (Intercept) -521.6985    16.8392  -30.98
>>   AgeGr         62.6786     0.9783   64.07
>>   PositionD    139.4682     7.8568   17.75
>>   PositionM    141.2227     7.7072   18.32
>>   PositionF    135.1241     8.1911   16.50
>>
>>   Correlation of Fixed Effects:
>>              (Intr) AgeGr  PostnD PostnM
>>   AgeGr     -0.910
>>   PositionD -0.359 -0.009
>>   PositionM -0.375  0.001  0.810
>>   PositionF -0.349 -0.003  0.756  0.782
>>
>>           model=lmer(YoYo~AgeGr+Position+(1|Id))
>>   summary(glht(model,linfct=mcp(Position="Tukey")))
>>
>>     Simultaneous Tests for General Linear Hypotheses
>>
>>   Multiple Comparisons of Means: Tukey Contrasts
>>
>>   Fit: lmer(formula = YoYo ~ AgeGr + Position + (1 | Id))
>>
>>   Linear Hypotheses:
>>                Estimate Std. Error z value Pr(>|z|)
>>   D - GK == 0  139.468      7.857  17.751   <1e-04 ***
>>   M - GK == 0  141.223      7.707  18.323   <1e-04 ***
>>   F - GK == 0  135.124      8.191  16.496   <1e-04 ***
>>   M - D == 0     1.754      4.799   0.366    0.983
>>   F - D == 0    -4.344      5.616  -0.774    0.862
>>   F - M == 0    -6.099      5.267  -1.158    0.645
>>   ---
>>   Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
>>   (Adjusted p values reported -- single-step method)
>>
>>            [[alternative HTML version deleted]]
>>
>>   ______________________________________________
>>   R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
>>
>>  https://stat.ethz.ch/mailman/listinfo/r-help [
>> https://stat.ethz.ch/mailman/listinfo/r-help] [
>> https://stat.ethz.ch/mailman/listinfo/r-help [
>> https://stat.ethz.ch/mailman/listinfo/r-help]]
>>   PLEASE do read the posting guide
>>   http://www.R-project.org/posting-guide.html [
>> http://www.R-project.org/posting-guide.html] [
>> http://www.R-project.org/posting-guide.html [
>> http://www.R-project.org/posting-guide.html]]
>>   and provide commented, minimal, self-contained, reproducible code.
>>
>>           [[alternative HTML version deleted]]
>>
>>   ______________________________________________
>>   R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
>>   https://stat.ethz.ch/mailman/listinfo/r-help [
>> https://stat.ethz.ch/mailman/listinfo/r-help] [
>> https://stat.ethz.ch/mailman/listinfo/r-help [
>> https://stat.ethz.ch/mailman/listinfo/r-help]]
>>   PLEASE do read the posting guide
>>   http://www.R-project.org/posting-guide.html [
>> http://www.R-project.org/posting-guide.html] [
>> http://www.R-project.org/posting-guide.html [
>> http://www.R-project.org/posting-guide.html]]
>>   and provide commented, minimal, self-contained, reproducible code.
>>
>>   ____________________________________________________________
>>   FREE ONLINE PHOTOSHARING - Share your photos online with your friends
>> and family!
>>   Visit http://www.inbox.com/photosharing [
>> http://www.inbox.com/photosharing] [http://www.inbox.com/photosharing [
>> http://www.inbox.com/photosharing]] to find out more!
>>
>>   ______________________________________________
>>   R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
>>   https://stat.ethz.ch/mailman/listinfo/r-help [
>> https://stat.ethz.ch/mailman/listinfo/r-help] [
>> https://stat.ethz.ch/mailman/listinfo/r-help [
>> https://stat.ethz.ch/mailman/listinfo/r-help]]
>>   PLEASE do read the posting guide
>> http://www.R-project.org/posting-guide.html [
>> http://www.R-project.org/posting-guide.html] [
>> http://www.R-project.org/posting-guide.html [
>> http://www.R-project.org/posting-guide.html]]
>>   and provide commented, minimal, self-contained, reproducible code.
>>
>>   --
>>   Michael
>>   http://www.dewey.myzen.co.uk/home.html [
>> http://www.dewey.myzen.co.uk/home.html] [
>> http://www.dewey.myzen.co.uk/home.html [
>> http://www.dewey.myzen.co.uk/home.html]]
>>
>>  ____________________________________________________________
>>  Can't remember your password? Do you need a strong and secure password?
>>  Use Password manager! It stores your passwords & protects your account.
>>  Check it out at http://mysecurelogon.com/manager [
>> http://mysecurelogon.com/manager]
>>
>> ____________________________________________________________
>> FREE 3D EARTH SCREENSAVER - Watch the Earth right on your desktop!
>> Check it out at http://www.inbox.com/earth
>>
>>
>>
>

	[[alternative HTML version deleted]]



More information about the R-help mailing list