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.
>
> On Aug 11, 2017, at 7:56 AM, Dennis F. Kahlbaum <kbomb at 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 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. >> > > ______________________________________________ > 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.David Winsemius Alameda, CA, USA 'Any technology distinguishable from magic is insufficiently advanced.' -Gehm's Corollary to Clarke's Third Law
Hi David: Thanks for your response and suggestions. --Dennis On 8/11/17 2:03 PM, David Winsemius wrote:>> On Aug 11, 2017, at 7:56 AM, Dennis F. Kahlbaum <kbomb at 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 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. >> ______________________________________________ >> 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. > David Winsemius > Alameda, CA, USA > > 'Any technology distinguishable from magic is insufficiently advanced.' -Gehm's Corollary to Clarke's Third Law > > > > >