[R] PROC MIXED RANDOM equivalence in R nlme

Dennis F. Kahlbaum kbomb at umich.edu
Fri Aug 11 16:56:52 CEST 2017


Hi Robert:

Thanks for your response, as well.

I'm sorry. As you've discovered, I made some capitalization errors when 
posting my R code.  The actual R code does use "lme", "data", and "random".

The dataframe is indeed named "emiss" and each item in the formula is a 
column in the dataframe. I used the following R code to read in the 
comma-delimited file (first line contains headings):

     emiss <- read.table("data.csv", header=TRUE, sep=",") # NA for 
missing data (default)

A review of the dataframe "emiss" shows that all items, except "STUDY", 
"VEHICLE" and "NEW", are of type "num". The items "STUDY", "VEHICLE", 
and "NEW" are of type "factor".

I will work on creating some "fake data". However, since I don't have 
SAS v6.12 anymore, I can't provide the so-called "correct results".

As for the R hang, your wording is correct: "R gets caught in a 
processing loop that produces no errors or warnings", even after 15 
minutes on an 8-core Mac Pro.  As I recall, the SAS code produced an 
answer in 15 seconds or less on a single core Mac II back in 1999.

Thanks again!

--Dennis



On 8/11/17 9:27 AM, Robert Baer wrote:
>
>
> On 8/10/2017 8:34 AM, Dennis F. Kahlbaum wrote:
>> -- snip --
> I don't have real help, but I'll remind you that R is case sensitive, 
> and it looks like that will be at least one problem in the solution 
> your are working on below:
> lme not LME
> data not DATA
> random = RANDOM
>> ------------------------------------------------------------------
>>
>> The R code I've devised for the PROC MIXED statement is shown below:
>>
>> ------------------------------------------------------------------
>> FitTHC <- LME(ln_thc ~ rv + t5 + t9 + ar + ol + ox + su + bz,
>>           DATA = emiss,
>>           RANDOM = ??????? )
>> ------------------------------------------------------------------
>>
>> As indicated, the problem I'm having is in constructing the 
>> equivalent code for the RANDOM and any remaining settings. I've tried
>>
>> RANDOM = ~1 + rv + t5 + t9 + ar + ol + ox + su + bz | new)
>>
>> but R hangs 
> Are the items in your random formula columns in a dataframe named 
> emiss?  Do they have data types?  Even if the data are proprietary 
> some fake data can make the problem more concrete.
> You are saying "gets caught in a processing loop that produces no 
> errors or warnings"???
>
>> and never produces a result. Therefore, what is the equivalent code 
>> for the SAS RANDOM?
>>
>> Thanks!
>>
>> ______________________________________________
>> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
>> 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.
>



More information about the R-help mailing list