Re: [R] regression and lmer

From: Mike Marchywka <marchywka_at_hotmail.com>
Date: Mon, 18 Apr 2011 08:14:16 -0400



Date: Mon, 18 Apr 2011 03:27:40 -0700
From: lamprianou_at_yahoo.com
To: r-help_at_r-project.org
Subject: Re: [R] regression and lmer

Dear all,

I hope this is the right place to ask this question.  
( hotmail not marking your text sorry I can';t find option to change that  )
opinions vary but mods seem to let this by and personally it seems appropriate to discuss general questions about what R can do.

( end my text )

I am reviewing a research where the analyst(s) are using a linear regression model. The dependent variable (DV) is a continuous measure. The independent variables (IVs) are a mixture of linear and categorical variables.

( my text )

No one wants to do homework or do your job for free but open free peer review should not be a problem philosophically.

( / my text, afraid to use less-than for inciting hotmail )

The
author investigates whether performance (DV - continuous linear) is a function of age (continuous IV1 - measured in years), previous performance (continuous IV2), country (categorical IV3 - six countries), the percentage of PhD graduates in each country (continuous IV4 - country level data - apparently only six different percentages since we have only six countries) and population of country (continuous IV5 - country level data - again only six numbers here, one for each country population).

My own opinion is that the lm function cannot be used with country level data as IVs (for example IV4 and IV5 cannot be entered into the model because they are country level data). If IV4 and IV5 are included in the model, it is possible that the model will not be able to be defined because we only have six countries and it is very likely that the levels of counties (IV3) may be confounding with IV4 and IV5. This also calls for multicollinearity issues, right? I would like to suggest to the analyst to use lmer using the IV3 as a random variable and IV4 and IV5 as IV at the second level of the two-level model.

The questions are: (a) Is it true that IV4 and IV5 cannot be entered in a one-level regression if we also have IV3?, (b) can I use an lm function to check for multicollinearity between IV3, IV4 and IV5? and (c) If we use a two-level regression model, does lmer cope well with only six coutnries as a random effect?

( my txt)

So you have presumably a large number of test subjects per country and a small number
( n~6 ) of countries. You could ask a number of questions such as, " do the mean performances
change from country to country by more than that expected given the observed distributions of performances within country?" You could also ask a question like, " if I try to describe performance as a function of country attriubte what fitting parameters minimize an error between fit and observation?" Apparently author tried to write an expression like

average_performance= a[country_index] + m1*some_attribute_of_country + m2* some_other_attribute_of_country + b

and then expected the fittring algoright to pick a[i], b,m1, and m2 in such a way as to minimize the resulting error. The reported fits hopefully minimize the error function but then you need to exmaine the second derivative in various directions, so you have to ask how the error varies as you change a[i],b,m1, and m2. ( Ignore b right now and assume it s included in a[i]). I guess if you can find a direction where the error can not change due to these contraints then it would seem to be impossible for the fit to come up with unique values. If you change each a[i] and m1 by some amounts, for example, can you pick those amounts to not change anything?

( /my text )

Thank you for your help

Jason

        [[alternative HTML version deleted]]



R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.                                                

R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Mon 18 Apr 2011 - 12:19:23 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Mon 18 Apr 2011 - 12:40:30 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive