> On Aug 11, 2017, at 7:56 AM, Dennis F. Kahlbaum <kb...@umich.edu> wrote:
>
> 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.
You really should start over. Compose a correct version of code and data and
resubmit to the place this questions _should_ have gone in the first place: the
R-SIG-Mixed-models mailing list (correct spelling and link should be on the
listinfo page.)
>
> 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@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.
>>
>
> ______________________________________________
> R-help@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.
David Winsemius
Alameda, CA, USA
'Any technology distinguishable from magic is insufficiently advanced.'
-Gehm's Corollary to Clarke's Third Law
______________________________________________
R-help@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.
